I took Jennifer's solution and without fail it solved the problem. Thanks for the help everyone!<div><br>Dan<br><br><div class="gmail_quote">On Mon, Feb 13, 2012 at 6:05 AM, Jennifer Adams <span dir="ltr"><<a href="mailto:jma@cola.iges.org">jma@cola.iges.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><br><div><div class="im"><div>On Feb 11, 2012, at 11:02 PM, Dan Leins wrote:</div><br>
<blockquote type="cite">All,<div><br></div><div>I have a GrADS script I use, among other things, to calculate run total precip over my model domain throughout my model run (output format is grib2). I create images of numerous fields on an hourly basis (say 72hrs worth for example), so I produce 72 runtotal precip images. As the script runs and iterates through a loop, my code to calculate the runtotal looks like this:</div>
<div><br></div><div>'define runtotal=sum(apcpsfc.1, t=2, t='dis_t')'</div><div>'display runtotal'</div><div><br></div><div>where dis_t is the hour I am currently working on.</div></blockquote><div>
<br></div></div><div>It sounds like this expression is inside a loop where dis_t starts at 2 and increases to the end of the time series. For each repetition of the loop, you are re-doing all the I/O for each time step from t=2 to t=dis_t. Try something like this instead:</div>
<div><br></div><div>'set gxout fwrite'</div><div>'define runtotal=apcpsfc.1(t=2)'</div><div>'d runtotal'</div><div>t=3</div><div>while (t<=dis_t)</div><div> 'define runtotal=runtotal+apacpsfc.1(t='t')'</div>
<div> 'd runtotal'</div><div> t=t+1</div><div>endwhile</div><div>'disable fwrite'</div><div><br></div><div>--Jennifer</div><div><br></div><br><blockquote type="cite"><div class="im"><div><br></div><div>
My script works fine and fast early on, but slows down over time. Naturally, I assume the amount of time it takes to 'sum' hours of precip will increase as I get further into my run. However, if I kill and restart my script at an arbitrary time (dis_t=60) for example, this calculation speeds right back up again. Likewise if I comment out this portion of my script, the rest of my script moves along at the same pace the entire time. I'm thinking some resources aren't being released properly while the script is running since the problem seems to go away if I kill/restart my script. </div>
<div><br></div><div>Is there a more efficient way to calculate a runtotal than the method shown above? I've been careful to undefine any variables after each iteration of my loop so I don't think that's a factor. </div>
<div><br></div><div>Any help would be great!</div><div>Thanks,</div><div>Dan Leins</div></div><div class="im">
_______________________________________________<br>gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br><a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</div></blockquote></div><br><div>
<span style="border-spacing:0px 0px;text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:12px;white-space:normal;font-family:Helvetica;word-spacing:0px"><span style="border-spacing:0px 0px;text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:12px;white-space:normal;font-family:Helvetica;word-spacing:0px"><span style="border-spacing:0px 0px;text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:12px;white-space:normal;font-family:Helvetica;word-spacing:0px"><div>
--</div><div>Jennifer M. Adams</div><div>IGES/COLA</div><div>4041 Powder Mill Road, Suite 302</div><div>Calverton, MD 20705</div><div><a href="mailto:jma@cola.iges.org" target="_blank">jma@cola.iges.org</a></div><div><br>
</div><br></span></span></span>
</div>
<br></div><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div>