Re: Bug in objective.jnl

[ansley@xxxxxxxxxxxxx (Ansley Manke)]

Hello Lev,

There's now a better way to do objective interpolation to a grid than the 
objective.jnl script that you're using.  Are you running Ferret version 4.91 or
5 ?  There's a newer objective analysis capability using Ferret's External 
Functions.

Take a look at objective_analysis_demo.jnl, also viewable as a demo from the
Ferret web page, under "On-line demonstrations", External functions examples:

http://ferret.wrc.noaa.gov/Ferret/Demos/objective_analysis_demo/objective_analysis_demo.html

Ansley Manke


Lev Tarasoff wrote:  Thu, 04 Nov 1999
> I've been trying to grid some xyz topo data for Scotland using the
> GO file objective.jnl, but the routine is shifting the data by what
> appears to be 3 degrees longitude to the east and what is worse, wrapping
> the data around the grid. How do I fix this (just shifting the input x
> coords does not halt the wrapping problem)?
> I'm running on a PC, linux RH 6.0, and I my script is as follows:
> 
> DEFINE AXIS/X=1:219:1 ia
> DEFINE AXIS/Y=0:122:0.5 ja
> DEFINE GRID/X=ia/Y=ja ija
> show grid ija
> FILE/SKIP=1/VAR="xpts ypts zpts"/grid=ija  scotmap4.dat
> GO objective xpts ypts zpts "352.75:359.95:0.05" "55.5:59.45:0.05" " " " "
> hdobj.unf
> 
> Thanks
> 
> Lev Tarasoff