RE: Problem using Objective Analysis

["Donald S Dunbar" <dsd@xxxxxxx>]

Ok, that seems to work.  I think the documentation is less than clear on
this point, however.

I'm actually using scatter2gridlaplace_xy and it works well--to a point.
There seems to be a limit on the number of scattered data points that can be
used; somewhere around 350.  Above this limit it still executes, but the
results are wrong.  Any further info on this?

-----Original Message-----
From: owner-ferret_users@ferret.wrc.noaa.gov
[mailto:owner-ferret_users@ferret.wrc.noaa.gov]On Behalf Of Ansley Manke
Sent: May 9, 2000 12:35 PM
To: Donald S Dunbar
Cc: ferret_users@ferret.wrc.noaa.gov
Subject: Re: Problem using Objective Analysis


Hi Donald,
I think I have the answer.  You need to get the numbers from your XAX  and
YAX axes into X and Y variables; so instead of:

     let h_grd=scatter2gridgauss_xy(lon,lat,h,xax,yax,2.)

the syntax is

     let h_grd=scatter2gridgauss_xy(lon,lat,h,x[gx=xax],y[gy=yax],2.)

Give it a try!

Ansley Manke


Donald S Dunbar wrote:

> I'm trying to use the objective analysis routines in ferret vs. 5.  Here
is
> my script:
>
> ========== 8< plot.jnl ============================
> !
> set window/aspect=`213/185`
> define axis/y=12.814:13.195:0.001799 YAX
> define axis/x=121.447:121.786:0.001847 XAX
> define grid/x=xax/y=yax g
> file/grid=g/var=lon,lat,h "./ferret.dat"
> set var/bad=0 h
> let h_grd=scatter2gridgauss_xy(lon,lat,h,xax,yax,2.)
> shade/lev=(0)(10)(25)(50)(100,600,50) h_grd,lon,lat
> contour/ov/lev=(0)(10)(25)(50)(100,600,50) h,lon,lat
> frame/file=plot.gif
> ========== 8< plot.jnl ============================
>
> When I try "GO PLOT" I get the following:
>
> yes? go plot
> !
> set window/aspect=`213/185`
>  !-> set window/aspect=1.1514
> define axis/y=12.814:13.195:0.001799 YAX
> define axis/x=121.447:121.786:0.001847 XAX
> define grid/x=xax/y=yax g
> file/grid=g/var=lon,lat,h "./ferret.dat"
> set var/bad=0 h
> let h_grd=scatter2gridgauss_xy(lon,lat,h,xax,yax,2.)
> shade/lev=(0)(10)(25)(50)(100,600,50) h_grd,lon,lat
>  **ERROR: variable unknown or not in data set: YAX
> shade/lev=(0)(10)(25)(50)(100,600,50) h_grd,lon,lat
> Command file, command group, or REPEAT execution aborted
> yes?
>
> Why doesn't this work????
>
> Don Dunbar
>
> PS
>
> Is there PDF formatted vs. 5 manual available?

--
Ansley Manke  Pacific Marine Environmental Laboratory  Seattle WA
(206)526-6246