[Thread Prev][Thread Next][Index]
[ferret_users] How to plot Monthly Global dataset of single variable?
- To: ferret_users@xxxxxxxx
- Subject: [ferret_users] How to plot Monthly Global dataset of single variable?
- From: Ashley Watson <ashleywats@xxxxxxxxx>
- Date: Thu, 06 Jul 2006 11:13:54 +0200
- Domainkey-signature: a=rsa-sha1; q=dns; c=nofws; s=beta; d=gmail.com;h=received:message-id:date:from:to:subject:mime-version:content-type;b=CyIIA5x8BWy57OYdc0ulRwoITCwCeAdlG6YEPP9UqbExM8++n1oqLMLion/UShjFg74M7Z5GLoLs+yXhyo75aSojmiG7rspTipBPUlsm2jNS6Qka9gFMX6ds6yDJ+kebQWyehbmMa+ONORQv4KpKsthml9b5ph6E/bTK5G9Xk/I=
- Sender: owner-ferret_users@xxxxxxxxxxxxx
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