Skip to content

interp_pts2profile with input time #215

@jmoralesFusion

Description

@jmoralesFusion

No correct interpolation when passing time array in input.

To reproduce

import tofu as tf
shot = 54719
multi = tf.imas2tofu.MultiIDSLoader(shot=shot, \
          ids=['wall', 'pulse_schedule', 'soft_x_rays', 'equilibrium'])
sxr_data = multi.to_Data('soft_x_rays', plot=False)
sxr_cam = sxr_data.lCam[0]
sxr_cam.set_dsino([2.4, 0])
sxr_pts = sxr_cam.dsino['pts']
sxr_plasma = multi.to_Plasma2D()
rhopn_all = sxr_plasma.interp_pts2profile(sxr_pts, quant='equilibrium.2drhopn')
rhopn = sxr_plasma.interp_pts2profile(sxr_pts, quant='equilibrium.2drhopn', t=np.r_[rhopn_all[1][10]])

rhopn_all differs from rhopn

Bug seems to be located in:

tofu/data/_comp.py line 431

Code :

                        ind = indt == indtu[ii]
                        val[ind,...] = mplTriLinInterp(mpltri,                        
                                                       vquant[indtq[ii],:],           
                                                       trifinder=trifind)(r,z).filled(fill_value) 

Metadata

Metadata

Assignees

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions