In [2]:

    

#Test igrfpy by plotting some igrf main field
import igrfpy
import datetime
import numpy as np
import matplotlib.pyplot as pp
%load_ext autoreload


    

In [17]:

    

%autoreload
nvals = 50
times = datetime.datetime(2000,6,17,12,0,0)
lats = [-90.+i*(180./(nvals-1)) for i in range(nvals)]
lons = 0.
alts = 1.
be,bn,bu = igrfpy.getmainfield(times,lats,lons,alts)


    

    




Running IGRF for 50 values

In [18]:

    

f = pp.figure()
a = f.add_subplot(1,1,1)
a.plot(lats,be,'b-',label='B_east')
a.plot(lats,bn,'r-',label='B_north')
a.plot(lats,bu,'g-',label='B_up')
a.legend()
pp.show()