<HTML>
<HEAD>
<META content="text/html; charset=iso-8859-1" http-equiv=Content-Type>
<META content="OPENWEBMAIL" name=GENERATOR>
</HEAD>
<BODY bgColor=#ffffff>
Dan,
<br />
<br />I haven't seen you had already undefined that variable in your original script.
<br />Well, if the processing time is really bothering you, you can try to code a more complex script.
<br />
<br />For example, a second try would be to fwrite a runtotal file (e.g., runtotal.gra), which would save the grib2 precip field for each iteration.
<br />In this case, the runtotal.gra file will change its content every time your script advance in time.
<br />
<br />For each iteration, your script will have to open a control file to access the runtotal.gra data (with a fix runtotal.ctl file), define a var with the sum of totalrunprecip field and the present grib2 precip field, close the runtotal.ctl file, and finally make a figure and fwrite a new runtotal.gra (over the older one). These 4 steps have to be repeated for each iteration, spending always the same time in each iteration. A cumbersome part will be to handle the time correctly to access the runtotal.gra data, since the runtotal.ctl file will have a different TDEF form from that of your grib2 file. An easy way to handle this, is to keep an runtotal.ctl with the same initial TDEF that that your grib2 files, and call the runtotal.gra data with something like: 'd precip.2(t=1) + precip.1', where precip.2(t=1) is the VAR from the runtotal.ctl file.
<br />
<br />Is it clear? Hope yes.
<br />
<br />Anyway, if you are interested in trying to code it, I can help you if necessary.
<br />Regards,
<br />Ricardo
<br />
<br /><font size="2"><b>On Sun, 12 Feb 2012 19:06:54 -0700, Dan Leins wrote</b>
<br />> I gave that a
shot and still no luck. I'm somewhat at a loss as to why the performance
degrades so much over time given this command.
<br />>
<br />> Dan
<br />>
<br />> On Sun, Feb 12, 2012 at 8:09 AM, Ricardo Hallak <span dir="ltr"><<a href="mailto:hallak@model.iag.usp.br">hallak@model.iag.usp.br</a>></span>
wrote:
<br />>
<blockquote style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;" class="gmail_quote">
<br />>
Hi
Dan,
<br />>
<br />> maybe
using:
<br />>
<br />> undefine
runtotal
<br />>
<br />> after each iteration may help in the
performance.
<br />>
<br />>
Regards,
<br />>
Ricardo
<br />>
<br />> <font><b>On Sat, 11 Feb 2012 21:02:18 -0700, Dan Leins
wrote</b>
<br />>
<br />> >
All,
<br />> >
<br />> > 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:
<br />> >
<br />> > 'define runtotal=sum(apcpsfc.1, t=2,
t='dis_t')'
<br />> > 'display
runtotal'
<br />> >
<br />> >
where dis_t is the hour I am currently working
on.
<br />> >
<br />> >
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.
<br />> >
<br />> > 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.
<br />> >
<br />> > Any help would be
great!
<br />> >
Thanks,
<br />> > Dan
Leins
<br />>
<br />> Ricardo Hallak
<br />>
<br />>
Department of Atmospheric Sciences (DCA)
<br />>
Institute of Astronomy, Geophysics and Atmospheric Sciences (IAG)
<br />>
University of São Paulo (USP)
<br />>
Rua do Matão, 1226 - Cidade Universitária - São Paulo SP 05508-900
<br />>
Phone.: <a target="_blank" value="+551130912850" href="tel:%2B55%20%2811%29%203091-2850">+55 (11) 3091-2850</a>; 3091-4703
<br />>
Fax: <a target="_blank" value="+551130914714" href="tel:%2B55%20%2811%29%203091-4714">+55 (11) 3091-4714</a>
<br />>
<br />>
</font>
<br />> _______________________________________________
<br />>
gradsusr mailing list
<br />>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>
<br />>
<a target="_blank" href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a>
<br />>
<br />> </blockquote>
<br />
<br />
<br />Ricardo Hallak
<br />
<br />
Department of Atmospheric Sciences (DCA)
<br />
Institute of Astronomy, Geophysics and Atmospheric Sciences (IAG)
<br />
University of São Paulo (USP)
<br />
Rua do Matão, 1226 - Cidade Universitária - São Paulo SP 05508-900
<br />
Phone.: +55 (11) 3091-2850; 3091-4703
<br />
Fax: +55 (11) 3091-4714
<br />
<br />
</font>
</BODY>
</HTML>