search

dstansby / pfsspy

Potential Field Source Surface model package for Python
https://pfsspy.readthedocs.io/
Other
40 stars 17 forks source link

The value of pfsspy.Output.br (or bc) is NaN #361

Closed rbluosolar closed 1 year ago

rbluosolar commented 1 year ago

I’m trying to obatin pfss field data with higher spatial resolution from the HMI synoptic map. So I make the larger size of pfsspy.input data, but I obtain the 'Nan' value of pfsspy.output.br (or bc). The figure I provided is the code about this question.

How can I solve this problem? Thank you!

image
dstansby commented 1 year ago

Have you double checked that hmi_map_pfss doesn't have any NaNs in the second case? Even though in theory your code should be getting rid of them, something might be going wrong in the nan_to_num or resample stage.

rbluosolar commented 1 year ago

Yes, I check agin. There are no NaN in hmi_map_pfss (Fig 1). It seams that there are 'NaN' when the size of pfss_in is larger than some sizes (Fig2, same pfss_in size, but different numbers of cells nrho in radial direction ).

Fig 1:

image

Fig 2:

image
dstansby commented 1 year ago

Thanks for opening this issue. I'm about to archive pfsspy as a project, so I'm afraid I won't be able to help you with it. I would recommend asking in https://community.openastronomy.org/c/sunpy/5 instead where you might be able to find some help.