<DIV>Good afternoon everyone,<BR><BR>I need to calculate the moisture flux and the moisture divergence of 1992 summer over the continental United States. The gs file I wrote is as following, but I am not sure if it is right, so anyone can look at this file and give me some comments? In addition, attached are the calculated result using this file.</DIV>
<DIV> </DIV>
<DIV>Best,</DIV>
<DIV>Zhongfang</DIV>
<DIV> </DIV>
<DIV>The gs file are:<BR>-------------------------------------------------<BR>'reinit'<BR>'set grid off'<BR>'set grads off'<BR>'set poli on'<BR>'set parea 1 10 1 8'<BR>'sdfopen shum.mon.mean.nc' <BR>'sdfopen uwnd.mon.mean.nc' <BR>'sdfopen vwnd.mon.mean.nc'<BR>'sdfopen pres.mon.mean.nc'<BR>'enable print sum_1992.gmf'<BR>'set lat 24.99 50.1'<BR>'set xlint 10'<BR>'set lon -130.01 -69.99'<BR>'set ylint 10'<BR>'set cint 20'<BR>'set t 1 600'<BR>'set font 0'<BR>'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6'<BR>'set mpdset cnworld mres'<BR>'set map 1 1 1'<BR>'draw map'<BR>'set gxout shaded'<BR>'set ylopts 1 4 0.18'<BR>'set xlopts 1 4 0.18'<BR>'set cthick 5'<BR>'define m=shum.1*uwnd.2'<BR>'define n=shum.1*vwnd.3'<BR>'define qu=vint(pres.4(z=1),m,700)/1000'<BR>'define qusum=ave(qu,t=534,t=536)'<BR>'define qv=vint(pres.4(z=1),n,700)/1000'<BR>'define qvsum=ave(qv,t=534,t=536)'<BR>'set t 6'<BR>'set gxout shaded'<BR>'d qusum;qvsum'<BR>'d hdivg(qusum,qvsum)'<BR>'set gxout vector'<BR>'set ccolor 1'<BR>'run cbar'<BR>'d qusum;qvsum'<BR>*'set cint 20'<BR>'draw string 1.9 7 Moisture transport'<BR>'draw string 9.5 7 JJA-1992'<BR>'print'<BR>'disable print'<BR>-------------------------------------------------<BR><BR><BR>Any help you can provide would be greatly appreciated, and thank you for your time.<BR><BR><BR><BR></DIV><BR><BR><SPAN title="neteasefooter">
<HR>
<A href="http://allyes.nie.163.com/main/adfclick?db=afanie&bid=1254,611,23&cid=145,4,1&sid=1351&show=ignore&url=http://tx2.163.com/fab.html">09年新晋3D主流网游《天下贰》，网易六年亿资打造</A> </SPAN><BR><BR><SPAN title="neteasefooter">
<HR>
<A href="http://allyes.nie.163.com/main/adfclick?db=afanie&bid=1258,613,23&cid=147,4,1&sid=1355&show=ignore&url=http://tx2.163.com/fab.html">"雪见"杨幂邀你共品3D国韵网游《天下贰》，海报免费领</A> </SPAN><br><br><span title="neteasefooter"/><hr/>
<a href="http://allyes.nie.163.com/main/adfclick?db=afanie&bid=1260,614,23&cid=148,4,1&sid=1357&show=ignore&url=http://tx2.163.com/fab.html">看陆川杨幂新片《琴棋书画》，品网易3D国韵网游《天下贰》</a>
</span>