Document Actions
RECT_TO_CURV
RECT_TO_CURV (V, lon_bounds_out, lat_bounds_out, missing_allowed) Uses Bilinear Interpolation to regrid data on a recilinear grid to a curvilinear grid.
|
Arguments: |
V |
Variable on the source rectilinear grid |
|
lon_bounds_out |
Destination grid longitudes, in degrees (2-D field of longitudes describing the curvilinear grid) |
|
|
lat_bounds_out |
Destination grid latitudes, in degrees (2-D field of latitudes describing the curvilinear grid) |
|
|
missing_allowed |
Number of missing values allowed among the 4 surrounding source cells: 0 to 3 |
|
|
Result Axes: |
X |
Inherited from lon_bounds_out |
|
Y |
Inherited from lat_bounds_out |
|
|
Z |
Inherited from V |
|
|
T |
Inherited from V |
Example: put a rectilinear data set on a curvilinear grid for comparison.
yes? use my_curvilinear_data.nc
yes? use levitus_climatology
yes? show data
currently SET data sets:
1> /.my_curvilinear_data.nc
name title I J K L
GEOLON geographic longitude of grid 1:180 1:173 ... ...
GEOLAT geographic latitude of grid 1:180 1:173 ... ...
SALT salinity 1:180 1:173 1:30 1:12
AIRT air temperature 1:180 1:173 1:30 1:12
2> /home/ja9/tmap/fer_dsets/data/levitus_climatology.cdf (default)
name title I J K L
TEMP TEMPERATURE 1:360 1:180 1:20 ...
SALT SALINITY 1:360 1:180 1:20 ...
yes? ! For convenience, define variables for arguments to RECT_TO_CURV
yes? let lonout = geolon[d=1]
yes? let latout = geolat[d=1]
yes? let a = rect_to_curv(temp[d=2,k=1], lonout, latout, 2)
yes? shade a, lonout, latout
yes? ! Compare to a variable on the curvilinear grid
yes? shade a-airt[d=1], lonout, latout
yes? ! save the new variable for all depths, with the curvilinear
yes? ! latitude and longitude variables
yes? let temp_curv = rect_to_curv(temp[d=2], lonout, latout, 2)
yes? save/file=tcurv.nc temp_curv, lonout, latout
LISTing to file tcurv.nc