[Met_help] [rt.rap.ucar.edu #90587] History for Problem with cyber corrected MET 8.1

John Halley Gotway via RT met_help at ucar.edu
Tue Jul 9 12:06:42 MDT 2019


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

John, I am getting different behavior in grid stat if I run the 8.1_beta1 version (not cyber secure) and the recently released cyber secure.  The command line I ran is:

/h/WXQC/met-8.1_beta1/bin/grid_stat /h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042300.0012.EPOSTOUT /h/data/global/WXQC/data/met/nc_obs/wwmca/201904231200_p5.nc  /h/data/global/WXQC/data/met/config/met_config/GridStatConfig_geps_tc_updated -outdir /h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC -v 3

Below is the listing for the 8.1_beta1 run.  The key area in in red.

DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Processing PROB(TCDC<20.000)/L0/PROB versus Cloud_PctSFC, for smoothing method NEAREST(1), over region FULL, using 258836 pairs.
DEBUG 2: Computing Probabilistic Statistics.
DEBUG 3: MetGrib2DataFile::data_plane() - Found exact match for 'PROB(TCDC>80.000)/L0/PROB' in GRIB2 record 170 field 1 of GRIB2 file '/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042300.0012.EPOSTOUT'
DEBUG 3: Rescaling probabilistic field from [0,100] to [0,1].
DEBUG 3: Reading forecast data for PROB(TCDC>80.000)/L0/PROB.
DEBUG 3: Reading observation data for Cloud_PctSFC.
WARNING:
WARNING: Resetting the regridding width from 2 to 1 for regridding method "NEAREST".
WARNING:
WARNING:
WARNING: Resetting the regridding width from 2 to 1 for regridding method "NEAREST".
WARNING:
WARNING:
WARNING: Resetting the regridding width from 2 to 1 for regridding method "NEAREST".
WARNING:
DEBUG 3: Applying conversion function.
DEBUG 1: Regridding climatology fmcf8100(*,*) to the verification grid.
DEBUG 3: For forecast valid at 20190423_120000, found 1 climatology field(s) valid at 20190423_120000.
DEBUG 3: Found 1 climatology fields.
DEBUG 3: Found 1 climatology mean and 0 climatology standard deviation field(s) for forecast PROB(TCDC>80.000)/L0/PROB.
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Processing PROB(TCDC>80.000)/L0/PROB versus Cloud_PctSFC, for smoothing method NEAREST(1), over region FULL, using 258836 pairs.
DEBUG 2: Computing Probabilistic Statistics.
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 1: Output file: /h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps_all_120000L_20190423_120000V.stat
DEBUG 1: Output file: /h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps_all_120000L_20190423_120000V_pct.txt
DEBUG 1: Output file: /h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps_all_120000L_20190423_120000V_pstd.txt
DEBUG 1: Output file: /h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps_all_120000L_20190423_120000V_pairs.nc

When I run the cyber secure version:

/h/WXQC/met-8.1/bin/grid_stat /h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042300.0012.EPOSTOUT /h/data/global/WXQC/data/met/nc_obs/wwmca/201904231200_p5.nc  /h/data/global/WXQC/data/met/config/met_config/GridStatConfig_geps_tc_updated -outdir /h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC -v 3

I get the following listing:

DEBUG 3: Use the forecast grid.
DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 720 Ny: 361 lat_ll: -90.000 lon_ll: -0.000 delta_lat: 0.500 delta_lon: 0.500
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=1667564125
DEBUG 1: Forecast File: /h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042300.0012.EPOSTOUT
DEBUG 1: Observation File: /h/data/global/WXQC/data/met/nc_obs/wwmca/201904231200_p5.nc
DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for 'PROB(TCDC<20.000)/L0/PROB' in GRIB2 record 169 field 1 of GRIB2 file '/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042300.0012.EPOSTOUT'
DEBUG 3: Rescaling probabilistic field from [0,100] to [0,1].
DEBUG 3: Reading forecast data for PROB(TCDC<20.000)/L0/PROB.
DEBUG 3: Reading observation data for Cloud_PctSFC.
WARNING:
WARNING: Resetting the regridding width from 2 to 1 for regridding method "NEAREST".
WARNING:
WARNING:
WARNING: Resetting the regridding width from 2 to 1 for regridding method "NEAREST".
WARNING:
WARNING: MetNcCFDataFile::collect_time_offsets(VarInfo &) -> Not found time out of 1.
ERROR  :
ERROR  : MetNcCFDataFile::collect_time_offsets(VarInfo &) -> variable "fmcf0020" has unknown error (999).
ERROR  :

The fcst, ob, and climo files will be sent in AMRDEC Safe to you.

Thanks
Bob Craig


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

Subject: Problem with cyber corrected MET 8.1
From: John Halley Gotway
Time: Mon Jun 17 12:19:11 2019

Hi Bob,

I was able to grab those files and reproduce the behavior you're
seeing.
The issue specifically is in the "porclds*.nc" files.  While you've
listed
it as being CF-compliant, there are a couple of issues there.

I understand that the first beta release of met-8.1 was able to read
these
files.  Subsequent development in that library led us to firm up the
requirements.  The immediate problem is that the data variable is
missing
the time dimension.  After adding the time dimension, it worked as it
had
before:

OLD:
        int fmcf8100(lon, lat) ;
NEW:
        int fmcf8100(time, lon, lat) ;

Please give that a shot if possible and let me know how it goes.

Thanks,
John

On Thu, Jun 13, 2019 at 9:19 AM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:

>
> Thu Jun 13 09:19:44 2019: Request 90587 was acted upon.
> Transaction: Ticket created by robert.craig.2 at us.af.mil
>        Queue: met_help
>      Subject: Problem with cyber corrected MET 8.1
>        Owner: Nobody
>   Requestors: robert.craig.2 at us.af.mil
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=90587 >
>
>
> John, I am getting different behavior in grid stat if I run the
8.1_beta1
> version (not cyber secure) and the recently released cyber secure.
The
> command line I ran is:
>
> /h/WXQC/met-8.1_beta1/bin/grid_stat
>
/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042300.0012.EPOSTOUT
> /h/data/global/WXQC/data/met/nc_obs/wwmca/201904231200_p5.nc
>
/h/data/global/WXQC/data/met/config/met_config/GridStatConfig_geps_tc_updated
> -outdir /h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC -v 3
>
> Below is the listing for the 8.1_beta1 run.  The key area in in red.
>
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Processing PROB(TCDC<20.000)/L0/PROB versus Cloud_PctSFC,
for
> smoothing method NEAREST(1), over region FULL, using 258836 pairs.
> DEBUG 2: Computing Probabilistic Statistics.
> DEBUG 3: MetGrib2DataFile::data_plane() - Found exact match for
> 'PROB(TCDC>80.000)/L0/PROB' in GRIB2 record 170 field 1 of GRIB2
file
>
'/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042300.0012.EPOSTOUT'
> DEBUG 3: Rescaling probabilistic field from [0,100] to [0,1].
> DEBUG 3: Reading forecast data for PROB(TCDC>80.000)/L0/PROB.
> DEBUG 3: Reading observation data for Cloud_PctSFC.
> WARNING:
> WARNING: Resetting the regridding width from 2 to 1 for regridding
method
> "NEAREST".
> WARNING:
> WARNING:
> WARNING: Resetting the regridding width from 2 to 1 for regridding
method
> "NEAREST".
> WARNING:
> WARNING:
> WARNING: Resetting the regridding width from 2 to 1 for regridding
method
> "NEAREST".
> WARNING:
> DEBUG 3: Applying conversion function.
> DEBUG 1: Regridding climatology fmcf8100(*,*) to the verification
grid.
> DEBUG 3: For forecast valid at 20190423_120000, found 1 climatology
> field(s) valid at 20190423_120000.
> DEBUG 3: Found 1 climatology fields.
> DEBUG 3: Found 1 climatology mean and 0 climatology standard
deviation
> field(s) for forecast PROB(TCDC>80.000)/L0/PROB.
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Processing PROB(TCDC>80.000)/L0/PROB versus Cloud_PctSFC,
for
> smoothing method NEAREST(1), over region FULL, using 258836 pairs.
> DEBUG 2: Computing Probabilistic Statistics.
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 1: Output file:
>
/h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps_all_120000L_20190423_120000V.stat
> DEBUG 1: Output file:
>
/h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps_all_120000L_20190423_120000V_pct.txt
> DEBUG 1: Output file:
>
/h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps_all_120000L_20190423_120000V_pstd.txt
> DEBUG 1: Output file:
>
/h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps_all_120000L_20190423_120000V_pairs.nc
>
> When I run the cyber secure version:
>
> /h/WXQC/met-8.1/bin/grid_stat
>
/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042300.0012.EPOSTOUT
> /h/data/global/WXQC/data/met/nc_obs/wwmca/201904231200_p5.nc
>
/h/data/global/WXQC/data/met/config/met_config/GridStatConfig_geps_tc_updated
> -outdir /h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC -v 3
>
> I get the following listing:
>
> DEBUG 3: Use the forecast grid.
> DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 720 Ny: 361
lat_ll:
> -90.000 lon_ll: -0.000 delta_lat: 0.500 delta_lon: 0.500
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=1667564125
> DEBUG 1: Forecast File:
>
/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042300.0012.EPOSTOUT
> DEBUG 1: Observation File:
/h/data/global/WXQC/data/met/nc_obs/wwmca/
> 201904231200_p5.nc
> DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> 'PROB(TCDC<20.000)/L0/PROB' in GRIB2 record 169 field 1 of GRIB2
file
>
'/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042300.0012.EPOSTOUT'
> DEBUG 3: Rescaling probabilistic field from [0,100] to [0,1].
> DEBUG 3: Reading forecast data for PROB(TCDC<20.000)/L0/PROB.
> DEBUG 3: Reading observation data for Cloud_PctSFC.
> WARNING:
> WARNING: Resetting the regridding width from 2 to 1 for regridding
method
> "NEAREST".
> WARNING:
> WARNING:
> WARNING: Resetting the regridding width from 2 to 1 for regridding
method
> "NEAREST".
> WARNING:
> WARNING: MetNcCFDataFile::collect_time_offsets(VarInfo &) -> Not
found
> time out of 1.
> ERROR  :
> ERROR  : MetNcCFDataFile::collect_time_offsets(VarInfo &) ->
variable
> "fmcf0020" has unknown error (999).
> ERROR  :
>
> The fcst, ob, and climo files will be sent in AMRDEC Safe to you.
>
> Thanks
> Bob Craig
>
>

------------------------------------------------
Subject: RE: [Non-DoD Source] Re: [rt.rap.ucar.edu #90587] Problem with cyber corrected MET 8.1
From: robert.craig.2 at us.af.mil
Time: Tue Jun 18 13:27:03 2019

John, what is the ncatted command to add a dimension?

Thanks
Bob

-----Original Message-----
From: John Halley Gotway via RT <met_help at ucar.edu>
Sent: Monday, June 17, 2019 1:19 PM
To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
Cc: SITTEL, MATTHEW C CTR USAF AFMC AFLCMC/HBAW-OL
<matthew.sittel.ctr at us.af.mil>; SEDLACEK, WILLIAM T GS-12 USAF ACC 16
WS/WXEV <william.sedlacek at us.af.mil>
Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #90587] Problem with
cyber corrected MET 8.1

Hi Bob,

I was able to grab those files and reproduce the behavior you're
seeing.
The issue specifically is in the "porclds*.nc" files.  While you've
listed it as being CF-compliant, there are a couple of issues there.

I understand that the first beta release of met-8.1 was able to read
these files.  Subsequent development in that library led us to firm up
the requirements.  The immediate problem is that the data variable is
missing the time dimension.  After adding the time dimension, it
worked as it had
before:

OLD:
        int fmcf8100(lon, lat) ;
NEW:
        int fmcf8100(time, lon, lat) ;

Please give that a shot if possible and let me know how it goes.

Thanks,
John

On Thu, Jun 13, 2019 at 9:19 AM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:

>
> Thu Jun 13 09:19:44 2019: Request 90587 was acted upon.
> Transaction: Ticket created by robert.craig.2 at us.af.mil
>        Queue: met_help
>      Subject: Problem with cyber corrected MET 8.1
>        Owner: Nobody
>   Requestors: robert.craig.2 at us.af.mil
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=90587
> >
>
>
> John, I am getting different behavior in grid stat if I run the
> 8.1_beta1 version (not cyber secure) and the recently released cyber
> secure.  The command line I ran is:
>
> /h/WXQC/met-8.1_beta1/bin/grid_stat
>
/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042
> 300.0012.EPOSTOUT
> /h/data/global/WXQC/data/met/nc_obs/wwmca/201904231200_p5.nc
>
/h/data/global/WXQC/data/met/config/met_config/GridStatConfig_geps_tc_
> updated -outdir
> /h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC -v 3
>
> Below is the listing for the 8.1_beta1 run.  The key area in in red.
>
> DEBUG 2:
>
----------------------------------------------------------------------
> ----------
> DEBUG 2:
> DEBUG 2: Processing PROB(TCDC<20.000)/L0/PROB versus Cloud_PctSFC,
for
> smoothing method NEAREST(1), over region FULL, using 258836 pairs.
> DEBUG 2: Computing Probabilistic Statistics.
> DEBUG 3: MetGrib2DataFile::data_plane() - Found exact match for
> 'PROB(TCDC>80.000)/L0/PROB' in GRIB2 record 170 field 1 of GRIB2
file
>
'/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042300.0012.EPOSTOUT'
> DEBUG 3: Rescaling probabilistic field from [0,100] to [0,1].
> DEBUG 3: Reading forecast data for PROB(TCDC>80.000)/L0/PROB.
> DEBUG 3: Reading observation data for Cloud_PctSFC.
> WARNING:
> WARNING: Resetting the regridding width from 2 to 1 for regridding
> method "NEAREST".
> WARNING:
> WARNING:
> WARNING: Resetting the regridding width from 2 to 1 for regridding
> method "NEAREST".
> WARNING:
> WARNING:
> WARNING: Resetting the regridding width from 2 to 1 for regridding
> method "NEAREST".
> WARNING:
> DEBUG 3: Applying conversion function.
> DEBUG 1: Regridding climatology fmcf8100(*,*) to the verification
grid.
> DEBUG 3: For forecast valid at 20190423_120000, found 1 climatology
> field(s) valid at 20190423_120000.
> DEBUG 3: Found 1 climatology fields.
> DEBUG 3: Found 1 climatology mean and 0 climatology standard
deviation
> field(s) for forecast PROB(TCDC>80.000)/L0/PROB.
> DEBUG 2:
> DEBUG 2:
>
----------------------------------------------------------------------
> ----------
> DEBUG 2:
> DEBUG 2: Processing PROB(TCDC>80.000)/L0/PROB versus Cloud_PctSFC,
for
> smoothing method NEAREST(1), over region FULL, using 258836 pairs.
> DEBUG 2: Computing Probabilistic Statistics.
> DEBUG 2:
> DEBUG 2:
>
----------------------------------------------------------------------
> ----------
> DEBUG 2:
> DEBUG 1: Output file:
>
/h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps
> _all_120000L_20190423_120000V.stat
> DEBUG 1: Output file:
>
/h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps
> _all_120000L_20190423_120000V_pct.txt
> DEBUG 1: Output file:
>
/h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps
> _all_120000L_20190423_120000V_pstd.txt
> DEBUG 1: Output file:
>
/h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps
> _all_120000L_20190423_120000V_pairs.nc
>
> When I run the cyber secure version:
>
> /h/WXQC/met-8.1/bin/grid_stat
>
/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042
> 300.0012.EPOSTOUT
> /h/data/global/WXQC/data/met/nc_obs/wwmca/201904231200_p5.nc
>
/h/data/global/WXQC/data/met/config/met_config/GridStatConfig_geps_tc_
> updated -outdir
> /h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC -v 3
>
> I get the following listing:
>
> DEBUG 3: Use the forecast grid.
> DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 720 Ny: 361
lat_ll:
> -90.000 lon_ll: -0.000 delta_lat: 0.500 delta_lon: 0.500
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=1667564125
> DEBUG 1: Forecast File:
>
/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042
> 300.0012.EPOSTOUT DEBUG 1: Observation File:
> /h/data/global/WXQC/data/met/nc_obs/wwmca/
> 201904231200_p5.nc
> DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> 'PROB(TCDC<20.000)/L0/PROB' in GRIB2 record 169 field 1 of GRIB2
file
>
'/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042300.0012.EPOSTOUT'
> DEBUG 3: Rescaling probabilistic field from [0,100] to [0,1].
> DEBUG 3: Reading forecast data for PROB(TCDC<20.000)/L0/PROB.
> DEBUG 3: Reading observation data for Cloud_PctSFC.
> WARNING:
> WARNING: Resetting the regridding width from 2 to 1 for regridding
> method "NEAREST".
> WARNING:
> WARNING:
> WARNING: Resetting the regridding width from 2 to 1 for regridding
> method "NEAREST".
> WARNING:
> WARNING: MetNcCFDataFile::collect_time_offsets(VarInfo &) -> Not
found
> time out of 1.
> ERROR  :
> ERROR  : MetNcCFDataFile::collect_time_offsets(VarInfo &) ->
variable
> "fmcf0020" has unknown error (999).
> ERROR  :
>
> The fcst, ob, and climo files will be sent in AMRDEC Safe to you.
>
> Thanks
> Bob Craig
>
>



------------------------------------------------
Subject: Problem with cyber corrected MET 8.1
From: John Halley Gotway
Time: Tue Jun 18 16:10:38 2019

Bob,

ncatted is a command line utility that's part of the NetCDF Operators
(NCO).  However, what's needed here is adding the time dimension to
the
existing variable, not an attribute.  So ncatted won't do the trick.

To test this myself, I used ncdump and ncgen:

ncdump porclds_04_12z_fmcf8100.nc > porclds_04_12z_fmcf8100.ncdump
... and then edit that ncdump file as follows...
int fmcf8100(time, lon, lat) ;
... and then recreate the NetCDF file ...
ncgen -o porclds_04_12z_fmcf8100_MOD.nc porclds_04_12z_fmcf8100.ncdump

And then met-8.1 can read that modified porclds_04_12z_fmcf8100_MOD.nc
file.

However, in my testing, I did also run into a separate problem.  I
fixed
this for inclusion in the met-8.1.1 bugfix release and it's described
here:
https://github.com/NCAR/MET/issues/1146

Thanks,
John

On Tue, Jun 18, 2019 at 1:27 PM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=90587 >
>
> John, what is the ncatted command to add a dimension?
>
> Thanks
> Bob
>
> -----Original Message-----
> From: John Halley Gotway via RT <met_help at ucar.edu>
> Sent: Monday, June 17, 2019 1:19 PM
> To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
> Cc: SITTEL, MATTHEW C CTR USAF AFMC AFLCMC/HBAW-OL <
> matthew.sittel.ctr at us.af.mil>; SEDLACEK, WILLIAM T GS-12 USAF ACC 16
> WS/WXEV <william.sedlacek at us.af.mil>
> Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #90587] Problem with
cyber
> corrected MET 8.1
>
> Hi Bob,
>
> I was able to grab those files and reproduce the behavior you're
seeing.
> The issue specifically is in the "porclds*.nc" files.  While you've
listed
> it as being CF-compliant, there are a couple of issues there.
>
> I understand that the first beta release of met-8.1 was able to read
these
> files.  Subsequent development in that library led us to firm up the
> requirements.  The immediate problem is that the data variable is
missing
> the time dimension.  After adding the time dimension, it worked as
it had
> before:
>
> OLD:
>         int fmcf8100(lon, lat) ;
> NEW:
>         int fmcf8100(time, lon, lat) ;
>
> Please give that a shot if possible and let me know how it goes.
>
> Thanks,
> John
>
> On Thu, Jun 13, 2019 at 9:19 AM robert.craig.2 at us.af.mil via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Thu Jun 13 09:19:44 2019: Request 90587 was acted upon.
> > Transaction: Ticket created by robert.craig.2 at us.af.mil
> >        Queue: met_help
> >      Subject: Problem with cyber corrected MET 8.1
> >        Owner: Nobody
> >   Requestors: robert.craig.2 at us.af.mil
> >       Status: new
> >  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=90587
> > >
> >
> >
> > John, I am getting different behavior in grid stat if I run the
> > 8.1_beta1 version (not cyber secure) and the recently released
cyber
> > secure.  The command line I ran is:
> >
> > /h/WXQC/met-8.1_beta1/bin/grid_stat
> >
/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042
> > 300.0012.EPOSTOUT
> > /h/data/global/WXQC/data/met/nc_obs/wwmca/201904231200_p5.nc
> >
/h/data/global/WXQC/data/met/config/met_config/GridStatConfig_geps_tc_
> > updated -outdir
> > /h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC -v 3
> >
> > Below is the listing for the 8.1_beta1 run.  The key area in in
red.
> >
> > DEBUG 2:
> >
----------------------------------------------------------------------
> > ----------
> > DEBUG 2:
> > DEBUG 2: Processing PROB(TCDC<20.000)/L0/PROB versus Cloud_PctSFC,
for
> > smoothing method NEAREST(1), over region FULL, using 258836 pairs.
> > DEBUG 2: Computing Probabilistic Statistics.
> > DEBUG 3: MetGrib2DataFile::data_plane() - Found exact match for
> > 'PROB(TCDC>80.000)/L0/PROB' in GRIB2 record 170 field 1 of GRIB2
file
> >
>
'/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042300.0012.EPOSTOUT'
> > DEBUG 3: Rescaling probabilistic field from [0,100] to [0,1].
> > DEBUG 3: Reading forecast data for PROB(TCDC>80.000)/L0/PROB.
> > DEBUG 3: Reading observation data for Cloud_PctSFC.
> > WARNING:
> > WARNING: Resetting the regridding width from 2 to 1 for regridding
> > method "NEAREST".
> > WARNING:
> > WARNING:
> > WARNING: Resetting the regridding width from 2 to 1 for regridding
> > method "NEAREST".
> > WARNING:
> > WARNING:
> > WARNING: Resetting the regridding width from 2 to 1 for regridding
> > method "NEAREST".
> > WARNING:
> > DEBUG 3: Applying conversion function.
> > DEBUG 1: Regridding climatology fmcf8100(*,*) to the verification
grid.
> > DEBUG 3: For forecast valid at 20190423_120000, found 1
climatology
> > field(s) valid at 20190423_120000.
> > DEBUG 3: Found 1 climatology fields.
> > DEBUG 3: Found 1 climatology mean and 0 climatology standard
deviation
> > field(s) for forecast PROB(TCDC>80.000)/L0/PROB.
> > DEBUG 2:
> > DEBUG 2:
> >
----------------------------------------------------------------------
> > ----------
> > DEBUG 2:
> > DEBUG 2: Processing PROB(TCDC>80.000)/L0/PROB versus Cloud_PctSFC,
for
> > smoothing method NEAREST(1), over region FULL, using 258836 pairs.
> > DEBUG 2: Computing Probabilistic Statistics.
> > DEBUG 2:
> > DEBUG 2:
> >
----------------------------------------------------------------------
> > ----------
> > DEBUG 2:
> > DEBUG 1: Output file:
> >
/h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps
> > _all_120000L_20190423_120000V.stat
> > DEBUG 1: Output file:
> >
/h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps
> > _all_120000L_20190423_120000V_pct.txt
> > DEBUG 1: Output file:
> >
/h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps
> > _all_120000L_20190423_120000V_pstd.txt
> > DEBUG 1: Output file:
> >
/h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC/grid_stat_geps
> > _all_120000L_20190423_120000V_pairs.nc
> >
> > When I run the cyber secure version:
> >
> > /h/WXQC/met-8.1/bin/grid_stat
> >
/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042
> > 300.0012.EPOSTOUT
> > /h/data/global/WXQC/data/met/nc_obs/wwmca/201904231200_p5.nc
> >
/h/data/global/WXQC/data/met/config/met_config/GridStatConfig_geps_tc_
> > updated -outdir
> > /h/data/global/WXQC/data/met/gdstat/mdlob_pairs/geps/TC -v 3
> >
> > I get the following listing:
> >
> > DEBUG 3: Use the forecast grid.
> > DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 720 Ny: 361
lat_ll:
> > -90.000 lon_ll: -0.000 delta_lat: 0.500 delta_lon: 0.500
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=1667564125
> > DEBUG 1: Forecast File:
> >
/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042
> > 300.0012.EPOSTOUT DEBUG 1: Observation File:
> > /h/data/global/WXQC/data/met/nc_obs/wwmca/
> > 201904231200_p5.nc
> > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> > 'PROB(TCDC<20.000)/L0/PROB' in GRIB2 record 169 field 1 of GRIB2
file
> >
>
'/h/data/global/ENSEMBLE/data/post/GEPS/global/2019042300/grib2.2019042300.0012.EPOSTOUT'
> > DEBUG 3: Rescaling probabilistic field from [0,100] to [0,1].
> > DEBUG 3: Reading forecast data for PROB(TCDC<20.000)/L0/PROB.
> > DEBUG 3: Reading observation data for Cloud_PctSFC.
> > WARNING:
> > WARNING: Resetting the regridding width from 2 to 1 for regridding
> > method "NEAREST".
> > WARNING:
> > WARNING:
> > WARNING: Resetting the regridding width from 2 to 1 for regridding
> > method "NEAREST".
> > WARNING:
> > WARNING: MetNcCFDataFile::collect_time_offsets(VarInfo &) -> Not
found
> > time out of 1.
> > ERROR  :
> > ERROR  : MetNcCFDataFile::collect_time_offsets(VarInfo &) ->
variable
> > "fmcf0020" has unknown error (999).
> > ERROR  :
> >
> > The fcst, ob, and climo files will be sent in AMRDEC Safe to you.
> >
> > Thanks
> > Bob Craig
> >
> >
>
>
>
>

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


More information about the Met_help mailing list