[gradsusr] script for Precipitable water
Davide Sacchetti
davide.sacchetti at arpal.gov.it
Thu Nov 21 05:27:44 EST 2013
I don't know why you are dividing by 9.8, but the problem is another
one:
vint(ps,q,300) is different from vint(1000,q,300)
ps is a vector, 1000 is a number
try:
vint(const(spfhprs,1000), spfhprs, 300)
it should work
bye bye
Davide
On Wed, 2013-11-20 at 15:25 +0000, Tidiane Diedhiou wrote:
> hi, i have used vint(ps,q,300) but it make an error. here is the
> result
>
> open Q??0_ECMWF_i2009083100.ctl
>
> LON set to -30 40
> LAT set to -10 45
> LEV set to 1000 1000
> Time values set: 2009:8:31:0 2009:8:31:Error from VINT: Incompatible
> grids. 0
> E set to 1 1
> ga-> q file
> File 1 : Q??0_ECMWF_i2009083100.grib2
> Descriptor: Q??0_ECMWF_i2009083100.ctl
> Binary: Q??0_ECMWF_i2009083100.grib2
> Type = Gridded
> Xsize = 71 Ysize = 56 Zsize = 5 Tsize = 41 Esize = 51
> Number of Variables = 1
> spfhprs 5 0,1,0 ** (1000 850 700 500 300) Specific Humidity
> [kg/kg]
> ga-> set t 1
> Time values set: 2009:8:31:0 2009:8:31:0
> ga-> set e 1
> E set to 1 1
> ga-> define PWA= vint(1000,spfhprs,300)/9.8
> Error from VINT: Incompatible grids.
> Define memory allocation size = 31808 bytes
> *** glibc detected *** grads: corrupted double-linked list: 0x0a15c358
> ***
> ======= Backtrace: =========
> /lib/i386-linux-gnu/libc.so.6(+0x6ebc2)[0xfa8bc2]
> /lib/i386-linux-gnu/libc.so.6(+0x6f91b)[0xfa991b]
> /lib/i386-linux-gnu/libc.so.6(cfree+0x6d)[0xfac94d]
> grads(gadef+0x17d1)[0x80ca0d1]
> grads(gacmd+0x1b54)[0x80e1ac4]
> grads(main+0x5d6)[0x80520b6]
> /lib/i386-linux-gnu/libc.so.6(__libc_start_main+0xf3)[0xf53113]
> grads[0x8052661]
> ======= Memory map: ========
> 00110000-001e4000 r-xp 00000000 08:05
> 4980994 /usr/lib/libnetcdf.so.6.0.0
> 001e4000-00221000 r--p 000d4000 08:05
> 4980994 /usr/lib/libnetcdf.so.6.0.0
> 00221000-00223000 rw-p 00111000 08:05
> 4980994 /usr/lib/libnetcdf.so.6.0.0
> 00223000-00228000 rw-p 00000000 00:00 0
> 00229000-00231000 r-xp 00000000 08:05
> 4981902 /usr/lib/libshp.so.1.0.1
> 00231000-00232000 r--p 00007000 08:05
> 4981902 /usr/lib/libshp.so.1.0.1
> 00232000-00233000 rw-p 00008000 08:05
> 4981902 /usr/lib/libshp.so.1.0.1
> 00233000-0025b000 r-xp 00000000 08:05
> 132035 /lib/i386-linux-gnu/libm-2.13.so
> 0025b000-0025c000 r--p 00028000 08:05
> 132035 /lib/i386-linux-gnu/libm-2.13.so
> 0025c000-0025d000 rw-p 00029000 08:05
> 132035 /lib/i386-linux-gnu/libm-2.13.so
> 0025d000-0027a000 r-xp 00000000 08:05
> 4986254 /usr/lib/i386-linux-gnu/libxcb.so.1.1.0
> 0027a000-0027b000 r--p 0001c000 08:05
> 4986254 /usr/lib/i386-linux-gnu/libxcb.so.1.1.0
> 0027b000-0027c000 rw-p 0001d000 08:05
> 4986254 /usr/lib/i386-linux-gnu/libxcb.so.1.1.0
> 0027c000-0027f000 r-xp 00000000 08:05
> 132016 /lib/i386-linux-gnu/libdl-2.13.so
> 0027f000-00280000 r--p 00002000 08:05
> 132016 /lib/i386-linux-gnu/libdl-2.13.so
> 00280000-00281000 rw-p 00003000 08:05
> 132016 /lib/i386-linux-gnu/libdl-2.13.so
> 00281000-002b3000 r-xp 00000000 08:05
> 4986059 /usr/lib/i386-linux-gnu/libfontconfig.so.1.4.4
> 002b3000-002b4000 ---p 00032000 08:05
> 4986059 /usr/lib/i386-linux-gnu/libfontconfig.so.1.4.4
> 002b4000-002b5000 r--p 00032000 08:05
> 4986059 /usr/lib/i386-linux-gnu/libfontconfig.so.1.4.4
> 002b5000-002b6000 rw-p 00033000 08:05
> 4986059 /usr/lib/i386-linux-gnu/libfontconfig.so.1.4.4
> 002b6000-00348000 r-xp 00000000 08:05
> 4983725 /usr/lib/i386-linux-gnu/libfreetype.so.6.6.2
> 00348000-0034c000 r--p 00091000 08:05
> 4983725 /usr/lib/i386-linux-gnu/libfreetype.so.6.6.2
> 0034c000-0034d000 rw-p 00095000 08:05
> 4983725 /usr/lib/i386-linux-gnu/libfreetype.so.6.6.2
> 0034d000-00350000 r-xp 00000000 08:05
> 132032 /lib/i386-linux-gnu/libgpg-error.so.0.8.0
> 00350000-00351000 r--p 00002000 08:05
> 132032 /lib/i386-linux-gnu/libgpg-error.so.0.8.0
> 00351000-00352000 rw-p 00003000 08:05
> 132032 /lib/i386-linux-gnu/libgpg-error.so.0.8.0
> 00352000-00353000 r-xp 00000000 00:00 0 [vdso]
> 00353000-00395000 r-xp 00000000 08:05
> 4986109 /usr/lib/i386-linux-gnu/libjasper.so.1.0.0
> 00395000-00396000 r--p 00041000 08:05
> 4986109 /usr/lib/i386-linux-gnu/libjasper.so.1.0.0
> 00396000-00399000 rw-p 00042000 08:05
> 4986109 /usr/lib/i386-linux-gnu/libjasper.so.1.0.0
> 00399000-0039f000 rw-p 00000000 00:00 0
> 0039f000-003b6000 r-xp 00000000 08:05
> 132065 /lib/i386-linux-gnu/libpthread-2.13.so
> 003b6000-003b7000 r--p 00016000 08:05
> 132065 /lib/i386-linux-gnu/libpthread-2.13.so
> 003b7000-003b8000 rw-p 00017000 08:05
> 132065 /lib/i386-linux-gnu/libpthread-2.13.so
> 003b8000-003ba000 rw-p 00000000 00:00 0
> 003ba000-003da000 r-xp 00000000 08:05
> 4980985 /usr/lib/libhdf5_hl.so.6.0.3
> 003da000-003db000 r--p 0001f000 08:05
> 4980985 /usr/lib/libhdf5_hl.so.6.0.3
> 003db000-003dc000 rw-p 00020000 08:05
> 4980985 /usr/lib/libhdf5_hl.so.6.0.3
> 003dc000-003ed000 rw-p 00000000 00:00 0
> 003ed000-003f2000 r-xp 00000000 08:05
> 4985966 /usr/lib/i386-linux-gnu/libXdmcp.so.6.0.0
> 003f2000-003f3000 r--p 00004000 08:05
> 4985966 /usr/lib/i386-linux-gnu/libXdmcp.so.6.0.0
> 003f3000-003f4000 rw-p 00005000 08:05
> 4985966 /usr/lib/i386-linux-gnu/libXdmcp.so.6.0.0
> 003f7000-00411000 r-xp 00000000 08:05
> 4983318 /usr/lib/libudunits2.so.0.1.0
> 00411000-00412000 r--p 0001a000 08:05
> 4983318 /usr/lib/libudunits2.so.0.1.0
> 00412000-00413000 rw-p 0001b000 08:05
> 4983318 /usr/lib/libudunits2.so.0.1.0
> 00413000-00467000 r-xp 00000000 08:05
> 4986035 /usr/lib/i386-linux-gnu/libcurl-gnutls.so.4.2.0
> 00467000-00468000 ---p 00054000 08:05
> 4986035 /usr/lib/i386-linux-gnu/libcurl-gnutls.so.4.2.0
> 00468000-00469000 r--p 00054000 08:05
> 4986035 /usr/lib/i386-linux-gnu/libcurl-gnutls.so.4.2.0
> 00469000-0046a000 rw-p 00055000 08:05
> 4986035 /usr/lib/i386-linux-gnu/libcurl-gnutls.so.4.2.0
> 0046a000-0049b000 r-xp 00000000 08:05
> 4986102 /usr/lib/i386-linux-gnu/libidn.so.11.6.5
> 0049b000-0049c000 r--p 00030000 08:05
> 4986102 /usr/lib/i386-linux-gnu/libidn.so.11.6.5
> 0049c000-0049d000 rw-p 00031000 08:05
> 4986102 /usr/lib/i386-linux-gnu/libidn.so.11.6.5
> 0049d000-004aa000 r-xp 00000000 08:05
> 4986123 /usr/lib/i386-linux-gnu/liblber-2.4.so.2.7.0
> 004aa000-004ab000 r--p 0000c000 08:05
> 4986123 /usr/lib/i386-linux-gnu/liblber-2.4.so.2.7.0
> 004ab000-004ac000 rw-p 0000d000 08:05
> 4986123 /usr/lib/i386-linux-gnu/liblber-2.4.so.2.7.0
> 004ac000-004b3000 r-xp 00000000 08:05
> 132069 /lib/i386-linux-gnu/librt-2.13.so
> 004b3000-004b4000 r--p 00006000 08:05
> 132069 /lib/i386-linux-gnu/librt-2.13.so
> 004b4000-004b5000 rw-p 00007000 08:05
> 132069 /lib/i386-linux-gnu/librt-2.13.so
> 004b5000-004c5000 r-xp 00000000 08:05
> 4986225 /usr/lib/i386-linux-gnu/libtasn1.so.3.1.11
> 004c5000-004c6000 r--p 0000f000 08:05
> 4986225 /usr/lib/i386-linux-gnu/libtasn1.so.3.1.11
> 004c6000-004c7000 rw-p 00010000 08:05
> 4986225 /usr/lib/i386-linux-gnu/libtasn1.so.3.1.11
> 004c8000-004ee000 r-xp 00000000 08:05
> 132021 /lib/i386-linux-gnu/libexpat.so.1.5.2
> 004ee000-004ef000 ---p 00026000 08:05
> 132021 /lib/i386-linux-gnu/libexpat.so.1.5.2
> 004ef000-004f1000 r--p 00026000 08:05
> 132021 /lib/i386-linux-gnu/libexpat.so.1.5.2
> 004f1000-004f2000 rw-p 00028000 08:05
> 132021 /lib/i386-linux-gnu/libexpat.so.1.5.2
> 004f2000-004f4000 r-xp 00000000 08:05
> 132034 /lib/i386-linux-gnu/libkeyutils.so.1.3
> 004f4000-004f5000 r--p 00001000 08:05
> 132034 /lib/i386-linux-gnu/libkeyutil
>
>
>
> 2013/11/20 Davide Sacchetti <davide.sacchetti at arpal.gov.it>
> see: vint at http://grads.iges.org/grads/gadoc/gadocindex.html
> bye bye
> Davide
>
> On Wed, 2013-11-20 at 10:45 +0000, Tidiane Diedhiou wrote:
> > I want to calculate the precipitable water from the specific
> humidity.
> >
> > could you please tell me how to do the script in grads?
> >
> > thank you
> >
> > Tidiane
> >
>
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
>
> --
> --> Attenzione cambio indirizzo: davide.sacchetti at arpal.gov.it
> --
> Davide Sacchetti
> Centro Funzionale Meteo Idrologico di Protezione Civile della
> Regione Liguria
> ARPAL Unità Tecnica Complessa di livello Regionale
> V.le Brigare Partigiane 2 16121 Genova (I)
> tel: +39 010 6437535 fax: +39 010 6437520
> mail: davide.sacchetti at arpal.gov.it web:
> www.meteoliguria.it
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
--
--> Attenzione cambio indirizzo: davide.sacchetti at arpal.gov.it
--
Davide Sacchetti
Centro Funzionale Meteo Idrologico di Protezione Civile della Regione Liguria
ARPAL Unità Tecnica Complessa di livello Regionale
V.le Brigare Partigiane 2 16121 Genova (I)
tel: +39 010 6437535 fax: +39 010 6437520
mail: davide.sacchetti at arpal.gov.it web: www.meteoliguria.it
More information about the gradsusr
mailing list