I'm playing around with some data that's on a cubed-sphere grid,
and am wondering if the description of the grid is causing the
usual "Y coord field lies in different plane from data to be
plotted" error when using the "shade ts,lon,lat" command.
The data of interest and the longitudes and latitudes are all
described as one-dimensional vectors, i.e., a 1-D array of N
points:
yes? use test.nc
yes? show data
currently SET data sets:
1> ./test.nc (default)
name title I J
K L M N
TS Surface temperature (radiative) 1:777602 ...
... 1:1 ... ...
LAT latitude 1:777602 ...
... ... ... ...
LON longitude 1:777602 ...
... ... ... …
and 'ncdump -h' shows:
prompt> ncdump -h test.nc | egrep
"(lon|lat|TS)"
float TS(time, ncol) ;
TS:units = "K" ;
TS:long_name =
"Surface temperature (radiative)" ;
TS:cell_methods =
"time: mean" ;
double lat(ncol) ;
lat:long_name =
"latitude" ;
lat:units =
"degrees_north" ;
lat:axis = "Y" ;
double lon(ncol) ;
lon:long_name =
"longitude" ;
lon:units =
"degrees_east" ;
lon:axis = "X" ;
Giving lon and lat the appropriate "axis" attributes doesn't
change anything about how FERRET sees the data.
Any suggestions appreciated!