[Met_help] [rt.rap.ucar.edu #80868] History for Where is my climo

John Halley Gotway via RT met_help at ucar.edu
Fri Jun 16 11:20:11 MDT 2017


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

John, I am so close to getting this working I can taste it - however one more hurdle.  When I run the command below, I get a the output below. I highlighted in red the concerning message.



/h/WXQC/met-6.0/bin/grid_stat h/data/global/WXQC/data/geps_verify/qcsave/2017060500/grib2.2017060500.0024 /h/data/global/WXQC/data/met/tmp/941504fa-76aa-41e7-a1b3-a0057a209172 /h/WXQC/met-6.0/data/config/GridStatConfig_geps_updated -outdir /h/data/global/WXQC/data/met/gdstat/geps -v 3



['/h/WXQC/met-6.0/bin/grid_stat', '/h/data/global/WXQC/data/geps_verify/qcsave/2017060500/grib2.2017060500.0024', '/h/data/global/WXQC/data/met/tmp/941504fa-76aa-41e7-a1b3-a0057a209172', '/h/WXQC/met-6.0/data/config/GridStatConfig_geps_updated', '-outdir', '/h/data/global/WXQC/data/met/gdstat/geps', '-v', '3']

DEBUG 1: Reading user-defined grib1 MET_GRIB_TABLES file: /h/data/global/WXQC/data/met/config/grib2/grib1_af_131_57.txt

DEBUG 1: Reading user-defined grib2 MET_GRIB_TABLES file: /h/data/global/WXQC/data/met/config/grib2/grib2_geps.txt

DEBUG 1: Default Config File: /h/WXQC/met-6.0/share/met/config/GridStatConfig_default

DEBUG 1: User Config File: /h/WXQC/met-6.0/data/config/GridStatConfig_geps_updated

DEBUG 3: Use the matching forecast and observation grids.

DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 360 Ny: 181 lat_ll: -90.000 lon_ll: -0.000 delta_lat: 1.000 delta_lon: 1.000

GSL_RNG_TYPE=mt19937

GSL_RNG_SEED=815022681

DEBUG 1: Forecast File: /h/data/global/WXQC/data/geps_verify/qcsave/2017060500/grib2.2017060500.0024

DEBUG 1: Observation File: /h/data/global/WXQC/data/met/tmp/941504fa-76aa-41e7-a1b3-a0057a209172

DEBUG 3: Reading forecast data for PWS15/Z10.

DEBUG 3: MetGrib2DataFile::data_plane() - Found exact match for 'PWS15/Z10' in GRIB2 record 242 field 1 of GRIB2 file '/h/data/global/WXQC/data/geps_verify/qcsave/2017060500/grib2.2017060500.0024'

DEBUG 3: Reading observation data for PWS15(*,*).

DEBUG 3: Found 0 climatology mean field(s) for forecast PWS15/Z10.

DEBUG 2: Processing masking regions.

DEBUG 3: Processing grid mask: FULL

DEBUG 2:

DEBUG 2: --------------------------------------------------------------------------------

DEBUG 2:

DEBUG 3: Smoothing field using the NEAREST(1) RectangularTemplate interpolation method.

DEBUG 3: Smoothing field using the NEAREST(1) RectangularTemplate interpolation method.

DEBUG 2: Processing PWS15/Z10 versus PWS15(*,*), for interpolation method NEAREST(1), over region FULL, using 65160 pairs.

DEBUG 2:

DEBUG 2: --------------------------------------------------------------------------------

DEBUG 2:

DEBUG 1: Output file: /h/data/global/WXQC/data/met/gdstat/geps/grid_stat_240000L_20170606_000000V.stat

DEBUG 1: Output file: /h/data/global/WXQC/data/met/gdstat/geps/grid_stat_240000L_20170606_000000V_pct.txt

DEBUG 1: Output file: /h/data/global/WXQC/data/met/gdstat/geps/grid_stat_240000L_20170606_000000V_pstd.txt

DEBUG 1: Output file: /h/data/global/WXQC/data/met/gdstat/geps/grid_stat_240000L_20170606_000000V_pairs.nc



I have three questions:



1)      why isn't the code finding the climatology mean field?  The netCDF climo file is present and has data so why doesn't the software use it?



2)      If it does find the climo file and process it, will the climo data end up in the pairs.nc as another netCDF variable?



3)      In the config file, I have pstd and pct turned on but when these files are created they are empty, why?






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

Subject: Where is my climo
From: John Halley Gotway
Time: Thu Jun 15 10:24:26 2017

Bob,

Looking in the config file you sent previously I see the following
setting:

climo_mean = {

   field     = "PWS15";

The "field" setting there should look the same as it is in the "fcst"
and
"obs" sections...

climo_mean = {

   field = [

      {

        name = "PWS15";

        level = [ "(*,*)" ];

...

Make sense?

Thanks,
John

On Thu, Jun 15, 2017 at 8:27 AM, robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:

>
> Thu Jun 15 08:27:59 2017: Request 80868 was acted upon.
> Transaction: Ticket created by robert.craig.2 at us.af.mil
>        Queue: met_help
>      Subject: Where is my climo
>        Owner: Nobody
>   Requestors: robert.craig.2 at us.af.mil
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=80868 >
>
>
> John, I am so close to getting this working I can taste it - however
one
> more hurdle.  When I run the command below, I get a the output
below. I
> highlighted in red the concerning message.
>
>
>
> /h/WXQC/met-6.0/bin/grid_stat
h/data/global/WXQC/data/geps_verify/qcsave/
> 2017060500/grib2.2017060500.0024 /h/data/global/WXQC/data/met/
> tmp/941504fa-76aa-41e7-a1b3-a0057a209172 /h/WXQC/met-
6.0/data/config/GridStatConfig_geps_updated
> -outdir /h/data/global/WXQC/data/met/gdstat/geps -v 3
>
>
>
> ['/h/WXQC/met-6.0/bin/grid_stat', '/h/data/global/WXQC/data/
> geps_verify/qcsave/2017060500/grib2.2017060500.0024',
> '/h/data/global/WXQC/data/met/tmp/941504fa-76aa-41e7-a1b3-
a0057a209172',
> '/h/WXQC/met-6.0/data/config/GridStatConfig_geps_updated', '-
outdir',
> '/h/data/global/WXQC/data/met/gdstat/geps', '-v', '3']
>
> DEBUG 1: Reading user-defined grib1 MET_GRIB_TABLES file:
> /h/data/global/WXQC/data/met/config/grib2/grib1_af_131_57.txt
>
> DEBUG 1: Reading user-defined grib2 MET_GRIB_TABLES file:
> /h/data/global/WXQC/data/met/config/grib2/grib2_geps.txt
>
> DEBUG 1: Default Config File: /h/WXQC/met-6.0/share/met/
> config/GridStatConfig_default
>
> DEBUG 1: User Config File: /h/WXQC/met-6.0/data/config/
> GridStatConfig_geps_updated
>
> DEBUG 3: Use the matching forecast and observation grids.
>
> DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 360 Ny: 181
lat_ll:
> -90.000 lon_ll: -0.000 delta_lat: 1.000 delta_lon: 1.000
>
> GSL_RNG_TYPE=mt19937
>
> GSL_RNG_SEED=815022681
>
> DEBUG 1: Forecast File: /h/data/global/WXQC/data/geps_
> verify/qcsave/2017060500/grib2.2017060500.0024
>
> DEBUG 1: Observation File: /h/data/global/WXQC/data/met/
> tmp/941504fa-76aa-41e7-a1b3-a0057a209172
>
> DEBUG 3: Reading forecast data for PWS15/Z10.
>
> DEBUG 3: MetGrib2DataFile::data_plane() - Found exact match for
> 'PWS15/Z10' in GRIB2 record 242 field 1 of GRIB2 file
> '/h/data/global/WXQC/data/geps_verify/qcsave/2017060500/
> grib2.2017060500.0024'
>
> DEBUG 3: Reading observation data for PWS15(*,*).
>
> DEBUG 3: Found 0 climatology mean field(s) for forecast PWS15/Z10.
>
> DEBUG 2: Processing masking regions.
>
> DEBUG 3: Processing grid mask: FULL
>
> DEBUG 2:
>
> DEBUG 2:
------------------------------------------------------------
> --------------------
>
> DEBUG 2:
>
> DEBUG 3: Smoothing field using the NEAREST(1) RectangularTemplate
> interpolation method.
>
> DEBUG 3: Smoothing field using the NEAREST(1) RectangularTemplate
> interpolation method.
>
> DEBUG 2: Processing PWS15/Z10 versus PWS15(*,*), for interpolation
method
> NEAREST(1), over region FULL, using 65160 pairs.
>
> DEBUG 2:
>
> DEBUG 2:
------------------------------------------------------------
> --------------------
>
> DEBUG 2:
>
> DEBUG 1: Output file: /h/data/global/WXQC/data/met/
> gdstat/geps/grid_stat_240000L_20170606_000000V.stat
>
> DEBUG 1: Output file: /h/data/global/WXQC/data/met/
> gdstat/geps/grid_stat_240000L_20170606_000000V_pct.txt
>
> DEBUG 1: Output file: /h/data/global/WXQC/data/met/
> gdstat/geps/grid_stat_240000L_20170606_000000V_pstd.txt
>
> DEBUG 1: Output file: /h/data/global/WXQC/data/met/
> gdstat/geps/grid_stat_240000L_20170606_000000V_pairs.nc
>
>
>
> I have three questions:
>
>
>
> 1)      why isn't the code finding the climatology mean field?  The
netCDF
> climo file is present and has data so why doesn't the software use
it?
>
>
>
> 2)      If it does find the climo file and process it, will the
climo data
> end up in the pairs.nc as another netCDF variable?
>
>
>
> 3)      In the config file, I have pstd and pct turned on but when
these
> files are created they are empty, why?
>
>
>
>
>
>

------------------------------------------------
Subject: RE: [Non-DoD Source] Re: [rt.rap.ucar.edu #80868] Where is my climo
From: robert.craig.2 at us.af.mil
Time: Thu Jun 15 11:37:25 2017

Thanks, that solved the lack of climo data but I still don't have any
output in the PSTD or PCT files.   Did you get output from the
datafiles I sent?

thanks

-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Thursday, June 15, 2017 11:24 AM
To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #80868] Where is my
climo

Bob,

Looking in the config file you sent previously I see the following
setting:

climo_mean = {

   field     = "PWS15";

The "field" setting there should look the same as it is in the "fcst"
and "obs" sections...

climo_mean = {

   field = [

      {

        name = "PWS15";

        level = [ "(*,*)" ];

...

Make sense?

Thanks,
John

On Thu, Jun 15, 2017 at 8:27 AM, robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:

>
> Thu Jun 15 08:27:59 2017: Request 80868 was acted upon.
> Transaction: Ticket created by robert.craig.2 at us.af.mil
>        Queue: met_help
>      Subject: Where is my climo
>        Owner: Nobody
>   Requestors: robert.craig.2 at us.af.mil
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=80868
> >
>
>
> John, I am so close to getting this working I can taste it - however
> one more hurdle.  When I run the command below, I get a the output
> below. I highlighted in red the concerning message.
>
>
>
> /h/WXQC/met-6.0/bin/grid_stat
> h/data/global/WXQC/data/geps_verify/qcsave/
> 2017060500/grib2.2017060500.0024 /h/data/global/WXQC/data/met/
> tmp/941504fa-76aa-41e7-a1b3-a0057a209172
> /h/WXQC/met-6.0/data/config/GridStatConfig_geps_updated
> -outdir /h/data/global/WXQC/data/met/gdstat/geps -v 3
>
>
>
> ['/h/WXQC/met-6.0/bin/grid_stat', '/h/data/global/WXQC/data/
> geps_verify/qcsave/2017060500/grib2.2017060500.0024',
> '/h/data/global/WXQC/data/met/tmp/941504fa-76aa-41e7-a1b3-
a0057a209172
> ', '/h/WXQC/met-6.0/data/config/GridStatConfig_geps_updated',
> '-outdir', '/h/data/global/WXQC/data/met/gdstat/geps', '-v', '3']
>
> DEBUG 1: Reading user-defined grib1 MET_GRIB_TABLES file:
> /h/data/global/WXQC/data/met/config/grib2/grib1_af_131_57.txt
>
> DEBUG 1: Reading user-defined grib2 MET_GRIB_TABLES file:
> /h/data/global/WXQC/data/met/config/grib2/grib2_geps.txt
>
> DEBUG 1: Default Config File: /h/WXQC/met-6.0/share/met/
> config/GridStatConfig_default
>
> DEBUG 1: User Config File: /h/WXQC/met-6.0/data/config/
> GridStatConfig_geps_updated
>
> DEBUG 3: Use the matching forecast and observation grids.
>
> DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 360 Ny: 181
lat_ll:
> -90.000 lon_ll: -0.000 delta_lat: 1.000 delta_lon: 1.000
>
> GSL_RNG_TYPE=mt19937
>
> GSL_RNG_SEED=815022681
>
> DEBUG 1: Forecast File: /h/data/global/WXQC/data/geps_
> verify/qcsave/2017060500/grib2.2017060500.0024
>
> DEBUG 1: Observation File: /h/data/global/WXQC/data/met/
> tmp/941504fa-76aa-41e7-a1b3-a0057a209172
>
> DEBUG 3: Reading forecast data for PWS15/Z10.
>
> DEBUG 3: MetGrib2DataFile::data_plane() - Found exact match for
> 'PWS15/Z10' in GRIB2 record 242 field 1 of GRIB2 file
> '/h/data/global/WXQC/data/geps_verify/qcsave/2017060500/
> grib2.2017060500.0024'
>
> DEBUG 3: Reading observation data for PWS15(*,*).
>
> DEBUG 3: Found 0 climatology mean field(s) for forecast PWS15/Z10.
>
> DEBUG 2: Processing masking regions.
>
> DEBUG 3: Processing grid mask: FULL
>
> DEBUG 2:
>
> DEBUG 2:
------------------------------------------------------------
> --------------------
>
> DEBUG 2:
>
> DEBUG 3: Smoothing field using the NEAREST(1) RectangularTemplate
> interpolation method.
>
> DEBUG 3: Smoothing field using the NEAREST(1) RectangularTemplate
> interpolation method.
>
> DEBUG 2: Processing PWS15/Z10 versus PWS15(*,*), for interpolation
> method NEAREST(1), over region FULL, using 65160 pairs.
>
> DEBUG 2:
>
> DEBUG 2:
------------------------------------------------------------
> --------------------
>
> DEBUG 2:
>
> DEBUG 1: Output file: /h/data/global/WXQC/data/met/
> gdstat/geps/grid_stat_240000L_20170606_000000V.stat
>
> DEBUG 1: Output file: /h/data/global/WXQC/data/met/
> gdstat/geps/grid_stat_240000L_20170606_000000V_pct.txt
>
> DEBUG 1: Output file: /h/data/global/WXQC/data/met/
> gdstat/geps/grid_stat_240000L_20170606_000000V_pstd.txt
>
> DEBUG 1: Output file: /h/data/global/WXQC/data/met/
> gdstat/geps/grid_stat_240000L_20170606_000000V_pairs.nc
>
>
>
> I have three questions:
>
>
>
> 1)      why isn't the code finding the climatology mean field?  The
netCDF
> climo file is present and has data so why doesn't the software use
it?
>
>
>
> 2)      If it does find the climo file and process it, will the
climo data
> end up in the pairs.nc as another netCDF variable?
>
>
>
> 3)      In the config file, I have pstd and pct turned on but when
these
> files are created they are empty, why?
>
>
>
>
>
>




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


More information about the Met_help mailing list