[Met_help] [rt.rap.ucar.edu #69314] History for Empty output from Grid_stat

John Halley Gotway via RT met_help at ucar.edu
Wed Oct 15 11:26:57 MDT 2014


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

Dear MET Help,

I'm testing Grid_stat to compare my WRF simulation outputs to CCSM4 outputs.

I created some polygon files that represent subregions of my WRF tropical
channel domain and they reside at
/glade/scratch/jyoo/met-5.0/share/met/poly/tropicalsubregion.

Then, I included those polygon information in the GridStatConfig_ccsm file.

When, I executed grid_stat as below, it seems that only one polygon (NIN)
area was incorporated into the calculation.
Moreover, the output files are still empty.
What am I missing in this test?
Please let me know.
Thank you.

Regards,

Jinwoong Yoo
UNM


[jyoo at yslogin4 analysis]$ ls gridstatout/[jyoo at yslogin4 analysis]$
/glade/p/ral/jnt/MET/MET_releases/met-5.0/bin/grid_stat
/glade/scratch/jyoo/DOMAINS/postprd/lgm/postprd/outcopygb/WRFPRS_d01.00_2latlon
/glade/scratch/jyoo/domain/mean/analysis/LGMCCSM/6hrly/
b40.lgm21ka.1deg.003M.cam2.h3.1870-01.nc GridStatConfig_ccsm -outdir
/glade/scratch/jyoo/domain/mean/analysis/gridstatout
DEBUG 1: Default Config File:
/glade/p/ral/jnt/MET/MET_releases/met-5.0/share/met/config/GridStatConfig_default
DEBUG 1: User Config File: GridStatConfig_ccsm
WARNING:
WARNING: NcCfFile::open() -> could not extract init time from file name
WARNING:
WARNING:
WARNING: NcCfFile::get_grid_from_dimensions() -> Found multiple variables
for latitude, using "lat".
WARNING:
WARNING:
WARNING: NcCfFile::get_grid_from_dimensions() -> Found multiple variables
for longitude, using "lon".
WARNING:
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=59061588
DEBUG 1: Forecast File:
/glade/scratch/jyoo/DOMAINS/postprd/lgm/postprd/outcopygb/WRFPRS_d01.00_2latlon
DEBUG 1: Observation File:
/glade/scratch/jyoo/domain/mean/analysis/LGMCCSM/6hrly/
b40.lgm21ka.1deg.003M.cam2.h3.1870-01.nc
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Processing PRMSL/R247 veDEBUG 2: Processing PRMSL/R247 versus
PSL(0,*,*), for interpolation method UW_MEAN(1), over region NIN, using
4260 pairs.
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for interpolation method
UW_MEAN(1), over region NIN, using 10404 pairs.
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for interpolation method
UW_MEAN(1), over region NIN, using 3264 pairs.
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for interpolation method
UW_MEAN(1), over region NIN, using 2880 pairs.
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for interpolation method
UW_MEAN(1), over region NIN, using 3216 pairs.
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for interpolation method
UW_MEAN(1), over region NIN, using 3264 pairs.
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for interpolation method
UW_MEAN(1), over region NIN, using 6480 pairs.
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
NetCDF: Numeric conversion not representable


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

Subject: Empty output from Grid_stat
From: John Halley Gotway
Time: Wed Oct 08 13:26:39 2014

Jinwoong,

You can specify multiple masking regions as a comma-separated list of
file
names in the config file, like this:
   poly = [
"/glade/scratch/jyoo/met-
5.0/share/met/poly/tropicalsubregion/NAT.poly",

"/glade/scratch/jyoo/met-
5.0/share/met/poly/tropicalsubregion/NEP.poly",

"/glade/scratch/jyoo/met-
5.0/share/met/poly/tropicalsubregion/NIN.poly" ];

I suspect this piece of it is working fine.  It just aborted when it
tried
to write the first set of output to the NetCDF matched pairs output
file.

I'll look into that and let you know.  As I mentioned, this is the
first
time I've seen that error.

Thanks,
John

On Wed, Oct 8, 2014 at 12:55 PM, Jinwoong Yoo via RT
<met_help at ucar.edu>
wrote:

>
> Wed Oct 08 12:55:15 2014: Request 69314 was acted upon.
> Transaction: Ticket created by jinwoong.yoo at gmail.com
>        Queue: met_help
>      Subject: Empty output from Grid_stat
>        Owner: Nobody
>   Requestors: jinwoong.yoo at gmail.com
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=69314 >
>
>
> Dear MET Help,
>
> I'm testing Grid_stat to compare my WRF simulation outputs to CCSM4
> outputs.
>
> I created some polygon files that represent subregions of my WRF
tropical
> channel domain and they reside at
> /glade/scratch/jyoo/met-5.0/share/met/poly/tropicalsubregion.
>
> Then, I included those polygon information in the
GridStatConfig_ccsm file.
>
> When, I executed grid_stat as below, it seems that only one polygon
(NIN)
> area was incorporated into the calculation.
> Moreover, the output files are still empty.
> What am I missing in this test?
> Please let me know.
> Thank you.
>
> Regards,
>
> Jinwoong Yoo
> UNM
>
>
> [jyoo at yslogin4 analysis]$ ls gridstatout/[jyoo at yslogin4 analysis]$
> /glade/p/ral/jnt/MET/MET_releases/met-5.0/bin/grid_stat
>
>
/glade/scratch/jyoo/DOMAINS/postprd/lgm/postprd/outcopygb/WRFPRS_d01.00_2latlon
> /glade/scratch/jyoo/domain/mean/analysis/LGMCCSM/6hrly/
> b40.lgm21ka.1deg.003M.cam2.h3.1870-01.nc GridStatConfig_ccsm -outdir
> /glade/scratch/jyoo/domain/mean/analysis/gridstatout
> DEBUG 1: Default Config File:
>
> /glade/p/ral/jnt/MET/MET_releases/met-
5.0/share/met/config/GridStatConfig_default
> DEBUG 1: User Config File: GridStatConfig_ccsm
> WARNING:
> WARNING: NcCfFile::open() -> could not extract init time from file
name
> WARNING:
> WARNING:
> WARNING: NcCfFile::get_grid_from_dimensions() -> Found multiple
variables
> for latitude, using "lat".
> WARNING:
> WARNING:
> WARNING: NcCfFile::get_grid_from_dimensions() -> Found multiple
variables
> for longitude, using "lon".
> WARNING:
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=59061588
> DEBUG 1: Forecast File:
>
>
/glade/scratch/jyoo/DOMAINS/postprd/lgm/postprd/outcopygb/WRFPRS_d01.00_2latlon
> DEBUG 1: Observation File:
> /glade/scratch/jyoo/domain/mean/analysis/LGMCCSM/6hrly/
> b40.lgm21ka.1deg.003M.cam2.h3.1870-01.nc
> DEBUG 2:
> DEBUG 2:
>
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Processing PRMSL/R247 veDEBUG 2: Processing PRMSL/R247
versus
> PSL(0,*,*), for interpolation method UW_MEAN(1), over region NIN,
using
> 4260 pairs.
> DEBUG 2: Computing Categorical Statistics.
> DEBUG 2: Computing Continuous Statistics.
> DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for interpolation
method
> UW_MEAN(1), over region NIN, using 10404 pairs.
> DEBUG 2: Computing Categorical Statistics.
> DEBUG 2: Computing Continuous Statistics.
> DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for interpolation
method
> UW_MEAN(1), over region NIN, using 3264 pairs.
> DEBUG 2: Computing Categorical Statistics.
> DEBUG 2: Computing Continuous Statistics.
> DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for interpolation
method
> UW_MEAN(1), over region NIN, using 2880 pairs.
> DEBUG 2: Computing Categorical Statistics.
> DEBUG 2: Computing Continuous Statistics.
> DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for interpolation
method
> UW_MEAN(1), over region NIN, using 3216 pairs.
> DEBUG 2: Computing Categorical Statistics.
> DEBUG 2: Computing Continuous Statistics.
> DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for interpolation
method
> UW_MEAN(1), over region NIN, using 3264 pairs.
> DEBUG 2: Computing Categorical Statistics.
> DEBUG 2: Computing Continuous Statistics.
> DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for interpolation
method
> UW_MEAN(1), over region NIN, using 6480 pairs.
> DEBUG 2: Computing Categorical Statistics.
> DEBUG 2: Computing Continuous Statistics.
> NetCDF: Numeric conversion not representable
>
>

------------------------------------------------
Subject: Empty output from Grid_stat
From: Jinwoong Yoo
Time: Wed Oct 08 13:35:24 2014

Dear John,

I noticed that I kept "NIN" in all of the ascii files of polygon as I
copied the first file for the next region's polygon. So I fixed the
problem.
However, the output files from Grid_stat are still all empty.

Thank you.

Jinwoong Yoo
UNM

On Wed, Oct 8, 2014 at 1:26 PM, John Halley Gotway via RT
<met_help at ucar.edu
> wrote:

> Jinwoong,
>
> You can specify multiple masking regions as a comma-separated list
of file
> names in the config file, like this:
>    poly = [
> "/glade/scratch/jyoo/met-
5.0/share/met/poly/tropicalsubregion/NAT.poly",
>
> "/glade/scratch/jyoo/met-
5.0/share/met/poly/tropicalsubregion/NEP.poly",
>
> "/glade/scratch/jyoo/met-
5.0/share/met/poly/tropicalsubregion/NIN.poly" ];
>
> I suspect this piece of it is working fine.  It just aborted when it
tried
> to write the first set of output to the NetCDF matched pairs output
file.
>
> I'll look into that and let you know.  As I mentioned, this is the
first
> time I've seen that error.
>
> Thanks,
> John
>
> On Wed, Oct 8, 2014 at 12:55 PM, Jinwoong Yoo via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > Wed Oct 08 12:55:15 2014: Request 69314 was acted upon.
> > Transaction: Ticket created by jinwoong.yoo at gmail.com
> >        Queue: met_help
> >      Subject: Empty output from Grid_stat
> >        Owner: Nobody
> >   Requestors: jinwoong.yoo at gmail.com
> >       Status: new
> >  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=69314 >
> >
> >
> > Dear MET Help,
> >
> > I'm testing Grid_stat to compare my WRF simulation outputs to
CCSM4
> > outputs.
> >
> > I created some polygon files that represent subregions of my WRF
tropical
> > channel domain and they reside at
> > /glade/scratch/jyoo/met-5.0/share/met/poly/tropicalsubregion.
> >
> > Then, I included those polygon information in the
GridStatConfig_ccsm
> file.
> >
> > When, I executed grid_stat as below, it seems that only one
polygon (NIN)
> > area was incorporated into the calculation.
> > Moreover, the output files are still empty.
> > What am I missing in this test?
> > Please let me know.
> > Thank you.
> >
> > Regards,
> >
> > Jinwoong Yoo
> > UNM
> >
> >
> > [jyoo at yslogin4 analysis]$ ls gridstatout/[jyoo at yslogin4 analysis]$
> > /glade/p/ral/jnt/MET/MET_releases/met-5.0/bin/grid_stat
> >
> >
>
/glade/scratch/jyoo/DOMAINS/postprd/lgm/postprd/outcopygb/WRFPRS_d01.00_2latlon
> > /glade/scratch/jyoo/domain/mean/analysis/LGMCCSM/6hrly/
> > b40.lgm21ka.1deg.003M.cam2.h3.1870-01.nc GridStatConfig_ccsm
-outdir
> > /glade/scratch/jyoo/domain/mean/analysis/gridstatout
> > DEBUG 1: Default Config File:
> >
> >
> /glade/p/ral/jnt/MET/MET_releases/met-
5.0/share/met/config/GridStatConfig_default
> > DEBUG 1: User Config File: GridStatConfig_ccsm
> > WARNING:
> > WARNING: NcCfFile::open() -> could not extract init time from file
name
> > WARNING:
> > WARNING:
> > WARNING: NcCfFile::get_grid_from_dimensions() -> Found multiple
variables
> > for latitude, using "lat".
> > WARNING:
> > WARNING:
> > WARNING: NcCfFile::get_grid_from_dimensions() -> Found multiple
variables
> > for longitude, using "lon".
> > WARNING:
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=59061588
> > DEBUG 1: Forecast File:
> >
> >
>
/glade/scratch/jyoo/DOMAINS/postprd/lgm/postprd/outcopygb/WRFPRS_d01.00_2latlon
> > DEBUG 1: Observation File:
> > /glade/scratch/jyoo/domain/mean/analysis/LGMCCSM/6hrly/
> > b40.lgm21ka.1deg.003M.cam2.h3.1870-01.nc
> > DEBUG 2:
> > DEBUG 2:
> >
> >
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Processing PRMSL/R247 veDEBUG 2: Processing PRMSL/R247
versus
> > PSL(0,*,*), for interpolation method UW_MEAN(1), over region NIN,
using
> > 4260 pairs.
> > DEBUG 2: Computing Categorical Statistics.
> > DEBUG 2: Computing Continuous Statistics.
> > DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for
interpolation
> method
> > UW_MEAN(1), over region NIN, using 10404 pairs.
> > DEBUG 2: Computing Categorical Statistics.
> > DEBUG 2: Computing Continuous Statistics.
> > DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for
interpolation
> method
> > UW_MEAN(1), over region NIN, using 3264 pairs.
> > DEBUG 2: Computing Categorical Statistics.
> > DEBUG 2: Computing Continuous Statistics.
> > DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for
interpolation
> method
> > UW_MEAN(1), over region NIN, using 2880 pairs.
> > DEBUG 2: Computing Categorical Statistics.
> > DEBUG 2: Computing Continuous Statistics.
> > DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for
interpolation
> method
> > UW_MEAN(1), over region NIN, using 3216 pairs.
> > DEBUG 2: Computing Categorical Statistics.
> > DEBUG 2: Computing Continuous Statistics.
> > DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for
interpolation
> method
> > UW_MEAN(1), over region NIN, using 3264 pairs.
> > DEBUG 2: Computing Categorical Statistics.
> > DEBUG 2: Computing Continuous Statistics.
> > DEBUG 2: Processing PRMSL/R247 versus PSL(0,*,*), for
interpolation
> method
> > UW_MEAN(1), over region NIN, using 6480 pairs.
> > DEBUG 2: Computing Categorical Statistics.
> > DEBUG 2: Computing Continuous Statistics.
> > NetCDF: Numeric conversion not representable
> >
> >
>
>

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


More information about the Met_help mailing list