John,<br>I was able to successfully run copygb with your help and now am looking to run pcp_combine, subtract, to get a one hour precip field from my WRF output files. I have isolated the two files (and actually uploaded them to the server described in previous emails), and used the following command to run pcp_combine assuming the following:<br>
<br> "accum1" indicates the accumulation interval to be used from in_file1 in HH format (required). -I set this as 12<br> "accum2" indicates the accumulation interval to be used from in_file2 in HH format (required). -I set this as 11<br>
<br>I got this:<br><br>/bin/pcp_combine -subtract WRFPRS_d02.012_00 12 WRFPRS_d02.011_00 11 difference.out<br>Reading input file: WRFPRS_d02.012_00<br>Reading input file: WRFPRS_d02.011_00<br>Performing subtraction command.<br>
Writing output file: difference.out<br><br>Figuring that pcp_combine took the difference between two GRIB files, the output should also be grip, no? So I tried to convert it to view using Grads, but got an error when I ran <a href="http://grib2ctl.pl">grib2ctl.pl</a>: <br>
"/<a href="http://grib2ctl.pl">grib2ctl.pl</a> difference.out > difference_ctl.ctl<br>found unidentified data <br>Big problem:<br> either difference.out is missing or not a grib file<br> or wgrib is not on your path<br>
<br>Please help. Thank you very much for your help so far.<br><br>Jeff Duda<br><br><div class="gmail_quote">On Fri, Dec 18, 2009 at 3:45 PM, John Halley Gotway <span dir="ltr"><<a href="mailto:johnhg@rap.ucar.edu">johnhg@rap.ucar.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Jeff,<br>
<br>
Well I'd first check to see what's actually in that GRIB file using the<br>
"wgrib" command. For example, try this:<br>
wgrib WRFPRS_d02.012_00 | grep APCP<br>
<br>
That'll dump out the contents of that GRIB file, looking only for the APCP<br>
records (accumulated precip). Look to see if that file really does<br>
contain a record for 0-12 hour accumulated precip.<br>
<br>
If you have more questions about this, it'd probably be easier for you to<br>
just send me some sample data to take a look at. If you'd like to do<br>
that, you could post it to our anonymous ftp site in the following way:<br>
<br>
cd "the directory that contains the files you want to send"<br>
ftp <a href="http://ftp.rap.ucar.edu" target="_blank">ftp.rap.ucar.edu</a><br>
username = anonymous<br>
password = "your email address"<br>
cd incoming/irap/met_help<br>
mkdir duda_data<br>
cd duda_data<br>
put WRFPRS_d02.012_00<br>
put WRFPRS_d02.011_00<br>
bye<br>
<br>
Hope that helps.<br>
<font color="#888888"><br>
John<br>
</font><div><div></div><div class="h5"><br>
> John,<br>
> I am Bill Gallus' student and I am trying to use MODE to verify a WRF<br>
> model<br>
> run using accumulated precip. To do this, I need to match Stage_IV data<br>
> to<br>
> my model grid. To do that, though, I need to get a difference in precip<br>
> file to compare with Stage_IV data (which has 1 hourly accumulated<br>
> precip).<br>
> My WRF files use run total accumulated precip. I have run WPP on the<br>
> NetCDF<br>
> output of the WRF files and now am trying to use pcp_combine -subtract to<br>
> basically subtract the total accumulated precip from one file from one an<br>
> hour before it (my NetCDF output interval is 15 minutes). When I tried to<br>
> run pcp_combine, though, I got an error. Can you help me with this,<br>
> please? Thank you.<br>
><br>
><br>
> Here's a screen grab of the command and error:<br>
><br>
> -bash-3.2$ /usr/local/METv2.0/bin/pcp_combine -subtract WRFPRS_d02.012_00<br>
> 12<br>
> WRFPRS_d02.011_00 11 20080719_18-17_precipdiff<br>
> Reading input file: WRFPRS_d02.012_00<br>
><br>
><br>
> ERROR: get_field() -> can't find grib code 61 with accumulation of 12<br>
> hours<br>
> in grib file: WRFPRS_d02.012_00<br>
><br>
> Jeff Duda<br>
> --<br>
> Jeff Duda<br>
> Iowa State University<br>
> Meteorology Graduate Student<br>
> 3134 Agronomy Hall<br>
> <a href="http://www.meteor.iastate.edu/%7Ejdduda" target="_blank">www.meteor.iastate.edu/~jdduda</a><br>
><br>
<br>
<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Jeff Duda<br>Iowa State University<br>Meteorology Graduate Student<br>3134 Agronomy Hall<br><a href="http://www.meteor.iastate.edu/~jdduda">www.meteor.iastate.edu/~jdduda</a><br>