<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Good Afternoon Users,</p>
<p><br>
</p>
<p>I have a grib1 file that contains 500mb heights from 1985 to 2014, using the NCEP/NCAR Reanalysis Dataset. I am trying to create an averaged climatology file to create monthly anomalies with. I've written a small GrADS script to try and create this climo
 file but, I'm not having much luck with it. The script outputs a .grib1 file but, no control file, as advertised.</p>
<p><br>
</p>
<p>What am I doing wrong?</p>
<p><br>
</p>
<p>Script Code:</p>
<p><br>
</p>
<p>function pull (args)<br>
<br>
'reinit'<br>
<br>
rc = gsfallow('on')<br>
<br>
currdate = subwrd(args,1)<br>
curr_run = subwrd(args,2)<br>
regions = subwrd(args,3)<br>
month = subwrd(args,4)<br>
date = sublin(currdate,1)<br>
mon = substr(date,5,2)<br>
day = substr(date,7,2)<br>
yr = substr(date,1,4)<br>
<br>
path1 = '/mnt/homes/perrya/'<br>
path2 = 'new.climate.atlas/'<br>
model = ncep_reanal<br>
<br>
var = hgt500<br>
<br>
'open 'path1''path2'hgt500.mon.1985.2014.grib.ctl'<br>
<br>
count = 0<br>
rec = 360<br>
incr = 1; t = 0<br>
fhr = 0<br>
inc = 1<br>
<br>
'set gxout fwrite'<br>
'set fwrite 'path1''path2'hgt500.1985.2014.mean.grib'<br>
<br>
while(count&lt;=rec)<br>
<br>
'set t 't<br>
<br>
say t<br>
<br>
'set lon -320 40' <br>
'set lat -30 90'<br>
<br>
'define hgt500 = (hgt500mb.1/100)'<br>
<br>
'define climo=ave(hgt500,t=0,t='t',1mo)'<br>
<br>
'd climo'<br>
<br>
count = count &#43; incr<br>
t = t &#43; incr<br>
fhr = fhr &#43; inc<br>
<br>
endwhile<br>
<br>
'quit'<br>
<br>
</p>
<p><br>
</p>
<div id="Signature">
<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:Calibri,Arial,Helvetica,sans-serif">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<div style="font-family:Tahoma; font-size:13px">
<div style="font-size:13px; font-family:Tahoma">
<div style="font-size:13px; font-family:Tahoma">
<div style="font-size:13px; font-family:Tahoma">
<div style="font-size:13px; font-family:Tahoma">
<div style="font-size:13px; font-family:Tahoma">
<div style="font-size:13px; font-family:Tahoma">
<div style="font-size:13px; font-family:Tahoma">
<div style="font-size:13px; font-family:Tahoma">
<div style="font-size:13px; font-family:Tahoma"><span style="font-family:Calibri,Arial,Helvetica,sans-serif"><font size="3"><font face="Arial"><font face="Tahoma" size="2"><font size="2"></font></font></font></font>
<div style=""><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:11pt; color:rgb(0,111,201)"><font style="color:rgb(0,111,201)"><span style="color:rgb(0,111,201)"><br>
Aaron D. Perry</span></font><font style="color:rgb(0,111,201)"><br>
<span style="font-family:Calibri,Arial,Helvetica,sans-serif"><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:11pt; color:rgb(0,111,201)"><font style="color:rgb(0,111,201)"><span style="color:rgb(0,111,201)">Class of 2016</span><span style="color:rgb(0,111,201)"></span></font><font style="color:rgb(0,111,201)"><span style="color:rgb(0,111,201)"></span></font></span></span><br style="color:rgb(0,111,201)">
<span style="color:rgb(0,111,201)">Lyndon State College</span><br>
Lyndonville, VT<span style="color:rgb(0,111,201)"></span></font><font style="color:rgb(0,111,201)"><br style="color:rgb(0,111,201)">
<span style="color:rgb(0,111,201)">Twitter: </span><a id="LPNoLP" href="https://twitter.com/AaronPerryLSC" tabindex="0"><span style="color:rgb(0,111,201); text-decoration:underline"></span></a><a id="LPNoLP" href="https://www.twitter.com/AaronPerryLSC"></a><a id="LPNoLP" href="https://twitter.com/AaronPerryLSC" tabindex="0"><span style="color:rgb(0,111,201); text-decoration:underline">@AaronPerryLSC</span></a><span style="color:rgb(0,111,201)"></span><font style="color:rgb(0,111,201)"><span style="color:rgb(0,111,201)"></span></font><br style="color:rgb(0,111,201)">
<span style="color:rgb(0,111,201)"></span></font></span></div>
<span style="font-size:11pt; color:rgb(0,111,201)"></span><span style="font-size:11pt; color:rgb(0,111,201)"></span><span style="font-size:11pt; color:rgb(0,111,201)"></span>
<div style=""><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:11pt; color:rgb(0,111,201)"><font style="color:rgb(0,111,201)"><span style="color:rgb(0,111,201)">Mobile: 617-780-4312<br>
<span style="font-family:Calibri,Arial,Helvetica,sans-serif"><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:11pt; color:rgb(0,111,201)"><font style="color:rgb(0,111,201)"><font style="color:rgb(0,111,201)"><span style="color:rgb(0,111,201)">Web:
</span><span style="color:rgb(0,111,201)"></span></font></font></span></span><a id="LPNoLP" href="http://shearanalytics.net/"><span style="font-family:Calibri,Arial,Helvetica,sans-serif"><span style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:11pt; color:rgb(0,111,201)"><font style="color:rgb(0,111,201)"><font style="color:rgb(0,111,201)"><span style="color:rgb(0,111,201)">shearanalytics.net</span></font></font></span></span></a><br>
Email: <a id="LPNoLP" href="http://Aaron.Perry@lyndonstate.edu">Aaron.Perry@lyndonstate.edu</a><br>
</span></font></span></div>
</span></div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>