Old IGRF-13 in pyamps/mlt_utils.py

[pacesm]

The currently latest version 1.7 seems to be still using embedded IGRF-13 coefficients failing for epochs > 2025.0

pyAMPS/pyamps/mlt_utils.py

Lines 39 to 52 in fbad688

	# first make arrays of IGRF dipole coefficients. This is used to make rotation matrix from geo to cd coords
	# these values are from https://www.ngdc.noaa.gov/IAGA/vmod/igrf12coeffs.txt
	time =[1900.0, 1905.0, 1910.0, 1915.0, 1920.0, 1925.0, 1930.0, 1935.0, 1940.0, 1945.0, 1950.0, 1955.0, 1960.0, 1965.0, 1970.0, 1975.0, 1980.0, 1985.0, 1990.0, 1995.0, 2000.0, 2005.0, 2010.0, 2015.0, 2020.0, 2025.0]
	g10 = [-31543, -31464, -31354, -31212, -31060, -30926, -30805, -30715, -30654, -30594, -30554, -30500, -30421, -30334, -30220, -30100, -29992, -29873, -29775, -29692, -29619.4, -29554.63, -29496.57, -29441.46, -29404.8]
	g11 = [ -2298, -2298, -2297, -2306, -2317, -2318, -2316, -2306, -2292, -2285, -2250, -2215, -2169, -2119, -2068, -2013, -1956, -1905, -1848, -1784, -1728.2, -1669.05, -1586.42, -1501.77, -1450.9]
	h11 = [ 5922, 5909, 5898, 5875, 5845, 5817, 5808, 5812, 5821, 5810, 5815, 5820, 5791, 5776, 5737, 5675, 5604, 5500, 5406, 5306, 5186.1, 5077.99, 4944.26, 4795.99, 4652.5]
	g10sv = 5.7 # secular variations
	g11sv = 7.4
	h11sv = -25.9
	g10.append(g10[-1] + g10sv * 5) # append 2025 values using secular variation
	g11.append(g11[-1] + g11sv * 5)
	h11.append(h11[-1] + h11sv * 5)
	igrf_dipole = pd.DataFrame({'g10':g10, 'g11':g11, 'h11':h11}, index = time)
	igrf_dipole['B0'] = np.sqrt(igrf_dipole.g10**2 + igrf_dipole.g11**2 + igrf_dipole.h11**2)