IPE Physics: density structures

[twfang]

[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 ***@***.***> 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 <#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 ***@***.***> 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 > < #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: ***@***.*** — You are receiving this because you were assigned. Reply to this email directly, view it on GitHub <#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 ***@***.***> 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 > < #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: ***@***.*** — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#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.