Re: [ferret_users] Problem regridding GFDL ECDA v3.1 ocean analysis to 1x1

[judy <jhz8181@xxxxxxxxx>]

Thanks Russ -- I wonder how long it would have taken for me to spot 
that? I feel like a kid at Christmas, I've got the top level of salinity 
and the northern hemisphere in rectilinear!! Now I need to get the other 
49 levels and find out where the southern hemisphere went to...

thanks again!
judy

On 8/10/16 8:02 PM, Russ Fiedler wrote:
> Hi Judy,
>
> You're not seeing the wood for the trees! You've misspelled "let" with 
> "lat" in your last command. It looks as if you are on the right track 
> with what you are doing.
>
>
>
> Russ
>
> On 11/08/16 05:45, judy wrote:
> > Hi Ferreters,
> >
> > I'm very new to curvilinear grids (and ferret!). I'm using an example 
> > I found in the users guide -> appendix a:functions -> CURV_TO_RECT  
> > to guide me in remapping GFDL to 1x1. Can anyone help with my error? 
> > Or better yet, is there a script someone might share from a previous 
> > successful effort?
> >
> > yes? use gfdl-salt-1980.nc
> > yes? show data
> >      currently SET data sets:
> >     1> ./gfdl-salt-1980.nc  (default)
> >  name     title                             I         J K L         
> > M         N
> >  GEOLON_T tracer longitude                 1:360     1:200 ...       
> > ...       ...       ...
> >  GEOLAT_T tracer latitude                  1:360     1:200 ...       
> > ...       ...       ...
> >  SALT     Salinity                         1:360     1:200 1:50      
> > 1:12      ...       ...
> >  AVERAGE_DT
> >           Length of average period         ...       ... ... 
> > 1:12      ...       ...
> >
> > yes? let lonin = geolon_t[d=1]
> > yes? let latin = geolat_t[d=1]
> > yes? define axis/x=0:360:1/units=degrees xax
> > yes? define axis/y=0:90:1/units=degrees yax
> > yes? let lonlatout = y[gy=yax] + x[gx=xax]
> > yes? lat map = curv_to_rect_map ( lonin,latin,lonlatout,2)
> >  **ERROR: unknown command: lat
> >
> > thanks!
> >
> > judy