[Thread Prev][Thread Next][Index]

[ferret_users] How to plot Monthly Global dataset of single variable?



Hi all,
      I have monthly global dataset  of  255 months  from  198112 to 200302 . ( file size 66096000 bytes ). I read the data with the following commands.
I think there is no problem in reading the data.

Set memory/size=80
DEFINE AXIS/X=1:`360*180`:1 binary_x
DEFINE AXIS/T=1:255:1 binary_t
DEFINE GRID/X=binary_x/T=binary_t binary_g
FILE/FORMAT=stream/G=binary_g/VAR=sst OISST_198112_200302

Then I worked to define the grid for data according to section 4.1.3 of Chapter 2.

DEFINE AXIS/MODULO/X=0.5:359.5:1 onedeg_x
DEFINE AXIS/Y=-89.5:89.5:1 onedeg_y
DEFINE AXIS/T=15-Dec-1981:15-Feb-2003:1/UNITS=month month_t
DEFINE GRID/X=onedeg_x/Y=onedeg_y/T=month_t onedeg_g
LET dummy = x[GX=onedeg_g] + y[GY=onedeg_g] + t[GT=onedeg_g]
LET field = RESHAPE(sst[i=1:`360*180`],dummy)
I got following message while defining time axis,

*** NOTE: /UNIT=MONTHS is ambiguous ... using 1/12 of 365.2425 days

I tried "shade sst[l=1] " but got the error
**ERROR: dimensions improperly specified: must be a 2D region

I dont understand how to view my data with shade/fill  commnad after these steps.
Looking forward for help.

Ashley




[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement