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