[Met_help] [rt.rap.ucar.edu #53230] History for Fwd: mode error

Paul Oldenburg via RT met_help at ucar.edu
Fri Feb 17 09:23:54 MST 2012


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




Dear Met help,

I tried to use MODE today, and a problem arise that seemed like Mode has difficulty to find the APCP variable.
The forecast grib file was generated from UPP, and here is what it looks like with wgrib command:

rec 312:11433654:date 2011052000 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0) grid=255 sfc 0-0hr acc:
 APCP=Total precipitation [kg/m^2]
 timerange 4 P1 0 P2 0 TimeU 1  nx 212 ny 167 GDS grid 3 num_in_ave 0 missing 0
 center 7 subcenter 0 process 89 Table 2 scan: WE:SN winds(grid)
 Lambert Conf: Lat1 21.734000 Lon1 -116.920000 Lov -98.941000
     Latin1 35.605000 Latin2 35.605000 LatSP 0.000000 LonSP 0.000000
     North Pole (212 x 167) Dx 18.000000 Dy 18.000000 scan 64 mode 137
 min/max data 0 0  num bits 0  BDS_Ref 0  DecScale 0 BinScale 0

the configure file I have for MODE is:
  fcst_field = "61/A3";

However, after I typed in the command, MODE said:
ERROR:: process_fcst_obs_files() -> APCP_03/A3 not found in file: /..../WRFPRS_d01.000

So I tried to use pcp_combine to convert grib to netcdf file.
However, there are also bunch of errors encountered:
ERROR: sum_grib_files() -> Cannot find a file with a valid time of 20110520_000000 and accumulation time of 030000 matching the regular expression ".*"
ERROR: get_field() -> can't find grib code 4 with accumulation of 030000 in GRIB file: /..../WRFPRS_d01.000

Here is config file.

Thanks!
Di





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

Subject: Re: [rt.rap.ucar.edu #53230] Fwd: mode error
From: Paul Oldenburg
Time: Wed Feb 01 14:42:27 2012

Di,

Notice at the end of the first line of wgrib output: "0-0hr acc".
This means that record contains accumulation data of
0 hour length.  You should look for a record with a three hour
accumulation, which might look like "0-3hr acc" for
example.  You should not have to use pcp_combine if the accumulation
interval that you are interested in is present in
the GRIB file.  Please let me know if you have any other questions.

Paul


On 02/01/2012 01:31 PM, Wu, Di[SCIENCE SYSTEMS APPLICATIONS] via RT
wrote:
>
> Wed Feb 01 13:31:22 2012: Request 53230 was acted upon.
> Transaction: Ticket created by di.wu at nasa.gov
>         Queue: met_help
>       Subject: Fwd: mode error
>         Owner: Nobody
>    Requestors: di.wu at nasa.gov
>        Status: new
>   Ticket<URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=53230>
>
>
>
>
>
> Dear Met help,
>
> I tried to use MODE today, and a problem arise that seemed like Mode
has difficulty to find the APCP variable.
> The forecast grib file was generated from UPP, and here is what it
looks like with wgrib command:
>
> rec 312:11433654:date 2011052000 APCP kpds5=61 kpds6=1 kpds7=0
levels=(0,0) grid=255 sfc 0-0hr acc:
>   APCP=Total precipitation [kg/m^2]
>   timerange 4 P1 0 P2 0 TimeU 1  nx 212 ny 167 GDS grid 3 num_in_ave
0 missing 0
>   center 7 subcenter 0 process 89 Table 2 scan: WE:SN winds(grid)
>   Lambert Conf: Lat1 21.734000 Lon1 -116.920000 Lov -98.941000
>       Latin1 35.605000 Latin2 35.605000 LatSP 0.000000 LonSP
0.000000
>       North Pole (212 x 167) Dx 18.000000 Dy 18.000000 scan 64 mode
137
>   min/max data 0 0  num bits 0  BDS_Ref 0  DecScale 0 BinScale 0
>
> the configure file I have for MODE is:
>    fcst_field = "61/A3";
>
> However, after I typed in the command, MODE said:
> ERROR:: process_fcst_obs_files() ->  APCP_03/A3 not found in file:
/..../WRFPRS_d01.000
>
> So I tried to use pcp_combine to convert grib to netcdf file.
> However, there are also bunch of errors encountered:
> ERROR: sum_grib_files() ->  Cannot find a file with a valid time of
20110520_000000 and accumulation time of 030000 matching the regular
expression ".*"
> ERROR: get_field() ->  can't find grib code 4 with accumulation of
030000 in GRIB file: /..../WRFPRS_d01.000
>
> Here is config file.
>
> Thanks!
> Di
>
>
>


------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #53230] Fwd: mode error
From: Wu, Di[SCIENCE SYSTEMS APPLICATIONS]
Time: Wed Feb 01 15:38:02 2012

Thanks for the suggestion, it works for MODE now.

However, I am still having trouble to get this pcp_combine working.  I
tried different commands and it wasn't working, so what the command
line should be if I am trying to combine this three hourly
precipitation files into one 3-hour acc. precip. file?

info for three files:
rec 1:0:date 2011051923 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=202 sfc 0-1hr acc: bitmap: 2377 undef
  APCP=Total precipitation [kg/m^2]
  timerange 4 P1 0 P2 1 TimeU 1  nx 65 ny 43 GDS grid 5 num_in_ave 0
missing 0
  center 7 subcenter 4 process 182 Table 2 scan: WE:SN winds(grid)
  polar stereo: Lat1 7.838000 Long1 -141.028000 Orient -105.000000
     north pole (65 x 43) Dx 190500 Dy 190500 scan 64 mode 9
  min/max data 0 2.96  num bits 9  BDS_Ref 0  DecScale 2 BinScale 0

rec 1:0:date 2011052000 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=202 sfc 0-1hr acc: bitmap: 2377 undef
  APCP=Total precipitation [kg/m^2]
  timerange 4 P1 0 P2 1 TimeU 1  nx 65 ny 43 GDS grid 5 num_in_ave 0
missing 0
  center 7 subcenter 4 process 182 Table 2 scan: WE:SN winds(grid)
  polar stereo: Lat1 7.838000 Long1 -141.028000 Orient -105.000000
     north pole (65 x 43) Dx 190500 Dy 190500 scan 64 mode 9
  min/max data 0 38.79  num bits 12  BDS_Ref 0  DecScale 2 BinScale 0

rec 1:0:date 2011052001 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=202 sfc 0-1hr acc: bitmap: 2377 undef
  APCP=Total precipitation [kg/m^2]
  timerange 4 P1 0 P2 1 TimeU 1  nx 65 ny 43 GDS grid 5 num_in_ave 0
missing 0
  center 7 subcenter 4 process 182 Table 2 scan: WE:SN winds(grid)
  polar stereo: Lat1 7.838000 Long1 -141.028000 Orient -105.000000
     north pole (65 x 43) Dx 190500 Dy 190500 scan 64 mode 9
  min/max data 0 24.6  num bits 12  BDS_Ref 0  DecScale 2 BinScale 0

Thanks,
Di


On Feb 1, 2012, at 4:42 PM, Paul Oldenburg via RT wrote:

> Di,
>
> Notice at the end of the first line of wgrib output: "0-0hr acc".
This means that record contains accumulation data of
> 0 hour length.  You should look for a record with a three hour
accumulation, which might look like "0-3hr acc" for
> example.  You should not have to use pcp_combine if the accumulation
interval that you are interested in is present in
> the GRIB file.  Please let me know if you have any other questions.
>
> Paul
>
>
> On 02/01/2012 01:31 PM, Wu, Di[SCIENCE SYSTEMS APPLICATIONS] via RT
wrote:
>>
>> Wed Feb 01 13:31:22 2012: Request 53230 was acted upon.
>> Transaction: Ticket created by di.wu at nasa.gov
>>        Queue: met_help
>>      Subject: Fwd: mode error
>>        Owner: Nobody
>>   Requestors: di.wu at nasa.gov
>>       Status: new
>>  Ticket<URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=53230>
>>
>>
>>
>>
>>
>> Dear Met help,
>>
>> I tried to use MODE today, and a problem arise that seemed like
Mode has difficulty to find the APCP variable.
>> The forecast grib file was generated from UPP, and here is what it
looks like with wgrib command:
>>
>> rec 312:11433654:date 2011052000 APCP kpds5=61 kpds6=1 kpds7=0
levels=(0,0) grid=255 sfc 0-0hr acc:
>>  APCP=Total precipitation [kg/m^2]
>>  timerange 4 P1 0 P2 0 TimeU 1  nx 212 ny 167 GDS grid 3 num_in_ave
0 missing 0
>>  center 7 subcenter 0 process 89 Table 2 scan: WE:SN winds(grid)
>>  Lambert Conf: Lat1 21.734000 Lon1 -116.920000 Lov -98.941000
>>      Latin1 35.605000 Latin2 35.605000 LatSP 0.000000 LonSP
0.000000
>>      North Pole (212 x 167) Dx 18.000000 Dy 18.000000 scan 64 mode
137
>>  min/max data 0 0  num bits 0  BDS_Ref 0  DecScale 0 BinScale 0
>>
>> the configure file I have for MODE is:
>>   fcst_field = "61/A3";
>>
>> However, after I typed in the command, MODE said:
>> ERROR:: process_fcst_obs_files() ->  APCP_03/A3 not found in file:
/..../WRFPRS_d01.000
>>
>> So I tried to use pcp_combine to convert grib to netcdf file.
>> However, there are also bunch of errors encountered:
>> ERROR: sum_grib_files() ->  Cannot find a file with a valid time of
20110520_000000 and accumulation time of 030000 matching the regular
expression ".*"
>> ERROR: get_field() ->  can't find grib code 4 with accumulation of
030000 in GRIB file: /..../WRFPRS_d01.000
>>
>> Here is config file.
>>
>> Thanks!
>> Di
>>
>>
>>
>
>


------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #53230] Fwd: mode error
From: Paul Oldenburg
Time: Wed Feb 01 15:56:23 2012

Di,

First of all, it appears that you are trying to accumulate 1hr precip
from three files with three different
initialization times (2011051923, 2011052000 and 2011052001).  I'm not
sure what your verification experiment is, but
this seems like a strange way to combine forecasted precip amounts.
If this data will be used as the observation data,
then it seems less peculiar.

Try this:

pcp_combine -add [file1] 1 [file2] 1 [file3] 1 [out_file] -v 3

Please let me know if you have any trouble.

Paul


On 02/01/2012 03:38 PM, Wu, Di[SCIENCE SYSTEMS APPLICATIONS] via RT
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=53230>
>
> Thanks for the suggestion, it works for MODE now.
>
> However, I am still having trouble to get this pcp_combine working.
I tried different commands and it wasn't working, so what the command
line should be if I am trying to combine this three hourly
precipitation files into one 3-hour acc. precip. file?
>
> info for three files:
> rec 1:0:date 2011051923 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=202 sfc 0-1hr acc: bitmap: 2377 undef
>    APCP=Total precipitation [kg/m^2]
>    timerange 4 P1 0 P2 1 TimeU 1  nx 65 ny 43 GDS grid 5 num_in_ave
0 missing 0
>    center 7 subcenter 4 process 182 Table 2 scan: WE:SN winds(grid)
>    polar stereo: Lat1 7.838000 Long1 -141.028000 Orient -105.000000
>       north pole (65 x 43) Dx 190500 Dy 190500 scan 64 mode 9
>    min/max data 0 2.96  num bits 9  BDS_Ref 0  DecScale 2 BinScale 0
>
> rec 1:0:date 2011052000 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=202 sfc 0-1hr acc: bitmap: 2377 undef
>    APCP=Total precipitation [kg/m^2]
>    timerange 4 P1 0 P2 1 TimeU 1  nx 65 ny 43 GDS grid 5 num_in_ave
0 missing 0
>    center 7 subcenter 4 process 182 Table 2 scan: WE:SN winds(grid)
>    polar stereo: Lat1 7.838000 Long1 -141.028000 Orient -105.000000
>       north pole (65 x 43) Dx 190500 Dy 190500 scan 64 mode 9
>    min/max data 0 38.79  num bits 12  BDS_Ref 0  DecScale 2 BinScale
0
>
> rec 1:0:date 2011052001 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=202 sfc 0-1hr acc: bitmap: 2377 undef
>    APCP=Total precipitation [kg/m^2]
>    timerange 4 P1 0 P2 1 TimeU 1  nx 65 ny 43 GDS grid 5 num_in_ave
0 missing 0
>    center 7 subcenter 4 process 182 Table 2 scan: WE:SN winds(grid)
>    polar stereo: Lat1 7.838000 Long1 -141.028000 Orient -105.000000
>       north pole (65 x 43) Dx 190500 Dy 190500 scan 64 mode 9
>    min/max data 0 24.6  num bits 12  BDS_Ref 0  DecScale 2 BinScale
0
>
> Thanks,
> Di
>
>
> On Feb 1, 2012, at 4:42 PM, Paul Oldenburg via RT wrote:
>
>> Di,
>>
>> Notice at the end of the first line of wgrib output: "0-0hr acc".
This means that record contains accumulation data of
>> 0 hour length.  You should look for a record with a three hour
accumulation, which might look like "0-3hr acc" for
>> example.  You should not have to use pcp_combine if the
accumulation interval that you are interested in is present in
>> the GRIB file.  Please let me know if you have any other questions.
>>
>> Paul
>>
>>
>> On 02/01/2012 01:31 PM, Wu, Di[SCIENCE SYSTEMS APPLICATIONS] via RT
wrote:
>>>
>>> Wed Feb 01 13:31:22 2012: Request 53230 was acted upon.
>>> Transaction: Ticket created by di.wu at nasa.gov
>>>         Queue: met_help
>>>       Subject: Fwd: mode error
>>>         Owner: Nobody
>>>    Requestors: di.wu at nasa.gov
>>>        Status: new
>>>   Ticket<URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=53230>
>>>
>>>
>>>
>>>
>>>
>>> Dear Met help,
>>>
>>> I tried to use MODE today, and a problem arise that seemed like
Mode has difficulty to find the APCP variable.
>>> The forecast grib file was generated from UPP, and here is what it
looks like with wgrib command:
>>>
>>> rec 312:11433654:date 2011052000 APCP kpds5=61 kpds6=1 kpds7=0
levels=(0,0) grid=255 sfc 0-0hr acc:
>>>   APCP=Total precipitation [kg/m^2]
>>>   timerange 4 P1 0 P2 0 TimeU 1  nx 212 ny 167 GDS grid 3
num_in_ave 0 missing 0
>>>   center 7 subcenter 0 process 89 Table 2 scan: WE:SN winds(grid)
>>>   Lambert Conf: Lat1 21.734000 Lon1 -116.920000 Lov -98.941000
>>>       Latin1 35.605000 Latin2 35.605000 LatSP 0.000000 LonSP
0.000000
>>>       North Pole (212 x 167) Dx 18.000000 Dy 18.000000 scan 64
mode 137
>>>   min/max data 0 0  num bits 0  BDS_Ref 0  DecScale 0 BinScale 0
>>>
>>> the configure file I have for MODE is:
>>>    fcst_field = "61/A3";
>>>
>>> However, after I typed in the command, MODE said:
>>> ERROR:: process_fcst_obs_files() ->   APCP_03/A3 not found in
file: /..../WRFPRS_d01.000
>>>
>>> So I tried to use pcp_combine to convert grib to netcdf file.
>>> However, there are also bunch of errors encountered:
>>> ERROR: sum_grib_files() ->   Cannot find a file with a valid time
of 20110520_000000 and accumulation time of 030000 matching the
regular expression ".*"
>>> ERROR: get_field() ->   can't find grib code 4 with accumulation
of 030000 in GRIB file: /..../WRFPRS_d01.000
>>>
>>> Here is config file.
>>>
>>> Thanks!
>>> Di
>>>
>>>
>>>
>>
>>
>


------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #53230] Fwd: mode error
From: John Halley Gotway
Time: Wed Feb 01 16:10:20 2012

Paul and Di,

The pcp_combine "-add" option should do the trick as Paul suggested.
Alternatively, the "-sum" option would work as well:

pcp_combine -sum 00000000_000000 01 20110522_000000 03 apcp_03h.nc
-pcpdir /path/to/directory/containing/the/grib files

This command will search through the files in the directory you
specify, looking for 3 1-hourly records that can be summed up into a
3-hourly accumulation that's valid at the time you specify
(20110522_000000).  The long string of zero's tells the tool to not
worry about keeping the init times constants across the files - since
you're working with observation files.

Please see the pcp_combine usage statement or the MET User's Guide for
more details.

Either way should work - it's up to you to use whichever is most
convenient.

Thanks,
John Halley Gotway

On 02/01/2012 03:56 PM, Paul Oldenburg via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=53230>
>
> Di,
>
> First of all, it appears that you are trying to accumulate 1hr
precip from three files with three different
> initialization times (2011051923, 2011052000 and 2011052001).  I'm
not sure what your verification experiment is, but
> this seems like a strange way to combine forecasted precip amounts.
If this data will be used as the observation data,
> then it seems less peculiar.
>
> Try this:
>
> pcp_combine -add [file1] 1 [file2] 1 [file3] 1 [out_file] -v 3
>
> Please let me know if you have any trouble.
>
> Paul
>
>
> On 02/01/2012 03:38 PM, Wu, Di[SCIENCE SYSTEMS APPLICATIONS] via RT
wrote:
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=53230>
>>
>> Thanks for the suggestion, it works for MODE now.
>>
>> However, I am still having trouble to get this pcp_combine working.
I tried different commands and it wasn't working, so what the command
line should be if I am trying to combine this three hourly
precipitation files into one 3-hour acc. precip. file?
>>
>> info for three files:
>> rec 1:0:date 2011051923 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=202 sfc 0-1hr acc: bitmap: 2377 undef
>>     APCP=Total precipitation [kg/m^2]
>>     timerange 4 P1 0 P2 1 TimeU 1  nx 65 ny 43 GDS grid 5
num_in_ave 0 missing 0
>>     center 7 subcenter 4 process 182 Table 2 scan: WE:SN
winds(grid)
>>     polar stereo: Lat1 7.838000 Long1 -141.028000 Orient
-105.000000
>>        north pole (65 x 43) Dx 190500 Dy 190500 scan 64 mode 9
>>     min/max data 0 2.96  num bits 9  BDS_Ref 0  DecScale 2 BinScale
0
>>
>> rec 1:0:date 2011052000 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=202 sfc 0-1hr acc: bitmap: 2377 undef
>>     APCP=Total precipitation [kg/m^2]
>>     timerange 4 P1 0 P2 1 TimeU 1  nx 65 ny 43 GDS grid 5
num_in_ave 0 missing 0
>>     center 7 subcenter 4 process 182 Table 2 scan: WE:SN
winds(grid)
>>     polar stereo: Lat1 7.838000 Long1 -141.028000 Orient
-105.000000
>>        north pole (65 x 43) Dx 190500 Dy 190500 scan 64 mode 9
>>     min/max data 0 38.79  num bits 12  BDS_Ref 0  DecScale 2
BinScale 0
>>
>> rec 1:0:date 2011052001 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=202 sfc 0-1hr acc: bitmap: 2377 undef
>>     APCP=Total precipitation [kg/m^2]
>>     timerange 4 P1 0 P2 1 TimeU 1  nx 65 ny 43 GDS grid 5
num_in_ave 0 missing 0
>>     center 7 subcenter 4 process 182 Table 2 scan: WE:SN
winds(grid)
>>     polar stereo: Lat1 7.838000 Long1 -141.028000 Orient
-105.000000
>>        north pole (65 x 43) Dx 190500 Dy 190500 scan 64 mode 9
>>     min/max data 0 24.6  num bits 12  BDS_Ref 0  DecScale 2
BinScale 0
>>
>> Thanks,
>> Di
>>
>>
>> On Feb 1, 2012, at 4:42 PM, Paul Oldenburg via RT wrote:
>>
>>> Di,
>>>
>>> Notice at the end of the first line of wgrib output: "0-0hr acc".
This means that record contains accumulation data of
>>> 0 hour length.  You should look for a record with a three hour
accumulation, which might look like "0-3hr acc" for
>>> example.  You should not have to use pcp_combine if the
accumulation interval that you are interested in is present in
>>> the GRIB file.  Please let me know if you have any other
questions.
>>>
>>> Paul
>>>
>>>
>>> On 02/01/2012 01:31 PM, Wu, Di[SCIENCE SYSTEMS APPLICATIONS] via
RT wrote:
>>>>
>>>> Wed Feb 01 13:31:22 2012: Request 53230 was acted upon.
>>>> Transaction: Ticket created by di.wu at nasa.gov
>>>>          Queue: met_help
>>>>        Subject: Fwd: mode error
>>>>          Owner: Nobody
>>>>     Requestors: di.wu at nasa.gov
>>>>         Status: new
>>>>    Ticket<URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=53230>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> Dear Met help,
>>>>
>>>> I tried to use MODE today, and a problem arise that seemed like
Mode has difficulty to find the APCP variable.
>>>> The forecast grib file was generated from UPP, and here is what
it looks like with wgrib command:
>>>>
>>>> rec 312:11433654:date 2011052000 APCP kpds5=61 kpds6=1 kpds7=0
levels=(0,0) grid=255 sfc 0-0hr acc:
>>>>    APCP=Total precipitation [kg/m^2]
>>>>    timerange 4 P1 0 P2 0 TimeU 1  nx 212 ny 167 GDS grid 3
num_in_ave 0 missing 0
>>>>    center 7 subcenter 0 process 89 Table 2 scan: WE:SN
winds(grid)
>>>>    Lambert Conf: Lat1 21.734000 Lon1 -116.920000 Lov -98.941000
>>>>        Latin1 35.605000 Latin2 35.605000 LatSP 0.000000 LonSP
0.000000
>>>>        North Pole (212 x 167) Dx 18.000000 Dy 18.000000 scan 64
mode 137
>>>>    min/max data 0 0  num bits 0  BDS_Ref 0  DecScale 0 BinScale 0
>>>>
>>>> the configure file I have for MODE is:
>>>>     fcst_field = "61/A3";
>>>>
>>>> However, after I typed in the command, MODE said:
>>>> ERROR:: process_fcst_obs_files() ->    APCP_03/A3 not found in
file: /..../WRFPRS_d01.000
>>>>
>>>> So I tried to use pcp_combine to convert grib to netcdf file.
>>>> However, there are also bunch of errors encountered:
>>>> ERROR: sum_grib_files() ->    Cannot find a file with a valid
time of 20110520_000000 and accumulation time of 030000 matching the
regular expression ".*"
>>>> ERROR: get_field() ->    can't find grib code 4 with accumulation
of 030000 in GRIB file: /..../WRFPRS_d01.000
>>>>
>>>> Here is config file.
>>>>
>>>> Thanks!
>>>> Di
>>>>
>>>>
>>>>
>>>
>>>
>>
>

------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #53230] Fwd: mode error
From: Wu, Di[SCIENCE SYSTEMS APPLICATIONS]
Time: Thu Feb 02 09:30:44 2012

Hi, Paul and John

I tried both methods and only the -add works at this time.  The -sum
gives me the same complain as it did before unfortunately.
I was trying to sum three level 4 data together.  The -sum command
seems looking for some particular strings that either I or the obs.
data didn't provide right.
Anyway, I am happy with the fact it is working, and thank you both for
the help!

Di

On Feb 1, 2012, at 6:10 PM, John Halley Gotway via RT wrote:

> Paul and Di,
>
> The pcp_combine "-add" option should do the trick as Paul suggested.
Alternatively, the "-sum" option would work as well:
>
> pcp_combine -sum 00000000_000000 01 20110522_000000 03 apcp_03h.nc
-pcpdir /path/to/directory/containing/the/grib files
>
> This command will search through the files in the directory you
specify, looking for 3 1-hourly records that can be summed up into a
3-hourly accumulation that's valid at the time you specify
> (20110522_000000).  The long string of zero's tells the tool to not
worry about keeping the init times constants across the files - since
you're working with observation files.
>
> Please see the pcp_combine usage statement or the MET User's Guide
for more details.
>
> Either way should work - it's up to you to use whichever is most
convenient.
>
> Thanks,
> John Halley Gotway
>
> On 02/01/2012 03:56 PM, Paul Oldenburg via RT wrote:
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=53230>
>>
>> Di,
>>
>> First of all, it appears that you are trying to accumulate 1hr
precip from three files with three different
>> initialization times (2011051923, 2011052000 and 2011052001).  I'm
not sure what your verification experiment is, but
>> this seems like a strange way to combine forecasted precip amounts.
If this data will be used as the observation data,
>> then it seems less peculiar.
>>
>> Try this:
>>
>> pcp_combine -add [file1] 1 [file2] 1 [file3] 1 [out_file] -v 3
>>
>> Please let me know if you have any trouble.
>>
>> Paul
>>
>>
>> On 02/01/2012 03:38 PM, Wu, Di[SCIENCE SYSTEMS APPLICATIONS] via RT
wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=53230>
>>>
>>> Thanks for the suggestion, it works for MODE now.
>>>
>>> However, I am still having trouble to get this pcp_combine
working.  I tried different commands and it wasn't working, so what
the command line should be if I am trying to combine this three hourly
precipitation files into one 3-hour acc. precip. file?
>>>
>>> info for three files:
>>> rec 1:0:date 2011051923 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=202 sfc 0-1hr acc: bitmap: 2377 undef
>>>    APCP=Total precipitation [kg/m^2]
>>>    timerange 4 P1 0 P2 1 TimeU 1  nx 65 ny 43 GDS grid 5
num_in_ave 0 missing 0
>>>    center 7 subcenter 4 process 182 Table 2 scan: WE:SN
winds(grid)
>>>    polar stereo: Lat1 7.838000 Long1 -141.028000 Orient
-105.000000
>>>       north pole (65 x 43) Dx 190500 Dy 190500 scan 64 mode 9
>>>    min/max data 0 2.96  num bits 9  BDS_Ref 0  DecScale 2 BinScale
0
>>>
>>> rec 1:0:date 2011052000 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=202 sfc 0-1hr acc: bitmap: 2377 undef
>>>    APCP=Total precipitation [kg/m^2]
>>>    timerange 4 P1 0 P2 1 TimeU 1  nx 65 ny 43 GDS grid 5
num_in_ave 0 missing 0
>>>    center 7 subcenter 4 process 182 Table 2 scan: WE:SN
winds(grid)
>>>    polar stereo: Lat1 7.838000 Long1 -141.028000 Orient
-105.000000
>>>       north pole (65 x 43) Dx 190500 Dy 190500 scan 64 mode 9
>>>    min/max data 0 38.79  num bits 12  BDS_Ref 0  DecScale 2
BinScale 0
>>>
>>> rec 1:0:date 2011052001 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=202 sfc 0-1hr acc: bitmap: 2377 undef
>>>    APCP=Total precipitation [kg/m^2]
>>>    timerange 4 P1 0 P2 1 TimeU 1  nx 65 ny 43 GDS grid 5
num_in_ave 0 missing 0
>>>    center 7 subcenter 4 process 182 Table 2 scan: WE:SN
winds(grid)
>>>    polar stereo: Lat1 7.838000 Long1 -141.028000 Orient
-105.000000
>>>       north pole (65 x 43) Dx 190500 Dy 190500 scan 64 mode 9
>>>    min/max data 0 24.6  num bits 12  BDS_Ref 0  DecScale 2
BinScale 0
>>>
>>> Thanks,
>>> Di
>>>
>>>
>>> On Feb 1, 2012, at 4:42 PM, Paul Oldenburg via RT wrote:
>>>
>>>> Di,
>>>>
>>>> Notice at the end of the first line of wgrib output: "0-0hr acc".
This means that record contains accumulation data of
>>>> 0 hour length.  You should look for a record with a three hour
accumulation, which might look like "0-3hr acc" for
>>>> example.  You should not have to use pcp_combine if the
accumulation interval that you are interested in is present in
>>>> the GRIB file.  Please let me know if you have any other
questions.
>>>>
>>>> Paul
>>>>
>>>>
>>>> On 02/01/2012 01:31 PM, Wu, Di[SCIENCE SYSTEMS APPLICATIONS] via
RT wrote:
>>>>>
>>>>> Wed Feb 01 13:31:22 2012: Request 53230 was acted upon.
>>>>> Transaction: Ticket created by di.wu at nasa.gov
>>>>>         Queue: met_help
>>>>>       Subject: Fwd: mode error
>>>>>         Owner: Nobody
>>>>>    Requestors: di.wu at nasa.gov
>>>>>        Status: new
>>>>>   Ticket<URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=53230>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Dear Met help,
>>>>>
>>>>> I tried to use MODE today, and a problem arise that seemed like
Mode has difficulty to find the APCP variable.
>>>>> The forecast grib file was generated from UPP, and here is what
it looks like with wgrib command:
>>>>>
>>>>> rec 312:11433654:date 2011052000 APCP kpds5=61 kpds6=1 kpds7=0
levels=(0,0) grid=255 sfc 0-0hr acc:
>>>>>   APCP=Total precipitation [kg/m^2]
>>>>>   timerange 4 P1 0 P2 0 TimeU 1  nx 212 ny 167 GDS grid 3
num_in_ave 0 missing 0
>>>>>   center 7 subcenter 0 process 89 Table 2 scan: WE:SN
winds(grid)
>>>>>   Lambert Conf: Lat1 21.734000 Lon1 -116.920000 Lov -98.941000
>>>>>       Latin1 35.605000 Latin2 35.605000 LatSP 0.000000 LonSP
0.000000
>>>>>       North Pole (212 x 167) Dx 18.000000 Dy 18.000000 scan 64
mode 137
>>>>>   min/max data 0 0  num bits 0  BDS_Ref 0  DecScale 0 BinScale 0
>>>>>
>>>>> the configure file I have for MODE is:
>>>>>    fcst_field = "61/A3";
>>>>>
>>>>> However, after I typed in the command, MODE said:
>>>>> ERROR:: process_fcst_obs_files() ->    APCP_03/A3 not found in
file: /..../WRFPRS_d01.000
>>>>>
>>>>> So I tried to use pcp_combine to convert grib to netcdf file.
>>>>> However, there are also bunch of errors encountered:
>>>>> ERROR: sum_grib_files() ->    Cannot find a file with a valid
time of 20110520_000000 and accumulation time of 030000 matching the
regular expression ".*"
>>>>> ERROR: get_field() ->    can't find grib code 4 with
accumulation of 030000 in GRIB file: /..../WRFPRS_d01.000
>>>>>
>>>>> Here is config file.
>>>>>
>>>>> Thanks!
>>>>> Di
>>>>>
>>>>>
>>>>>
>>>>
>>>>
>>>
>>
>


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


More information about the Met_help mailing list