[gradsusr] writing data for radius-height plot

sim.aberson sim.aberson at noaa.gov
Fri Dec 29 14:03:41 EST 2017


I have a set of data, azimuthal averages of 54 variables by height and 
radius, and I would like to plot them using GrADS.  I have 39 height 
levels and 135 radii.  I am writing the binary file with

       do l=1,54
         do j=1,39
           write(13)(rdat(l,j,i),i=1,135)
         end do
       end do

[I note that I also wrote the data to an ascii file, and it looks fine. 
The values are between -10 and 10.]

The control file I am using is:

  DSET ^fort.12
  TITLE Wavenumbers
  OPTIONS sequential
  UNDEF    9.99999956E+24
  XDEF 135 LINEAR 0 2
  YDEF 1 LINEAR 1 1
  ZDEF 39 LEVELS 0.035 0.1 0.2 0.4 0.6 0.8 1. 1.2
                 1.5 1.8 2.1 2.4 2.7 3. 3.4 3.8
                 4.2 4.6 5. 5.5 6. 6.5 7. 7.5 8.
                 8.5 9. 9.75 10.5 11.25 12. 13.
                 14. 15. 16. 17. 18. 19. 20.
  TDEF 1 LINEAR 1jan1950 1yr
  VARS 54
  uu1r  39  99 radial wind wavenumber 1 real
...
  ENDVARS

When I try to plot:

162_> grads

Grid Analysis and Display System (GrADS) Version 2.0.a7.1
Copyright (c) 1988-2008 by Brian Doty and the
Institute for Global Environment and Society (IGES)
GrADS comes with ABSOLUTELY NO WARRANTY
See file COPYRIGHT for more information

Config: v2.0.a7.1 little-endian readline printim grib2 netcdf hdf4-sds 
hdf5 geotiff
Issue 'q config' command for more information.
GX Package Initialization: Size = 11 8.5
ga-> open fourier.ctl
Scanning description file:  fourier.ctl
Data file fort.12 is open as file 1
LON set to 0 268
LAT set to 1 1
LEV set to 0.035 20
Time values set: 1950:1:1:0 1950:1:1:0
E set to 1 1
ga-> set lev 0.1 20
LEV set to 0.1 20
ga-> d uu1r
Contouring: 3e+29 to 3e+30 interval 3e+29
ga-> printim test.png png
ga-> quit
No hardcopy metafile open
GX package terminated

I tried all sorts of things (removed sequential, different variables, 
real versus integer values in control file, etc.) to no avail.  I'm on a 
mac running Sierra Version 10.12.6.  Any help greatly appreciated.

Thanks,
Sim
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.png
Type: image/png
Size: 10175 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20171229/2002c303/attachment.png 


More information about the gradsusr mailing list