<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV>Dear Sushant Puranik,</DIV>
<DIV> </DIV>
<DIV> Thank you !</DIV>
<DIV> </DIV>
<DIV> My purpose of writing a GrADS script using vint() is to obtain the vertically integrated moisture transport vector fields as shown in the paper by Professor Krishnamurti (attached). However, by adapting your method, I notice that the magnitude of vector scale is 1000 instead of 0.2 etc (see spmtest.doc). What is the unit of the output ? Something is missing and hope you can enlighten me. Should I use the rotational motion fields from the streamfunctions instead of the pure wind vectors to obtain the mositure transport ? I am also interested to have script to compute the dry and moist static energies, if available.</DIV>
<DIV> </DIV>
<DIV> Best regards.</DIV>
<DIV> </DIV>
<DIV>(Mr) OOI See-hai</DIV>
<DIV> </DIV>
<DIV><BR><BR>--- On <B>Wed, 25/8/10, sushant puranik <I><sushantpuranik@gmail.com></I></B> wrote:<BR></DIV>
<BLOCKQUOTE style="BORDER-LEFT: rgb(16,16,255) 2px solid; PADDING-LEFT: 5px; MARGIN-LEFT: 5px"><BR>From: sushant puranik <sushantpuranik@gmail.com><BR>Subject: Re: [gradsusr] use of vint for an area<BR>To: "GrADS Users Forum" <gradsusr@gradsusr.org><BR>Date: Wednesday, 25 August, 2010, 1:59 PM<BR><BR>
<DIV id=yiv1806871567>Hello Ooi See Hai<BR>I am using grads on Linux (FC 8) OS. Earlier while trying VINT() i also came across the same problem.<BR>After that i adopted following procedure<BR><BR>[Sushant@atmmtech1 ~]$ gradsnc<BR><BR>Grid Analysis and Display System (GrADS) Version 1.8SL11<BR>Copyright (c) 1988-2001 by Brian Doty <BR>Center for Ocean-Land-Atmosphere Studies<BR>Institute for Global Environment and Society <BR>All Rights Reserved<BR><BR>Config: v1.8SL11 32-bit little-endian readline sdf/xdf netcdf lats printim<BR><BR>Issue 'q config' command for more information.<BR><BR>Landscape mode? (no for portrait): <BR>GX Package Initialization: Size = 11 8.5 <BR>ga-> c<BR>ga-> sdfopen /home/Sushant/Desktop/NcepSpecifichumidity/<A href="http://shum.1989.nc/" rel=nofollow target=_blank>shum.1989.nc</A> <BR>Scanning self-describing file: /home/Sushant/Desktop/NcepSpecifichumidity/<A href="http://shum.1989.nc/" rel=nofollow
target=_blank>shum.1989.nc</A><BR>Found 1 data variables in SDF file.<BR>Found displayable variable shum with 8 levels in SDF file.<BR>SDF file /home/Sushant/Desktop/NcepSpecifichumidity/<A href="http://shum.1989.nc/" rel=nofollow target=_blank>shum.1989.nc</A> is open as file 1<BR>LON set to 0 360 <BR>LAT set to -90 90 <BR>LEV set to 1000 1000 <BR>Time values set: 1989:1:1:0 1989:1:1:0 <BR>ga-> sdfopen /home/Sushant/Desktop/Ncepmeriwind/data/<A href="http://vwnd.1989.nc/" rel=nofollow target=_blank>vwnd.1989.nc</A> <BR>Scanning self-describing file: /home/Sushant/Desktop/Ncepmeriwind/data/<A href="http://vwnd.1989.nc/" rel=nofollow target=_blank>vwnd.1989.nc</A><BR>Found 1 data variables in SDF file.<BR>Found displayable variable vwnd with 17 levels in SDF file.<BR>SDF file /home/Sushant/Desktop/Ncepmeriwind/data/<A href="http://vwnd.1989.nc/" rel=nofollow target=_blank>vwnd.1989.nc</A> is open as file 2<BR>ga-> set gxout
shaded<BR>ga-> set lev 1000<BR>LEV set to 1000 1000 <BR>ga-> define a=(shum*vwnd.2)<BR>Define memory allocation size = 42340 bytes<BR>ga-> define b=vint(lev,a,300)<BR>Define memory allocation size = 42340 bytes<BR>ga-> d b<BR>Contouring: -1500 to 1200 interval 300 <BR><BR>and the output of it i have attached here with this mail.<BR><BR>regards <BR><BR>Sushant Puranik<BR>Dept. of Atmospheric & Space Sciences,<BR>University of Pune,<BR>Pune-07,<BR>India.<BR><BR><BR><BR>
<DIV class=yiv1806871567gmail_quote>On Tue, Aug 24, 2010 at 11:14 AM, See Hai Ooi <SPAN dir=ltr><<A href="http://aa.mc310.mail.yahoo.com/mc/compose?to=axl419@yahoo.com" rel=nofollow target=_blank ymailto="mailto:axl419@yahoo.com">axl419@yahoo.com</A>></SPAN> wrote:<BR>
<BLOCKQUOTE style="BORDER-LEFT: rgb(204,204,204) 1px solid; MARGIN: 0pt 0pt 0pt 0.8ex; PADDING-LEFT: 1ex" class=yiv1806871567gmail_quote>
<TABLE border=0 cellSpacing=0 cellPadding=0>
<TBODY>
<TR>
<TD style="FONT-FAMILY: inherit; font-size-adjust: inherit; font-stretch: inherit" vAlign=top>
<DIV>Dear Mr Arlindo,</DIV>
<DIV> </DIV>
<DIV> It appears to me that the vint function is only valid for a point. When I employed it for an area, I encountered the following problem (I use window version of GrADS under VISTA OS) :-</DIV>
<DIV> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">C:\>grads</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Starting X server under \Resources\Xming</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Starting grads under c:\OPENGR~1\Contents\Cygwin\Versions\20A7OG~1.3\i686 ...</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3 face="Times New Roman"> </FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Copyright (c) 1988-2008 by Brian Doty and the</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Institute for Global Environment and Society (IGES)</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">GrADS comes with ABSOLUTELY NO WARRANTY</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">See file COPYRIGHT for more information</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3 face="Times New Roman"> </FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf hdf4-sds hdf5</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">opendap-grids,stn athena geotiff</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Issue 'q config' command for more information.</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Loading User Defined Extensions table </bin/gex/udxt> ... ok.</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Landscape mode? ('n' for portrait):</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">GX Package Initialization: Size = 11 8.5</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">ga-> sdfopen i:\laptop\cdc\shum\<A href="http://shum.2010.nc/" rel=nofollow target=_blank>shum.2010.nc</A></FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Scanning self-describing file:<SPAN> </SPAN>i:\laptop\cdc\shum\<A href="http://shum.2010.nc/" rel=nofollow target=_blank>shum.2010.nc</A></FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">SDF file i:\laptop\cdc\shum\<A href="http://shum.2010.nc/" rel=nofollow target=_blank>shum.2010.nc</A> is open as file 1</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">LON set to 0 360</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">LAT set to -90 90</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">LEV set to 1000 1000</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Time values set: 2010:1:1:0 2010:1:1:0</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">E set to 1 1</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">ga-> a=vint(1000,shum,500)</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Error from VINT:<SPAN> </SPAN>Incompatible grids.</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Define memory allocation size = 84680 bytes</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman"><SPAN> </SPAN>6 [main] grads 3100 _cygtls::handle_exceptions: Exception: STATUS_ACCESS_V</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">IOLATION</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Exception: STATUS_ACCESS_VIOLATION at eip=610C3938</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><SPAN lang=ES><FONT size=3><FONT face="Times New Roman">eax=5F676E69 ebx=11A74040 ecx=00000026 edx=7373696D esi=00000069 edi=11A73E48</FONT></FONT></SPAN> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">ebp=00000328 esp=0022A9D0 program=C:\OpenGrADS\Contents\Cygwin\Versions\2.0.a7.o</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">ga.3\i686\grads.exe, pid 3100, thread main</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">cs=001B ds=0023 es=0023 fs=003B gs=0000 ss=0023</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Stack trace:</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">Frame<SPAN> </SPAN>Function<SPAN> </SPAN>Args</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman"><SPAN> </SPAN>6910 [main] grads 3100 _cygtls::handle_exceptions: Exception: STATUS_ACCESS_V</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman">IOLATION</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3><FONT face="Times New Roman"><SPAN> </SPAN>8471 [main] grads 3100 _cygtls::handle_exceptions: Error while dumping state</FONT></FONT> </DIV>
<P style="MARGIN: 0cm 0cm 0pt" class=yiv1806871567MsoNormal><FONT size=3 face="Times New Roman">(probably corrupted stack)</FONT> </DIV>
<DIV><BR> Can you help to resolve my problem ?</DIV>
<DIV> </DIV>
<DIV> Thank you and best regards.</DIV>
<DIV> </DIV>
<DIV>Sincerely,</DIV>
<DIV> </DIV>
<DIV>Ooi See Hai</DIV>
<DIV> </DIV></TD></TR></TBODY></TABLE><BR></BLOCKQUOTE></DIV></DIV></BLOCKQUOTE></td></tr></table><br>