Hi Christelle,
I don't see the variable "SST" in your netCDF CDL below, but
assuming that it has just been truncated, the issue appears to be
that the dataset in question has curvilinear coordinates. What you
are seeing on the X and Y axes of your plot is likely the index
values (if this is right the axis tics will not be formatted in
degrees. GO land is trying to overlay lat/long positions on a
coordinate system of the index ranges.)
If so you need to be using the curvilinear (3-argument) syntax of
the SHADE command
yes? SHADE/i=1:224/j=1:224/l=1 SST, XLON, XLAT;
GO land
See
http://ferret.pmel.noaa.gov/Ferret/documentation/users-guide/customizing-plots/MAP-PROJECTIONS-AND-CURVILINEAR-COORDINATES#_VPINDEXENTRY_875
- Steve
On 2/8/2013 12:30 AM, CASTET Christelle
wrote:
Hello I
am a fairly new user of FERRET and trying to plot data that are in
the netcdf format. My issue is when I want to map the land
contour with GO LAND, the land appears all squeezed in the bottom
of my plot:
SHADE/i=1:224/j=1:224/l=1 SST;GO land
or SHADE/x=1:224/y=1:224/l=1 SST;GO land
I can't figure out if there is something wrong with my commands or
the way my netcdf file is written, or if I have to use a
particular map projection (I have looked at some demo scripts
without finding any answer to my problem).
Thank you in advance for any help you can give me.
Christelle Castet
*Here is the output of show data :
* currently SET data sets:
1> ./wrfout_d01_2011-03-12_00:00:00.nc (default)
name title I J
K L
TIMES ... ...
... 1:24
LU_INDEX 1:224 1:224
... 1:24
ZNU 1:39 ...
... 1:24
ZNW 1:40 ...
... 1:24
ZS 1:4 ...
... 1:24
DZS 1:4 ...
... 1:24
VAR_SSO 1:224 1:224
... 1:24
LAP_HGT 1:224 1:224
... 1:24
U 1:225 1:224
1:39 1:24
V 1:224 1:225
1:39 1:24
W 1:224 1:224
1:40 1:24
*And from ncdump -h :
*netcdf wrfout_d01_2011-03-12_00\:00\:00 {
dimensions:
Time = UNLIMITED ; // (24 currently)
DateStrLen = 19 ;
west_east = 224 ;
south_north = 224 ;
bottom_top = 39 ;
bottom_top_stag = 40 ;
soil_layers_stag = 4 ;
west_east_stag = 225 ;
south_north_stag = 225 ;
variables:
char Times(Time, DateStrLen) ;
float LU_INDEX(Time, south_north, west_east) ;
LU_INDEX:FieldType = 104 ;
LU_INDEX:MemoryOrder = "XY " ;
LU_INDEX:description = "LAND USE CATEGORY" ;
LU_INDEX:units = "" ;
LU_INDEX:stagger = "" ;
LU_INDEX:coordinates = "XLONG XLAT" ;
float XLAT(Time, south_north, west_east) ;
XLAT:FieldType = 104 ;
XLAT:MemoryOrder = "XY " ;
XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;
XLAT:units = "degree_north" ;
XLAT:stagger = "" ;
float XLONG(Time, south_north, west_east) ;
XLONG:FieldType = 104 ;
XLONG:MemoryOrder = "XY " ;
XLONG:description = "LONGITUDE, WEST IS NEGATIVE" ;
XLONG:units = "degree_east" ;
XLONG:stagger = "" ;
float XLAT_U(Time, south_north, west_east_stag) ;
XLAT_U:FieldType = 104 ;
XLAT_U:MemoryOrder = "XY " ;
XLAT_U:description = "LATITUDE, SOUTH IS NEGATIVE" ;
XLAT_U:units = "degree_north" ;
XLAT_U:stagger = "X" ;
XLAT_U:coordinates = "XLONG_U XLAT_U" ;
float XLONG_U(Time, south_north, west_east_stag) ;
XLONG_U:FieldType = 104 ;
XLONG_U:MemoryOrder = "XY " ;
XLONG_U:description = "LONGITUDE, WEST IS NEGATIVE"
;
XLONG_U:units = "degree_east" ;
XLONG_U:stagger = "X" ;
XLONG_U:coordinates = "XLONG_U XLAT_U" ;
float XLAT_V(Time, south_north_stag, west_east) ;
XLAT_V:FieldType = 104 ;
XLAT_V:MemoryOrder = "XY " ;
XLAT_V:description = "LATITUDE, SOUTH IS NEGATIVE" ;
XLAT_V:units = "degree_north" ;
XLAT_V:stagger = "Y" ;
XLAT_V:coordinates = "XLONG_V XLAT_V" ;
float XLONG_V(Time, south_north_stag, west_east) ;
XLONG_V:FieldType = 104 ;
XLONG_V:MemoryOrder = "XY " ;
XLONG_V:description = "LONGITUDE, WEST IS NEGATIVE"
;
XLONG_V:units = "degree_east" ;
XLONG_V:stagger = "Y" ;
XLONG_V:coordinates = "XLONG_V XLAT_V" ;
|