<div>Hi OOI </div>
<div> </div>
<div> </div>
<div>I don't know exactly but your question about the units can be solved if you divide the variable 'shum' by 1000 to make it kg/kg as in NCEP data it given as grams/kg.</div>
<div> </div>
<div> </div>
<div>hope it works</div>
<div> </div>
<div>Cheers</div>
<div>Rashed</div>
<div><br><br> </div>
<div class="gmail_quote">On Wed, Aug 25, 2010 at 3:19 PM, See Hai Ooi <span dir="ltr"><<a href="mailto:axl419@yahoo.com">axl419@yahoo.com</a>></span> wrote:<br>
<blockquote style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">
<table border="0" cellspacing="0" cellpadding="0">
<tbody>
<tr>
<td valign="top">
<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><<a href="mailto:sushantpuranik@gmail.com" target="_blank">sushantpuranik@gmail.com</a>></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 <<a href="mailto:sushantpuranik@gmail.com" target="_blank">sushantpuranik@gmail.com</a>><br>Subject: Re: [gradsusr] use of vint for an area<br>
To: "GrADS Users Forum" <<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a>><br>Date: Wednesday, 25 August, 2010, 1:59 PM
<div>
<div></div>
<div class="h5"><br><br>
<div>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>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">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">
<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"><font size="3"><font face="Times New Roman">C:\>grads</font></font> </p></td></tr></tbody></table></blockquote></div>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">Starting X server under \Resources\Xming</font></font> </p></div>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">Starting grads under c:\OPENGR~1\Contents\Cygwin\Versions\20A7OG~1.3\i686 ...</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3" face="Times New Roman"> </font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">Copyright (c) 1988-2008 by Brian Doty and the</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">Institute for Global Environment and Society (IGES)</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">GrADS comes with ABSOLUTELY NO WARRANTY</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">See file COPYRIGHT for more information</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3" face="Times New Roman"> </font>
<p style="MARGIN: 0cm 0cm 0pt"><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>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">opendap-grids,stn athena geotiff</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">Issue 'q config' command for more information.</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">Loading User Defined Extensions table </bin/gex/udxt> ... ok.</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">Landscape mode? ('n' for portrait):</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">GX Package Initialization: Size = 11 8.5</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><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>
<p style="MARGIN: 0cm 0cm 0pt"><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>
<p style="MARGIN: 0cm 0cm 0pt"><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>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">LON set to 0 360</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">LAT set to -90 90</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">LEV set to 1000 1000</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">Time values set: 2010:1:1:0 2010:1:1:0</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">E set to 1 1</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">ga-> a=vint(1000,shum,500)</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">Error from VINT:<span> </span>Incompatible grids.</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">Define memory allocation size = 84680 bytes</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman"><span> </span>6 [main] grads 3100 _cygtls::handle_exceptions: Exception: STATUS_ACCESS_V</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">IOLATION</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">Exception: STATUS_ACCESS_VIOLATION at eip=610C3938</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><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>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">ebp=00000328 esp=0022A9D0 program=C:\OpenGrADS\Contents\Cygwin\Versions\2.0.a7.o</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">ga.3\i686\grads.exe, pid 3100, thread main</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">cs=001B ds=0023 es=0023 fs=003B gs=0000 ss=0023</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">Stack trace:</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">Frame<span> </span>Function<span> </span>Args</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman"><span> </span>6910 [main] grads 3100 _cygtls::handle_exceptions: Exception: STATUS_ACCESS_V</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman">IOLATION</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3"><font face="Times New Roman"><span> </span>8471 [main] grads 3100 _cygtls::handle_exceptions: Error while dumping state</font></font>
<p style="MARGIN: 0cm 0cm 0pt"><font size="3" face="Times New Roman">(probably corrupted stack)</font>
<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>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></div></div></blockquote></td></tr></tbody></table><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></blockquote></div><br>