Precipitation 3/6 h

James Burn j.burn at GEOS.COM
Tue Feb 3 06:10:42 EST 2009


Hi Fil

When looping (with _count) to plot a set of GFS 3hrly precip charts I use:

if (math_fmod(_count,2))
     'define APCPsfc6hr=APCPsfc'
     'set t '_count-1
     'define APCPsfcnew=APCPsfc6hr-APCPsfc'

     'd APCPsfcnew'
     'set t '_count
else
    'd APCPsfc'
endif

Which works nicely. I've not sorted the 3 hourly precip in a meteogram yet,
but I imagine something similar would work.

Good luck
James

On Fri, 23 Jan 2009 15:00:06 +0100, Fil <fil at METEOPT.COM> wrote:

>Hi everyone,
>
>I'm trying to plot a precipitation meteogram from GFS 3h data but, for
>example, +33h timestep has precipitation from +30h to +33h but +36h
timestep
>also has precipitation from +30h to +33h and +33h to +36h. So what I wanted
>to know is how can I subtract precipitation like 30h = 30h - 27h, 36h = 36h
>- 33h, 42h = 42h - 39h, etc, so that +30h/+36h/+42h/etc timestep doesn't
>include precipitation from +27h/+33h/+39h/etc timestep.
>
>Also how can I sum precipitation from +3h to 180h? I'm trying to draw the
>total precipitation above the bars graph and I tried something like:
>
>totalprec=sum(APCPsfc,t=2,t=180)
>'draw string x y 'totalprec
>
>But it doesn't work at all.
>
>Hope someone can help me.
>
>Cheers.



More information about the gradsusr mailing list