pyglow icon indicating copy to clipboard operation
pyglow copied to clipboard

Getting nan's for 06-2018

Open shaharkadmiel opened this issue 5 years ago • 2 comments

Point.ap, Point.ap_daily, and Point.apmsis are all nan's for 06-2018. I have updated the indices and the problem persists.

shaharkadmiel avatar Aug 19 '19 13:08 shaharkadmiel

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()

bharding512 avatar Aug 19 '19 17:08 bharding512

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.

2018_2019.zip

FasilGibdaw avatar Aug 31 '19 12:08 FasilGibdaw