[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