[gradsusr] script for Precipitable water
Tidiane Diedhiou
tige83diedhiou at gmail.com
Thu Nov 21 08:37:54 EST 2013
hi Davide
I'm very gratefull about the script you gave to me.
thank.
Tidiane
2013/11/21 Davide Sacchetti <davide.sacchetti at arpal.gov.it>
> 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
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20131121/048f6144/attachment-0001.html
More information about the gradsusr
mailing list