newbie: trying to print ann average vertical profile numbers at specific lat, lon for monthly gridded binary file
Steve Piper
scpiper at UCSD.EDU
Fri Sep 26 17:43:23 EDT 2008
Hi,
I just started using Grads. Great program!
I would like to write out in ascii, the annual average vertical
profile at
all 9 levels at a specified location (lat and lon) in
the grid for variable bioa.
I am reading a gridded data file that was written with Fortran
statement declaration
real x(72,44,9,12,26) ! lon, lat, levels, months, variables
My ctl file is:
DSET /Users/piper/Grads-test/CSU.gurney.L1.mmean.le.3D
* fortran sequential file c(72,44,9,12,26)
* lon lat lev month tracer
* sequential means fortran binary with record length integers
in record
* otherwise would use stream or direct access
* SUN is big-endian (had to swap to make it so with zz.f)
* OPTIONS big_endian for SUN, mac powerpc?
OPTIONS sequential
* input file is from little_endian machine
OPTIONS little_endian
TITLE CSU Gurney 3D monthly output for T3L1
UNDEF 1.e36
XDEF 72 linear -180 5
YDEF 44 linear -86 4
ZDEF 9 levels 1000 925 850 700 500 400 300 200 100
* would like to use middle of every 365 calendar month - this is
close
TDEF 12 linear 16jan94 1mo
VARS 26
fos90 9 -1,20 fos90
* -1 = non-standard binary, 20 = var and time reordered
fos95 9 -1,20 fos95
bioa 9 -1,20 bioa
dum 9 -1,20 dum
bio1 9 -1,20 bio1
bio2 9 -1,20 bio2
bio3 9 -1,20 bio3
bio4 9 -1,20 bio4
bio5 9 -1,20 bio5
bio6 9 -1,20 bio6
bio7 9 -1,20 bio7
bio8 9 -1,20 bio8
bio9 9 -1,20 bio9
bio10 9 -1,20 bio10
bio11 9 -1,20 bio11
oce1 9 -1,20 oce1
oce2 9 -1,20 oce2
oce3 9 -1,20 oce3
oce4 9 -1,20 oce4
oce5 9 -1,20 oce5
oce6 9 -1,20 oce6
oce7 9 -1,20 oce7
oce8 9 -1,20 oce8
oce9 9 -1,20 oce9
oce10 9 -1,20 oce10
oce11 9 -1,20 oce11
ENDVARS
set lon -180 180 # default is 0 360 when I issue q dim
set lat -86 86 # this was the default
set lev 1000 100 # default was 1000
set t 1 12 # default was t 1
q dim
Default file number is: 1
X is varying Lon = -180 to 180 X = 1 to 73
Y is varying Lat = -86 to 86 Y = 1 to 44
Z is varying Lev = 1000 to 100 Z = 1 to 9
T is varying Time = 00Z16JAN1994 to 00Z16DEC1994 T = 1 to 12
E is fixed Ens = 1 E = 1
Not sure how to produce annual average array newbioa(72,44,9) from
bioa(72,44,9,12)
Tried this:
define newbioa = ave(bioa,t=1,t=12)
It gives 108 lines (=9 x 12) of output stating
Averaging. dim = 3, start = 1, end =12
ga-> q dim
Default file number is: 1
X is varying Lon = -180 to 180 X = 1 to 73
Y is varying Lat = -86 to 86 Y = 1 to 44
Z is varying Lev = 1000 to 100 Z = 1 to 9
T is varying Time = 00Z16JAN1994 to 00Z16DEC1994 T = 1 to 12
E is fixed Ens = 1 E = 1
d gr2stn(newbioa,60.75,89.38)
Operation error: Invalid dimension environment
Too many varying dimensions
set t 1
ga-> q dim
Default file number is: 1
X is varying Lon = -180 to 180 X = 1 to 73
Y is varying Lat = -86 to 86 Y = 1 to 44
Z is varying Lev = 1000 to 100 Z = 1 to 9
T is fixed Time = 00Z16JAN1994 T = 1
E is fixed Ens = 1 E = 1
d gr2stn(newbioa,60.75,89.38)
Operation error: Invalid dimension environment
Looping dimension does not vary
Not sure why the gr2stn statements do not work. Help?
Am I on the right track?
Thanks!
Steve
More information about the gradsusr
mailing list