pinga-lab / magnetic-ellipsoid Goto Github PK

The coordinate transformation matrices should be defined directly by using the "geological angles" strike, dip and rake, like Clark et al. (1986, p. 182) did. By using their approach, however, the matrix for oblate ellipsoids must be modified. I will try to implement the equations presented by Clark et al. (1986) and modify the equations for oblate ellipsoids.

• Clark, D., Saul, S., and Emerson, D.: Magnetic and gravity anomalies of a triaxial ellipsoid, Exploration Geophysics, 17, 189–200, 1986.

This is great stuff. Do you guys have an implementation of this for python3, and using your new software libraries (verde, etc.)?

The three test_*_ignore_missing_prop() tests are failing for me using numpy 1.11.
Do they pass for you under Windows @birocoles?

The fields produced by ellipsoids with approximatelly the same orientation and shape should be close to each other. These numerical test fail and I guess the problem is the orientation matrices. I need to verify the position and orientation of the ellipsoids.