[Thread Prev][Thread Next][Index]

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



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


[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement