[Met_help] MODE inputs

John Halley Gotway johnhg at rap.ucar.edu
Tue Feb 24 13:59:20 MST 2009


Please read the usage statement for the pcp_combine "-subtract" command a
little closer.

For the "-add" and "-subtract" commands, you're supposed to pass the two
actual filenames you'd like to to use, each followed by the accumulation
interval you'd like to select from the file.  You don't pass anytime
initial or valid time information for the add and subtract commands.

So in your case, you'd run something like:

cd /groups/wrf/firewx/output/22may2007_mp6.cu3.bl1.icbcenkf0/outer/postprd
pcp_combine -subtract WRFPRS_d01.15 15 WRFPRS_d01.14 14

This will subtract: 15 hour accumulation in WRFPRS_d01.15 minus 14 hour
accumulation in WRFPRS_d01.14.

The add and subtract commands don't do anything fancy.  You run it once
for each file operation you'd like to perform and tell it which two files
to use.  So once you set it up the way you'd like, I'd suggest running it
with a script.

Hope that helps.


> Mr. Gotway,
> I am glad a solution already exists!  Unfortunately I get this...
> [soholmbe at u101 pcp_combine]$ ./pcp_combine -subtract 2007-05-22_23:00:00
> 14 2007-05-23_00:00:00 15 wrf_00z23may2007_1hr.nc -pcpdir
> /groups/wrf/firewx/output/22may2007_mp6.cu3.bl1.icbcenkf0/outer/postprd
> Reading input file: 2007-05-22_23:00:00
>   GribFile::open(const char *) -> unable to open grib file
> 2007-05-22_23:00:00
> [soholmbe at u101 pcp_combine]$
> My wrfout is hourly from 09z to 06z, I have 22 separate wrfout files.  I
> have WRFPRS_d01.00 through WRFPRS_d01.21 from the WPP in this location but
> they are not being read.  Have you seen this error before?  Here are the
> top two lines of the WPP file for the May 23 at 00z.
> [soholmbe at u101 postprd]$ wgrib WRFPRS_d01.15
> 1:0:d=07052209:PRMSL:kpds5=2:kpds6=102:kpds7=0:TR=0:P1=15:P2=0:TimeU=1:MSL:15hr
> fcst:NAve=0
> 2:80896:d=07052209:PRES:kpds5=1:kpds6=109:kpds7=1:TR=0:P1=15:P2=0:TimeU=1:hybrid
> lev 1:15hr fcst:NAve=0
> Does this look right to you, is something wrong with the WPP?  Did I get
> the accum1 and accum2 correct?
> Thank you!
> Shelley Holmberg
> -----Original Message-----
> From: John Halley Gotway [mailto:johnhg at rap.ucar.edu]
> Sent: Thu 2/19/2009 2:50 PM
> To: Holmberg, Shelley; met_help at ucar.edu
> Subject: Re: [Met_help] MODE inputs
> Shelley,
> OK, I understand your question now.  You'll want to use a combination of
> WPP and the PCP-Combine tool in MET.
> Suppose you're running your WRF-ARW model and creating an output file
> every 3 hours out to 24 hours.  So it'll write files every 3 hours that be
> default contain precipitation accumulated up until that
> time.  So the 12 hour file will contain precip accumulated over 12 hours.
> And the 15 hour file will contain precip accumulated over 15 hours.  When
> you run these files through WPP, corresponding GRIB
> files will be created using the same accumulation intervals.
> Then you can run the PCP-Combine tool on those GRIB files using the
> "-subtract" option.  That will allow you to "subtract" the 12 hour
> accumulation from the 15 hour accumulation, yielding a 3 hour
> accumulation.  You can refer to the MET User's Guide for details on
> running the "-subtract" option.  We put in the "-subtract" option for the
> very reason your asking about.
> However, you won't be able to verify 1 hour accumulations unless you're
> writing a WRF-out file every hour.
> Hope that helps.
> John
> Holmberg, Shelley wrote:
>> Mr. Gotway,
>> Sorry if I was not clear before.  My problem is getting the correct
>> accumulation period for the total accumulated precipitation in the
>> wrfout file.  Right now, the precipitation accumulates over 22 hours and
>> I would like to have 1 hour (or 3 hour) accumulations for each of the 22
>> hours.  The variable "Accumulated total precipitation," Grib 61 will
>> work prefectly but is it possible to get the right accumulation period
>> with the WPP?
>> Thank you!
>> Shelley
>> -----Original Message-----
>> From: John Halley Gotway [mailto:johnhg at rap.ucar.edu]
>> Sent: Wed 2/18/2009 6:08 PM
>> To: Holmberg, Shelley
>> Cc: met_help at ucar.edu
>> Subject: Re: [Met_help] MODE inputs
>> Shelley,
>> Let me see if I understand your situation correctly.  You are running
>> the
>> WRF-ARW model and are generating WRF NetCDF output files.  You're
>> running
>> those NetCDF output files through the WRF Post-Processor (WPP) and
>> generating GRIB files.  Then you trying to use the PCP-Combine tool in
>> MET
>> to combine the precipitation fields.
>> Please note that PCP-Combine currently only reads GRIB files.  So you'd
>> need to run your WRF output through WPP prior to running MET.
>> You say that you'd like to combine RAINC and RAINNC from the model.  Now
>> I'm not sure when you run your data through WPP, how does it write those
>> fields out?  Perhaps they're written as large scale precip (NCPCP) and
>> convective precip (ACPCP).
>> And your question is this: Is there a way to get PCP-Combine to add
>> together NCPCP and ACPCP?  The short answer is no, not currently.
>> But if I understand correctly, I'm wondering why don't you just use the
>> accumulated precipitation field, APCP (i.e. GRIB code 61)?  The StageIV
>> analysis, against which you're verifying, is accumulated precip (code
>> 61).
>>  Is there a reason that you'd like to combine RAINC and RAINNC, rather
>> than just using accumulated precip directly?
>> If there is some reason why combining these two precip types makes
>> sense,
>> we could probably add support for it to the PCP-Combine tool.
>> Do I understand the situation correctly?  Does verifying APCP accomplish
>> what you're trying to do?
>> Thanks,
>> John Halley Gotway
>> johnhg at ucar.edu
>>> Hi!
>>> I have been working with the MODE tool within MET off and on for a year
>>> now.  I have to tool properly configured as it will work with the
>>> canned
>>> data and some of my own data.  My main issue is that I am trying to
>>> compare 1 or 3 hour surface accumulations of rainfall from WRF output
>>> (convective and grid scale) and Stage IV observations.  I have
>>> successfully used copygb and pcp_combine to format the Stage IV obs.
>>> but
>>> with the modeled data, I haven't been able to: add the rainc and rainnc
>>> variables; get the proper accumulation period; and then write out a
>>> netcdf
>>> or GRIB file.  Is there a program available that will properly process
>>> my
>>> modeled data?
>>> Thank you!
>>> Shelley
>>> _______________________________________________
>>> Met_help mailing list
>>> Met_help at mailman.ucar.edu
>>> http://mailman.ucar.edu/mailman/listinfo/met_help

More information about the Met_help mailing list