pyglow
pyglow copied to clipboard
Getting nan's for 06-2018
Point.ap
, Point.ap_daily
, and Point.apmsis
are all nan's for 06-2018. I have updated the indices and the problem persists.
This is the same issue as #105. Our source of Kp/Ap was NOAA/NGDC, but as of mid-2018 they stopped providing this information. We should use GFZ Potsdam.
Unfortunately there is no easy workaround until this problem is fixed. If you know the values already you can specify them using user_ind=True:
pt = pyglow.Point(t, lat, lon, alt, user_ind=True)
pt.ap = 10.
pt.ap_daily = 10.
pt.apmsis = [10., 10., 10., 10., 10., 10.]
pt.run_msis()
A quick fix for @shaharkadmiel is here. I use the Kp/Ap fom GFZ Potsdam and F10.7 from OMNIWeb to reproduce kpap as the same format as NOAA/NGDC manually. Unzip the attached file, copy and paste the two files to your installed ../pyglow/kpap folder (not in the folder you downloaded) and it works fine.