We present a general framework for two-dimensional finite difference modeling of magnetotelluric data in the presence of general anisotropy. Our approach is modular, allowing differential operators for a range of formulations of the governing equations, defined on several possible discrete grids, to be constructed from a basic set of first...