grads question

[JEFFREY S GALL]

I am attempting to create a .dat file and .ctl file for use in creating radius-height (r-z) plots in Grads.  I have 4 different times and the dimensions of the data are 32 points in the radial direction and 30 points in the vertical direction.  Suppose that I want to create a r-z plot of temperature.  This is what I have done:

DO T=1,4  !T=time counter


       WRITE(541,REC=C) (R(I,T),I=1,IR) !R=radius, IR=32
       WRITE(541,REC=C) (Z(K,T),K=1,KX) !Z=height, KX=30       WRITE(541,REC=C) ((T_GRADS(K,I,C),K=1,KX),I=1,IR) !T_GRADS = temperature END DO Now, my .ctl file looks like:DSET    0.3sf_grads.dat
byteswapped
TITLE        Sample
UNDEF           -9999
XDEF               30  linear  1  1
YDEF               32   linear 1  1
ZDEF                1   linear 1  1
TDEF             4  linear   12:00Z15SEP1999   180MN
vars              3
z      0 99  height
r      0 99  radius
t_az   0 99  TEMP AZ. AVG.
endvars
When I attempt to create my r-z plot, I get a field is constant, value=0 error.  What am I doing wrong here.  Any help/advice in fixing this problem would be most appreciated.  Thanks,
 Jeff





-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070206/6ab55c53/attachment.html