Moisture Convergence

[Charles Seman]

Hi Michael and See Hai Ooi,

I would prefer not to be addressed as "Dr";  "Chuck", "Charles", or "Hey
There!" are fine...

Thanks Michael, for helping with this!  Here is the "mc.doc" attachment:

----------------------------------------------------------------------------------
* To obtain average moisture convergence at different sectorial points
* by Ooi See Hai on 11 June 2005
*
'reinit'
*
'open d:\n2\ctl\tmp850.ctl'
'open d:\n2\ctl\rh850.ctl'
'open d:\n2\ctl\ugrd850.ctl'
'open d:\n2\ctl\vgrd850.ctl'
*
 'set gxout fwrite'
 'set fwrite d:\n2\mc.dat'
*
 m=1
 while (m <= 732)
  'set t 'm
*
  'tc = (tmp-273.16)'
  'td = tc-((14.55+0.114*tc)*(1-0.01*rh.2) +
pow((2.5+0.007*tc)*(1-0.01*rh.2),3) +
(15.9+0.117*tc)*pow((1-0.01*rh.2),14))'
  'vapr = 6.112*exp((17.67*td)/(td+243.5))'
  'e = vapr*1.001+(850-100)/900*0.0034'
  'mixr = 0.62197*(e/(850-e))*1000'
  'mconv = (-1)*hdivg((ugrd.3)*mixr,(vgrd.4)*mixr)*1e4'
*
* at 17.5N 112.5E
*
  'mc=aave(mconv,lon=110,lon=115,lat=15,lat=20)'
  'd mc'
*
* at 12.N 112.5E
*
  'mc=aave(mconv,lon=110,lon=115,lat=10,lat=15)'
  'd mc'
*
* at 7.5N 107.5E
*
  'mc=aave(mconv,lon=105,lon=110,lat=5,lat=10)'
  'd mc'
*
* at 2.5N 107.5E
*
  'mc=aave(mconv,lon=105,lon=110,lat=0,lat=5)'
  'd mc'
*
*
  m=m+2
 endwhile
*
 'disable fwrite'
*
'reinit'





dset  d:\n2\mc.dat
undef -9.99e+33
title 850 hPa moisture convergence [*1e4/sec]
xdef   1 linear 1 1
ydef   1 linear 1 1
zdef   1 linear 1 1
tdef   366 linear 00z01jan2000 1dy
vars 4
mc1   1 99 1  * moisture convergence at 17.5N 112.5E
mc2   1 99 1  * moisture convergence at 12.5N 112.5E
mc3   1 99 1  * moisture convergence at  7.5N 107.5E
mc4   1 99 1  * moisture convergence at  2.5N 107.5E
endvars




ga-> open d:\n2\mc
Scanning description file:  d:\n2\mc
Data file d:\n2\mc.dat is open as file 1
LON set to 1 1
LAT set to 1 1
LEV set to 1 1
Time values set: 2000:1:1:0 2000:1:1:0
ga-> set grads off
ga-> set parea 1.0 10.0 1.0 7.0
ga-> set lon 1
LON set to 1 1
ga-> set lat 1
LAT set to 1 1
ga-> set t 275 366
Time values set: 2000:10:1:0 2000:12:31:0
ga-> set xflip on
ga-> set vrange -1.2 1.2
1-D axis limits set: -1.2 to 1.2
ga-> set yflip on
ga-> d mc1
ga-> d mc2
ga-> d mc3
ga-> d mc4
ga-> define zl=0
Define memory allocation size = 368 bytes
ga-> set cmark 0
cmark = 0
ga-> d zl
ga-> set annot 1 3
SET ANNOT values:  color = 1  thickness = 3
ga-> draw xlab Date
ga-> draw ylab <- Convergence/Divergence ->
ga-> draw title Moisture Convergence (*1e4 /sec)\ at 17.5E 112.5N (open
circle), 12.5E 112.5N (closed circle)\ 7.5N 107.5E (open square) & 2.5N
107.5E (closed square)
ga-> enable print d:\n2\mc.gmf
ga-> print
ga-> disable print
Hardcopy output file is closed
ga->



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20050614/774d0ef9/attachment.html