[gradsusr] I cannot create a gridded data file to be read by grads
Dan_fauba
dbarrera at agro.uba.ar
Mon Nov 17 11:40:54 EST 2014
HI,
I am a new grads user. I cannot create properly a gridded data file to
be read by grads.
The list of messages:
$ grads
grads: Symbol `_ZTVN6libdap15BaseTypeFactoryE' has different size in
shared object, consider re-linking
Grid Analysis and Display System (GrADS) Version 2.0.1
Copyright (c) 1988-2011 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.1 little-endian readline printim grib2 netcdf hdf4-sds hdf5
opendap-grids,stn geotiff shapefile
Issue 'q config' command for more detailed configuration information
Landscape mode? ('n' for portrait): q config n
GX Package Initialization: Size = 11 8.5
ga-> run racc
No hardcopy metafile open
All files closed; all defined objects released;
All GrADS attributes have been reinitialized
Most GrADS attributes have been reset
Cannot plot color bar: No shading information
Hardcopy output file is closed
ga-> open racc.ctl
Scanning description file: racc.ctl
Data file amsub_20140304_1838bin.dat is open as file 2
ga-> query file
File 1 : rainfall - ch4
Descriptor: racc.ctl
Binary: amsub_20140304_1838bin.dat
Type = Gridded
Xsize = 315 Ysize = 250 Zsize = 1 Tsize = 1 Esize = 1
Number of Variables = 1
rr 0 -9,40,2 ** rr (mm/h)
ga-> d rr
Constant field. Value = 0
ga-> quit
No hardcopy metafile open
GX package terminated
I created a binary gridded data file with fortran 90, as follows:
irecord= nfil1*ncol1 !250*315*4 = 78750 * 4 = 315000
open(10,file=s183_02_bin,access='direct',form='unformatted',recl=irecord,status='unknown')
write(10,rec=1) ((rg1(i,j),i=noes,nest),j=nsur,nnor)
close(10)
where rg1(i,j) is real*4.
Thanks in advance for your help.
Daniel
grads user member: dbarrera at agro.uba.ar
-------------- next part --------------
$ grads
grads: Symbol `_ZTVN6libdap15BaseTypeFactoryE' has different size in shared object, consider re-linking
Grid Analysis and Display System (GrADS) Version 2.0.1
Copyright (c) 1988-2011 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.1 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn geotiff shapefile
Issue 'q config' command for more detailed configuration information
Landscape mode? ('n' for portrait): q config n
GX Package Initialization: Size = 11 8.5
ga-> run racc
No hardcopy metafile open
All files closed; all defined objects released;
All GrADS attributes have been reinitialized
Most GrADS attributes have been reset
Cannot plot color bar: No shading information
Hardcopy output file is closed
ga-> open racc.ctl
Scanning description file: racc.ctl
Data file amsub_20140304_1838bin.dat is open as file 2
ga-> query file
File 1 : rainfall - ch4
Descriptor: racc.ctl
Binary: amsub_20140304_1838bin.dat
Type = Gridded
Xsize = 315 Ysize = 250 Zsize = 1 Tsize = 1 Esize = 1
Number of Variables = 1
rr 0 -9,40,2 ** rr (mm/h)
ga-> d rr
Constant field. Value = 0
ga-> quit
No hardcopy metafile open
GX package terminated
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 140924.0638.bin.gif
Type: image/gif
Size: 10575 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20141117/2da2949e/attachment-0001.gif
-------------- next part --------------
ÿÿø*4!÷ÿ
-------------- next part --------------
A non-text attachment was scrubbed...
Name: amsub_20140304_1838bin.dat
Type: application/x-ns-proxy-autoconfig
Size: 315000 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20141117/2da2949e/attachment-0002.bin
-------------- next part --------------
A non-text attachment was scrubbed...
Name: amsub_20140304_1838.dat
Type: application/x-ns-proxy-autoconfig
Size: 732488 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20141117/2da2949e/attachment-0003.bin
-------------- next part --------------
'reinit'
'reset'
'open racc.ctl'
'enable print out.gms'
*'set grads off'
*TITULOS
'set string 2 bc'
'set strsiz 0.15'
'set font 5'
'draw string 5.5 8.2 INTENSIDAD DE PRECIPITACION (mm/h)'
'set string 1 bc'
'draw string 5.5 7.9 01/10/2014 13:38 UTC'
*TEXTO INFERIOR
'set strsiz 0.14'
'draw string 5.5 0.2 Tecnica 183WSL corrida en CONAE a partir de imagenes AMSU-B/NOAA19'
*DIBUJO DE LA VARIABLE
*'set mpdset /home/hidroestimador/sa'
'set gxout shaded'
*'set cint 3'
*'set ccolor 1'
'set clevs 1 2 4 6 8 10 12 15 20 25 30 35 40'
'set ccols 0 9 14 4 11 5 13 3 10 7 12 8 2 6'
'd rr'
*'d rr/10'
'run cbarn 0.8 1 9.5 3.8'
*'run /home/hidroestimador/cbarn 0.8 1 9.5 3.8'
*'printim /home/hidroestimador/HYDRO/imagenes/rainestimation/INSTANTANEA/141001.1338.bin.gif x750 white'
'printim 140924.0638.bin.gif x750 white'
*grafico='/home/hidroestimador/HYDRO/imagenes/rainestimation/INSTANTANEA/141001.1338.bin.gif'
'disable print'
*'quit'
-------------- next part --------------
dset amsub_20140304_1838bin.dat
undef -9.0
OPTIONS little_endian
*options yrev
*options yrev byteswapped
title rainfall - ch4
ydef 250 linear -44.95 0.1
xdef 315 linear -74.95 0.1
tdef 1 linear 17:45Z23Aug2013 30mn
*tdef 1 linear 13:38Z01oct14 30mn
*tdef 1 linear 11:37Z01oct2014 30mn
zdef 1 levels 1000
vars 1
*rr 0 99 goes data
rr 0 -9,40,2 ** rr (mm/h)
ENDVARS
More information about the gradsusr
mailing list