Python load_model_shc Examples

Example #1

File: model_loaders.py Project: shineusn/pyMagMod

class TestComposedModelDiffConstrained(TestCase, ComposedModelTestMixIn):
    parameters = ("time", "location")
    options = {"scale": [1, 1, -1]}
    components = [
        (load_model_shc_combined(CHAOS6_STATIC,
                                 CHAOS6_CORE_LATEST,
                                 to_mjd2000=decimal_year_to_mjd2000_simple),
         1.0, {
             "max_degree": 15
         }),
        (load_model_shc(CHAOS6_CORE_LATEST,
                        to_mjd2000=decimal_year_to_mjd2000_simple), -1.0, {
                            "min_degree": 5,
                            "max_degree": 15
                        }),
    ]
    reference_values = (6201.125, (30.0, 40.0, 6400.0),
                        (31512.85733779784, 2660.412538490814,
                         30567.851101678763))
    validity = (-1058.4945, 7195.49805)

Example #2

File: field_lines.py Project: shineusn/pyMagMod

 def model(self):
     if not hasattr(self, "_model"):
         self._model = load_model_shc(IGRF12,
                                      interpolate_in_decimal_years=True)
     return self._model

Example #3

File: model_loaders.py Project: shineusn/pyMagMod

 def load(self):
     return load_model_shc(CHAOS6_CORE_LATEST,
                           validity_start=2000.0,
                           validity_end=2018.0)

Example #4

File: model_loaders.py Project: shineusn/pyMagMod

 def load(self):
     return load_model_shc(CHAOS6_CORE_LATEST)

Example #5

File: model_loaders.py Project: shineusn/pyMagMod

 def load(self):
     return load_model_shc(CHAOS6_STATIC)

Example #6

File: model_loaders.py Project: shineusn/pyMagMod

 def load(self):
     return load_model_shc(MF7)

Example #7

File: model_loaders.py Project: shineusn/pyMagMod

 def load(self):
     return load_model_shc(LCS1)

Example #8

File: model_loaders.py Project: shineusn/pyMagMod

 def load(self):
     return load_model_shc(SIFM)

Example #9

File: model_loaders.py Project: shineusn/pyMagMod

 def load(self):
     return load_model_shc(IGRF12, interpolate_in_decimal_years=True)

Example #10

File: earth_mag_mods.py Project: zhixin-xue/geoist

)
from magmod.util import datetime_to_decimal_year, vnorm

from magmod.magnetic_model.parser_mio import parse_swarm_mio_file
from magmod.magnetic_model.tests.data import SWARM_MIO_SHA_2_TEST_DATA
import magmod._pymm as pymm

print(SWARM_MIO_SHA_2_TEST_DATA) #DIFI4 is a type of MIO SHA model
print(mjd2000_to_decimal_year(7305))
print(mjd2000_to_decimal_year([5479., 7305., 6392.0]))
d1 = dt.datetime(2015,1,1)   # import time , location(lat, lon)
d11 = datetime_to_decimal_year(d1) # datetime to decimal year
loc = (30.0, 40.0, 1000.0)
wmm2015 = load_model_wmm(WMM_2015)  #load wmm2015 model
igrf11 = load_model_igrf(IGRF11)    #load igrf11 model
igrf12 = load_model_shc(IGRF12, interpolate_in_decimal_years=True)    #load igrf12 model

igrf13 = load_model_igrf(IGRF13)
wmm2020 = load_model_wmm(WMM_2020)

emm = load_model_emm(EMM_2010_STATIC, EMM_2010_SECVAR)  #load emm model
options = {"scale": [1, 1, -1]}   #-1 is Z direction

wmm2020.eval(decimal_year_to_mjd2000(d11), loc, 0, 0, **options) 
igrf13.eval(decimal_year_to_mjd2000(d11), loc, 0, 0, **options)

wmm2015.eval(decimal_year_to_mjd2000(d11), loc, 0, 0, **options) # 0,0 mean input,output using GEODETIC_ABOVE_WGS84
igrf11.eval(decimal_year_to_mjd2000(d11), loc, 0, 0, **options)
igrf12.eval(decimal_year_to_mjd2000(d11), loc, 0, 0, **options) #North-X, East-Y, Vertical-Z
vnorm(igrf12.eval(decimal_year_to_mjd2000(d11), (30.0, 40.0, 1000.0), 0, 0)) #Total intensity
mjd2000_to_decimal_year([5479., 7305., 6392.0])