Include Weimer 2005 into dynamo solver

[twfang]

The electrodynamics in IPE currently includes the Heelis high latitude empirical model. We need to get the Weimer 2005 into the dynamo solver as well.

[twfang]

A new branch called ipe_refactor_new_weimer has been created for the development. This branch includes weimer 2005 in WAM, solar density that required for weimer 2005, and the latest IPE_refactor code. Solar wind density has been brought in to IPE_refactor in this branch as well.

[naomimaruyama]

weimer 2005 has now been implemented in the dynamo solver. In other words, weimer2005 replaces the heelis potential model.

here is a preliminary result. (top) input to the solver: (middle) output from the solver (bottom) EXB drift

solar wind condition: angle=180 bt=5 swvel=500 swden=0.01

Further fine tune is required to get it to work properly.

[twfang]

@IonospherePlasmasphereElectrodynamics Do you have any updates on this implementation, either for empirical or for the self-consistent package? Are you still using ipe_refactor_devwam_weimer05? If you have any codes that you want to commit to the branch, please do so. Otherwise, I would suggest you open a new branch from our latest branch called development_ipe_refactor_merage, which will make the merge simpler later.

[naomimaruyama]

I am testing the new code with the new weimer 2005 with the branch ipe_refactor_devwam_weimer05. I can easily transfer the code modification to a new branch from development_ipe_refactor_merage.

[twfang]

Thanks! Then, will it be ok to delete the ipe_refactor_devwam_weimer05 branch? I will then create a branch from development_ipe_refactor_merge for you to work on?

[naomimaruyama]

sure! thanks!

On Thu, Sep 19, 2019 at 1:00 PM twfang notifications@github.com wrote:

Thanks! Then, will it be ok to delete the ipe_refactor_devwam_weimer05 branch? I will then create a branch from development_ipe_refactor_merge for you to work on?

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/SWPC-IPE/WAM-IPE/issues/309?email_source=notifications&email_token=ADIZZYX7UT3GGFTPE3RN3ALQKPD4HA5CNFSM4IFZ7JPKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOD7EPUQQ#issuecomment-533264962, or mute the thread https://github.com/notifications/unsubscribe-auth/ADIZZYTBOQKEB5HJK47XV43QKPD4HANCNFSM4IFZ7JPA .

--

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]

I have created a new branch called development_weimer05 out of development_ipe_refactor_merge for Weimer05 implementation in IPE.

[naomimaruyama]

Thank you very much! I am going to move the changes over there.

On Mon, Sep 23, 2019 at 1:45 PM twfang notifications@github.com wrote:

I have created a new branch called development_weimer05 out of development_ipe_refactor_merge for Weimer05 implementation in IPE.

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/SWPC-IPE/WAM-IPE/issues/309?email_source=notifications&email_token=ADIZZYUY2WCIWNJAMHDMFB3QLEMFNA5CNFSM4IFZ7JPKYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOD7MBHSQ#issuecomment-534254538, or mute the thread https://github.com/notifications/unsubscribe-auth/ADIZZYXXNDRHF5ZT66GDDOLQLEMFNANCNFSM4IFZ7JPA .

--

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]

@IonospherePlasmasphereElectrodynamics Could you please check if you get potentials in the southern hemisphere from your outputs? Thank you very much!

[twfang]

@IonospherePlasmasphereElectrodynamics Naomi, when you get a chance, could you please plot your results to see if you get potentials in the southern hemisphere? Thank you!

[twfang]

Weimer has been included in the dynamo solver. However, there are still some problems in the latest package that I am currently debugging. Below shows the potentials from the two different dynamo packages.

Heelis package

Weimer package

[twfang]

While exploring the problem, I did a test run which set high-latitude potential to be 0. Once the high-latitude potential is turned off. The model with self-consistent electrodynamics became unstable. It looks like the field-line integrals are failing at some locations, which leads to bad electric potential and electric fields.

Above figure shows the field-line integrals 1 (top) and 3 (bottom) from the 1st to 4th time steps (left to right). Bad field-line integrals start to show up at 2nd time-step near lp=138 (mlat~ 7.1 degree) and become worse after that.

The problem then shows up in potentials (top), Ed1 (middle), and Ed2 (bottom).

Line plots show the problem tubes.

The bad field line integrals are related to the NaN values in those tubes. I will merge George's changes to see if the problem can be solved.

[twfang]

The weimer_ipe_fixes branch has been updated - George is still testing/looking at results but this branch has Weimer + electrodynamics and appears to be working smoothly (only tried standalone so far)

The latest NaNs traced to a problem with the transport (in the lp direction). There was an assumption built-in that convection would always be in the lp-1, lp, lp+1 range

• which is not always true - and that was throwing some of the inner (low-latitude) flux-tubes out of whack!

George has run the 3 days - 16th 17th 18th with no problems (well, at least none that I've found)