Hi Judy,
Yes, 2 would be a reasonable value. A larger radius will give you
a smoother result, a smaller one would in general capture
variations in the original data better, but at the cost of perhaps
some cells not being filled in.
-Ansley
On 8/11/2016 6:52 AM, judy wrote:
Russ, got it! Thanks. Also just to
confirm... In the example I see this line:
yes? let map =
curv_to_rect_map ( lonin,latin,lonlatout,10)
The CURV_TO_RECT docs define a radius argument, which "should
be slightly larger than the size of the input curvilinear grid
cells so that any output grid location will have some input data
contributing to its value." The example output grid is 5x5, so
the '10' is 10-degree radius? In my case I have an output grid
of 1x1 so I used '2' as a 2-degree radius as a reasonable start
point. Do I have this right?
On 8/10/16 9:35 PM, Russ Fiedler wrote:
Hi
Judy,
yes? define axis/y=0:90:1/units=degrees yax
This only defines the Northern hemisphere (0N to 90N).
I'll leave the global definition as an exercise for the
reader...
Russ
On 11/08/16 10:57, judy wrote:
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
|