Hourly precipitation plots for hourly WRF output

Hi,

I'm a new NCL user and am having some trouble figuring out the best way
to generate plots of hourly precipitation from my hourly wrfout_* files.

There is a sample script online (
http://www.mmm.ucar.edu/wrf/users/graphics/NCL/Examples/Precip/wrf_Precip.ncl)
that almost gets the job done, but I'm having trouble because I have a
list of individual files with hourly output at one time, instead of one
file with several times in it. [As an aside, I also tried using NCO to
concatenate all of the files into one huge file, and running the example
script on it, but I then got errors referencing issues with the
WRFUserARW.ncl file; an example is pasted below.]
I have been trying to use addfiles to generate my list but then have
lots of subsequent issues with functions like wrf_user_list_times,
wrf_map, wrf_user_getvar, etc, that require a single file instead of a
variable of type "list" as is generated with addfiles.

Ideally, I would love to have a script that loops over each time/file in
my model run and generates a separate plot of hourly precipitation for
that time (I can do this for total precipitation easily, as that only
needs the time in the file itself, whereas hourly precip needs to be
calculated by subtracting off the total precip from the time before).

If anyone has a similar script that they wouldn't mind sharing, or any
advice on how this would be most easily done, I would be very appreciative!

Thanks so much,
Kelly Mahoney

---------------------------------------------------------
errors when trying to run sample script with one large, NCO-concatenated
file:
> Coordinates:
> fatal:Number of subscripts do not match number of dimensions of
> variable,(2) Subscripts used, (3) Subscripts expected
> fatal:Execute: Error occurred at or near line 1373 in file
> $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Feb 23 14:52:05 2010