search

NOAA-SWPC / IPE

Ionosphere Plasmasphere Electrodynamics Model
GNU General Public License v3.0
6 stars 10 forks source link

Plasma Depletion in IPE #25

Closed twfang closed 3 years ago

twfang commented 4 years ago

The Kp dependent plasma depletion is included in IPE (plamsa_depletion branch). To test the code, 2015 St. Pactirk Day's storm simulation is carried out and several plots are documented in this issue.

twfang commented 4 years ago

Here is the TEC from 2015 03/18 during the recovery phase of the storm. A slight decrease can be seen in TEC with the depletion code. I will work on plots showing an equatorial plane from the native IPE grid.

TEC from default run 20150318_ori

TEC with depletion code 20150318_depleted

twfang commented 4 years ago

Here are some O+ results in the native IPE grid. The plots are in log scale, both at 3UT on 2015 3/18, which is 3 hours after the peak DST. Not sure if this is a good time to look into. But you can decrease density once the new code is implemented.

Original 201503180300_ori

With depletion code 201503180300_depleted

twfang commented 4 years ago

Here is the plot with and without the plasma depletion code. They are from the same time 2015/03/18 3UT. The values come from the apex of each tube (equatorial plane).

Plots on top are the H+ density while on the bottom are the He+. mp=1 is where the magnetic longitude=0 (i.e. American sector, 22LT) while mp=40 is where geographic longitude ~ 100E (~10LT). X-axis is the magnetic latitude of the foot-print each tube.

density_equ

Please let me know if this is the correct plot for seeing the effect. Thank you!

twfang commented 3 years ago

More work is required for this task. We will open a new issue once the new development is initiated.