[Met_help] [rt.rap.ucar.edu #49377] History for Request for preparing gridded observations in grib1 format-need respond

Tue Sep 13 14:41:50 MDT 2011

----------------------------------------------------------------
  Initial Request
----------------------------------------------------------------

Dear Tressa
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,
jam
***********************************************************************************************************
/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).
***********************************************************************************************************

On Mon, Jun 13, 2011 at 9:50 PM, krishna kishore <osurikishore at gmail.com> wrote:
>
> Hi John,
> Thank you for your help in instaling METV3.0 and running point-stat.
> Now i want to use grid-stat for rainfall validation for Indian Domain.
> For that, first i should prepare the rainfall gridded data in grib1 format.
> Can you help me in that.
> If you have any utilitities to convert TRMM or (CMORPHS etc) any rainfall data to the grib1 format (that MET expects), Please share with me.  That will help for my PhD thesis.
> Hope the positive response,
> Thanks,
> Krishna.
>
> --
> ###############################
> O. Krishna Kishore
> Project Associate,
> Room No: 411, Block - 6,
> Centre for Atmospheric Sciences,
> IIT Delhi, Hauz Khas - 110016
> New Delhi.
> ###############################
>

--
###############################
O. Krishna Kishore
Project Associate,
Room No: 411, Block - 6,
Centre for Atmospheric Sciences,
IIT Delhi, Hauz Khas - 110016
New Delhi.
###############################

--
###############################
O. Krishna Kishore
Project Associate,
Room No: 411, Block - 6,
Centre for Atmospheric Sciences,
IIT Delhi, Hauz Khas - 110016
New Delhi.
###############################

----------------------------------------------------------------
  Complete Ticket History
----------------------------------------------------------------

Subject: Re: [rt.rap.ucar.edu #49377] Request for preparing gridded observations in grib1 format-need respond
From: John Halley Gotway
Time: Thu Sep 01 15:15:52 2011

Jam,

This looks like the same question we recently received from Krishna
Kishore.  Here's the response I just sent:

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 03:12 PM, RAL HelpDesk {for ucharuso at purdue.edu} wrote:
>
> Thu Sep 01 15:12:31 2011: Request 49377 was acted upon.
> Transaction: Ticket created by ucharuso at purdue.edu
>        Queue: met_help
>      Subject: Request for preparing gridded observations in grib1
format-need respond
>        Owner: Nobody
>   Requestors: ucharuso at purdue.edu
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=49377 >
>
>
> Dear Tressa
> 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,
> jam
>
***********************************************************************************************************
> /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).
>
***********************************************************************************************************
>
>
> On Mon, Jun 13, 2011 at 9:50 PM, krishna kishore
<osurikishore at gmail.com> wrote:
>>
>> Hi John,
>> Thank you for your help in instaling METV3.0 and running point-
stat.
>> Now i want to use grid-stat for rainfall validation for Indian
Domain.
>> For that, first i should prepare the rainfall gridded data in grib1
format.
>> Can you help me in that.
>> If you have any utilitities to convert TRMM or (CMORPHS etc) any
rainfall data to the grib1 format (that MET expects), Please share
with me.  That will help for my PhD thesis.
>> Hope the positive response,
>> Thanks,
>> Krishna.
>>
>> --
>> ###############################
>> O. Krishna Kishore
>> Project Associate,
>> Room No: 411, Block - 6,
>> Centre for Atmospheric Sciences,
>> IIT Delhi, Hauz Khas - 110016
>> New Delhi.
>> ###############################
>>
>
>
>
> --
> ###############################
> O. Krishna Kishore
> Project Associate,
> Room No: 411, Block - 6,
> Centre for Atmospheric Sciences,
> IIT Delhi, Hauz Khas - 110016
> New Delhi.
> ###############################
>
>
>
>
> --
> ###############################
> O. Krishna Kishore
> Project Associate,
> Room No: 411, Block - 6,
> Centre for Atmospheric Sciences,
> IIT Delhi, Hauz Khas - 110016
> New Delhi.
> ###############################

------------------------------------------------