[Met_help] [rt.rap.ucar.edu #51478] History for Matching grids for MODE
Paul Oldenburg via RT
met_help at ucar.edu
Mon Nov 21 14:41:38 MST 2011
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello Met Help,
I am having problems with matching grids when running MODE.
My observational data is re-gridded to match the forecast data, yet I still
get an error of no matching grid. The observational data I am using is the
NCEP stage 4 gridded precipitation acquired via the link provided on the
met users page.
Here is a wgrib -V dump of the forecast data (4km gridsize):
rec 1:0:date 2005091500 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0) grid=255
> sfc 0-6hr acc:
> APCP=Total precipitation [kg/m^2]
> timerange 4 P1 0 P2 6 TimeU 1 nx 441 ny 288 GDS grid 3 num_in_ave 0
> missing 0
> center 7 subcenter 0 process 125 Table 2 scan: WE:SN winds(grid)
> Lambert Conf: Lat1 40.338000 Lon1 -93.247000 Lov -100.008000
> Latin1 50.000000 Latin2 40.000000 LatSP 0.000000 LonSP 0.000000
> North Pole (441 x 288) Dx 4.000000 Dy 4.000000 scan 64 mode 136
> min/max data 0 21.729 num bits 15 BDS_Ref 0 DecScale 3 BinScale 0
>
and Here is a wgrib -V dump of the observational data (regridded with
copygb to 4km):
rec 1:0:date 2005091412 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0) grid=255
> sfc 12-18hr acc: bitmap: 11739 undef
> APCP=Total precipitation [kg/m^2]
> timerange 4 P1 12 P2 18 TimeU 1 nx 441 ny 228 GDS grid 3 num_in_ave 0
> missing 0
> center 7 subcenter 4 process 182 Table 2 scan: WE:SN winds(grid)
> Lambert Conf: Lat1 40.338000 Lon1 -93.247000 Lov -100.008000
> Latin1 50.000000 Latin2 40.000000 LatSP 0.000000 LonSP 0.000000
> North Pole (441 x 228) Dx 4.000000 Dy 4.000000 scan 64 mode 8
> min/max data 0 39.11 num bits 12 BDS_Ref 0 DecScale 2 BinScale 0
>
MODE works for my other domain for which I have regridded the observational
data in the same way. For reference I have put the wgrib -V dumps of both
those files as well below:
wgrib of forecast 36km (this is the working one)
rec 1:0:date 2005091500 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0) grid=255
> sfc 0-6hr acc:
> APCP=Total precipitation [kg/m^2]
> timerange 4 P1 0 P2 0 TimeU 1 nx 181 ny 111 GDS grid 3 num_in_ave 0
> missing 0
> center 7 subcenter 0 process 125 Table 2 scan: WE:SN winds(grid)
> Lambert Conf: Lat1 17.952000 Lon1 -128.359000 Lov -100.008000
> Latin1 50.000000 Latin2 40.000000 LatSP 0.000000 LonSP 0.000000
> North Pole (181 x 111) Dx 36.000000 Dy 36.000000 scan 64 mode 136
> min/max data 0 0 num bits 0 BDS_Ref 0 DecScale 0 BinScale 0
>
wgrib of observation 36km (this is the working one, matching the one above
properly)
rec 1:0:date 2005091412 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0) grid=255
> sfc 12-18hr acc: bitmap: 18488 undef
> APCP=Total precipitation [kg/m^2]
> timerange 4 P1 12 P2 18 TimeU 1 nx 181 ny 111 GDS grid 3 num_in_ave 0
> missing 0
> center 7 subcenter 4 process 182 Table 2 scan: WE:SN winds(grid)
> Lambert Conf: Lat1 17.952000 Lon1 -128.359000 Lov -100.008000
> Latin1 50.000000 Latin2 40.000000 LatSP 0.000000 LonSP 0.000000
> North Pole (181 x 111) Dx 36.000000 Dy 36.000000 scan 64 mode 8
> min/max data 0 20.79 num bits 12 BDS_Ref 0 DecScale 2 BinScale 0
>
while some of the numbers in the wgrib dump do not match, they seem to be
the same ones that don't match in the data that is running. I can't seem to
work out why the other domain is not working. My only thought is that the
grid size could be the issue and the re-gridding of the stage-4 is not
working properly.
Any help would be great!
Thanks much,
-Gavin Chensue
--
Gavin D. Chensue
University of Michigan - College of Engineering
Dept. of Atmospheric Oceanic and Space Sciences
(734) 546-0725
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #51478] Matching grids for MODE
From: Paul Oldenburg
Time: Tue Nov 15 08:58:56 2011
Gavin,
The grids have differing numbers of points in the y-direction
(forecast: ny 288, obs: ny 228). If you fix your copygb
command to make that adjustment, I'm pretty sure it should work, since
everything else looks identical. Please let us
know if you have any questions.
Thanks,
Paul
On 11/15/2011 08:47 AM, Gavin Chensue via RT wrote:
>
> Tue Nov 15 08:47:03 2011: Request 51478 was acted upon.
> Transaction: Ticket created by gchensue at umich.edu
> Queue: met_help
> Subject: Matching grids for MODE
> Owner: Nobody
> Requestors: gchensue at umich.edu
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51478 >
>
>
> Hello Met Help,
>
> I am having problems with matching grids when running MODE.
> My observational data is re-gridded to match the forecast data, yet
I still
> get an error of no matching grid. The observational data I am using
is the
> NCEP stage 4 gridded precipitation acquired via the link provided on
the
> met users page.
>
> Here is a wgrib -V dump of the forecast data (4km gridsize):
>
> rec 1:0:date 2005091500 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=255
>> sfc 0-6hr acc:
>> APCP=Total precipitation [kg/m^2]
>> timerange 4 P1 0 P2 6 TimeU 1 nx 441 ny 288 GDS grid 3
num_in_ave 0
>> missing 0
>> center 7 subcenter 0 process 125 Table 2 scan: WE:SN winds(grid)
>> Lambert Conf: Lat1 40.338000 Lon1 -93.247000 Lov -100.008000
>> Latin1 50.000000 Latin2 40.000000 LatSP 0.000000 LonSP
0.000000
>> North Pole (441 x 288) Dx 4.000000 Dy 4.000000 scan 64 mode
136
>> min/max data 0 21.729 num bits 15 BDS_Ref 0 DecScale 3
BinScale 0
>>
>
> and Here is a wgrib -V dump of the observational data (regridded
with
> copygb to 4km):
>
> rec 1:0:date 2005091412 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=255
>> sfc 12-18hr acc: bitmap: 11739 undef
>> APCP=Total precipitation [kg/m^2]
>> timerange 4 P1 12 P2 18 TimeU 1 nx 441 ny 228 GDS grid 3
num_in_ave 0
>> missing 0
>> center 7 subcenter 4 process 182 Table 2 scan: WE:SN winds(grid)
>> Lambert Conf: Lat1 40.338000 Lon1 -93.247000 Lov -100.008000
>> Latin1 50.000000 Latin2 40.000000 LatSP 0.000000 LonSP
0.000000
>> North Pole (441 x 228) Dx 4.000000 Dy 4.000000 scan 64 mode 8
>> min/max data 0 39.11 num bits 12 BDS_Ref 0 DecScale 2 BinScale
0
>>
>
> MODE works for my other domain for which I have regridded the
observational
> data in the same way. For reference I have put the wgrib -V dumps of
both
> those files as well below:
>
> wgrib of forecast 36km (this is the working one)
>
> rec 1:0:date 2005091500 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=255
>> sfc 0-6hr acc:
>> APCP=Total precipitation [kg/m^2]
>> timerange 4 P1 0 P2 0 TimeU 1 nx 181 ny 111 GDS grid 3
num_in_ave 0
>> missing 0
>> center 7 subcenter 0 process 125 Table 2 scan: WE:SN winds(grid)
>> Lambert Conf: Lat1 17.952000 Lon1 -128.359000 Lov -100.008000
>> Latin1 50.000000 Latin2 40.000000 LatSP 0.000000 LonSP
0.000000
>> North Pole (181 x 111) Dx 36.000000 Dy 36.000000 scan 64 mode
136
>> min/max data 0 0 num bits 0 BDS_Ref 0 DecScale 0 BinScale 0
>>
>
> wgrib of observation 36km (this is the working one, matching the one
above
> properly)
>
> rec 1:0:date 2005091412 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
grid=255
>> sfc 12-18hr acc: bitmap: 18488 undef
>> APCP=Total precipitation [kg/m^2]
>> timerange 4 P1 12 P2 18 TimeU 1 nx 181 ny 111 GDS grid 3
num_in_ave 0
>> missing 0
>> center 7 subcenter 4 process 182 Table 2 scan: WE:SN winds(grid)
>> Lambert Conf: Lat1 17.952000 Lon1 -128.359000 Lov -100.008000
>> Latin1 50.000000 Latin2 40.000000 LatSP 0.000000 LonSP
0.000000
>> North Pole (181 x 111) Dx 36.000000 Dy 36.000000 scan 64 mode
8
>> min/max data 0 20.79 num bits 12 BDS_Ref 0 DecScale 2 BinScale
0
>>
>
> while some of the numbers in the wgrib dump do not match, they seem
to be
> the same ones that don't match in the data that is running. I can't
seem to
> work out why the other domain is not working. My only thought is
that the
> grid size could be the issue and the re-gridding of the stage-4 is
not
> working properly.
>
> Any help would be great!
>
> Thanks much,
> -Gavin Chensue
>
>
------------------------------------------------
Subject: Matching grids for MODE
From: Gavin Chensue
Time: Tue Nov 15 09:01:04 2011
my goodness, thank you, I have been looking at so many numbers this
morning
I missed that quite obvious error.
thank you!
-Gavin
On Tue, Nov 15, 2011 at 10:58 AM, Paul Oldenburg via RT
<met_help at ucar.edu>wrote:
> Gavin,
>
> The grids have differing numbers of points in the y-direction
(forecast:
> ny 288, obs: ny 228). If you fix your copygb
> command to make that adjustment, I'm pretty sure it should work,
since
> everything else looks identical. Please let us
> know if you have any questions.
>
> Thanks,
>
> Paul
>
>
> On 11/15/2011 08:47 AM, Gavin Chensue via RT wrote:
> >
> > Tue Nov 15 08:47:03 2011: Request 51478 was acted upon.
> > Transaction: Ticket created by gchensue at umich.edu
> > Queue: met_help
> > Subject: Matching grids for MODE
> > Owner: Nobody
> > Requestors: gchensue at umich.edu
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=51478 >
> >
> >
> > Hello Met Help,
> >
> > I am having problems with matching grids when running MODE.
> > My observational data is re-gridded to match the forecast data,
yet I
> still
> > get an error of no matching grid. The observational data I am
using is
> the
> > NCEP stage 4 gridded precipitation acquired via the link provided
on the
> > met users page.
> >
> > Here is a wgrib -V dump of the forecast data (4km gridsize):
> >
> > rec 1:0:date 2005091500 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
> grid=255
> >> sfc 0-6hr acc:
> >> APCP=Total precipitation [kg/m^2]
> >> timerange 4 P1 0 P2 6 TimeU 1 nx 441 ny 288 GDS grid 3
num_in_ave 0
> >> missing 0
> >> center 7 subcenter 0 process 125 Table 2 scan: WE:SN
winds(grid)
> >> Lambert Conf: Lat1 40.338000 Lon1 -93.247000 Lov -100.008000
> >> Latin1 50.000000 Latin2 40.000000 LatSP 0.000000 LonSP
0.000000
> >> North Pole (441 x 288) Dx 4.000000 Dy 4.000000 scan 64 mode
136
> >> min/max data 0 21.729 num bits 15 BDS_Ref 0 DecScale 3
BinScale 0
> >>
> >
> > and Here is a wgrib -V dump of the observational data (regridded
with
> > copygb to 4km):
> >
> > rec 1:0:date 2005091412 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
> grid=255
> >> sfc 12-18hr acc: bitmap: 11739 undef
> >> APCP=Total precipitation [kg/m^2]
> >> timerange 4 P1 12 P2 18 TimeU 1 nx 441 ny 228 GDS grid 3
num_in_ave 0
> >> missing 0
> >> center 7 subcenter 4 process 182 Table 2 scan: WE:SN
winds(grid)
> >> Lambert Conf: Lat1 40.338000 Lon1 -93.247000 Lov -100.008000
> >> Latin1 50.000000 Latin2 40.000000 LatSP 0.000000 LonSP
0.000000
> >> North Pole (441 x 228) Dx 4.000000 Dy 4.000000 scan 64 mode
8
> >> min/max data 0 39.11 num bits 12 BDS_Ref 0 DecScale 2
BinScale 0
> >>
> >
> > MODE works for my other domain for which I have regridded the
> observational
> > data in the same way. For reference I have put the wgrib -V dumps
of both
> > those files as well below:
> >
> > wgrib of forecast 36km (this is the working one)
> >
> > rec 1:0:date 2005091500 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
> grid=255
> >> sfc 0-6hr acc:
> >> APCP=Total precipitation [kg/m^2]
> >> timerange 4 P1 0 P2 0 TimeU 1 nx 181 ny 111 GDS grid 3
num_in_ave 0
> >> missing 0
> >> center 7 subcenter 0 process 125 Table 2 scan: WE:SN
winds(grid)
> >> Lambert Conf: Lat1 17.952000 Lon1 -128.359000 Lov -100.008000
> >> Latin1 50.000000 Latin2 40.000000 LatSP 0.000000 LonSP
0.000000
> >> North Pole (181 x 111) Dx 36.000000 Dy 36.000000 scan 64
mode 136
> >> min/max data 0 0 num bits 0 BDS_Ref 0 DecScale 0 BinScale 0
> >>
> >
> > wgrib of observation 36km (this is the working one, matching the
one
> above
> > properly)
> >
> > rec 1:0:date 2005091412 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
> grid=255
> >> sfc 12-18hr acc: bitmap: 18488 undef
> >> APCP=Total precipitation [kg/m^2]
> >> timerange 4 P1 12 P2 18 TimeU 1 nx 181 ny 111 GDS grid 3
num_in_ave 0
> >> missing 0
> >> center 7 subcenter 4 process 182 Table 2 scan: WE:SN
winds(grid)
> >> Lambert Conf: Lat1 17.952000 Lon1 -128.359000 Lov -100.008000
> >> Latin1 50.000000 Latin2 40.000000 LatSP 0.000000 LonSP
0.000000
> >> North Pole (181 x 111) Dx 36.000000 Dy 36.000000 scan 64
mode 8
> >> min/max data 0 20.79 num bits 12 BDS_Ref 0 DecScale 2
BinScale 0
> >>
> >
> > while some of the numbers in the wgrib dump do not match, they
seem to be
> > the same ones that don't match in the data that is running. I
can't seem
> to
> > work out why the other domain is not working. My only thought is
that the
> > grid size could be the issue and the re-gridding of the stage-4 is
not
> > working properly.
> >
> > Any help would be great!
> >
> > Thanks much,
> > -Gavin Chensue
> >
> >
>
>
>
--
Gavin D. Chensue
University of Michigan - College of Engineering
Dept. of Atmospheric Oceanic and Space Sciences
(734) 546-0725
------------------------------------------------
More information about the Met_help
mailing list