Re: Problem using Objective Analysis

[Ansley Manke <ansley@xxxxxxxxxxxxx>]

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