search

NOAA-SWPC / GSMWAM-IPE

Coupled NEMS app for WAM-IPE
5 stars 11 forks source link

Fixes for NmF2 and hmF2 in SWIO #30

Closed twfang closed 4 years ago

twfang commented 4 years ago

For NmF2 and hmF2, we need to search from the top of the profile to find the density peak at F-region. The current method simply finds the maximum value resulting in a peak at E-region at some locations. Let's discuss how to correct the calculation here.

twfang commented 4 years ago

I have looked into the updated outputs. There are two issues to hmF2.

  1. Model behavior at 200 km can be a problem and largely comes from the minor species
  2. Because the problem 1 tends to create a larger density at around 200km, it leads to a peak at around 200km instead of F-region or negative values associated with the quadratic fit.

201303160100

TEC_NmF2_hmF2_03160100

Ne density profiles for a negative point

neprf_negetive_03160100

Look into a neighboring point where hmF2 is near 200km

neprf_low_03160100
twfang commented 4 years ago

I have used Raffaele's code to extend his run to the 2nd day. Just to see how things behave. Here are the results from 201303170100. The problem for the odd point seems to come from O+ instead of minor ions.

TEC_NmF2_hmF2_03170100 neprf_odd_03170100
twfang commented 4 years ago

The current method searches for peaks over a range of 200 and 1000 km and picks the largest density peak for NmF2 and hmF2.