[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