<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
<title></title>
</head>
<body bgcolor="#ffffff" text="#000000">
Hi Michael and See Hai Ooi,<br>
<br>
I would prefer not to be addressed as "Dr"; "Chuck", "Charles", or
"Hey There!" are fine...<br>
<br>
Thanks Michael, for helping with this! Here is the "mc.doc" attachment:<br>
<br>
----------------------------------------------------------------------------------<br>
* To obtain average moisture convergence at different sectorial points<br>
* by Ooi See Hai on 11 June 2005<br>
*<br>
'reinit'<br>
*<br>
'open d:\n2\ctl\tmp850.ctl'<br>
'open d:\n2\ctl\rh850.ctl'<br>
'open d:\n2\ctl\ugrd850.ctl'<br>
'open d:\n2\ctl\vgrd850.ctl'<br>
*<br>
'set gxout fwrite'<br>
'set fwrite d:\n2\mc.dat'<br>
*<br>
m=1<br>
while (m <= 732)<br>
'set t 'm<br>
*<br>
'tc = (tmp-273.16)'<br>
'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))'<br>
'vapr = 6.112*exp((17.67*td)/(td+243.5))'<br>
'e = vapr*1.001+(850-100)/900*0.0034'<br>
'mixr = 0.62197*(e/(850-e))*1000'<br>
'mconv = (-1)*hdivg((ugrd.3)*mixr,(vgrd.4)*mixr)*1e4'<br>
*<br>
* at 17.5N 112.5E<br>
*<br>
'mc=aave(mconv,lon=110,lon=115,lat=15,lat=20)'<br>
'd mc'<br>
*<br>
* at 12.N 112.5E<br>
*<br>
'mc=aave(mconv,lon=110,lon=115,lat=10,lat=15)'<br>
'd mc'<br>
*<br>
* at 7.5N 107.5E<br>
*<br>
'mc=aave(mconv,lon=105,lon=110,lat=5,lat=10)'<br>
'd mc'<br>
*<br>
* at 2.5N 107.5E<br>
*<br>
'mc=aave(mconv,lon=105,lon=110,lat=0,lat=5)'<br>
'd mc'<br>
*<br>
*<br>
m=m+2<br>
endwhile<br>
*<br>
'disable fwrite'<br>
*<br>
'reinit'<br>
<br>
<br>
<br>
<br>
<br>
dset d:\n2\mc.dat<br>
undef -9.99e+33<br>
title 850 hPa moisture convergence [*1e4/sec]<br>
xdef 1 linear 1 1<br>
ydef 1 linear 1 1<br>
zdef 1 linear 1 1<br>
tdef 366 linear 00z01jan2000 1dy<br>
vars 4<br>
mc1 1 99 1 * moisture convergence at 17.5N 112.5E<br>
mc2 1 99 1 * moisture convergence at 12.5N 112.5E<br>
mc3 1 99 1 * moisture convergence at 7.5N 107.5E<br>
mc4 1 99 1 * moisture convergence at 2.5N 107.5E<br>
endvars<br>
<br>
<br>
<br>
<br>
ga-> open d:\n2\mc<br>
Scanning description file: d:\n2\mc<br>
Data file d:\n2\mc.dat is open as file 1<br>
LON set to 1 1<br>
LAT set to 1 1<br>
LEV set to 1 1<br>
Time values set: 2000:1:1:0 2000:1:1:0<br>
ga-> set grads off<br>
ga-> set parea 1.0 10.0 1.0 7.0<br>
ga-> set lon 1<br>
LON set to 1 1<br>
ga-> set lat 1<br>
LAT set to 1 1<br>
ga-> set t 275 366<br>
Time values set: 2000:10:1:0 2000:12:31:0<br>
ga-> set xflip on<br>
ga-> set vrange -1.2 1.2<br>
1-D axis limits set: -1.2 to 1.2<br>
ga-> set yflip on<br>
ga-> d mc1<br>
ga-> d mc2<br>
ga-> d mc3<br>
ga-> d mc4<br>
ga-> define zl=0<br>
Define memory allocation size = 368 bytes<br>
ga-> set cmark 0<br>
cmark = 0<br>
ga-> d zl<br>
ga-> set annot 1 3<br>
SET ANNOT values: color = 1 thickness = 3<br>
ga-> draw xlab Date<br>
ga-> draw ylab <- Convergence/Divergence -><br>
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)<br>
ga-> enable print d:\n2\mc.gmf<br>
ga-> print<br>
ga-> disable print<br>
Hardcopy output file is closed<br>
ga-><br>
<br>
<br>
<br>