[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