|
|
@ -3,7 +3,7 @@ |
|
|
|
generate field tables from IGRF12 |
|
|
|
generate field tables from IGRF12 |
|
|
|
''' |
|
|
|
''' |
|
|
|
|
|
|
|
|
|
|
|
import pyigrf12 as igrf |
|
|
|
import igrf12 as igrf |
|
|
|
import numpy as np |
|
|
|
import numpy as np |
|
|
|
import datetime |
|
|
|
import datetime |
|
|
|
from pathlib import Path |
|
|
|
from pathlib import Path |
|
|
@ -54,7 +54,7 @@ if __name__ == '__main__': |
|
|
|
|
|
|
|
|
|
|
|
for i,lat in enumerate(lats): |
|
|
|
for i,lat in enumerate(lats): |
|
|
|
for j,lon in enumerate(lons): |
|
|
|
for j,lon in enumerate(lons): |
|
|
|
mag = igrf.igrf(date=date, glat=lat, glon=lon, alt=0., isv=0, itype=1) |
|
|
|
mag = igrf.igrf(date, glat=lat, glon=lon, alt_km=0., isv=0, itype=1) |
|
|
|
intensity_table[i][j] = mag.total/1e5 |
|
|
|
intensity_table[i][j] = mag.total/1e5 |
|
|
|
inclination_table[i][j] = mag.incl |
|
|
|
inclination_table[i][j] = mag.incl |
|
|
|
declination_table[i][j] = mag.decl |
|
|
|
declination_table[i][j] = mag.decl |
|
|
|