grads question

[Lemmy]

The record length of R, Z (1D) and T_grads (2D) are not consistent


JEFFREY S GALL <jsg229 at PSU.EDU> wrote:
  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







---------------------------------
 Now you can scan emails quickly with a reading pane. Get the new Yahoo! Mail.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070206/e16ff26e/attachment.html