Hi All,

If you are using GrADS, this info may help.

A student asked:
>
> I am trying to give grads a go with the 2d plots.  if
> I try to plot an x vs z plot, I get nothing.  BUT if I
> do an x vs y plot (changing the z values to y values), I
> think I get something similar to fovell's plots.  Its
> hard to debug without being able to do a x-z plot.  Is
> is possible to have only 2D data in grads (or do I need to
> pretend y is z?)  here is my code for writing out
> th(i,k) and pprt(i,k) and the control file to view it.
> Do you see anything that I am doing wrong here?

Comments below!

>
> writing out grads data
> ! open output file
>
> OPEN(10,FILE='output.gds',FORM='UNFORMATTED',ACCESS='DIRECT',
> RECL=4*NX)
>  irec=1

This is fine.

>
> ! Write output file
>  do k=2,nz-1
>    WRITE(10,REC=IREC) (th(i,k),i=2,nx-1)
>    irec = irec+1
>  enddo
>  do k=2,nz-1
>    WRITE(10,REC=IREC) (pprt(i,k),i=2,nx-1)
>    irec = irec+1
>  enddo

Take care here because you used NX in RECL when you opened the file, but you
are actually only writing from 2 to NX-1 (i.e. a total of nx-2 values).

>
> --------------------------------------------------
>
> CTL file:
> dset   ./output.gds
> title  results from simple model
> undef   -99999
> xdef   25 linear 0 .1
> ydef    1 linear 1  1
> zdef   25 linear 0 .1
> tdef   1 linear 00Z01JAN2000 5mn
> vars   2
> theta    1 99 arbitrary variable
> ppert    1 99 arbitrary variable
> endvars

This explains the plotting issue in GrADS.

The first number on the variable lines is the number of levels on which the
variable exists.  So, replace with:

vars   2
theta    25 99 arbitrary variable
ppert    25 99 arbitrary variable
endvars

That should do it!

All the best,
Matt

****************************************************************************