<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV>Dear Mr Charles Seman,</DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Great !&nbsp; You have resolved my problem.&nbsp; Grateful if you could kindly modify my final script (as attached) and delete/add any comments to ensure that the methodology is accurate.&nbsp; Again, thank you very much indeed.</DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Let me answer your questions and inform you certain approaches taken by me to finalise my script:-</DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.&nbsp; The vertical levels are 1000, 925, 850, 700, 600, 500, 400 and 300 hPa.</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.&nbsp; In spmtest2, I set it as 1000 hPa, i.e. the default level.</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3.&nbsp; In model.ctl, all the required elements (ps, u, v, hum etc) are in one grib file.&nbsp; For my case, I need to open an additional file of mean sea level pressure.&nbsp; I am not sure how would the vint() &nbsp;function read this mean sea level pressure file, uwnd and vwnd files of standard levels.&nbsp;&nbsp;My "lev" does not vary across the lon/lat domain but it does vary vertically.</DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;In accordance with your suggestions,&nbsp;my corrected analysis is similar to that by Professor Krishnamurthi (He also used NCEP/NCAR Reanalysis Data.).&nbsp; Regretably, he did not mention the unit of the vector scale (just 0.1).&nbsp; In GrADS, the scale factor (f) is set to 100.&nbsp; Am I right to divide my final result by 100&nbsp; (to get 4 kg/ms instead of 0.1 ??) ?&nbsp;&nbsp;I think I understand why GrADS set the top level as 275.&nbsp; Is it because the next standard level is 250 hPa ?&nbsp; In fact, I found out that it appears to have no impact on the output by using either 300 or 275.&nbsp; Can I just leave the standard levels as they are ?</DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I believe I can proceed to write my script to obtain vertically integrated moist static energy.&nbsp; I have one question to ask.&nbsp; For specific humidity, the topmost level is 300 hPa.&nbsp; For others, like winds, the levels extend to 50 hPa.&nbsp; As you know, the moist static energy consists of internal energy, potential energy and latent energy.&nbsp; Can I sum all these components (say, specific humidity to 300 hPa while winds, to 100 hPa) with differing numbers of standard levels ?</DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Best regards.</DIV>
<DIV>&nbsp;</DIV>
<DIV>Ooi</DIV>
<DIV><BR><BR>--- On <B>Sat, 28/8/10, Charles Seman <I>&lt;Charles.Seman@noaa.gov&gt;</I></B> wrote:<BR></DIV>
<BLOCKQUOTE style="BORDER-LEFT: rgb(16,16,255) 2px solid; PADDING-LEFT: 5px; MARGIN-LEFT: 5px"><BR>From: Charles Seman &lt;Charles.Seman@noaa.gov&gt;<BR>Subject: Re: [gradsusr] use of vint for an area<BR>To: "See Hai Ooi" &lt;axl419@yahoo.com&gt;<BR>Date: Saturday, 28 August, 2010, 2:13 AM<BR><BR>
<DIV id=yiv1148383546>Dear See Hai Ooi,<BR><BR>Regarding setting a vertical pressure range for defining "um" and "vm":&nbsp; do this first, then go back to set to a single level before defining "ium" and "ivm" (you can try it both ways, with a varying vertical dimension range and when the vertical dimension is set to a single level)... I tried it for a local dataset and found the vertically integrated specific humidity fields were the same for both methods:&nbsp; when the vertical level is set to a single level a 2D (x,y) field is generated by "vint" and when the vertical levels are set to a range then a 3D (x,y,p) field is generated by "vint".&nbsp; I did a difference of "isphum2(x,y,p)" and "isphum1(x,y)" and the result was zero.&nbsp; Note, by defining at a single level "isphum1(x,y)" becomes a "wild card" variable (<A class=yiv1148383546moz-txt-link-freetext href="http://grads.iges.org/grads/gadoc/variable.html#new" rel=nofollow
 target=_blank>http://grads.iges.org/grads/gadoc/variable.html#new</A>) and can be used if and when the vertical level is set outside the range used to define "isphum2(x,y,p)"... please see a copy of a GrADS session below where "isphum1(x,y)" and "isphum2(x,y,p)" are defined as shown:<BR>---<BR>ga-&gt; q dims<BR>Default file number is: 1 <BR>X is varying&nbsp;&nbsp; Lon = 178.75 to 188.75&nbsp;&nbsp; X = 72 to 76<BR>Y is varying&nbsp;&nbsp; Lat = -1 to 9&nbsp;&nbsp; Y = 45 to 50<BR>Z is fixed&nbsp;&nbsp;&nbsp;&nbsp; Lev = 1000&nbsp; Z = 1<BR>T is fixed&nbsp;&nbsp;&nbsp;&nbsp; Time = 12Z15JAN1981&nbsp; T = 1<BR>ga-&gt; isphum1= vint(ps,sphum,125)<BR>Define memory allocation size = 120 bytes<BR>Name already DEFINEd:&nbsp; isphum1.&nbsp;&nbsp; Will be deleted and replaced.<BR>ga-&gt; set z 1 13<BR>LEV set to 1000 100 <BR>ga-&gt; isphum2= vint(ps,sphum,125)<BR>Define memory allocation size = 1560 bytes<BR>Name already DEFINEd:&nbsp; isphum2.&nbsp;&nbsp; Will
 be deleted and replaced.<BR>ga-&gt; set z 14<BR>LEV set to 70 70 <BR>ga-&gt; d isphum1<BR>Contouring: 17300 to 18300 interval 100 <BR>ga-&gt; d isphum2<BR>Cannot contour grid - all undefined values <BR>ga-&gt; d isphum2(z=13) <BR>Contouring: 17300 to 18300 interval 100 <BR>ga-&gt; d isphum2(z=13)-isphum1<BR>Constant field.&nbsp; Value = 0<BR>ga-&gt; d isphum2(z=12)-isphum1<BR>Constant field.&nbsp; Value = 0<BR>ga-&gt; d isphum2(z=11)-isphum1<BR>Constant field.&nbsp; Value = 0<BR>ga-&gt; d isphum2(z=10)-isphum1<BR>Constant field.&nbsp; Value = 0<BR>ga-&gt; d isphum2(z=9)-isphum1<BR>Constant field.&nbsp; Value = 0<BR>ga-&gt; d isphum2(z=8)-isphum1<BR>Constant field.&nbsp; Value = 0<BR>ga-&gt; d isphum2(z=7)-isphum1<BR>Constant field.&nbsp; Value = 0<BR>ga-&gt; d isphum2(z=6)-isphum1<BR>Constant field.&nbsp; Value = 0<BR>ga-&gt; d isphum2(z=5)-isphum1<BR>Constant field.&nbsp; Value = 0<BR>ga-&gt; d isphum2(z=4)-isphum1<BR>Constant field.&nbsp; Value =
 0<BR>ga-&gt; d isphum2(z=3)-isphum1<BR>Constant field.&nbsp; Value = 0<BR>ga-&gt; d isphum2(z=2)-isphum1<BR>Constant field.&nbsp; Value = 0<BR>ga-&gt; d isphum2(z=1)-isphum1<BR>Constant field.&nbsp; Value = 0<BR><BR>---<BR><BR>I hope this helps,<BR>Chuck<BR><BR>Charles Seman wrote: 
<BLOCKQUOTE type="cite">Dear See Hai Ooi,<BR><BR>I took a look at your "spmtest2.doc" including the question regarding "Usage Note No 4 in GrADS document", copied here again from <A class=yiv1148383546moz-txt-link-freetext href="http://grads.iges.org/grads/gadoc/gradfuncvint.html" rel=nofollow target=_blank>http://grads.iges.org/grads/gadoc/gradfuncvint.html</A><BR>"It is usually a good idea to make the <CODE><I>top</I></CODE> pressure value to be at the top of a layer, which is midway between grid points. For example, if the default file (and the data) have pressure levels of ...,500,400,300,250,... then a good value for <CODE><I>top</I></CODE> might be 275, the value at the top of the layer that extends from 350 to 275 mb. "<BR><BR>My understanding of putting the top at 275 is this:&nbsp; given data pressure levels at 500, 400, 300, 250,... then the midpoint levels would be at 450, 350, 275 which is the halfway distance between immediate lower and
 upper levels... the layer for the data at 400: 450-350, and the layer for the data at 300: 350-275&nbsp; &lt;- note, this top layer is not symmetric about the data level...<BR><BR>Some questions and comments about the code in "spmtest2.doc"...<BR>1) what are the vertical pressure levels in your datasets?<BR>2) what is the range of vertical pressure levels set to before you define "um" and "vm"? if it is only set for a single level, try setting the vertical levels to an appropriate range (taking into account "Usage Note No 4 in GrADS document" for guidance)<BR>3) instead of using "lev" for the surface pressure (see "Example 1" in GrADS document: vint(ps,q,275)), I would try bringing in a surface pressure field appropriate for the data (may need to access another dataset?)... say "ps(x,y)" (does "lev" vary across the lon/lat domain?)... you could also try "Example 2" in GrADS document with your data to check on your code...<BR><BR>I hope this
 helps,<BR>Chuck<BR><BR><BR>********************************************************************<BR>Charles Seman&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <A class=yiv1148383546moz-txt-link-abbreviated href="http://aa.mc310.mail.yahoo.com/mc/compose?to=Charles.Seman@noaa.gov" rel=nofollow target=_blank ymailto="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</A><BR>U.S. Department of Commerce / NOAA / OAR<BR>Geophysical Fluid Dynamics Laboratory&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; voice: (609) 452-6547<BR>201 Forrestal Road&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fax: (609) 987-5063<BR>Princeton, NJ&nbsp;
 08540-6649&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <A class=yiv1148383546moz-txt-link-freetext href="http://www.gfdl.noaa.gov/~cjs/" rel=nofollow target=_blank>http://www.gfdl.noaa.gov/~cjs/</A><BR>********************************************************************<BR><BR></BLOCKQUOTE></DIV></BLOCKQUOTE></td></tr></table><br>