[Met_help] [rt.rap.ucar.edu #53051] History for Cloud File Mismatch

John Halley Gotway via RT met_help at ucar.edu
Thu Jan 26 13:45:14 MST 2012


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

Hi guys, when using MODE to compare a cloud analysis to cloud forecast, I am getting the following error: 

 

ERROR: process_fcst_obs_files() -> The forecast and observation grids do not match.

 

The grib header of the dcf cloud forecast is:  

 

rec 1:0:date 2012012600 TCDC kpds5=71 kpds6=200 kpds7=0 levels=(0,0) grid=0 atmos col 0hr fcst:

  TCDC=Total cloud cover [%]

  timerange 1 P1 0 P2 0 TimeU 1  nx 199 ny 149 GDS grid 3 num_in_ave 0 missing 0

  center 57 subcenter 2 process 36 Table 2 scan: WE:SN winds(grid)

  Lambert Conf: Lat1 0.888000 Lon1 72.538000 Lov 105.000000

      Latin1 60.000000 Latin2 30.000000 LatSP -90.000000 LonSP 0.000000

      North Pole (199 x 149) Dx 45.000000 Dy 45.000000 scan 64 mode 8

  min/max data 0 93  num bits 7  BDS_Ref 0  DecScale 0 BinScale 0

 

The grib header for the cloud analysis file (wwmca) is:

 

rec 1:0:date 2012012600 TCDC kpds5=71 kpds6=200 kpds7=0 levels=(0,0) grid=255 atmos col anl:

  TCDC=Total cloud cover [%]

  timerange 0 P1 0 P2 0 TimeU 1  nx 199 ny 149 GDS grid 3 num_in_ave 0 missing 0

  center 57 subcenter 2 process 32 Table 128 scan: WE:SN winds(grid)

  Lambert Conf: Lat1 0.888000 Lon1 72.537000 Lov 105.000000

      Latin1 60.000000 Latin2 30.000000 LatSP 0.000000 LonSP 0.000000

      North Pole (199 x 149) Dx 45.000000 Dy 45.000000 scan 64 mode 8

  min/max data 0 100  num bits 7  BDS_Ref 0  DecScale 0 BinScale 0

 

The only difference I noticed was in the thousands field of the Lon1.  Would that cause this error.

 

The command line I used is:

 

/home/qcteam/METv3.0.1/bin/mode /h/data/global/WXQC/data/dcf/t06a/fcsts_t06a_2012012600_000.GR1 /h/data/global/WXQC/data/temp/met/cloud_obs_file_71 /home/qcteam/METv3.0.1/data/config/DCFModeConfig_45 -outdir /h/data/global/WXQC/data/met/cloud/t06/

 

 

I have attached the to grib files and config file.

 

Thanks

 

Bob Craig

 

 



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

Subject: Re: [rt.rap.ucar.edu #53051] Cloud File Mismatch
From: John Halley Gotway
Time: Thu Jan 26 13:32:59 2012

Bob,

Yes, the difference in the thousanths place of the Lon1 parameter
would cause this error to occur.  When checking for grid equality, MET
is checking each of the components that define the grid.  To
resolve this, I'd suggest running one of the GRIB files through coypgb
to regrid it to the other's domain.  An example of using copygb to
regrid a lambert conformal grid can be found on the MET
website in the online tutorial:
    http://www.dtcenter.org/met/users/support/online_tutorial/METv3.0/copygb/run3.php

Alternatively, if you decide that you don't want MET to do this strict
of grid checking, you could modify the source code to use different
logic.  Checking for grid equality is done in the file
"grid_base.cc" in this function:
    bool operator==(const Grid & g1, const Grid & g2)

In earlier versions of MET, we simply checked that Nx and Ny were
same, but we ended up missing small differences like the one you have
in your data!

Hope that helps clarify.

John Halley Gotway
met_help at ucar.edu

On 01/26/2012 01:22 PM, Craig, Robert J Civ USAF AFWA 16 WS/WXN via RT
wrote:
>
> Thu Jan 26 13:22:23 2012: Request 53051 was acted upon.
> Transaction: Ticket created by Robert.Craig at offutt.af.mil
>         Queue: met_help
>       Subject: Cloud File Mismatch
>         Owner: Nobody
>    Requestors: Robert.Craig at offutt.af.mil
>        Status: new
>   Ticket<URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=53051>
>
>
> Hi guys, when using MODE to compare a cloud analysis to cloud
forecast, I am getting the following error:
>
>
>
> ERROR: process_fcst_obs_files() ->  The forecast and observation
grids do not match.
>
>
>
> The grib header of the dcf cloud forecast is:
>
>
>
> rec 1:0:date 2012012600 TCDC kpds5=71 kpds6=200 kpds7=0 levels=(0,0)
grid=0 atmos col 0hr fcst:
>
>    TCDC=Total cloud cover [%]
>
>    timerange 1 P1 0 P2 0 TimeU 1  nx 199 ny 149 GDS grid 3
num_in_ave 0 missing 0
>
>    center 57 subcenter 2 process 36 Table 2 scan: WE:SN winds(grid)
>
>    Lambert Conf: Lat1 0.888000 Lon1 72.538000 Lov 105.000000
>
>        Latin1 60.000000 Latin2 30.000000 LatSP -90.000000 LonSP
0.000000
>
>        North Pole (199 x 149) Dx 45.000000 Dy 45.000000 scan 64 mode
8
>
>    min/max data 0 93  num bits 7  BDS_Ref 0  DecScale 0 BinScale 0
>
>
>
> The grib header for the cloud analysis file (wwmca) is:
>
>
>
> rec 1:0:date 2012012600 TCDC kpds5=71 kpds6=200 kpds7=0 levels=(0,0)
grid=255 atmos col anl:
>
>    TCDC=Total cloud cover [%]
>
>    timerange 0 P1 0 P2 0 TimeU 1  nx 199 ny 149 GDS grid 3
num_in_ave 0 missing 0
>
>    center 57 subcenter 2 process 32 Table 128 scan: WE:SN
winds(grid)
>
>    Lambert Conf: Lat1 0.888000 Lon1 72.537000 Lov 105.000000
>
>        Latin1 60.000000 Latin2 30.000000 LatSP 0.000000 LonSP
0.000000
>
>        North Pole (199 x 149) Dx 45.000000 Dy 45.000000 scan 64 mode
8
>
>    min/max data 0 100  num bits 7  BDS_Ref 0  DecScale 0 BinScale 0
>
>
>
> The only difference I noticed was in the thousands field of the
Lon1.  Would that cause this error.
>
>
>
> The command line I used is:
>
>
>
> /home/qcteam/METv3.0.1/bin/mode
/h/data/global/WXQC/data/dcf/t06a/fcsts_t06a_2012012600_000.GR1
/h/data/global/WXQC/data/temp/met/cloud_obs_file_71
/home/qcteam/METv3.0.1/data/config/DCFModeConfig_45 -outdir
/h/data/global/WXQC/data/met/cloud/t06/
>
>
>
>
>
> I have attached the to grib files and config file.
>
>
>
> Thanks
>
>
>
> Bob Craig
>
>
>
>
>

------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #53051] Cloud File Mismatch
From: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Time: Thu Jan 26 13:40:55 2012

I was using copygb to convert but my config file had the 1000ths
difference in longitude.  That was the issue.

thanks
-----Original Message-----
From: John Halley Gotway via RT
[mailto:met_help at ucar.edu] 
Sent: Thursday, January 26, 2012 2:33 PM
To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Cc: Sittel, Matthew C CTR
USAF AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu #53051] Cloud File
Mismatch

Bob,

Yes, the difference in the thousanths place of the
Lon1 parameter would cause this error to occur.  When checking for
grid equality, MET is checking each of the components that define the
grid.  To 
resolve this, I'd suggest running one of the GRIB files
through coypgb to regrid it to the other's domain.  An example of
using copygb to regrid a lambert conformal grid can be found on the
MET 
website in the online tutorial:
http://www.dtcenter.org/met/users/support/online_tutorial/METv3.0/copygb/run3.php
Alternatively, if you decide that you don't want MET to do this strict
of grid checking, you could modify the source code to use different
logic.  Checking for grid equality is done in the file
"grid_base.cc" in this function:
    bool operator==(const Grid & g1,
const Grid & g2)

In earlier versions of MET, we simply checked that
Nx and Ny were same, but we ended up missing small differences like
the one you have in your data!

Hope that helps clarify.

John
Halley Gotway
met_help at ucar.edu

On 01/26/2012 01:22 PM, Craig,
Robert J Civ USAF AFWA 16 WS/WXN via RT wrote:
>
> Thu Jan 26
13:22:23 2012: Request 53051 was acted upon.
> Transaction: Ticket
created by Robert.Craig at offutt.af.mil
>         Queue: met_help
>
Subject: Cloud File Mismatch
>         Owner: Nobody
>
Requestors: Robert.Craig at offutt.af.mil
>        Status: new
>
Ticket<URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=53051>
>
>
> Hi guys, when using MODE to compare a cloud analysis to cloud
forecast, I am getting the following error:
>
>
>
> ERROR:
process_fcst_obs_files() ->  The forecast and observation grids do not
match.
>
>
>
> The grib header of the dcf cloud forecast is:
>
>
>
> rec 1:0:date 2012012600 TCDC kpds5=71 kpds6=200 kpds7=0
levels=(0,0) grid=0 atmos col 0hr fcst:
>
>    TCDC=Total cloud
cover [%]
>
>    timerange 1 P1 0 P2 0 TimeU 1  nx 199 ny 149 GDS
grid 3 num_in_ave 0 missing 0
>
>    center 57 subcenter 2 process
36 Table 2 scan: WE:SN winds(grid)
>
>    Lambert Conf: Lat1
0.888000 Lon1 72.538000 Lov 105.000000
>
>        Latin1 60.000000
Latin2 30.000000 LatSP -90.000000 LonSP 0.000000
>
>        North
Pole (199 x 149) Dx 45.000000 Dy 45.000000 scan 64 mode 8
>
>
min/max data 0 93  num bits 7  BDS_Ref 0  DecScale 0 BinScale 0
>
>
>
> The grib header for the cloud analysis file (wwmca) is:
>
>
>
> rec 1:0:date 2012012600 TCDC kpds5=71 kpds6=200 kpds7=0 levels=(0,0)
grid=255 atmos col anl:
>
>    TCDC=Total cloud cover [%]
>
>
timerange 0 P1 0 P2 0 TimeU 1  nx 199 ny 149 GDS grid 3 num_in_ave 0
missing 0
>
>    center 57 subcenter 2 process 32 Table 128 scan:
WE:SN winds(grid)
>
>    Lambert Conf: Lat1 0.888000 Lon1 72.537000
Lov 105.000000
>
>        Latin1 60.000000 Latin2 30.000000 LatSP
0.000000 LonSP 0.000000
>
>        North Pole (199 x 149) Dx
45.000000 Dy 45.000000 scan 64 mode 8
>
>    min/max data 0 100  num
bits 7  BDS_Ref 0  DecScale 0 BinScale 0
>
>
>
> The only
difference I noticed was in the thousands field of the Lon1.  Would
that cause this error.
>
>
>
> The command line I used is:
>
>
>
> /home/qcteam/METv3.0.1/bin/mode
/h/data/global/WXQC/data/dcf/t06a/fcsts_t06a_2012012600_000.GR1
/h/data/global/WXQC/data/temp/met/cloud_obs_file_71
/home/qcteam/METv3.0.1/data/config/DCFModeConfig_45 -outdir
/h/data/global/WXQC/data/met/cloud/t06/
>
>
>
>
>
> I have
attached the to grib files and config file.
>
>
>
> Thanks
>
>
>
> Bob Craig
>
>
>
>
>

------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #53051] Cloud File Mismatch
From: John Halley Gotway
Time: Thu Jan 26 13:44:48 2012

Bob,

Glad you were able to figure it out.

I'll resolve the ticket.

Thanks,
John

On 01/26/2012 01:40 PM, Craig, Robert J Civ USAF AFWA 16 WS/WXN via RT
wrote:
> I was using copygb to convert but my config file had the 1000ths
difference in longitude.  That was the issue.

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


More information about the Met_help mailing list