[Thread Prev][Thread Next][Index]
Re: scat2gridgauss input error?
Hello Ken,
If the file test.dat starts with the header line "time pres temp"
then you need to add /SKIP=1 to the SET DATA/EZ command.
Then the data can be read in properly.
You also might need to experiment with the gridding parameters, and
perhaps also try SCAT2GRIDLAPLACE_XY which uses a different
gridding algorithm. Your data points are quite spread out in the region
so you will likely need to make the scale and cutoff parameters larger.
The command:
yes? plot/vs time, pres
will show you where the scattered points are on the plot
Ansley Manke
"Kwok,Ken [PYR]" wrote:
> Hello Ferreters,
>
> I have a problem reading in an ascii file that looks something like this:
>
> time pres temp
> 17.000 980 14
> 17.200 973 13
> 17.300 966 10
> 17.500 951 9
> 17.660 942 8
> 17.830 928 5
> 18.000 913 4
> 18.100 900 4
> 18.450 895 3
> 18.560 900 4
> 18.750 930 6
> 18.890 955 7
> 19.000 968 8
>
> I wanted to plot the data with time on x-axis, and a contour of temp with
> height (height on y-axis). As advised, I am using the function
> SCAT2GRIDGAUSS since the time and pres data are irregular:
>
> yes? DEFINE AXIS/X=17:20:0.1 xax
> yes? DEFINE AXIS/Y=980:880:20 yax
> yes? SET DATA/EZ/VARIABLES="time,pres,temp" test.dat
> yes? LET out = SCAT2GRIDGAUSS_XY(time,pres,temp,x[gx=xax],y[gy=yax],2,2,2,2)
> yes? CONTOUR out
> : Success
> **TMAP ERR: Last or next-to-last record read:
> 17.000 980 14
> Data set: /usr/local/ferret/examples/test.dat
> Data file: /usr/local/ferret/examples/test.dat
>
> I keep getting the same error at the end and no plot is created. Can anyone
> please help me?
>
> Thanks,
> Ken
> ---------------------------------
> Kenneth Kwok
> Meteorological Services of Canada
> Pacific and Yukon Region
> Telephone: 604.664.9179
--
Ansley Manke Pacific Marine Environmental Laboratory Seattle WA (206)526-6246
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement