IPE Physics: density structures

[twfang]

Running George's no nan branch with the empirical electric field and standalone IPE with fixed drivers. So, days don't mean much in these plots.

0316

0317

0318

[twfang]

A few IPE standalone runs with different settings (with chemical equilibrium or diffusion). It seems like with diffusive solution turns on for both species, the TEC has fewer structures especially in the nighttime and early morning. In the older version of IPE (before refactored), both species were set to be the diffusive solution. Therefore, I suggest we set bot IHEPLS and INPLS to 1 in IPE_Plasma_Class.F90.

[twfang]

The high density bends at the high-latitude region are caused by O+ density. The problem became less severe when the transport is turned off. The figure shown here is from the 3rd day standalone IPE run with diffusive solution of N+/He+ and fixed drivers.

O+ at 350 km

O+ at 150 km

[twfang]

With setting the transport starting from lp=30 instead of going all the way to the pole, IPE looks quite reasonable. Still need to figure out what is wrong in high latitude transport.

[twfang]

Geroge is looking into the ksi_fac in the transport subroutine. It appears that the values look odd at high-latitudes and leads to strange transport for all ions. If we set the ksi_fac to 1 for lp 1 to 30 (the original one use lp from 1 to 17), the problem seems to be disappeared.

Below figures shows the WAM-IPE coupled run from 0320. This is results from the 5th-day run.

No-dynamo

self-consistent dynamo

[gmillward]

ksi_fac(alt) - ksi_fac(202km), as a function of altitude, for a few flux-tubes (around lp=20). Needs investigating:

[naomimaruyama]

George, this plot is very interesting. What do the profiles look like for Lp>30? thanks. naomi

On Mon, Feb 3, 2020 at 11:29 AM gmillward notifications@github.com wrote:

ksi_fac(alt) - ksi_fac(202km), as a function of altitude, for a few flux-tubes (around lp=20). Needs investigating: [image: ksi] https://user-images.githubusercontent.com/24465054/73679807-6273d680-4678-11ea-9f13-f22f5f43a8f8.png

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/NOAA-SWPC/IPE/issues/16?email_source=notifications&email_token=ADIZZYXN2XCDJN7LIF7PQQ3RBBO7VA5CNFSM4KHXTXMKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEKU4LTQ#issuecomment-581551566, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADIZZYQHIRBM3P4KMSJYPW3RBBO7VANCNFSM4KHXTXMA .

--

Naomi Maruyama CIRES, Univ. of Colorado Boulder NOAA Space Weather Prediction Center phone: +1 303.497.4857 <%28303%29%20497-4857> email: naomi.maruyama@noaa.gov

[gmillward]

Naomi

Not sure - will look.

Could you try plotting this out with your version? I'm looking here at ksi_fac for mp=40 and lp = 20 (and 21, 22, 23, 24) as a function of height.

ksi_fac should be a smooth function of height - not jagged - just wondering what your original code shows - this is just for the first time step on the 20th - but I think any run can be used

Cheers

George.

On Mon, Feb 3, 2020 at 11:37 AM Naomi Maruyama notifications@github.com wrote:

George, this plot is very interesting. What do the profiles look like for Lp>30? thanks. naomi

On Mon, Feb 3, 2020 at 11:29 AM gmillward notifications@github.com wrote:

ksi_fac(alt) - ksi_fac(202km), as a function of altitude, for a few flux-tubes (around lp=20). Needs investigating: [image: ksi] < https://user-images.githubusercontent.com/24465054/73679807-6273d680-4678-11ea-9f13-f22f5f43a8f8.png

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub < https://github.com/NOAA-SWPC/IPE/issues/16?email_source=notifications&email_token=ADIZZYXN2XCDJN7LIF7PQQ3RBBO7VA5CNFSM4KHXTXMKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEKU4LTQ#issuecomment-581551566 , or unsubscribe < https://github.com/notifications/unsubscribe-auth/ADIZZYQHIRBM3P4KMSJYPW3RBBO7VANCNFSM4KHXTXMA

.

--

Naomi Maruyama CIRES, Univ. of Colorado Boulder NOAA Space Weather Prediction Center phone: +1 303.497.4857 <%28303%29%20497-4857> email: naomi.maruyama@noaa.gov

— You are receiving this because you were assigned. Reply to this email directly, view it on GitHub https://github.com/NOAA-SWPC/IPE/issues/16?email_source=notifications&email_token=AF2U5HXUG4HDLAPQSEWJWFTRBBQAFA5CNFSM4KHXTXMKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEKU5I5Q#issuecomment-581555318, or unsubscribe https://github.com/notifications/unsubscribe-auth/AF2U5HTIRNLL3GLQDS5QSSLRBBQAFANCNFSM4KHXTXMA .

--

Dr. George Millward Cooperative Institute for Research in Environmental Sciences (CIRES) University of Colorado, Boulder

Visiting scientist at NOAA Space Weather Prediction Center Boulder, CO

george.millward@noaa.gov

[naomimaruyama]

George,

sure, let me have a look. will send you the plot shortly

On Mon, Feb 3, 2020 at 11:37 AM Naomi Maruyama notifications@github.com wrote:

George, this plot is very interesting. What do the profiles look like for Lp>30? thanks. naomi

On Mon, Feb 3, 2020 at 11:29 AM gmillward notifications@github.com wrote:

ksi_fac(alt) - ksi_fac(202km), as a function of altitude, for a few flux-tubes (around lp=20). Needs investigating: [image: ksi] < https://user-images.githubusercontent.com/24465054/73679807-6273d680-4678-11ea-9f13-f22f5f43a8f8.png

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub < https://github.com/NOAA-SWPC/IPE/issues/16?email_source=notifications&email_token=ADIZZYXN2XCDJN7LIF7PQQ3RBBO7VA5CNFSM4KHXTXMKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEKU4LTQ#issuecomment-581551566 , or unsubscribe < https://github.com/notifications/unsubscribe-auth/ADIZZYQHIRBM3P4KMSJYPW3RBBO7VANCNFSM4KHXTXMA

.

--

Naomi Maruyama CIRES, Univ. of Colorado Boulder NOAA Space Weather Prediction Center phone: +1 303.497.4857 <%28303%29%20497-4857> email: naomi.maruyama@noaa.gov

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/NOAA-SWPC/IPE/issues/16?email_source=notifications&email_token=ADIZZYR5PQN7DLPJ6MTN7T3RBBQAFA5CNFSM4KHXTXMKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEKU5I5Q#issuecomment-581555318, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADIZZYQMDHIWN3MKVNOAPDLRBBQAFANCNFSM4KHXTXMA .

--

Naomi Maruyama CIRES, Univ. of Colorado Boulder NOAA Space Weather Prediction Center phone: +1 303.497.4857 <%28303%29%20497-4857> email: naomi.maruyama@noaa.gov

[gmillward]

Naomi

Yeh, problem is still there (for me) with the tubes > 30 (if anything it's even larger)

[twfang]

After fixing N+/He+ movement and set transport for lp<30 to 1, the model appears to be very stable. Here is the TEC of March 29. This is the coupled WAM-IPE with self-consistent electrodynamics and time-varying drivers. The model run started on March 16.

[naomimaruyama]

George,

I checked the ksi_factor in my old simulation. they are almost constant for the height range between 200 and 340km: ~ 0.999991 The values get slightly smaller toward the apex altitude: 0.999978 I do not bother to make a plot because it is nearly constant. you can see the actual output file: hera:/scratch2/NCEPDEV/stmp1/Naomi.Maruyama/save/mpi20160330v6/run15/1580777321_ipe_hera_intel_serial2/fort.9000 the 2nd column is ht[meter], and the 3rd is factor_ksi for the flux tube mp=40, lp=20, as requested.

Remember the factor_ksi is the ratio of the magnetic field strength of the flux tube of (mp, lp) to an imaginary flux tube at t=t-1(so called stepback issue), thus it depends on the EXB drift and the time step. My values are very close to 1, because the time step is very small, 5sec in this case. what are your time step and EXB drift? we should use the same values for better comparison. Thanks.

On Mon, Feb 3, 2020 at 11:51 AM gmillward notifications@github.com wrote:

Naomi

Not sure - will look.

Could you try plotting this out with your version? I'm looking here at ksi_fac for mp=40 and lp = 20 (and 21, 22, 23, 24) as a function of height.

ksi_fac should be a smooth function of height - not jagged - just wondering what your original code shows - this is just for the first time step on the 20th - but I think any run can be used

Cheers

George.

On Mon, Feb 3, 2020 at 11:37 AM Naomi Maruyama notifications@github.com wrote:

George, this plot is very interesting. What do the profiles look like for Lp>30? thanks. naomi

On Mon, Feb 3, 2020 at 11:29 AM gmillward notifications@github.com wrote:

ksi_fac(alt) - ksi_fac(202km), as a function of altitude, for a few flux-tubes (around lp=20). Needs investigating: [image: ksi] <

https://user-images.githubusercontent.com/24465054/73679807-6273d680-4678-11ea-9f13-f22f5f43a8f8.png

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <

https://github.com/NOAA-SWPC/IPE/issues/16?email_source=notifications&email_token=ADIZZYXN2XCDJN7LIF7PQQ3RBBO7VA5CNFSM4KHXTXMKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEKU4LTQ#issuecomment-581551566

, or unsubscribe <

https://github.com/notifications/unsubscribe-auth/ADIZZYQHIRBM3P4KMSJYPW3RBBO7VANCNFSM4KHXTXMA

.

--

Naomi Maruyama CIRES, Univ. of Colorado Boulder NOAA Space Weather Prediction Center phone: +1 303.497.4857 <%28303%29%20497-4857> email: naomi.maruyama@noaa.gov

— You are receiving this because you were assigned. Reply to this email directly, view it on GitHub < https://github.com/NOAA-SWPC/IPE/issues/16?email_source=notifications&email_token=AF2U5HXUG4HDLAPQSEWJWFTRBBQAFA5CNFSM4KHXTXMKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEKU5I5Q#issuecomment-581555318 , or unsubscribe < https://github.com/notifications/unsubscribe-auth/AF2U5HTIRNLL3GLQDS5QSSLRBBQAFANCNFSM4KHXTXMA

.

--

Dr. George Millward Cooperative Institute for Research in Environmental Sciences (CIRES) University of Colorado, Boulder

Visiting scientist at NOAA Space Weather Prediction Center Boulder, CO

george.millward@noaa.gov

— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/NOAA-SWPC/IPE/issues/16?email_source=notifications&email_token=ADIZZYRNJFP4DPDAY45MXKLRBBRRZA5CNFSM4KHXTXMKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEKU6TDQ#issuecomment-581560718, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADIZZYUEKKY5LRH7W2UXYHDRBBRRZANCNFSM4KHXTXMA .

--

Naomi Maruyama CIRES, Univ. of Colorado Boulder NOAA Space Weather Prediction Center phone: +1 303.497.4857 <%28303%29%20497-4857> email: naomi.maruyama@noaa.gov

[twfang]

The run for Apr 14 with time-varying drivers and self-consistent electrodynamics. Ksi is set to 1 for lp<30. The run starting from March 16 and appears to be stable.

[gmillward]

[twfang]

Further discussion is now moved to #17

[twfang]