[ferret_users] Script crashed -- scat2gridgauss_xy

[Szymon Roziewski <szymon.roziewski@xxxxxxxxx>]

Dear Ferreters,
I have such a script:
set win/size=6.0/aspect=`1.1`
SET VIEWPORT FULL
CANCEL MODE LOGO
DEFINE AXIS/X=-5.666666667:5.25:0.083333333/unit=degree xlon
DEFINE AXIS/Y=-2.0:10.0:0.083333333/unit=degree ylat
DEFINE GRID/x=xlon/y=ylat grid_xlon_xlat
FILE/VARIABLES=swh_grid/COLUMNS=132/GRID=grid_xlon_xlat swh20101018090000.txt
set var/bad=-1 swh_grid
set var/units=meters swh_grid
set var/title="Significant Wave Height" swh_grid
DEFINE AXIS/X=-5.666666667:5.25:0.010416667/unit=degree xfine
DEFINE AXIS/Y=-2.0:10.0:0.010416667/unit=degree yfine
DEFINE GRID/x=xfine/y=yfine fine_grid
let xpts = xsequence(0*swh_grid+x[g=swh_grid])
let ypts = xsequence(0*swh_grid+y[g=swh_grid])
let datapts = xsequence(swh_grid)
let xscale=0.1
let yscale=0.1
let cutoff=2
let swh_smooth = scat2gridgauss_xy(xpts,ypts,datapts,x[g=fine_grid],y[g=fine_grid],xscale,yscale,cutoff,0)
list swh_smooth

wchich should do interpolation data, but it doesn't. The script follows up to the line with list swh_smooth and as reached as is iddle.
I think the line with scat2gridgauss_xy function failed. But I don't know the reason.
Why is that I don't know. I have very similar script which copes with cdf data and do it well.
But in this case previous approach failed. (When I try using txt data with flag bad=-1.)
I attach the data file.


Kind regards,
Szymon Roziewski

Attachment: swh20101018090000.tar.bz
Description: application/bzip