irfpy.moon.anomaly

Magnetic anomaly model.

Purucker’s model is implemented in AnomalyModelPME.

class irfpy.moon.anomaly.AnomalyModelPME[source]

Bases: object

Purucker’s Model.

Refer to M. E. Purucker 2010 in JGR planets. http://core2.gsfc.nasa.gov/research/purucker/moon_2010/

>>> am = AnomalyModelPME()
>>> print(am.br.shape)
(361, 361)
>>> print(am.lonlist[0:3])
[0. 1. 2.]
>>> print(am.latlist[0:3])
[90.  89.5 89. ]
>>> print(am.lonlist[130])
130.0
>>> print(am.latlist[75])
52.5
>>> print(am.br[130, 75])
0.69
>>> print(am.bt[130, 75])
-0.53
>>> print(am.bp[130, 75])
-0.49

Todo

Better user interface. Probably by giving longitude and latitude, user may get a vector.

Load a data file.

The data is obtained from http://core2.gsfc.nasa.gov/research/purucker/moon_2010/

The corrective model with file name brtps_d1_170_alt30_har14tol15_dip6cpadpr_phase30_correlation.xy is supported here. Download it, and place it to your prefered folder and specify by .irfpyrc, which should be [moon] section purucker_anomaly_data entry.

irfpy.moon.anomaly.doctests()[source]