inconsistent comparison results

[shihikoo]

Hi

I'm testing this python geopack package with the IDL routine.

I first tested with the parameters used in the test_geopack1.md as following: Time/dipole tilt: 2001-01-01/02:03:04 UT/ps = -0.533585131 rad Position: (xgsm,ygsm,zgsm) = (-5.1,0.3,2.8) Model parameter: T89: iopt = 2 T96, T01, and T04: par = [2., -87., 2., -5., 0., 0., ps, xgsm,ygsm,zgsm].

t89 | idl | 20.772 | -0.647 | -15.072   | python | 20.772 | -0.647 | -15.072 t96 | idl | 61.179 | -1.461 | -40.451   | python | 61.179 | -1.461 | -40.451 t01 | idl | 46.973 | 1.544 | -31.355   | python | 46.354 | 1.440 | -32.000 t04 | idl | -8.080 | -1.512 | -0.722   | python | 12.010 | -2.635 | -11.964

This result is similar to what is in test_geopack1.md except t04 for IDL routine.

I then tested for one of the MMS data, and the results are much different datetime | 2016-01-01/12:13:20 iopt / kp | 3 / 2 (xgsm,ygsm,zgsm) = | 6.96, 0.49, -0.66 parmod | 0.976, -50.20, 6.59, -5.512, 0., 0., tilt, x,y,z

t89 | idl | -3.657066785 | 1.085712758 | 13.86057416   | python | -2.37842015 | 1.038191716 | 14.38198039

t96 | idl | 22.97810632 | -1.106448782 | -4.147864164   | python | 50.00624944 | 6.397909797 | -147.0215389

t01 | idl | 5.001211425 | 0.02607965332 | 3.264103442   | python | 6.230733691 | 4.283414418 | -2.58591528

t04 | idl | -15.25291193 | -8.650237035 | -10.97633666   | python | 7.905332371 | -1.161844446 | -7.933264503

Could you help me with it?

Thanks

[tsssss]

Hello,

It looks like you got two questions here:

1. The IDL results for t04 are not the same as those in the example.
2. In the MMS example, IDL and python results do not compare in essentially all models.

Let's solve #1 first. Are you running the example with par = [2,-87,2,-5,0,0, tilt, x,y,z], where the tilt is in radian?

Thanks, Sheng

On Fri, 25 Jun 2021 at 11:45, Jing @.***> wrote:

Hi

I'm testing this python geopack package with the IDL routine.

I first tested with the parameters used in the test_geopack1.md as following: Time/dipole tilt: 2001-01-01/02:03:04 UT/ps = -0.533585131 rad Position: (xgsm,ygsm,zgsm) = (-5.1,0.3,2.8) Model parameter: T89: iopt = 2 T96, T01, and T04: par = [2., -87., 2., -5., 0., 0., ps, xgsm,ygsm,zgsm].

t89 | idl | 20.772 | -0.647 | -15.072 | python | 20.772 | -0.647 | -15.072 t96 | idl | 61.179 | -1.461 | -40.451 | python | 61.179 | -1.461 | -40.451 t01 | idl | 46.973 | 1.544 | -31.355 | python | 46.354 | 1.440 | -32.000 t04 | idl | -8.080 | -1.512 | -0.722 | python | 12.010 | -2.635 | -11.964

This result is similar to what is in test_geopack1.md except t04 for IDL routine.

I then tested for one of the MMS data, and the results are much different datetime | 2016-01-01/12:13:20 iopt / kp | 3 / 2 (xgsm,ygsm,zgsm) = | 6.96, 0.49, -0.66 parmod | 0.976, -50.20, 6.59, -5.512, 0., 0., tilt, x,y,z

t89 | idl | -3.657066785 | 1.085712758 | 13.86057416 | python | -2.37842015 | 1.038191716 | 14.38198039

t96 | idl | 22.97810632 | -1.106448782 | -4.147864164 | python | 50.00624944 | 6.397909797 | -147.0215389

t01 | idl | 5.001211425 | 0.02607965332 | 3.264103442 | python | 6.230733691 | 4.283414418 | -2.58591528

t04 | idl | -15.25291193 | -8.650237035 | -10.97633666 | python | 7.905332371 | -1.161844446 | -7.933264503

Could you help me with it?

Thanks

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/tsssss/geopack/issues/7, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADBMLBRBDN2CXTWVIPOBVYTTUSQDHANCNFSM47KDZOEQ .

[tsssss]

Hello,

I'll close this issue because I've tested that if tilt in radian is used, the IDL results for t04 are the same as those in the example. For Bx,y,z, I got 12.009596 -2.6345863 -11.964033

Please feel free to test your MMS example and contact me if there's still inconsistency.

Thanks, Sheng