[Met_help] [rt.rap.ucar.edu #49376] History for Request for preparing gridded observations in grib1 format
Tressa Fowler via RT
met_help at ucar.edu
Wed Oct 19 10:57:16 MDT 2011
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
MET-Help,
Thank you for your continuous help,
I am running Grid stat (for rainfall verification). My ARW model covers
IHOP region with Lambert conformal projection. My output is in 1hour
interval. and I want to verify the 3-hr accumulation rainfall.
My initial condition is 20020602_00 and forecast up to 72 hours.
I am verifying the rainfall accumulation of 49-51 hour forecast.
For that i am using PCP_COMBINE script. But it is giving error as follows:
The command I used is:
*#!/bin/sh
echo
echo "*** Running PCP-Combine to combine 1-hourly APCP accumulation
forecasts into a single 3 hour accumulation forecast ***"
../bin/pcp_combine 20020604_010000 1 20020604_030000 3 ../jam.nc -pcpdir
../model*
Please help what may be wrong and how to overcome....
Thanks,
Krishna.
***********************************************************************************************************
/home/ba01/u135/kosuri/met/met3/jam>./jam_pcp_combine.sh
*** Running PCP-Combine to combine 1-hourly APCP accumulation forecasts into
a single 3 hour accumulation forecast ***
Performing sum command: Init/In_Accum/Valid/Out_Accum Times =
20020604_010000/010000/20020604_030000/030000
ERROR: do_sum_command() -> The output accumulation time (030000) cannot be
greater than the lead time (020000).
***********************************************************************************************************
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #49376] Request for preparing gridded observations in grib1 format
From: John Halley Gotway
Time: Thu Sep 01 15:00:26 2011
Krishna,
Sorry for the delay in getting back to you. In the future, please
send questions about MET directly to MET-Help. I went ahead and
created a MET-Help ticket for your question.
It looks like you're trying sum up 3 1-hourly accumulations into a
single 3-hourly accumulation using the "-sum" option of pcp_combine.
The short answer to your question is to try running the
following command instead:
../bin/pcp_combine 20020604_000000 1 20020604_030000 3 ../jam.nc
-pcpdir ../model*
All I've done here is change the first date from 20020604_010000 to
20020604_000000. If you look at the usage statement for pcp_combine,
you'll see that the arguments for the "-sum" option are:
SUM_ARGS:
init_time
in_accum
valid_time
out_accum
out_file
[-pcpdir path]
[-pcprx reg_exp]
where "init_time" indicates the initialization time
of the input GRIB files in YYYYMMDD[_HH[MMSS]] format (required).
"in_accum" indicates the accumulation interval
of the input GRIB files in HH[MMSS] format (required).
"valid_time" indicates the desired valid time
in YYYYMMDD[_HH[MMSS]] format (required).
"out_accum" indicates the desired accumulation
interval for the output NetCDF file in HH[MMSS] format (required).
"out_file" indicates the name of the output
NetCDF file to be written consisting of the sum of the accumulation
intervals (required).
"-pcpdir path" overrides the default
precipitation directory (.) (optional).
"-pcprx reg_exp" overrides the default regular
expression for precipitation file naming convention (.*) (optional).
Note: Set init_time to 00000000_000000 when
summing observation files.
The error message you got stated that the "lead_time" can't be less
than the accumulation interval. The "lead_time" is computed as the
"valid_time" minus the "init_time". So in your case, since
you're computing a 3-hour accumulation from 1-hourly files, the
difference between the valid and init times needs to be 3 hours.
Please try that out and let me know if it doesn't work.
Also, another option you have is to run pcp_combine using the "-add"
option, where you explicitly list the file names and accumulation
intervals you want to add together. It would do the same thing
as the "-sum" option in this case, but with different arguments. Feel
free to use whichever you find most convenient. However, there was a
bug discovered using the "-add" option in the last 6
months. So be sure you're using the latest set of posted patches.
Hope that helps. Please let me know if your issues persist.
Thanks,
John Halley Gotway
met_help at ucar.edu
On 09/01/2011 02:50 PM, RAL HelpDesk {for John Halley Gotway} wrote:
>
> Thu Sep 01 14:50:07 2011: Request 49376 was acted upon.
> Transaction: Ticket created by johnhg
> Queue: met_help
> Subject: Request for preparing gridded observations in grib1
format
> Owner: johnhg
> Requestors: osurikishore at gmail.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49376 >
>
>
> MET-Help,
>
> Thank you for your continuous help,
> I am running Grid stat (for rainfall verification). My ARW model
covers
> IHOP region with Lambert conformal projection. My output is in
1hour
> interval. and I want to verify the 3-hr accumulation rainfall.
>
> My initial condition is 20020602_00 and forecast up to 72 hours.
> I am verifying the rainfall accumulation of 49-51 hour forecast.
>
> For that i am using PCP_COMBINE script. But it is giving error as
follows:
> The command I used is:
> *#!/bin/sh
> echo
> echo "*** Running PCP-Combine to combine 1-hourly APCP accumulation
> forecasts into a single 3 hour accumulation forecast ***"
> ../bin/pcp_combine 20020604_010000 1 20020604_030000 3 ../jam.nc
-pcpdir
> ../model*
>
> Please help what may be wrong and how to overcome....
> Thanks,
> Krishna.
>
***********************************************************************************************************
> /home/ba01/u135/kosuri/met/met3/jam>./jam_pcp_combine.sh
>
> *** Running PCP-Combine to combine 1-hourly APCP accumulation
forecasts into
> a single 3 hour accumulation forecast ***
> Performing sum command: Init/In_Accum/Valid/Out_Accum Times =
> 20020604_010000/010000/20020604_030000/030000
>
>
> ERROR: do_sum_command() -> The output accumulation time (030000)
cannot be
> greater than the lead time (020000).
>
***********************************************************************************************************
>
>
------------------------------------------------
Subject: Request for preparing gridded observations in grib1 format
From: krishna kishore
Time: Sun Sep 04 21:14:43 2011
Thanks John,
Still it is not working after I tried with your tips
I will give some more details about my data...
I prepared WRFPRS_d01.000 grib files for ARW model forecast using WPP
post
processing system. Right now, i am using these grib files for
PCP_COMBINE.
Will this step is correct?
And for your information, I am pasting error executable command I used
for
PCP_COMBINE:
***************************************************************
../bin/pcp_combine \
*20020604_000000* 1 20020604_030000 3 \
../jam.nc \
-pcpdir ./model -v 2
***************************************************************
../bin/pcp_combine \
*20020602_480000* 1 20020602_510000 3 \
../jam.nc \
-pcpdir ./model -v 2
***************************************************************
My doubts are 1. How *should i provide the init_time? is it
experiment
initial time or the time from which i want to add. *
*For more clarity for my doubt...I have used 2002-06-02_00 as initial
condition and had forecast up to 72 hours i.e., **
2002-06-05_00..............**I want to accumulate the rainfall between
48
hour and 51 hour forecasts....*
*What time should i give in "init_time" place...Is it
**20020602_480000 or
**20020604_000000*
I used in both ways...But same error as follows
====================================================================
/home/ba01/u135/kosuri/met/met3/jam>./test_pcp_combine.sh
*** Running PCP-Combine to combine 1-hourly APCP accumulation
forecasts into
a single 3 hour accumulation forecast ***
Performing sum command: Init/In_Accum/Valid/Out_Accum Times =
20020604_000000/010000/20020604_030000/030000
Searching for 3 files with accumulation times of 010000 to sum to a
total
accumulation time of 030000.
ERROR: sum_grib_files() -> Cannot find a file with a valid time of
20020604_030000 and accumulation time of 010000 matching the regular
expression ".*"
====================================================================
Could you help me please,
Thanks once again for your continuous help,
Krishna.
On Fri, Sep 2, 2011 at 2:30 AM, RAL HelpDesk {for John Halley Gotway}
<
met_help at ucar.edu> wrote:
> Krishna,
>
> Sorry for the delay in getting back to you. In the future, please
send
> questions about MET directly to MET-Help. I went ahead and created
a
> MET-Help ticket for your question.
>
> It looks like you're trying sum up 3 1-hourly accumulations into a
single
> 3-hourly accumulation using the "-sum" option of pcp_combine. The
short
> answer to your question is to try running the
> following command instead:
>
> ../bin/pcp_combine 20020604_000000 1 20020604_030000 3 ../jam.nc
-pcpdir
> ../model*
>
> All I've done here is change the first date from 20020604_010000 to
> 20020604_000000. If you look at the usage statement for
pcp_combine, you'll
> see that the arguments for the "-sum" option are:
>
> SUM_ARGS:
> init_time
> in_accum
> valid_time
> out_accum
> out_file
> [-pcpdir path]
> [-pcprx reg_exp]
>
> where "init_time" indicates the initialization time
of the
> input GRIB files in YYYYMMDD[_HH[MMSS]] format (required).
> "in_accum" indicates the accumulation
interval of
> the input GRIB files in HH[MMSS] format (required).
> "valid_time" indicates the desired valid time
in
> YYYYMMDD[_HH[MMSS]] format (required).
> "out_accum" indicates the desired
accumulation
> interval for the output NetCDF file in HH[MMSS] format (required).
> "out_file" indicates the name of the output
NetCDF
> file to be written consisting of the sum of the accumulation
intervals
> (required).
> "-pcpdir path" overrides the default
precipitation
> directory (.) (optional).
> "-pcprx reg_exp" overrides the default
regular
> expression for precipitation file naming convention (.*) (optional).
>
> Note: Set init_time to 00000000_000000 when
summing
> observation files.
>
> The error message you got stated that the "lead_time" can't be less
than
> the accumulation interval. The "lead_time" is computed as the
"valid_time"
> minus the "init_time". So in your case, since
> you're computing a 3-hour accumulation from 1-hourly files, the
difference
> between the valid and init times needs to be 3 hours. Please try
that out
> and let me know if it doesn't work.
>
> Also, another option you have is to run pcp_combine using the "-add"
> option, where you explicitly list the file names and accumulation
intervals
> you want to add together. It would do the same thing
> as the "-sum" option in this case, but with different arguments.
Feel free
> to use whichever you find most convenient. However, there was a bug
> discovered using the "-add" option in the last 6
> months. So be sure you're using the latest set of posted patches.
>
> Hope that helps. Please let me know if your issues persist.
>
> Thanks,
> John Halley Gotway
> met_help at ucar.edu
>
>
> On 09/01/2011 02:50 PM, RAL HelpDesk {for John Halley Gotway} wrote:
> >
> > Thu Sep 01 14:50:07 2011: Request 49376 was acted upon.
> > Transaction: Ticket created by johnhg
> > Queue: met_help
> > Subject: Request for preparing gridded observations in grib1
format
> > Owner: johnhg
> > Requestors: osurikishore at gmail.com
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49376 >
> >
> >
> > MET-Help,
> >
> > Thank you for your continuous help,
> > I am running Grid stat (for rainfall verification). My ARW model
covers
> > IHOP region with Lambert conformal projection. My output is in
1hour
> > interval. and I want to verify the 3-hr accumulation rainfall.
> >
> > My initial condition is 20020602_00 and forecast up to 72 hours.
> > I am verifying the rainfall accumulation of 49-51 hour forecast.
> >
> > For that i am using PCP_COMBINE script. But it is giving error as
> follows:
> > The command I used is:
> > *#!/bin/sh
> > echo
> > echo "*** Running PCP-Combine to combine 1-hourly APCP
accumulation
> > forecasts into a single 3 hour accumulation forecast ***"
> > ../bin/pcp_combine 20020604_010000 1 20020604_030000 3 ../jam.nc
-pcpdir
> > ../model*
> >
> > Please help what may be wrong and how to overcome....
> > Thanks,
> > Krishna.
> >
>
***********************************************************************************************************
> > /home/ba01/u135/kosuri/met/met3/jam>./jam_pcp_combine.sh
> >
> > *** Running PCP-Combine to combine 1-hourly APCP accumulation
forecasts
> into
> > a single 3 hour accumulation forecast ***
> > Performing sum command: Init/In_Accum/Valid/Out_Accum Times =
> > 20020604_010000/010000/20020604_030000/030000
> >
> >
> > ERROR: do_sum_command() -> The output accumulation time (030000)
cannot
> be
> > greater than the lead time (020000).
> >
>
***********************************************************************************************************
> >
> >
>
>
--
###############################
O. Krishna Kishore
Project Associate,
Room No: 411, Block - 6,
Centre for Atmospheric Sciences,
IIT Delhi, Hauz Khas - 110016
New Delhi.
###############################
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #49376] Request for preparing gridded observations in grib1 format
From: John Halley Gotway
Time: Tue Sep 06 15:21:45 2011
Krishna,
The "initialization" time is the experiment initialization time - and
NOT the time from which you want to start adding. Each GRIB file
contains 2 times - a model initialization time and a forecast
lead time. For your run that was initialized at 20020604_000000, each
GRIB output file from WPP will have the same model init time,
20020604_000000. However, the forecast times will differ for each
forecast lead time. If you run "wgrib -V" on these files, you'll see
the timing information listed out.
Now, I'm not sure why the first command isn't working. I looks like
it should work:
../bin/pcp_combine \
20020604_000000 1 20020604_030000 3 \
../jam.nc \
-pcpdir ./model -v 2
I'm assuming that your "model" sub-directory really does contain the
GRIB files you want to sum up. You could also try setting the logging
level higher "-v 3".
If you'd like, you could load up some sample GRIB data to our ftp
site, and I could run it through PCP-Combine to see if there's
anything odd going on.
If you'd like to go that route, you can find directions for posting
data to our anonymous FTP site here:
http://www.dtcenter.org/met/users/support/met_help.php#ftp
Thanks,
John Halley Gotway
On 09/04/2011 09:14 PM, RAL HelpDesk {for krishna kishore} wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49376 >
>
> Thanks John,
> Still it is not working after I tried with your tips
> I will give some more details about my data...
> I prepared WRFPRS_d01.000 grib files for ARW model forecast using
WPP post
> processing system. Right now, i am using these grib files for
PCP_COMBINE.
> Will this step is correct?
>
> And for your information, I am pasting error executable command I
used for
> PCP_COMBINE:
> ***************************************************************
> ../bin/pcp_combine \
> *20020604_000000* 1 20020604_030000 3 \
> ../jam.nc \
> -pcpdir ./model -v 2
> ***************************************************************
> ../bin/pcp_combine \
> *20020602_480000* 1 20020602_510000 3 \
> ../jam.nc \
> -pcpdir ./model -v 2
> ***************************************************************
> My doubts are 1. How *should i provide the init_time? is it
experiment
> initial time or the time from which i want to add. *
> *For more clarity for my doubt...I have used 2002-06-02_00 as
initial
> condition and had forecast up to 72 hours i.e., **
> 2002-06-05_00..............**I want to accumulate the rainfall
between 48
> hour and 51 hour forecasts....*
> *What time should i give in "init_time" place...Is it
**20020602_480000 or
> **20020604_000000*
>
> I used in both ways...But same error as follows
> ====================================================================
> /home/ba01/u135/kosuri/met/met3/jam>./test_pcp_combine.sh
>
> *** Running PCP-Combine to combine 1-hourly APCP accumulation
forecasts into
> a single 3 hour accumulation forecast ***
> Performing sum command: Init/In_Accum/Valid/Out_Accum Times =
> 20020604_000000/010000/20020604_030000/030000
> Searching for 3 files with accumulation times of 010000 to sum to a
total
> accumulation time of 030000.
>
>
> ERROR: sum_grib_files() -> Cannot find a file with a valid time of
> 20020604_030000 and accumulation time of 010000 matching the regular
> expression ".*"
> ====================================================================
>
> Could you help me please,
> Thanks once again for your continuous help,
> Krishna.
>
> On Fri, Sep 2, 2011 at 2:30 AM, RAL HelpDesk {for John Halley
Gotway} <
> met_help at ucar.edu> wrote:
>
>> Krishna,
>>
>> Sorry for the delay in getting back to you. In the future, please
send
>> questions about MET directly to MET-Help. I went ahead and created
a
>> MET-Help ticket for your question.
>>
>> It looks like you're trying sum up 3 1-hourly accumulations into a
single
>> 3-hourly accumulation using the "-sum" option of pcp_combine. The
short
>> answer to your question is to try running the
>> following command instead:
>>
>> ../bin/pcp_combine 20020604_000000 1 20020604_030000 3 ../jam.nc
-pcpdir
>> ../model*
>>
>> All I've done here is change the first date from 20020604_010000 to
>> 20020604_000000. If you look at the usage statement for
pcp_combine, you'll
>> see that the arguments for the "-sum" option are:
>>
>> SUM_ARGS:
>> init_time
>> in_accum
>> valid_time
>> out_accum
>> out_file
>> [-pcpdir path]
>> [-pcprx reg_exp]
>>
>> where "init_time" indicates the initialization
time of the
>> input GRIB files in YYYYMMDD[_HH[MMSS]] format (required).
>> "in_accum" indicates the accumulation
interval of
>> the input GRIB files in HH[MMSS] format (required).
>> "valid_time" indicates the desired valid
time in
>> YYYYMMDD[_HH[MMSS]] format (required).
>> "out_accum" indicates the desired
accumulation
>> interval for the output NetCDF file in HH[MMSS] format (required).
>> "out_file" indicates the name of the output
NetCDF
>> file to be written consisting of the sum of the accumulation
intervals
>> (required).
>> "-pcpdir path" overrides the default
precipitation
>> directory (.) (optional).
>> "-pcprx reg_exp" overrides the default
regular
>> expression for precipitation file naming convention (.*)
(optional).
>>
>> Note: Set init_time to 00000000_000000 when
summing
>> observation files.
>>
>> The error message you got stated that the "lead_time" can't be less
than
>> the accumulation interval. The "lead_time" is computed as the
"valid_time"
>> minus the "init_time". So in your case, since
>> you're computing a 3-hour accumulation from 1-hourly files, the
difference
>> between the valid and init times needs to be 3 hours. Please try
that out
>> and let me know if it doesn't work.
>>
>> Also, another option you have is to run pcp_combine using the "-
add"
>> option, where you explicitly list the file names and accumulation
intervals
>> you want to add together. It would do the same thing
>> as the "-sum" option in this case, but with different arguments.
Feel free
>> to use whichever you find most convenient. However, there was a
bug
>> discovered using the "-add" option in the last 6
>> months. So be sure you're using the latest set of posted patches.
>>
>> Hope that helps. Please let me know if your issues persist.
>>
>> Thanks,
>> John Halley Gotway
>> met_help at ucar.edu
>>
>>
>> On 09/01/2011 02:50 PM, RAL HelpDesk {for John Halley Gotway}
wrote:
>>>
>>> Thu Sep 01 14:50:07 2011: Request 49376 was acted upon.
>>> Transaction: Ticket created by johnhg
>>> Queue: met_help
>>> Subject: Request for preparing gridded observations in grib1
format
>>> Owner: johnhg
>>> Requestors: osurikishore at gmail.com
>>> Status: new
>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49376 >
>>>
>>>
>>> MET-Help,
>>>
>>> Thank you for your continuous help,
>>> I am running Grid stat (for rainfall verification). My ARW model
covers
>>> IHOP region with Lambert conformal projection. My output is in
1hour
>>> interval. and I want to verify the 3-hr accumulation rainfall.
>>>
>>> My initial condition is 20020602_00 and forecast up to 72 hours.
>>> I am verifying the rainfall accumulation of 49-51 hour forecast.
>>>
>>> For that i am using PCP_COMBINE script. But it is giving error as
>> follows:
>>> The command I used is:
>>> *#!/bin/sh
>>> echo
>>> echo "*** Running PCP-Combine to combine 1-hourly APCP
accumulation
>>> forecasts into a single 3 hour accumulation forecast ***"
>>> ../bin/pcp_combine 20020604_010000 1 20020604_030000 3 ../jam.nc
-pcpdir
>>> ../model*
>>>
>>> Please help what may be wrong and how to overcome....
>>> Thanks,
>>> Krishna.
>>>
>>
***********************************************************************************************************
>>> /home/ba01/u135/kosuri/met/met3/jam>./jam_pcp_combine.sh
>>>
>>> *** Running PCP-Combine to combine 1-hourly APCP accumulation
forecasts
>> into
>>> a single 3 hour accumulation forecast ***
>>> Performing sum command: Init/In_Accum/Valid/Out_Accum Times =
>>> 20020604_010000/010000/20020604_030000/030000
>>>
>>>
>>> ERROR: do_sum_command() -> The output accumulation time (030000)
cannot
>> be
>>> greater than the lead time (020000).
>>>
>>
***********************************************************************************************************
>>>
>>>
>>
>>
>
>
------------------------------------------------
Subject: Request for preparing gridded observations in grib1 format
From: krishna kishore
Time: Fri Sep 09 11:38:59 2011
Hi John,
Thank you for your mail.
I tried as you told...But still, it is not working.
So, I uploaded to your MET_help directory with my name "krishna_data"
Please help me in this regard and give details on the error
With best regards,
Krishna.
On Tue, Sep 6, 2011 at 5:21 PM, RAL HelpDesk {for John Halley Gotway}
<
met_help at ucar.edu> wrote:
> Krishna,
>
> The "initialization" time is the experiment initialization time -
and NOT
> the time from which you want to start adding. Each GRIB file
contains 2
> times - a model initialization time and a forecast
> lead time. For your run that was initialized at 20020604_000000,
each GRIB
> output file from WPP will have the same model init time,
20020604_000000.
> However, the forecast times will differ for each
> forecast lead time. If you run "wgrib -V" on these files, you'll
see the
> timing information listed out.
>
> Now, I'm not sure why the first command isn't working. I looks like
it
> should work:
>
> ../bin/pcp_combine \
> 20020604_000000 1 20020604_030000 3 \
> ../jam.nc \
> -pcpdir ./model -v 2
>
> I'm assuming that your "model" sub-directory really does contain the
GRIB
> files you want to sum up. You could also try setting the logging
level
> higher "-v 3".
>
> If you'd like, you could load up some sample GRIB data to our ftp
site, and
> I could run it through PCP-Combine to see if there's anything odd
going on.
>
> If you'd like to go that route, you can find directions for posting
data to
> our anonymous FTP site here:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> Thanks,
> John Halley Gotway
>
>
> On 09/04/2011 09:14 PM, RAL HelpDesk {for krishna kishore} wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49376 >
> >
> > Thanks John,
> > Still it is not working after I tried with your tips
> > I will give some more details about my data...
> > I prepared WRFPRS_d01.000 grib files for ARW model forecast using
WPP
> post
> > processing system. Right now, i am using these grib files for
> PCP_COMBINE.
> > Will this step is correct?
> >
> > And for your information, I am pasting error executable command I
used
> for
> > PCP_COMBINE:
> > ***************************************************************
> > ../bin/pcp_combine \
> > *20020604_000000* 1 20020604_030000 3 \
> > ../jam.nc \
> > -pcpdir ./model -v 2
> > ***************************************************************
> > ../bin/pcp_combine \
> > *20020602_480000* 1 20020602_510000 3 \
> > ../jam.nc \
> > -pcpdir ./model -v 2
> > ***************************************************************
> > My doubts are 1. How *should i provide the init_time? is it
> experiment
> > initial time or the time from which i want to add. *
> > *For more clarity for my doubt...I have used 2002-06-02_00 as
initial
> > condition and had forecast up to 72 hours i.e., **
> > 2002-06-05_00..............**I want to accumulate the rainfall
between 48
> > hour and 51 hour forecasts....*
> > *What time should i give in "init_time" place...Is it
**20020602_480000
> or
> > **20020604_000000*
> >
> > I used in both ways...But same error as follows
> >
====================================================================
> > /home/ba01/u135/kosuri/met/met3/jam>./test_pcp_combine.sh
> >
> > *** Running PCP-Combine to combine 1-hourly APCP accumulation
forecasts
> into
> > a single 3 hour accumulation forecast ***
> > Performing sum command: Init/In_Accum/Valid/Out_Accum Times =
> > 20020604_000000/010000/20020604_030000/030000
> > Searching for 3 files with accumulation times of 010000 to sum to
a total
> > accumulation time of 030000.
> >
> >
> > ERROR: sum_grib_files() -> Cannot find a file with a valid time of
> > 20020604_030000 and accumulation time of 010000 matching the
regular
> > expression ".*"
> >
====================================================================
> >
> > Could you help me please,
> > Thanks once again for your continuous help,
> > Krishna.
> >
> > On Fri, Sep 2, 2011 at 2:30 AM, RAL HelpDesk {for John Halley
Gotway} <
> > met_help at ucar.edu> wrote:
> >
> >> Krishna,
> >>
> >> Sorry for the delay in getting back to you. In the future,
please send
> >> questions about MET directly to MET-Help. I went ahead and
created a
> >> MET-Help ticket for your question.
> >>
> >> It looks like you're trying sum up 3 1-hourly accumulations into
a
> single
> >> 3-hourly accumulation using the "-sum" option of pcp_combine.
The short
> >> answer to your question is to try running the
> >> following command instead:
> >>
> >> ../bin/pcp_combine 20020604_000000 1 20020604_030000 3 ../jam.nc-
pcpdir
> >> ../model*
> >>
> >> All I've done here is change the first date from 20020604_010000
to
> >> 20020604_000000. If you look at the usage statement for
pcp_combine,
> you'll
> >> see that the arguments for the "-sum" option are:
> >>
> >> SUM_ARGS:
> >> init_time
> >> in_accum
> >> valid_time
> >> out_accum
> >> out_file
> >> [-pcpdir path]
> >> [-pcprx reg_exp]
> >>
> >> where "init_time" indicates the initialization
time of
> the
> >> input GRIB files in YYYYMMDD[_HH[MMSS]] format (required).
> >> "in_accum" indicates the accumulation
interval of
> >> the input GRIB files in HH[MMSS] format (required).
> >> "valid_time" indicates the desired valid
time in
> >> YYYYMMDD[_HH[MMSS]] format (required).
> >> "out_accum" indicates the desired
accumulation
> >> interval for the output NetCDF file in HH[MMSS] format
(required).
> >> "out_file" indicates the name of the
output
> NetCDF
> >> file to be written consisting of the sum of the accumulation
intervals
> >> (required).
> >> "-pcpdir path" overrides the default
> precipitation
> >> directory (.) (optional).
> >> "-pcprx reg_exp" overrides the default
regular
> >> expression for precipitation file naming convention (.*)
(optional).
> >>
> >> Note: Set init_time to 00000000_000000
when
> summing
> >> observation files.
> >>
> >> The error message you got stated that the "lead_time" can't be
less than
> >> the accumulation interval. The "lead_time" is computed as the
> "valid_time"
> >> minus the "init_time". So in your case, since
> >> you're computing a 3-hour accumulation from 1-hourly files, the
> difference
> >> between the valid and init times needs to be 3 hours. Please try
that
> out
> >> and let me know if it doesn't work.
> >>
> >> Also, another option you have is to run pcp_combine using the "-
add"
> >> option, where you explicitly list the file names and accumulation
> intervals
> >> you want to add together. It would do the same thing
> >> as the "-sum" option in this case, but with different arguments.
Feel
> free
> >> to use whichever you find most convenient. However, there was a
bug
> >> discovered using the "-add" option in the last 6
> >> months. So be sure you're using the latest set of posted
patches.
> >>
> >> Hope that helps. Please let me know if your issues persist.
> >>
> >> Thanks,
> >> John Halley Gotway
> >> met_help at ucar.edu
> >>
> >>
> >> On 09/01/2011 02:50 PM, RAL HelpDesk {for John Halley Gotway}
wrote:
> >>>
> >>> Thu Sep 01 14:50:07 2011: Request 49376 was acted upon.
> >>> Transaction: Ticket created by johnhg
> >>> Queue: met_help
> >>> Subject: Request for preparing gridded observations in
grib1
> format
> >>> Owner: johnhg
> >>> Requestors: osurikishore at gmail.com
> >>> Status: new
> >>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49376>
> >>>
> >>>
> >>> MET-Help,
> >>>
> >>> Thank you for your continuous help,
> >>> I am running Grid stat (for rainfall verification). My ARW
model
> covers
> >>> IHOP region with Lambert conformal projection. My output is in
1hour
> >>> interval. and I want to verify the 3-hr accumulation rainfall.
> >>>
> >>> My initial condition is 20020602_00 and forecast up to 72 hours.
> >>> I am verifying the rainfall accumulation of 49-51 hour forecast.
> >>>
> >>> For that i am using PCP_COMBINE script. But it is giving error
as
> >> follows:
> >>> The command I used is:
> >>> *#!/bin/sh
> >>> echo
> >>> echo "*** Running PCP-Combine to combine 1-hourly APCP
accumulation
> >>> forecasts into a single 3 hour accumulation forecast ***"
> >>> ../bin/pcp_combine 20020604_010000 1 20020604_030000 3
../jam.nc-pcpdir
> >>> ../model*
> >>>
> >>> Please help what may be wrong and how to overcome....
> >>> Thanks,
> >>> Krishna.
> >>>
> >>
>
***********************************************************************************************************
> >>> /home/ba01/u135/kosuri/met/met3/jam>./jam_pcp_combine.sh
> >>>
> >>> *** Running PCP-Combine to combine 1-hourly APCP accumulation
forecasts
> >> into
> >>> a single 3 hour accumulation forecast ***
> >>> Performing sum command: Init/In_Accum/Valid/Out_Accum Times =
> >>> 20020604_010000/010000/20020604_030000/030000
> >>>
> >>>
> >>> ERROR: do_sum_command() -> The output accumulation time (030000)
cannot
> >> be
> >>> greater than the lead time (020000).
> >>>
> >>
>
***********************************************************************************************************
> >>>
> >>>
> >>
> >>
> >
> >
>
>
--
###############################
O. Krishna Kishore
Project Associate,
Room No: 411, Block - 6,
Centre for Atmospheric Sciences,
IIT Delhi, Hauz Khas - 110016
New Delhi.
###############################
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #49376] Request for preparing gridded observations in grib1 format
From: John Halley Gotway
Time: Fri Sep 09 13:29:00 2011
Krishna,
Good news. I see the problem and it's easy to fix. You were assuming
that your GRIB files contained 1-hour accumulated precipitation, but
they do not. Instead, they contain a runtime accumulation -
meaning that the output file for forecast hour "n" contains
precipitation from hour 0 through hour n. At the end of this message,
I've listed out the wgrib output that demonstrates that you have
runtime accumulations rather than hourly ones. This is actually the
default when running the WRF-ARW core, so it's not a problem.
You can easily run PCP-Combine with different arguments to get the
accumulations you want. You'll just use the "-subtract" option. For
example, to get the 3-hour accumulation between hours 48 and
51, you'd run the following command:
pcp_combine -subtract WRFPRS_d01.051 51 WRFPRS_d01.048 48
WRFPRS_d01_48_51.nc
This subtracts the 48 accumulation in the second file from the 51 hour
accumulation in the first one. The result is a 3-hour accumulation
between hours 48 and 51. You can look at the usage statement
for pcp_combine to see more details about the "-subtract" option.
By the way, if it's important to you for some reason, you can
reconfigure your WRF run to dump out hourly accumulations if you'd
like. Take a look in the WRF User's Guide for more details, but I
believe you'd just set "prec_acc_dt" to 60 in the namelist file. That
tells WRF to "dump" the precipitation accumulation bucket every 60
minutes. So the resulting accumulations would be 1-hour. See
page 5-56 of the WRF-ARW user's guide for more info
(http://www.mmm.ucar.edu/wrf/users/docs/user_guide_V3/ARWUsersGuideV3.pdf).
If you were to do this, you'd revert back to using the "-sum" option
for PCP-Combine.
Hope that gets you going.
Thanks,
John
wgrib output:
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=0:TimeU=1:sfc:0-
0hr acc:NAve=0
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=1:TimeU=1:sfc:0-
1hr acc:NAve=0
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=2:TimeU=1:sfc:0-
2hr acc:NAve=0
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:TimeU=1:sfc:0-
3hr acc:NAve=0
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=4:TimeU=1:sfc:0-
4hr acc:NAve=0
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=5:TimeU=1:sfc:0-
5hr acc:NAve=0
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=6:TimeU=1:sfc:0-
6hr acc:NAve=0
103:8192888:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=7:TimeU=1:sfc:0-
7hr acc:NAve=0
103:8210588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=8:TimeU=1:sfc:0-
8hr acc:NAve=0
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=9:TimeU=1:sfc:0-
9hr acc:NAve=0
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=10:TimeU=1:sfc:0-
10hr acc:NAve=0
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=11:TimeU=1:sfc:0-
11hr acc:NAve=0
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=12:TimeU=1:sfc:0-
12hr acc:NAve=0
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=13:TimeU=1:sfc:0-
13hr acc:NAve=0
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=14:TimeU=1:sfc:0-
14hr acc:NAve=0
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=15:TimeU=1:sfc:0-
15hr acc:NAve=0
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=16:TimeU=1:sfc:0-
16hr acc:NAve=0
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=17:TimeU=1:sfc:0-
17hr acc:NAve=0
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=18:TimeU=1:sfc:0-
18hr acc:NAve=0
103:8080788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=19:TimeU=1:sfc:0-
19hr acc:NAve=0
103:8098488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=20:TimeU=1:sfc:0-
20hr acc:NAve=0
103:8098488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=21:TimeU=1:sfc:0-
21hr acc:NAve=0
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=22:TimeU=1:sfc:0-
22hr acc:NAve=0
103:8092588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=23:TimeU=1:sfc:0-
23hr acc:NAve=0
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=24:TimeU=1:sfc:0-
24hr acc:NAve=0
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=25:TimeU=1:sfc:0-
25hr acc:NAve=0
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=26:TimeU=1:sfc:0-
26hr acc:NAve=0
103:8086688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=27:TimeU=1:sfc:0-
27hr acc:NAve=0
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=28:TimeU=1:sfc:0-
28hr acc:NAve=0
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=29:TimeU=1:sfc:0-
29hr acc:NAve=0
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=30:TimeU=1:sfc:0-
30hr acc:NAve=0
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=31:TimeU=1:sfc:0-
31hr acc:NAve=0
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=32:TimeU=1:sfc:0-
32hr acc:NAve=0
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=33:TimeU=1:sfc:0-
33hr acc:NAve=0
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=34:TimeU=1:sfc:0-
34hr acc:NAve=0
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=35:TimeU=1:sfc:0-
35hr acc:NAve=0
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=36:TimeU=1:sfc:0-
36hr acc:NAve=0
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=37:TimeU=1:sfc:0-
37hr acc:NAve=0
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=38:TimeU=1:sfc:0-
38hr acc:NAve=0
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=39:TimeU=1:sfc:0-
39hr acc:NAve=0
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=40:TimeU=1:sfc:0-
40hr acc:NAve=0
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=41:TimeU=1:sfc:0-
41hr acc:NAve=0
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=42:TimeU=1:sfc:0-
42hr acc:NAve=0
103:8092588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=43:TimeU=1:sfc:0-
43hr acc:NAve=0
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=44:TimeU=1:sfc:0-
44hr acc:NAve=0
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=45:TimeU=1:sfc:0-
45hr acc:NAve=0
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=46:TimeU=1:sfc:0-
46hr acc:NAve=0
103:8122088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=47:TimeU=1:sfc:0-
47hr acc:NAve=0
103:8133888:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=48:TimeU=1:sfc:0-
48hr acc:NAve=0
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=49:TimeU=1:sfc:0-
49hr acc:NAve=0
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=50:TimeU=1:sfc:0-
50hr acc:NAve=0
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=51:TimeU=1:sfc:0-
51hr acc:NAve=0
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=52:TimeU=1:sfc:0-
52hr acc:NAve=0
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=53:TimeU=1:sfc:0-
53hr acc:NAve=0
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=54:TimeU=1:sfc:0-
54hr acc:NAve=0
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=55:TimeU=1:sfc:0-
55hr acc:NAve=0
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=56:TimeU=1:sfc:0-
56hr acc:NAve=0
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=57:TimeU=1:sfc:0-
57hr acc:NAve=0
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=58:TimeU=1:sfc:0-
58hr acc:NAve=0
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=59:TimeU=1:sfc:0-
59hr acc:NAve=0
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=60:TimeU=1:sfc:0-
60hr acc:NAve=0
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=61:TimeU=1:sfc:0-
61hr acc:NAve=0
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=62:TimeU=1:sfc:0-
62hr acc:NAve=0
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=63:TimeU=1:sfc:0-
63hr acc:NAve=0
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=64:TimeU=1:sfc:0-
64hr acc:NAve=0
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=65:TimeU=1:sfc:0-
65hr acc:NAve=0
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=66:TimeU=1:sfc:0-
66hr acc:NAve=0
On 09/09/2011 11:38 AM, RAL HelpDesk {for krishna kishore} wrote:
> Please help me in this regard and give details on the error
------------------------------------------------
Subject: Request for preparing gridded observations in grib1 format
From: krishna kishore
Time: Mon Sep 12 11:12:32 2011
Dear John,
Thank you very much for your help,
I tried with "subtract" option for WRF output and combine option for
OBSERVATIONS.
I was successful.
When running the Grid_stat, it is giving the grid mismatch error.
Then we
came to know that the dimension of WRF data and observations are not
in the
same grid.
But I dont know how to bring those on the same grid.
I uploaded all the data along with the wrf and observation dimension
file in
the ftp site */incoming/irap/met_help/krishna_data*.
I also uploaded the test_pcp_combine file and Gridstat_config file
also.
Can you please help to overcome this..
Krishna.
On Fri, Sep 9, 2011 at 3:29 PM, RAL HelpDesk {for John Halley Gotway}
<
met_help at ucar.edu> wrote:
> Krishna,
>
> Good news. I see the problem and it's easy to fix. You were
assuming that
> your GRIB files contained 1-hour accumulated precipitation, but they
do not.
> Instead, they contain a runtime accumulation -
> meaning that the output file for forecast hour "n" contains
precipitation
> from hour 0 through hour n. At the end of this message, I've listed
out the
> wgrib output that demonstrates that you have
> runtime accumulations rather than hourly ones. This is actually the
> default when running the WRF-ARW core, so it's not a problem.
>
> You can easily run PCP-Combine with different arguments to get the
> accumulations you want. You'll just use the "-subtract" option.
For
> example, to get the 3-hour accumulation between hours 48 and
> 51, you'd run the following command:
>
> pcp_combine -subtract WRFPRS_d01.051 51 WRFPRS_d01.048 48
> WRFPRS_d01_48_51.nc
>
> This subtracts the 48 accumulation in the second file from the 51
hour
> accumulation in the first one. The result is a 3-hour accumulation
between
> hours 48 and 51. You can look at the usage statement
> for pcp_combine to see more details about the "-subtract" option.
>
> By the way, if it's important to you for some reason, you can
reconfigure
> your WRF run to dump out hourly accumulations if you'd like. Take a
look in
> the WRF User's Guide for more details, but I
> believe you'd just set "prec_acc_dt" to 60 in the namelist file.
That
> tells WRF to "dump" the precipitation accumulation bucket every 60
minutes.
> So the resulting accumulations would be 1-hour. See
> page 5-56 of the WRF-ARW user's guide for more info (
>
http://www.mmm.ucar.edu/wrf/users/docs/user_guide_V3/ARWUsersGuideV3.pdf).
> If you were to do this, you'd revert back to using the "-sum"
option
> for PCP-Combine.
>
> Hope that gets you going.
>
> Thanks,
> John
>
> wgrib output:
>
>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=0:TimeU=1:sfc:0-
0hr
> acc:NAve=0
>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=1:TimeU=1:sfc:0-
1hr
> acc:NAve=0
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=2:TimeU=1:sfc:0-
2hr
> acc:NAve=0
>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:TimeU=1:sfc:0-
3hr
> acc:NAve=0
>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=4:TimeU=1:sfc:0-
4hr
> acc:NAve=0
>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=5:TimeU=1:sfc:0-
5hr
> acc:NAve=0
>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=6:TimeU=1:sfc:0-
6hr
> acc:NAve=0
>
103:8192888:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=7:TimeU=1:sfc:0-
7hr
> acc:NAve=0
>
103:8210588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=8:TimeU=1:sfc:0-
8hr
> acc:NAve=0
>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=9:TimeU=1:sfc:0-
9hr
> acc:NAve=0
>
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=10:TimeU=1:sfc:0-
10hr
> acc:NAve=0
>
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=11:TimeU=1:sfc:0-
11hr
> acc:NAve=0
>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=12:TimeU=1:sfc:0-
12hr
> acc:NAve=0
>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=13:TimeU=1:sfc:0-
13hr
> acc:NAve=0
>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=14:TimeU=1:sfc:0-
14hr
> acc:NAve=0
>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=15:TimeU=1:sfc:0-
15hr
> acc:NAve=0
>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=16:TimeU=1:sfc:0-
16hr
> acc:NAve=0
>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=17:TimeU=1:sfc:0-
17hr
> acc:NAve=0
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=18:TimeU=1:sfc:0-
18hr
> acc:NAve=0
>
103:8080788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=19:TimeU=1:sfc:0-
19hr
> acc:NAve=0
>
103:8098488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=20:TimeU=1:sfc:0-
20hr
> acc:NAve=0
>
103:8098488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=21:TimeU=1:sfc:0-
21hr
> acc:NAve=0
>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=22:TimeU=1:sfc:0-
22hr
> acc:NAve=0
>
103:8092588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=23:TimeU=1:sfc:0-
23hr
> acc:NAve=0
>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=24:TimeU=1:sfc:0-
24hr
> acc:NAve=0
>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=25:TimeU=1:sfc:0-
25hr
> acc:NAve=0
>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=26:TimeU=1:sfc:0-
26hr
> acc:NAve=0
>
103:8086688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=27:TimeU=1:sfc:0-
27hr
> acc:NAve=0
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=28:TimeU=1:sfc:0-
28hr
> acc:NAve=0
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=29:TimeU=1:sfc:0-
29hr
> acc:NAve=0
>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=30:TimeU=1:sfc:0-
30hr
> acc:NAve=0
>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=31:TimeU=1:sfc:0-
31hr
> acc:NAve=0
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=32:TimeU=1:sfc:0-
32hr
> acc:NAve=0
>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=33:TimeU=1:sfc:0-
33hr
> acc:NAve=0
>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=34:TimeU=1:sfc:0-
34hr
> acc:NAve=0
>
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=35:TimeU=1:sfc:0-
35hr
> acc:NAve=0
>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=36:TimeU=1:sfc:0-
36hr
> acc:NAve=0
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=37:TimeU=1:sfc:0-
37hr
> acc:NAve=0
>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=38:TimeU=1:sfc:0-
38hr
> acc:NAve=0
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=39:TimeU=1:sfc:0-
39hr
> acc:NAve=0
>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=40:TimeU=1:sfc:0-
40hr
> acc:NAve=0
>
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=41:TimeU=1:sfc:0-
41hr
> acc:NAve=0
>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=42:TimeU=1:sfc:0-
42hr
> acc:NAve=0
>
103:8092588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=43:TimeU=1:sfc:0-
43hr
> acc:NAve=0
>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=44:TimeU=1:sfc:0-
44hr
> acc:NAve=0
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=45:TimeU=1:sfc:0-
45hr
> acc:NAve=0
>
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=46:TimeU=1:sfc:0-
46hr
> acc:NAve=0
>
103:8122088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=47:TimeU=1:sfc:0-
47hr
> acc:NAve=0
>
103:8133888:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=48:TimeU=1:sfc:0-
48hr
> acc:NAve=0
>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=49:TimeU=1:sfc:0-
49hr
> acc:NAve=0
>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=50:TimeU=1:sfc:0-
50hr
> acc:NAve=0
>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=51:TimeU=1:sfc:0-
51hr
> acc:NAve=0
>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=52:TimeU=1:sfc:0-
52hr
> acc:NAve=0
>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=53:TimeU=1:sfc:0-
53hr
> acc:NAve=0
>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=54:TimeU=1:sfc:0-
54hr
> acc:NAve=0
>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=55:TimeU=1:sfc:0-
55hr
> acc:NAve=0
>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=56:TimeU=1:sfc:0-
56hr
> acc:NAve=0
>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=57:TimeU=1:sfc:0-
57hr
> acc:NAve=0
>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=58:TimeU=1:sfc:0-
58hr
> acc:NAve=0
>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=59:TimeU=1:sfc:0-
59hr
> acc:NAve=0
>
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=60:TimeU=1:sfc:0-
60hr
> acc:NAve=0
>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=61:TimeU=1:sfc:0-
61hr
> acc:NAve=0
>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=62:TimeU=1:sfc:0-
62hr
> acc:NAve=0
>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=63:TimeU=1:sfc:0-
63hr
> acc:NAve=0
>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=64:TimeU=1:sfc:0-
64hr
> acc:NAve=0
>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=65:TimeU=1:sfc:0-
65hr
> acc:NAve=0
>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=66:TimeU=1:sfc:0-
66hr
> acc:NAve=0
>
>
> On 09/09/2011 11:38 AM, RAL HelpDesk {for krishna kishore} wrote:
> > Please help me in this regard and give details on the error
>
>
--
###############################
O. Krishna Kishore
Project Associate,
Room No: 411, Block - 6,
Centre for Atmospheric Sciences,
IIT Delhi, Hauz Khas - 110016
New Delhi.
###############################
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #49376] Request for preparing gridded observations in grib1 format
From: John Halley Gotway
Time: Mon Sep 12 11:49:39 2011
Krishna,
Yes, in order to perform gridded verification, you must first put your
forecast and observation files on the same grid. In your case, you
have GRIB files for both your forecast (GRIB output of the
post processor) and observation (GRIB StageIV analysis). Since you
have GRIB data, you can use the copygb utility to regrid it. Copygb
is distributed as part of the WRF-PostProcessor and Unified
PostProcessor. In addition, it's available as a stand-alone package
here:
http://www.cpc.ncep.noaa.gov/products/wesley/copygb.html
Next, you need to decide if you want to regrid your observations onto
your forecast domain or vice-versa. I'll assume that you want to put
the StageIV obs onto your model domain. The next step is to
figure out which command line arguments you should use for copygb.
You can follow the instructions we have posted here for regridding to
a Lambert Conformal projection:
http://www.dtcenter.org/met/users/support/online_tutorial/METv3.0/copygb/run3.php
Here's your forecast grid information:
:Projection = "Lambert Conformal" ;
:scale_lat_1 = "60.000000" ;
:scale_lat_2 = "30.000000" ;
:lat_pin = "28.223000" ;
:lon_pin = "-110.190000" ;
:x_pin = "0.000000" ;
:y_pin = "0.000000" ;
:lon_orient = "-99.450000" ;
:d_km = "9.000000" ;
:r_km = "6371.200000" ;
:nx = "236" ;
:ny = "200 grid_points" ;
Following the instructions on that page, this translates to the
following grid definition:
255 3 236 200 28223 -110190 8 -99450 9000 9000 0 64 60000 30000
Prior to using the StageIV data you should run each file through
copygb as follows:
copygb -xg"255 3 236 200 28223 -110190 8 -99450 9000 9000 0 64
60000 30000" ST4.2002060300.01h ST4.2002060300.01h_regrid
Hopefully that'll get you going.
Thanks,
John
On 09/12/2011 11:12 AM, RAL HelpDesk {for krishna kishore} wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49376 >
>
> Dear John,
> Thank you very much for your help,
> I tried with "subtract" option for WRF output and combine option for
> OBSERVATIONS.
> I was successful.
> When running the Grid_stat, it is giving the grid mismatch error.
Then we
> came to know that the dimension of WRF data and observations are not
in the
> same grid.
> But I dont know how to bring those on the same grid.
>
> I uploaded all the data along with the wrf and observation dimension
file in
> the ftp site */incoming/irap/met_help/krishna_data*.
> I also uploaded the test_pcp_combine file and Gridstat_config file
also.
>
> Can you please help to overcome this..
> Krishna.
>
>
> On Fri, Sep 9, 2011 at 3:29 PM, RAL HelpDesk {for John Halley
Gotway} <
> met_help at ucar.edu> wrote:
>
>> Krishna,
>>
>> Good news. I see the problem and it's easy to fix. You were
assuming that
>> your GRIB files contained 1-hour accumulated precipitation, but
they do not.
>> Instead, they contain a runtime accumulation -
>> meaning that the output file for forecast hour "n" contains
precipitation
>> from hour 0 through hour n. At the end of this message, I've
listed out the
>> wgrib output that demonstrates that you have
>> runtime accumulations rather than hourly ones. This is actually
the
>> default when running the WRF-ARW core, so it's not a problem.
>>
>> You can easily run PCP-Combine with different arguments to get the
>> accumulations you want. You'll just use the "-subtract" option.
For
>> example, to get the 3-hour accumulation between hours 48 and
>> 51, you'd run the following command:
>>
>> pcp_combine -subtract WRFPRS_d01.051 51 WRFPRS_d01.048 48
>> WRFPRS_d01_48_51.nc
>>
>> This subtracts the 48 accumulation in the second file from the 51
hour
>> accumulation in the first one. The result is a 3-hour accumulation
between
>> hours 48 and 51. You can look at the usage statement
>> for pcp_combine to see more details about the "-subtract" option.
>>
>> By the way, if it's important to you for some reason, you can
reconfigure
>> your WRF run to dump out hourly accumulations if you'd like. Take
a look in
>> the WRF User's Guide for more details, but I
>> believe you'd just set "prec_acc_dt" to 60 in the namelist file.
That
>> tells WRF to "dump" the precipitation accumulation bucket every 60
minutes.
>> So the resulting accumulations would be 1-hour. See
>> page 5-56 of the WRF-ARW user's guide for more info (
>>
http://www.mmm.ucar.edu/wrf/users/docs/user_guide_V3/ARWUsersGuideV3.pdf).
>> If you were to do this, you'd revert back to using the "-sum"
option
>> for PCP-Combine.
>>
>> Hope that gets you going.
>>
>> Thanks,
>> John
>>
>> wgrib output:
>>
>>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=0:TimeU=1:sfc:0-
0hr
>> acc:NAve=0
>>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=1:TimeU=1:sfc:0-
1hr
>> acc:NAve=0
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=2:TimeU=1:sfc:0-
2hr
>> acc:NAve=0
>>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:TimeU=1:sfc:0-
3hr
>> acc:NAve=0
>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=4:TimeU=1:sfc:0-
4hr
>> acc:NAve=0
>>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=5:TimeU=1:sfc:0-
5hr
>> acc:NAve=0
>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=6:TimeU=1:sfc:0-
6hr
>> acc:NAve=0
>>
103:8192888:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=7:TimeU=1:sfc:0-
7hr
>> acc:NAve=0
>>
103:8210588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=8:TimeU=1:sfc:0-
8hr
>> acc:NAve=0
>>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=9:TimeU=1:sfc:0-
9hr
>> acc:NAve=0
>>
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=10:TimeU=1:sfc:0-
10hr
>> acc:NAve=0
>>
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=11:TimeU=1:sfc:0-
11hr
>> acc:NAve=0
>>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=12:TimeU=1:sfc:0-
12hr
>> acc:NAve=0
>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=13:TimeU=1:sfc:0-
13hr
>> acc:NAve=0
>>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=14:TimeU=1:sfc:0-
14hr
>> acc:NAve=0
>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=15:TimeU=1:sfc:0-
15hr
>> acc:NAve=0
>>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=16:TimeU=1:sfc:0-
16hr
>> acc:NAve=0
>>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=17:TimeU=1:sfc:0-
17hr
>> acc:NAve=0
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=18:TimeU=1:sfc:0-
18hr
>> acc:NAve=0
>>
103:8080788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=19:TimeU=1:sfc:0-
19hr
>> acc:NAve=0
>>
103:8098488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=20:TimeU=1:sfc:0-
20hr
>> acc:NAve=0
>>
103:8098488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=21:TimeU=1:sfc:0-
21hr
>> acc:NAve=0
>>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=22:TimeU=1:sfc:0-
22hr
>> acc:NAve=0
>>
103:8092588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=23:TimeU=1:sfc:0-
23hr
>> acc:NAve=0
>>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=24:TimeU=1:sfc:0-
24hr
>> acc:NAve=0
>>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=25:TimeU=1:sfc:0-
25hr
>> acc:NAve=0
>>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=26:TimeU=1:sfc:0-
26hr
>> acc:NAve=0
>>
103:8086688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=27:TimeU=1:sfc:0-
27hr
>> acc:NAve=0
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=28:TimeU=1:sfc:0-
28hr
>> acc:NAve=0
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=29:TimeU=1:sfc:0-
29hr
>> acc:NAve=0
>>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=30:TimeU=1:sfc:0-
30hr
>> acc:NAve=0
>>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=31:TimeU=1:sfc:0-
31hr
>> acc:NAve=0
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=32:TimeU=1:sfc:0-
32hr
>> acc:NAve=0
>>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=33:TimeU=1:sfc:0-
33hr
>> acc:NAve=0
>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=34:TimeU=1:sfc:0-
34hr
>> acc:NAve=0
>>
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=35:TimeU=1:sfc:0-
35hr
>> acc:NAve=0
>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=36:TimeU=1:sfc:0-
36hr
>> acc:NAve=0
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=37:TimeU=1:sfc:0-
37hr
>> acc:NAve=0
>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=38:TimeU=1:sfc:0-
38hr
>> acc:NAve=0
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=39:TimeU=1:sfc:0-
39hr
>> acc:NAve=0
>>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=40:TimeU=1:sfc:0-
40hr
>> acc:NAve=0
>>
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=41:TimeU=1:sfc:0-
41hr
>> acc:NAve=0
>>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=42:TimeU=1:sfc:0-
42hr
>> acc:NAve=0
>>
103:8092588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=43:TimeU=1:sfc:0-
43hr
>> acc:NAve=0
>>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=44:TimeU=1:sfc:0-
44hr
>> acc:NAve=0
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=45:TimeU=1:sfc:0-
45hr
>> acc:NAve=0
>>
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=46:TimeU=1:sfc:0-
46hr
>> acc:NAve=0
>>
103:8122088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=47:TimeU=1:sfc:0-
47hr
>> acc:NAve=0
>>
103:8133888:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=48:TimeU=1:sfc:0-
48hr
>> acc:NAve=0
>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=49:TimeU=1:sfc:0-
49hr
>> acc:NAve=0
>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=50:TimeU=1:sfc:0-
50hr
>> acc:NAve=0
>>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=51:TimeU=1:sfc:0-
51hr
>> acc:NAve=0
>>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=52:TimeU=1:sfc:0-
52hr
>> acc:NAve=0
>>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=53:TimeU=1:sfc:0-
53hr
>> acc:NAve=0
>>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=54:TimeU=1:sfc:0-
54hr
>> acc:NAve=0
>>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=55:TimeU=1:sfc:0-
55hr
>> acc:NAve=0
>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=56:TimeU=1:sfc:0-
56hr
>> acc:NAve=0
>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=57:TimeU=1:sfc:0-
57hr
>> acc:NAve=0
>>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=58:TimeU=1:sfc:0-
58hr
>> acc:NAve=0
>>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=59:TimeU=1:sfc:0-
59hr
>> acc:NAve=0
>>
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=60:TimeU=1:sfc:0-
60hr
>> acc:NAve=0
>>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=61:TimeU=1:sfc:0-
61hr
>> acc:NAve=0
>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=62:TimeU=1:sfc:0-
62hr
>> acc:NAve=0
>>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=63:TimeU=1:sfc:0-
63hr
>> acc:NAve=0
>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=64:TimeU=1:sfc:0-
64hr
>> acc:NAve=0
>>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=65:TimeU=1:sfc:0-
65hr
>> acc:NAve=0
>>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=66:TimeU=1:sfc:0-
66hr
>> acc:NAve=0
>>
>>
>> On 09/09/2011 11:38 AM, RAL HelpDesk {for krishna kishore} wrote:
>>> Please help me in this regard and give details on the error
>>
>>
>
>
------------------------------------------------
Subject: Request for preparing gridded observations in grib1 format
From: krishna kishore
Time: Mon Sep 12 16:50:47 2011
Dear John,
Thank you very much for kind and continuous help,
Now, I could able to run the grid stat with your help.
Now our statistics are looking good.
Here, we did match the observations grids (field) onto the model
forecast
grids. I know that, copygb will not interpolate.
My doubt is, whether the observed rainfall data at the grid points
will be
the same before and after the regridding using copygb or not?
And According to you, what would be the best way to do-------1.
whether
observation grids onto model grids or 2. Model grids onto
observations
grids....and What may be the differences in the statistics between
these two
methods.
I think, You understand my question/doubt?
Once again thank you very much for your help.
with best regards,
Krishna.
On Mon, Sep 12, 2011 at 11:19 PM, RAL HelpDesk {for John Halley
Gotway} <
met_help at ucar.edu> wrote:
> Krishna,
>
> Yes, in order to perform gridded verification, you must first put
your
> forecast and observation files on the same grid. In your case, you
have
> GRIB files for both your forecast (GRIB output of the
> post processor) and observation (GRIB StageIV analysis). Since you
have
> GRIB data, you can use the copygb utility to regrid it. Copygb is
> distributed as part of the WRF-PostProcessor and Unified
> PostProcessor. In addition, it's available as a stand-alone package
here:
> http://www.cpc.ncep.noaa.gov/products/wesley/copygb.html
>
> Next, you need to decide if you want to regrid your observations
onto your
> forecast domain or vice-versa. I'll assume that you want to put the
StageIV
> obs onto your model domain. The next step is to
> figure out which command line arguments you should use for copygb.
You can
> follow the instructions we have posted here for regridding to a
Lambert
> Conformal projection:
>
>
http://www.dtcenter.org/met/users/support/online_tutorial/METv3.0/copygb/run3.php
>
> Here's your forecast grid information:
> :Projection = "Lambert Conformal" ;
> :scale_lat_1 = "60.000000" ;
> :scale_lat_2 = "30.000000" ;
> :lat_pin = "28.223000" ;
> :lon_pin = "-110.190000" ;
> :x_pin = "0.000000" ;
> :y_pin = "0.000000" ;
> :lon_orient = "-99.450000" ;
> :d_km = "9.000000" ;
> :r_km = "6371.200000" ;
> :nx = "236" ;
> :ny = "200 grid_points" ;
>
> Following the instructions on that page, this translates to the
following
> grid definition:
> 255 3 236 200 28223 -110190 8 -99450 9000 9000 0 64 60000 30000
>
> Prior to using the StageIV data you should run each file through
copygb as
> follows:
> copygb -xg"255 3 236 200 28223 -110190 8 -99450 9000 9000 0 64
60000
> 30000" ST4.2002060300.01h ST4.2002060300.01h_regrid
>
> Hopefully that'll get you going.
>
> Thanks,
> John
>
>
> On 09/12/2011 11:12 AM, RAL HelpDesk {for krishna kishore} wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49376 >
> >
> > Dear John,
> > Thank you very much for your help,
> > I tried with "subtract" option for WRF output and combine option
for
> > OBSERVATIONS.
> > I was successful.
> > When running the Grid_stat, it is giving the grid mismatch error.
Then
> we
> > came to know that the dimension of WRF data and observations are
not in
> the
> > same grid.
> > But I dont know how to bring those on the same grid.
> >
> > I uploaded all the data along with the wrf and observation
dimension file
> in
> > the ftp site */incoming/irap/met_help/krishna_data*.
> > I also uploaded the test_pcp_combine file and Gridstat_config file
also.
> >
> > Can you please help to overcome this..
> > Krishna.
> >
> >
> > On Fri, Sep 9, 2011 at 3:29 PM, RAL HelpDesk {for John Halley
Gotway} <
> > met_help at ucar.edu> wrote:
> >
> >> Krishna,
> >>
> >> Good news. I see the problem and it's easy to fix. You were
assuming
> that
> >> your GRIB files contained 1-hour accumulated precipitation, but
they do
> not.
> >> Instead, they contain a runtime accumulation -
> >> meaning that the output file for forecast hour "n" contains
> precipitation
> >> from hour 0 through hour n. At the end of this message, I've
listed out
> the
> >> wgrib output that demonstrates that you have
> >> runtime accumulations rather than hourly ones. This is actually
the
> >> default when running the WRF-ARW core, so it's not a problem.
> >>
> >> You can easily run PCP-Combine with different arguments to get
the
> >> accumulations you want. You'll just use the "-subtract" option.
For
> >> example, to get the 3-hour accumulation between hours 48 and
> >> 51, you'd run the following command:
> >>
> >> pcp_combine -subtract WRFPRS_d01.051 51 WRFPRS_d01.048 48
> >> WRFPRS_d01_48_51.nc
> >>
> >> This subtracts the 48 accumulation in the second file from the 51
hour
> >> accumulation in the first one. The result is a 3-hour
accumulation
> between
> >> hours 48 and 51. You can look at the usage statement
> >> for pcp_combine to see more details about the "-subtract" option.
> >>
> >> By the way, if it's important to you for some reason, you can
> reconfigure
> >> your WRF run to dump out hourly accumulations if you'd like.
Take a
> look in
> >> the WRF User's Guide for more details, but I
> >> believe you'd just set "prec_acc_dt" to 60 in the namelist file.
That
> >> tells WRF to "dump" the precipitation accumulation bucket every
60
> minutes.
> >> So the resulting accumulations would be 1-hour. See
> >> page 5-56 of the WRF-ARW user's guide for more info (
> >>
>
http://www.mmm.ucar.edu/wrf/users/docs/user_guide_V3/ARWUsersGuideV3.pdf).
> >> If you were to do this, you'd revert back to using the "-sum"
option
> >> for PCP-Combine.
> >>
> >> Hope that gets you going.
> >>
> >> Thanks,
> >> John
> >>
> >> wgrib output:
> >>
> >>
>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=0:TimeU=1:sfc:0-
0hr
> >> acc:NAve=0
> >>
>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=1:TimeU=1:sfc:0-
1hr
> >> acc:NAve=0
> >>
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=2:TimeU=1:sfc:0-
2hr
> >> acc:NAve=0
> >>
>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:TimeU=1:sfc:0-
3hr
> >> acc:NAve=0
> >>
>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=4:TimeU=1:sfc:0-
4hr
> >> acc:NAve=0
> >>
>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=5:TimeU=1:sfc:0-
5hr
> >> acc:NAve=0
> >>
>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=6:TimeU=1:sfc:0-
6hr
> >> acc:NAve=0
> >>
>
103:8192888:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=7:TimeU=1:sfc:0-
7hr
> >> acc:NAve=0
> >>
>
103:8210588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=8:TimeU=1:sfc:0-
8hr
> >> acc:NAve=0
> >>
>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=9:TimeU=1:sfc:0-
9hr
> >> acc:NAve=0
> >>
>
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=10:TimeU=1:sfc:0-
10hr
> >> acc:NAve=0
> >>
>
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=11:TimeU=1:sfc:0-
11hr
> >> acc:NAve=0
> >>
>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=12:TimeU=1:sfc:0-
12hr
> >> acc:NAve=0
> >>
>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=13:TimeU=1:sfc:0-
13hr
> >> acc:NAve=0
> >>
>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=14:TimeU=1:sfc:0-
14hr
> >> acc:NAve=0
> >>
>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=15:TimeU=1:sfc:0-
15hr
> >> acc:NAve=0
> >>
>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=16:TimeU=1:sfc:0-
16hr
> >> acc:NAve=0
> >>
>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=17:TimeU=1:sfc:0-
17hr
> >> acc:NAve=0
> >>
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=18:TimeU=1:sfc:0-
18hr
> >> acc:NAve=0
> >>
>
103:8080788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=19:TimeU=1:sfc:0-
19hr
> >> acc:NAve=0
> >>
>
103:8098488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=20:TimeU=1:sfc:0-
20hr
> >> acc:NAve=0
> >>
>
103:8098488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=21:TimeU=1:sfc:0-
21hr
> >> acc:NAve=0
> >>
>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=22:TimeU=1:sfc:0-
22hr
> >> acc:NAve=0
> >>
>
103:8092588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=23:TimeU=1:sfc:0-
23hr
> >> acc:NAve=0
> >>
>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=24:TimeU=1:sfc:0-
24hr
> >> acc:NAve=0
> >>
>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=25:TimeU=1:sfc:0-
25hr
> >> acc:NAve=0
> >>
>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=26:TimeU=1:sfc:0-
26hr
> >> acc:NAve=0
> >>
>
103:8086688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=27:TimeU=1:sfc:0-
27hr
> >> acc:NAve=0
> >>
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=28:TimeU=1:sfc:0-
28hr
> >> acc:NAve=0
> >>
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=29:TimeU=1:sfc:0-
29hr
> >> acc:NAve=0
> >>
>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=30:TimeU=1:sfc:0-
30hr
> >> acc:NAve=0
> >>
>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=31:TimeU=1:sfc:0-
31hr
> >> acc:NAve=0
> >>
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=32:TimeU=1:sfc:0-
32hr
> >> acc:NAve=0
> >>
>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=33:TimeU=1:sfc:0-
33hr
> >> acc:NAve=0
> >>
>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=34:TimeU=1:sfc:0-
34hr
> >> acc:NAve=0
> >>
>
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=35:TimeU=1:sfc:0-
35hr
> >> acc:NAve=0
> >>
>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=36:TimeU=1:sfc:0-
36hr
> >> acc:NAve=0
> >>
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=37:TimeU=1:sfc:0-
37hr
> >> acc:NAve=0
> >>
>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=38:TimeU=1:sfc:0-
38hr
> >> acc:NAve=0
> >>
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=39:TimeU=1:sfc:0-
39hr
> >> acc:NAve=0
> >>
>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=40:TimeU=1:sfc:0-
40hr
> >> acc:NAve=0
> >>
>
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=41:TimeU=1:sfc:0-
41hr
> >> acc:NAve=0
> >>
>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=42:TimeU=1:sfc:0-
42hr
> >> acc:NAve=0
> >>
>
103:8092588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=43:TimeU=1:sfc:0-
43hr
> >> acc:NAve=0
> >>
>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=44:TimeU=1:sfc:0-
44hr
> >> acc:NAve=0
> >>
>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=45:TimeU=1:sfc:0-
45hr
> >> acc:NAve=0
> >>
>
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=46:TimeU=1:sfc:0-
46hr
> >> acc:NAve=0
> >>
>
103:8122088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=47:TimeU=1:sfc:0-
47hr
> >> acc:NAve=0
> >>
>
103:8133888:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=48:TimeU=1:sfc:0-
48hr
> >> acc:NAve=0
> >>
>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=49:TimeU=1:sfc:0-
49hr
> >> acc:NAve=0
> >>
>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=50:TimeU=1:sfc:0-
50hr
> >> acc:NAve=0
> >>
>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=51:TimeU=1:sfc:0-
51hr
> >> acc:NAve=0
> >>
>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=52:TimeU=1:sfc:0-
52hr
> >> acc:NAve=0
> >>
>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=53:TimeU=1:sfc:0-
53hr
> >> acc:NAve=0
> >>
>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=54:TimeU=1:sfc:0-
54hr
> >> acc:NAve=0
> >>
>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=55:TimeU=1:sfc:0-
55hr
> >> acc:NAve=0
> >>
>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=56:TimeU=1:sfc:0-
56hr
> >> acc:NAve=0
> >>
>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=57:TimeU=1:sfc:0-
57hr
> >> acc:NAve=0
> >>
>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=58:TimeU=1:sfc:0-
58hr
> >> acc:NAve=0
> >>
>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=59:TimeU=1:sfc:0-
59hr
> >> acc:NAve=0
> >>
>
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=60:TimeU=1:sfc:0-
60hr
> >> acc:NAve=0
> >>
>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=61:TimeU=1:sfc:0-
61hr
> >> acc:NAve=0
> >>
>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=62:TimeU=1:sfc:0-
62hr
> >> acc:NAve=0
> >>
>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=63:TimeU=1:sfc:0-
63hr
> >> acc:NAve=0
> >>
>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=64:TimeU=1:sfc:0-
64hr
> >> acc:NAve=0
> >>
>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=65:TimeU=1:sfc:0-
65hr
> >> acc:NAve=0
> >>
>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=66:TimeU=1:sfc:0-
66hr
> >> acc:NAve=0
> >>
> >>
> >> On 09/09/2011 11:38 AM, RAL HelpDesk {for krishna kishore} wrote:
> >>> Please help me in this regard and give details on the error
> >>
> >>
> >
> >
>
>
--
###############################
O. Krishna Kishore
Project Associate,
Room No: 411, Block - 6,
Centre for Atmospheric Sciences,
IIT Delhi, Hauz Khas - 110016
New Delhi.
###############################
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #49376] Request for preparing gridded observations in grib1 format
From: John Halley Gotway
Time: Tue Sep 13 09:32:32 2011
Krisha,
I'm glad you got the software working how you'd like. As for what's
the best way scientifically to do the interpolation, I'll defer to a
statistician here. One thing I can tell you is that you
should probably use the "budget interpolation" option for copygb which
makes sense for precip. You'd do this by passing the "-i3" command
line option when you call copygb. More information about all
the copygb options can be found here:
http://www.dtcenter.org/met/users/support/online_tutorial/METv3.0/copygb/copygb.txt
I'll go ahead and reassign this support ticket to a statistician. She
can give you more information about whether to interpolate to your
model or the observation domain.
John
On 09/12/2011 04:50 PM, RAL HelpDesk {for krishna kishore} wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49376 >
>
> Dear John,
> Thank you very much for kind and continuous help,
> Now, I could able to run the grid stat with your help.
> Now our statistics are looking good.
> Here, we did match the observations grids (field) onto the model
forecast
> grids. I know that, copygb will not interpolate.
>
> My doubt is, whether the observed rainfall data at the grid points
will be
> the same before and after the regridding using copygb or not?
>
> And According to you, what would be the best way to do-------1.
whether
> observation grids onto model grids or 2. Model grids onto
observations
> grids....and What may be the differences in the statistics between
these two
> methods.
>
> I think, You understand my question/doubt?
>
> Once again thank you very much for your help.
> with best regards,
> Krishna.
>
>
> On Mon, Sep 12, 2011 at 11:19 PM, RAL HelpDesk {for John Halley
Gotway} <
> met_help at ucar.edu> wrote:
>
>> Krishna,
>>
>> Yes, in order to perform gridded verification, you must first put
your
>> forecast and observation files on the same grid. In your case, you
have
>> GRIB files for both your forecast (GRIB output of the
>> post processor) and observation (GRIB StageIV analysis). Since you
have
>> GRIB data, you can use the copygb utility to regrid it. Copygb is
>> distributed as part of the WRF-PostProcessor and Unified
>> PostProcessor. In addition, it's available as a stand-alone
package here:
>> http://www.cpc.ncep.noaa.gov/products/wesley/copygb.html
>>
>> Next, you need to decide if you want to regrid your observations
onto your
>> forecast domain or vice-versa. I'll assume that you want to put
the StageIV
>> obs onto your model domain. The next step is to
>> figure out which command line arguments you should use for copygb.
You can
>> follow the instructions we have posted here for regridding to a
Lambert
>> Conformal projection:
>>
>>
http://www.dtcenter.org/met/users/support/online_tutorial/METv3.0/copygb/run3.php
>>
>> Here's your forecast grid information:
>> :Projection = "Lambert Conformal" ;
>> :scale_lat_1 = "60.000000" ;
>> :scale_lat_2 = "30.000000" ;
>> :lat_pin = "28.223000" ;
>> :lon_pin = "-110.190000" ;
>> :x_pin = "0.000000" ;
>> :y_pin = "0.000000" ;
>> :lon_orient = "-99.450000" ;
>> :d_km = "9.000000" ;
>> :r_km = "6371.200000" ;
>> :nx = "236" ;
>> :ny = "200 grid_points" ;
>>
>> Following the instructions on that page, this translates to the
following
>> grid definition:
>> 255 3 236 200 28223 -110190 8 -99450 9000 9000 0 64 60000 30000
>>
>> Prior to using the StageIV data you should run each file through
copygb as
>> follows:
>> copygb -xg"255 3 236 200 28223 -110190 8 -99450 9000 9000 0 64
60000
>> 30000" ST4.2002060300.01h ST4.2002060300.01h_regrid
>>
>> Hopefully that'll get you going.
>>
>> Thanks,
>> John
>>
>>
>> On 09/12/2011 11:12 AM, RAL HelpDesk {for krishna kishore} wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49376 >
>>>
>>> Dear John,
>>> Thank you very much for your help,
>>> I tried with "subtract" option for WRF output and combine option
for
>>> OBSERVATIONS.
>>> I was successful.
>>> When running the Grid_stat, it is giving the grid mismatch error.
Then
>> we
>>> came to know that the dimension of WRF data and observations are
not in
>> the
>>> same grid.
>>> But I dont know how to bring those on the same grid.
>>>
>>> I uploaded all the data along with the wrf and observation
dimension file
>> in
>>> the ftp site */incoming/irap/met_help/krishna_data*.
>>> I also uploaded the test_pcp_combine file and Gridstat_config file
also.
>>>
>>> Can you please help to overcome this..
>>> Krishna.
>>>
>>>
>>> On Fri, Sep 9, 2011 at 3:29 PM, RAL HelpDesk {for John Halley
Gotway} <
>>> met_help at ucar.edu> wrote:
>>>
>>>> Krishna,
>>>>
>>>> Good news. I see the problem and it's easy to fix. You were
assuming
>> that
>>>> your GRIB files contained 1-hour accumulated precipitation, but
they do
>> not.
>>>> Instead, they contain a runtime accumulation -
>>>> meaning that the output file for forecast hour "n" contains
>> precipitation
>>>> from hour 0 through hour n. At the end of this message, I've
listed out
>> the
>>>> wgrib output that demonstrates that you have
>>>> runtime accumulations rather than hourly ones. This is actually
the
>>>> default when running the WRF-ARW core, so it's not a problem.
>>>>
>>>> You can easily run PCP-Combine with different arguments to get
the
>>>> accumulations you want. You'll just use the "-subtract" option.
For
>>>> example, to get the 3-hour accumulation between hours 48 and
>>>> 51, you'd run the following command:
>>>>
>>>> pcp_combine -subtract WRFPRS_d01.051 51 WRFPRS_d01.048 48
>>>> WRFPRS_d01_48_51.nc
>>>>
>>>> This subtracts the 48 accumulation in the second file from the 51
hour
>>>> accumulation in the first one. The result is a 3-hour
accumulation
>> between
>>>> hours 48 and 51. You can look at the usage statement
>>>> for pcp_combine to see more details about the "-subtract" option.
>>>>
>>>> By the way, if it's important to you for some reason, you can
>> reconfigure
>>>> your WRF run to dump out hourly accumulations if you'd like.
Take a
>> look in
>>>> the WRF User's Guide for more details, but I
>>>> believe you'd just set "prec_acc_dt" to 60 in the namelist file.
That
>>>> tells WRF to "dump" the precipitation accumulation bucket every
60
>> minutes.
>>>> So the resulting accumulations would be 1-hour. See
>>>> page 5-56 of the WRF-ARW user's guide for more info (
>>>>
>>
http://www.mmm.ucar.edu/wrf/users/docs/user_guide_V3/ARWUsersGuideV3.pdf).
>>>> If you were to do this, you'd revert back to using the "-sum"
option
>>>> for PCP-Combine.
>>>>
>>>> Hope that gets you going.
>>>>
>>>> Thanks,
>>>> John
>>>>
>>>> wgrib output:
>>>>
>>>>
>>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=0:TimeU=1:sfc:0-
0hr
>>>> acc:NAve=0
>>>>
>>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=1:TimeU=1:sfc:0-
1hr
>>>> acc:NAve=0
>>>>
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=2:TimeU=1:sfc:0-
2hr
>>>> acc:NAve=0
>>>>
>>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:TimeU=1:sfc:0-
3hr
>>>> acc:NAve=0
>>>>
>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=4:TimeU=1:sfc:0-
4hr
>>>> acc:NAve=0
>>>>
>>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=5:TimeU=1:sfc:0-
5hr
>>>> acc:NAve=0
>>>>
>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=6:TimeU=1:sfc:0-
6hr
>>>> acc:NAve=0
>>>>
>>
103:8192888:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=7:TimeU=1:sfc:0-
7hr
>>>> acc:NAve=0
>>>>
>>
103:8210588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=8:TimeU=1:sfc:0-
8hr
>>>> acc:NAve=0
>>>>
>>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=9:TimeU=1:sfc:0-
9hr
>>>> acc:NAve=0
>>>>
>>
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=10:TimeU=1:sfc:0-
10hr
>>>> acc:NAve=0
>>>>
>>
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=11:TimeU=1:sfc:0-
11hr
>>>> acc:NAve=0
>>>>
>>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=12:TimeU=1:sfc:0-
12hr
>>>> acc:NAve=0
>>>>
>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=13:TimeU=1:sfc:0-
13hr
>>>> acc:NAve=0
>>>>
>>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=14:TimeU=1:sfc:0-
14hr
>>>> acc:NAve=0
>>>>
>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=15:TimeU=1:sfc:0-
15hr
>>>> acc:NAve=0
>>>>
>>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=16:TimeU=1:sfc:0-
16hr
>>>> acc:NAve=0
>>>>
>>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=17:TimeU=1:sfc:0-
17hr
>>>> acc:NAve=0
>>>>
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=18:TimeU=1:sfc:0-
18hr
>>>> acc:NAve=0
>>>>
>>
103:8080788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=19:TimeU=1:sfc:0-
19hr
>>>> acc:NAve=0
>>>>
>>
103:8098488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=20:TimeU=1:sfc:0-
20hr
>>>> acc:NAve=0
>>>>
>>
103:8098488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=21:TimeU=1:sfc:0-
21hr
>>>> acc:NAve=0
>>>>
>>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=22:TimeU=1:sfc:0-
22hr
>>>> acc:NAve=0
>>>>
>>
103:8092588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=23:TimeU=1:sfc:0-
23hr
>>>> acc:NAve=0
>>>>
>>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=24:TimeU=1:sfc:0-
24hr
>>>> acc:NAve=0
>>>>
>>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=25:TimeU=1:sfc:0-
25hr
>>>> acc:NAve=0
>>>>
>>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=26:TimeU=1:sfc:0-
26hr
>>>> acc:NAve=0
>>>>
>>
103:8086688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=27:TimeU=1:sfc:0-
27hr
>>>> acc:NAve=0
>>>>
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=28:TimeU=1:sfc:0-
28hr
>>>> acc:NAve=0
>>>>
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=29:TimeU=1:sfc:0-
29hr
>>>> acc:NAve=0
>>>>
>>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=30:TimeU=1:sfc:0-
30hr
>>>> acc:NAve=0
>>>>
>>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=31:TimeU=1:sfc:0-
31hr
>>>> acc:NAve=0
>>>>
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=32:TimeU=1:sfc:0-
32hr
>>>> acc:NAve=0
>>>>
>>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=33:TimeU=1:sfc:0-
33hr
>>>> acc:NAve=0
>>>>
>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=34:TimeU=1:sfc:0-
34hr
>>>> acc:NAve=0
>>>>
>>
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=35:TimeU=1:sfc:0-
35hr
>>>> acc:NAve=0
>>>>
>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=36:TimeU=1:sfc:0-
36hr
>>>> acc:NAve=0
>>>>
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=37:TimeU=1:sfc:0-
37hr
>>>> acc:NAve=0
>>>>
>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=38:TimeU=1:sfc:0-
38hr
>>>> acc:NAve=0
>>>>
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=39:TimeU=1:sfc:0-
39hr
>>>> acc:NAve=0
>>>>
>>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=40:TimeU=1:sfc:0-
40hr
>>>> acc:NAve=0
>>>>
>>
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=41:TimeU=1:sfc:0-
41hr
>>>> acc:NAve=0
>>>>
>>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=42:TimeU=1:sfc:0-
42hr
>>>> acc:NAve=0
>>>>
>>
103:8092588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=43:TimeU=1:sfc:0-
43hr
>>>> acc:NAve=0
>>>>
>>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=44:TimeU=1:sfc:0-
44hr
>>>> acc:NAve=0
>>>>
>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=45:TimeU=1:sfc:0-
45hr
>>>> acc:NAve=0
>>>>
>>
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=46:TimeU=1:sfc:0-
46hr
>>>> acc:NAve=0
>>>>
>>
103:8122088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=47:TimeU=1:sfc:0-
47hr
>>>> acc:NAve=0
>>>>
>>
103:8133888:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=48:TimeU=1:sfc:0-
48hr
>>>> acc:NAve=0
>>>>
>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=49:TimeU=1:sfc:0-
49hr
>>>> acc:NAve=0
>>>>
>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=50:TimeU=1:sfc:0-
50hr
>>>> acc:NAve=0
>>>>
>>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=51:TimeU=1:sfc:0-
51hr
>>>> acc:NAve=0
>>>>
>>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=52:TimeU=1:sfc:0-
52hr
>>>> acc:NAve=0
>>>>
>>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=53:TimeU=1:sfc:0-
53hr
>>>> acc:NAve=0
>>>>
>>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=54:TimeU=1:sfc:0-
54hr
>>>> acc:NAve=0
>>>>
>>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=55:TimeU=1:sfc:0-
55hr
>>>> acc:NAve=0
>>>>
>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=56:TimeU=1:sfc:0-
56hr
>>>> acc:NAve=0
>>>>
>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=57:TimeU=1:sfc:0-
57hr
>>>> acc:NAve=0
>>>>
>>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=58:TimeU=1:sfc:0-
58hr
>>>> acc:NAve=0
>>>>
>>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=59:TimeU=1:sfc:0-
59hr
>>>> acc:NAve=0
>>>>
>>
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=60:TimeU=1:sfc:0-
60hr
>>>> acc:NAve=0
>>>>
>>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=61:TimeU=1:sfc:0-
61hr
>>>> acc:NAve=0
>>>>
>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=62:TimeU=1:sfc:0-
62hr
>>>> acc:NAve=0
>>>>
>>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=63:TimeU=1:sfc:0-
63hr
>>>> acc:NAve=0
>>>>
>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=64:TimeU=1:sfc:0-
64hr
>>>> acc:NAve=0
>>>>
>>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=65:TimeU=1:sfc:0-
65hr
>>>> acc:NAve=0
>>>>
>>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=66:TimeU=1:sfc:0-
66hr
>>>> acc:NAve=0
>>>>
>>>>
>>>> On 09/09/2011 11:38 AM, RAL HelpDesk {for krishna kishore} wrote:
>>>>> Please help me in this regard and give details on the error
>>>>
>>>>
>>>
>>>
>>
>>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #49376] Request for preparing gridded observations in grib1 format
From: Tressa Fowler
Time: Tue Sep 13 12:40:36 2011
Dear Krisha,
The interpolation question is a difficult one without a perfect
answer. Generally, it is better to move the finer resolution product
to the coarser scale rather than moving from coarse to fine. However,
it is often recommended to keep the observational field intact and
adjust the models. When comparing several forecasts, the most
important consideration is to put everything on the same scale, which
will then usually be the coarsest of the group.
I am sorry this is not a definitive answer, but I hope you can refine
these recommendations to suit your analyses. If you have further
questions, please let me know.
Tressa
On Sep 13, 2011, at 9:32 AM, RAL HelpDesk {for John Halley Gotway}
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49376 >
>
> Krisha,
>
> I'm glad you got the software working how you'd like. As for what's
the best way scientifically to do the interpolation, I'll defer to a
statistician here. One thing I can tell you is that you
> should probably use the "budget interpolation" option for copygb
which makes sense for precip. You'd do this by passing the "-i3"
command line option when you call copygb. More information about all
> the copygb options can be found here:
>
http://www.dtcenter.org/met/users/support/online_tutorial/METv3.0/copygb/copygb.txt
>
> I'll go ahead and reassign this support ticket to a statistician.
She can give you more information about whether to interpolate to your
model or the observation domain.
>
> John
>
> On 09/12/2011 04:50 PM, RAL HelpDesk {for krishna kishore} wrote:
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49376 >
>>
>> Dear John,
>> Thank you very much for kind and continuous help,
>> Now, I could able to run the grid stat with your help.
>> Now our statistics are looking good.
>> Here, we did match the observations grids (field) onto the model
forecast
>> grids. I know that, copygb will not interpolate.
>>
>> My doubt is, whether the observed rainfall data at the grid points
will be
>> the same before and after the regridding using copygb or not?
>>
>> And According to you, what would be the best way to do-------1.
whether
>> observation grids onto model grids or 2. Model grids onto
observations
>> grids....and What may be the differences in the statistics between
these two
>> methods.
>>
>> I think, You understand my question/doubt?
>>
>> Once again thank you very much for your help.
>> with best regards,
>> Krishna.
>>
>>
>> On Mon, Sep 12, 2011 at 11:19 PM, RAL HelpDesk {for John Halley
Gotway} <
>> met_help at ucar.edu> wrote:
>>
>>> Krishna,
>>>
>>> Yes, in order to perform gridded verification, you must first put
your
>>> forecast and observation files on the same grid. In your case,
you have
>>> GRIB files for both your forecast (GRIB output of the
>>> post processor) and observation (GRIB StageIV analysis). Since
you have
>>> GRIB data, you can use the copygb utility to regrid it. Copygb is
>>> distributed as part of the WRF-PostProcessor and Unified
>>> PostProcessor. In addition, it's available as a stand-alone
package here:
>>> http://www.cpc.ncep.noaa.gov/products/wesley/copygb.html
>>>
>>> Next, you need to decide if you want to regrid your observations
onto your
>>> forecast domain or vice-versa. I'll assume that you want to put
the StageIV
>>> obs onto your model domain. The next step is to
>>> figure out which command line arguments you should use for copygb.
You can
>>> follow the instructions we have posted here for regridding to a
Lambert
>>> Conformal projection:
>>>
>>>
http://www.dtcenter.org/met/users/support/online_tutorial/METv3.0/copygb/run3.php
>>>
>>> Here's your forecast grid information:
>>> :Projection = "Lambert Conformal" ;
>>> :scale_lat_1 = "60.000000" ;
>>> :scale_lat_2 = "30.000000" ;
>>> :lat_pin = "28.223000" ;
>>> :lon_pin = "-110.190000" ;
>>> :x_pin = "0.000000" ;
>>> :y_pin = "0.000000" ;
>>> :lon_orient = "-99.450000" ;
>>> :d_km = "9.000000" ;
>>> :r_km = "6371.200000" ;
>>> :nx = "236" ;
>>> :ny = "200 grid_points" ;
>>>
>>> Following the instructions on that page, this translates to the
following
>>> grid definition:
>>> 255 3 236 200 28223 -110190 8 -99450 9000 9000 0 64 60000 30000
>>>
>>> Prior to using the StageIV data you should run each file through
copygb as
>>> follows:
>>> copygb -xg"255 3 236 200 28223 -110190 8 -99450 9000 9000 0 64
60000
>>> 30000" ST4.2002060300.01h ST4.2002060300.01h_regrid
>>>
>>> Hopefully that'll get you going.
>>>
>>> Thanks,
>>> John
>>>
>>>
>>> On 09/12/2011 11:12 AM, RAL HelpDesk {for krishna kishore} wrote:
>>>>
>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49376 >
>>>>
>>>> Dear John,
>>>> Thank you very much for your help,
>>>> I tried with "subtract" option for WRF output and combine option
for
>>>> OBSERVATIONS.
>>>> I was successful.
>>>> When running the Grid_stat, it is giving the grid mismatch error.
Then
>>> we
>>>> came to know that the dimension of WRF data and observations are
not in
>>> the
>>>> same grid.
>>>> But I dont know how to bring those on the same grid.
>>>>
>>>> I uploaded all the data along with the wrf and observation
dimension file
>>> in
>>>> the ftp site */incoming/irap/met_help/krishna_data*.
>>>> I also uploaded the test_pcp_combine file and Gridstat_config
file also.
>>>>
>>>> Can you please help to overcome this..
>>>> Krishna.
>>>>
>>>>
>>>> On Fri, Sep 9, 2011 at 3:29 PM, RAL HelpDesk {for John Halley
Gotway} <
>>>> met_help at ucar.edu> wrote:
>>>>
>>>>> Krishna,
>>>>>
>>>>> Good news. I see the problem and it's easy to fix. You were
assuming
>>> that
>>>>> your GRIB files contained 1-hour accumulated precipitation, but
they do
>>> not.
>>>>> Instead, they contain a runtime accumulation -
>>>>> meaning that the output file for forecast hour "n" contains
>>> precipitation
>>>>> from hour 0 through hour n. At the end of this message, I've
listed out
>>> the
>>>>> wgrib output that demonstrates that you have
>>>>> runtime accumulations rather than hourly ones. This is actually
the
>>>>> default when running the WRF-ARW core, so it's not a problem.
>>>>>
>>>>> You can easily run PCP-Combine with different arguments to get
the
>>>>> accumulations you want. You'll just use the "-subtract" option.
For
>>>>> example, to get the 3-hour accumulation between hours 48 and
>>>>> 51, you'd run the following command:
>>>>>
>>>>> pcp_combine -subtract WRFPRS_d01.051 51 WRFPRS_d01.048 48
>>>>> WRFPRS_d01_48_51.nc
>>>>>
>>>>> This subtracts the 48 accumulation in the second file from the
51 hour
>>>>> accumulation in the first one. The result is a 3-hour
accumulation
>>> between
>>>>> hours 48 and 51. You can look at the usage statement
>>>>> for pcp_combine to see more details about the "-subtract"
option.
>>>>>
>>>>> By the way, if it's important to you for some reason, you can
>>> reconfigure
>>>>> your WRF run to dump out hourly accumulations if you'd like.
Take a
>>> look in
>>>>> the WRF User's Guide for more details, but I
>>>>> believe you'd just set "prec_acc_dt" to 60 in the namelist file.
That
>>>>> tells WRF to "dump" the precipitation accumulation bucket every
60
>>> minutes.
>>>>> So the resulting accumulations would be 1-hour. See
>>>>> page 5-56 of the WRF-ARW user's guide for more info (
>>>>>
>>>
http://www.mmm.ucar.edu/wrf/users/docs/user_guide_V3/ARWUsersGuideV3.pdf).
>>>>> If you were to do this, you'd revert back to using the "-sum"
option
>>>>> for PCP-Combine.
>>>>>
>>>>> Hope that gets you going.
>>>>>
>>>>> Thanks,
>>>>> John
>>>>>
>>>>> wgrib output:
>>>>>
>>>>>
>>>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=0:TimeU=1:sfc:0-
0hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=1:TimeU=1:sfc:0-
1hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=2:TimeU=1:sfc:0-
2hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=3:TimeU=1:sfc:0-
3hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=4:TimeU=1:sfc:0-
4hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=5:TimeU=1:sfc:0-
5hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=6:TimeU=1:sfc:0-
6hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8192888:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=7:TimeU=1:sfc:0-
7hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8210588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=8:TimeU=1:sfc:0-
8hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=9:TimeU=1:sfc:0-
9hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=10:TimeU=1:sfc:0-
10hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=11:TimeU=1:sfc:0-
11hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=12:TimeU=1:sfc:0-
12hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=13:TimeU=1:sfc:0-
13hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=14:TimeU=1:sfc:0-
14hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=15:TimeU=1:sfc:0-
15hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=16:TimeU=1:sfc:0-
16hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=17:TimeU=1:sfc:0-
17hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=18:TimeU=1:sfc:0-
18hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8080788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=19:TimeU=1:sfc:0-
19hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8098488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=20:TimeU=1:sfc:0-
20hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8098488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=21:TimeU=1:sfc:0-
21hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=22:TimeU=1:sfc:0-
22hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8092588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=23:TimeU=1:sfc:0-
23hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=24:TimeU=1:sfc:0-
24hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=25:TimeU=1:sfc:0-
25hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=26:TimeU=1:sfc:0-
26hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8086688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=27:TimeU=1:sfc:0-
27hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=28:TimeU=1:sfc:0-
28hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=29:TimeU=1:sfc:0-
29hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=30:TimeU=1:sfc:0-
30hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8104388:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=31:TimeU=1:sfc:0-
31hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=32:TimeU=1:sfc:0-
32hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=33:TimeU=1:sfc:0-
33hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=34:TimeU=1:sfc:0-
34hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=35:TimeU=1:sfc:0-
35hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=36:TimeU=1:sfc:0-
36hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=37:TimeU=1:sfc:0-
37hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=38:TimeU=1:sfc:0-
38hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=39:TimeU=1:sfc:0-
39hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=40:TimeU=1:sfc:0-
40hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=41:TimeU=1:sfc:0-
41hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=42:TimeU=1:sfc:0-
42hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8092588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=43:TimeU=1:sfc:0-
43hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8110288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=44:TimeU=1:sfc:0-
44hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8116188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=45:TimeU=1:sfc:0-
45hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8127988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=46:TimeU=1:sfc:0-
46hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8122088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=47:TimeU=1:sfc:0-
47hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8133888:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=48:TimeU=1:sfc:0-
48hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=49:TimeU=1:sfc:0-
49hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=50:TimeU=1:sfc:0-
50hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8181088:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=51:TimeU=1:sfc:0-
51hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=52:TimeU=1:sfc:0-
52hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=53:TimeU=1:sfc:0-
53hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=54:TimeU=1:sfc:0-
54hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=55:TimeU=1:sfc:0-
55hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=56:TimeU=1:sfc:0-
56hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=57:TimeU=1:sfc:0-
57hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=58:TimeU=1:sfc:0-
58hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=59:TimeU=1:sfc:0-
59hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8175188:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=60:TimeU=1:sfc:0-
60hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8145688:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=61:TimeU=1:sfc:0-
61hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8139788:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=62:TimeU=1:sfc:0-
62hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8186988:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=63:TimeU=1:sfc:0-
63hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8157488:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=64:TimeU=1:sfc:0-
64hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8151588:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=65:TimeU=1:sfc:0-
65hr
>>>>> acc:NAve=0
>>>>>
>>>
103:8169288:d=02060200:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=0:P2=66:TimeU=1:sfc:0-
66hr
>>>>> acc:NAve=0
>>>>>
>>>>>
>>>>> On 09/09/2011 11:38 AM, RAL HelpDesk {for krishna kishore}
wrote:
>>>>>> Please help me in this regard and give details on the error
>>>>>
>>>>>
>>>>
>>>>
>>>
>>>
>>
>>
------------------------------------------------
More information about the Met_help
mailing list