!=======================================
use err672_samplexy_curv.nc
shade lon

! note the discontinuity
pause

let tabi=lon*0+i

!=======================================
file/var=t1,t2 15S.txt
shade/modulo/hlim=-180:180 tabi[d=1], lon[d=1], lat[d=1]
plot/vs/over/nolab t1,t2

!=======================================
list SAMPLEXY_CURV_NRST(tabi[d=1],lon[d=1],lat[d=1],t1,t2)

! all indices found are 53 --> WRONG
pause

!=======================================
let lon2=if i ge 53 then lon else lon-360
set var/units="degrees_east" lon2
shade/modulo/hlim=-180:180 tabi[d=1], lon2[d=1], lat[d=1]
plot/vs/over/nolab t1,t2

! sampling is CORRECT
list SAMPLEXY_CURV_NRST(tabi[d=1],lon2[d=1],lat[d=1],t1,t2)