Thank you for suggestion. Actually i want to calculate following integral. <br><br> <meta http-equiv="CONTENT-TYPE" content="text/html; charset=utf-8"><title></title><meta name="GENERATOR" content="OpenOffice.org 2.3 (Linux)">
<style type="text/css">
<!--
@page { size: 8.5in 11in; margin: 0.79in }
P { margin-bottom: 0.08in }
-->
</style>
<p style="margin-bottom: 0in;"> <font face="Liberation Serif, serif">
1000</font></p>
<p style="margin-bottom: 0in;"><font face="Liberation Serif, serif">Q
= (1/g) ∫ (V q)dp.
</font></p>
<p style="margin-bottom: 0in;"> <font face="Liberation Serif, serif">400</font></p>
<br>I used this kind of command in grads<br><br>'define ps=1000'<br>
'define a=vint(ps,vwnd*shum,400)'<br>'d a/9.81'<br>with this command i am able to make a plot.<br><br>whether it is correct or not ?<br><br>Thanking you<br><br><br>Sushant<br><br><div class="gmail_quote">
On Fri, Jul 23, 2010 at 10:38 AM, MY <span dir="ltr"><<a href="mailto:funnytyphoon@gmail.com">funnytyphoon@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi,<br>
<br>
please try<br>
<br>
'define ps=1000'<br>
'd vint(ps,vwnd*shum,400)'<br>
<br>
best wishes,<br>
++++++++++++++++++++++++++++++++++++++++++++++++++++<br>
YING, Ming<br>
Shanghai Typhoon Institute of China Meteorological Administration<br>
166 Puxi Road, Xuhui District, Shanghai, 200030<br>
PR China<br>
++++++++++++++++++++++++++++++++++++++++++++++++++++<br>
<br>
<br>
<br>
<br>
<br>
2010/7/23 sushant puranik <<a href="mailto:sushantpuranik@gmail.com">sushantpuranik@gmail.com</a>>:<br>
<div><div></div><div class="h5">> Hi<br>
> I am trying to find vertical integral in GrADS using "vint" command. For<br>
> that purpose i opened vwnd and specific humidity file as first and second<br>
> respectively. The q file command for both the variable are as<br>
><br>
> ga-> q file<br>
> File 1 : mean daily NMC reanalysis (2008)<br>
> Descriptor: /home/Sushant/Desktop/<a href="http://10may08v.nc" target="_blank">10may08v.nc</a><br>
> Binary: /home/Sushant/Desktop/<a href="http://10may08v.nc" target="_blank">10may08v.nc</a><br>
> Type = Gridded<br>
> Xsize = 29 Ysize = 25 Zsize = 8 Tsize = 1<br>
> Number of Variables = 1<br>
> vwnd 8 -999 mean Daily V wind<br>
><br>
> ga-> q file 2<br>
> File 2 : mean daily NMC reanalysis (2008)<br>
> Descriptor: /home/Sushant/Desktop/<a href="http://10may08q.nc" target="_blank">10may08q.nc</a><br>
> Binary: /home/Sushant/Desktop/<a href="http://10may08q.nc" target="_blank">10may08q.nc</a><br>
> Type = Gridded<br>
> Xsize = 29 Ysize = 25 Zsize = 8 Tsize = 1<br>
> Number of Variables = 1<br>
> shum 8 -999 mean Daily specific humidity<br>
><br>
> after this i used<br>
><br>
> ga-> d vint(1000,vwnd*shum.2,400)<br>
> gxstart=0, gxstop=29, gystart=0, gystop=25.<br>
> gzstart=0, gzstop=1, gtstart=0, gtstop=1.<br>
> gxstart=0, gxstop=29, gystart=0, gystop=25.<br>
> gzstart=0, gzstop=1, gtstart=0, gtstop=1.<br>
> Error from VINT: Incompatible grids.<br>
> Contouring: -0 to 2.4e+36 interval 3e+35<br>
> Segmentation fault<br>
><br>
> after showing this message, the grads window gets closed.<br>
><br>
> Can anyone has solution for this perticular problem.<br>
><br>
> Thanks in advance<br>
><br>
><br>
> --<br>
> Sushant<br>
><br>
><br>
</div></div><div><div></div><div class="h5">> _______________________________________________<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>
_______________________________________________<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><br clear="all"><br><br>