Poker Flat Incoherent Scatter Radar (PFISR) Observations of E-region Neutral Winds

This file contains Poker Flat Incoherent Scatter Radar (PFISR) E-region Neutral Winds Data.  This is associated with a submission to the Journal of Geophysical Research and once the article is accepted, I will put in the DOI for that article.

These data correspond to monthly data files that include the E-region neutral winds and other parameters for the from March 2013-June 2019.

These data were processed from the following files found in: https://data.amisr.com/database/tmp/Kaeppler/winds/ and https://data.amisr.com/database/tmp/Kaeppler/missing_IPY.tar.gz Please note that the error on the line of sight velocities may have been overestimated in these data and we scaled them by a eVLOS/sqrt(10).  Interested persons should contact Ashton Reimer or Roger Varney at SRI International for more information about these data, please see amisr.com

If you plan to work with the data, I would appreciate it if you would contact me via email.  We are very interested in collaborating.

File Documentation:

Please see the change log:

v0.5.1 - 10-23-2019: changed chi square to 0.01 for lower boundary

v0.5.3 - 12-02-2019: Added in that now passing in the Chi2 and Fitcode filtering by Config file
            Bigger change that I am scaling the AC dVlos by some sort of factor while Ashton figures this out.
            We decided that a conversative scaling would be to reduce the dVLOS by 1/sqrt(10).
            The chi square produced in the data Ashton sent me typically was around 0.01, so the uncertaintiies on the LOS velocities
            may be over estimated.  So we are just changing this as a temporary fix while Ashton fixes the uncertainty estimation.

v0.5.5 02 01 2020 - Added in calculation of Coriolis, Centrifugal, and Lorentz forcing
v0.5.5 02 10 2020 - Added a correction to qvert so that way I can calculate the lorentz term.
                    Found an error where qvert = 0 in the if statement goes to false.

v0.5.5 02 15 2020 - Put in  nuInscaler into the main program, scaling ALL kappas by the scaler number

v0.5.6 02 28 2020 -- Added some more vlos diagostics and the calculation of the scale height. Added Altitude offset

v0.5.6.2020.03.12_nuin_fracoff - testing putting in the Brekke formula for ion neutral collision frequency and took out frac

v0.5.7.2020.04.10 - Put in Ashton's revised ion neutral collision frequency formulas into IO.
                    Also wrote a testscript and at least for the file I used was only different by 2.5%.
                    Revised where the mag data is being pulled from since the URL is deprecated
                    Added in Kappa which is now being interpolate - plan to see where kappa =1 is located for the paper.
                    commented out nuin scaler just so I am not chasing my tail

                    test v0.5.7.2020.04.13_org commented back in original ion neutral collision frequency method
                    possible mistake that not summing up properly.

                    test v0.5.7.2020.04.13_newnuin_orgsum_noTr800 - new formula for nuin except took off Tr>800.
                    I expect this should be almost the same as before since the formulas are basically the same.
                    did the original sum using frac[0] and frac[1] want to see if I am underestimating

                    v0.5.7.2020.04.13_newnuin_orgsum_yesTr800 - same as above except now including Tr>800.

                    'v0.5.7.2020.04.13_newnuin_newsum_noTr800' - using the new sum now and new col freq

                    v0.5.7.2020.04.13_updatedorg - updated original uses original method but including the NO term

v0.6.0.2020.04.15 -- Now think I have the new ion neutral collision frequency working and validated.
                    Found a mistake in how I was calculating the ion neutral collision frequency that
                    the fraction weight I was using only included the O+ and O2+ terms and not NO+
                    Turns out I was basically weighting by about 0.5, so I was effectively reducing the
                    ion neutral collision frequency by about a factor of 0.5 or less...
                    From this point forward need to start using any results from > v0.6
                    This revision has changed previous results signi

v0.6.0.2020.04.21 -- updated to now include the temperature correction for the O2+

v0.6.1.2020.04.23 -- made a number of changes to the geomagnetic files and reprocessed from CDAweb.
                     Wrote new code to be able to process the files from CDAweb in the new format.
                     Also changed the geomagnetic data files

v0.6.1.2020.06.07 -- changed the generation of Monthly files to hopefully be in order now
                    Added in missingIPY files given to me by ashton, maybe improve data covarege
                    Some work going to need to be done to make sure that all of the 10, 15, and 20 minute data are there.