tsssss/geopack

Difference in output of t96 and t01 model between IDL and Python geopacks!

Closed this issue · 4 comments

Hello,
I am testing and comparing the outputs of this geopack for two different models with those of IDL and I am getting very different results.

I have attached the figures which I plotted for this purpose. For each panel in the figure I have plotted the difference between the output from IDL and Python (b_out_idl - b_out_py). Difference from T96 are plotted in the first row whereas the difference between T01 model are plotted in the second row. The value of x_gsm is at the top left position whereas the value of DST used for the code is at top right

For both T96 and T01 model far from the earth, into the tail, the difference between two outputs is minimal, almost always less than 2nT, and thus can be ignored in most cases (Figure 1, for x_gsm = 9.98 R_E).

However, there seem to be significant difference in output when we move close to the Earth or when we are on the dayside. Though the difference is smaller for T01 outputs compared to those of T96, they are still significant, specially far from the Earth on the dayside (Figure 2).

I wonder if anyone else have had similar issues.

The parameters I used for the codes are as follows:

time = 2015-01-01 00:00:00 GMT
par = [5, 0, 1, 1, 0, 0, 0, 0, 0, 0]
x_gsm = np.linspace(-15.1, 15, 61)
y_gsm = np.linspace(-15.1, 15, 61)
z_gsm = np.linspace(-15.1, 15, 61)

I can provide the full python and IDL codes if that will help with the reproduction of these differences.

Figure 1
Figure_3

Figure 2
Figure_2

Figure 3
Figure_1

Hey Sheng,
Yes I was able to reproduce the example which are in Readme file. They are exactly same.

I have uploaded the codes in my repository and have linked them here:

idl_code
python_code
idl_data_edit
plotting_code

I have tried to make the code consistent so that you can run them on your system with minimal edits. Also, please note that the codes are listed in the order they are run!

Hey Sheng,
Sorry for the delay in response. But yes it worked after I updated it with the latest version which is here on GitHub. Thank you so much for your help in debugging this!