[Met_help] [rt.rap.ucar.edu #90873] History for Defining decimal height levels
John Halley Gotway via RT
met_help at ucar.edu
Tue Jul 2 12:24:02 MDT 2019
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello,
For the GFS we verify the moisture and temperature in the top soil layer
(0-10 cm). For the grib1 files in MET the variable information would be
defined, name="TSOIL"; level="Z0-10";...for example,
*regrid_data_plane
/gpfs/hps3/emc/global/noscrub/Fanglin.Yang/stat/gfs/pgbf24.gfs.2019060100
G002 test.nc <http://test.nc> -field 'name="TSOIL"; level="Z0-10";'*DEBUG
1: Reading data file:
/gpfs/hps3/emc/global/noscrub/Fanglin.Yang/stat/gfs/pgbf24.gfs.2019060100
DEBUG 2: Input grid: Projection: Lat/Lon Nx: 360 Ny: 181 lat_ll: -90.000
lon_ll: -0.000 delta_lat: 1.000 delta_lon: 1.000
DEBUG 2: Output grid: Projection: Lat/Lon Nx: 144 Ny: 73 lat_ll: -90.000
lon_ll: -0.000 delta_lat: 2.500 delta_lon: 2.500
DEBUG 2: Interpolation options: method = NEAREST, width = 1, shape =
SQUARE, vld_thresh = 0.5
DEBUG 2: Range of input data (name="TSOIL"; level="Z0-10";) is 205.948 to
317.748.
DEBUG 2: Range of regridded data (name="TSOIL"; level="Z0-10";) is 206.148
to 316.748.
DEBUG 1: Writing output file: test.nc
For the next version of the GFS the output is being saved as GRIB2 files.
In the GRIB2 files the top soil layer is now defined in meters, so 0-0.1m,
but I can't seem to get this variable to be recongized by MET. Can the
height layer be defined as decimals?
*regrid_data_plane
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
G002 test.nc <http://test.nc> -field 'name="TSOIL"; level="Z0-0.1";'*
DEBUG 1: Reading data file:
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
DEBUG 2: Input grid: Projection: Lat/Lon Nx: 360 Ny: 181 lat_ll: -90.000
lon_ll: -0.000 delta_lat: 1.000 delta_lon: 1.000
DEBUG 2: Output grid: Projection: Lat/Lon Nx: 144 Ny: 73 lat_ll: -90.000
lon_ll: -0.000 delta_lat: 2.500 delta_lon: 2.500
DEBUG 2: Interpolation options: method = NEAREST, width = 1, shape =
SQUARE, vld_thresh = 0.5
WARNING:
WARNING: MetGrib2DataFile::data_plane() -> No matching record found for
'TSOIL/Z0.1-0'
WARNING:
ERROR :
ERROR : process_data_file() -> trouble getting field "name="TSOIL";
level="Z0-0.1";" from file
"/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2"
Though the above examples are using the MET calls directly, I'm actually
using METplus to do the verification. I was just messing around with
regrid_data_plane to try to see/figure out the variable level information
before setting in the METplus config file and running.
Mallory
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Defining decimal height levels
From: John Halley Gotway
Time: Tue Jul 02 11:00:59 2019
Hi Mallory,
The short answer is... just replace the 'Z' with an 'L' in your level
setting and it'll work. Listed below are examples of the
plot_data_plane
commands I ran to test this.
I see that the MET GRIB2 library code isn't interpreting this level
type
value as being a "vertical level type" (in which case the 'Z' would
work).
The level types are defined in this table:
https://www.nco.ncep.noaa.gov/pmb/docs/grib2/grib2_doc/grib2_table4-
5.shtml
In met-8.1, 'Z' matches level types 1, 101, 102, and 103.
I propose that 'Z' should also match...
*106: Depth Below Land Surface*
*117: Mixed Layer Depth*
*160: Depth Below Sea Level*
*161: Depth Below Water Surface*
Do you agree?
Thanks,
John
Yes, met-8.1 should read those levels just fine. To confirm I logged
on to
WCOSS and checked:
[John.H.Gotway at t10a2 ~]$ module use
/global/noscrub/Julie.Prestopnik/modulefiles
[John.H.Gotway at t10a2 ~]$ module load met/8.1
[John.H.Gotway at t10a2 ~]$ module load grib_util
[John.H.Gotway at t10a2 ~]$ $WGRIB2
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
| grep TSOIL
491:24142266:d=2019062412:TSOIL:0-0.1 m below ground:42 hour fcst:
493:24196608:d=2019062412:TSOIL:0.1-0.4 m below ground:42 hour fcst:
495:24250540:d=2019062412:TSOIL:0.4-1 m below ground:42 hour fcst:
497:24302377:d=2019062412:TSOIL:1-2 m below ground:42 hour fcst:
The TSOIL data we want is in records 491, 493, 495, and 497. You can
plot
them like this:
plot_data_plane
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
~/TSOIL_0-0.1.ps 'name="TSOIL"; level="L0-0.1";' -v 3
plot_data_plane
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
~/TSOIL_0.1-0.4.ps 'name="TSOIL"; level="L0.1-0.4";' -v 3
plot_data_plane
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
~/TSOIL_0.4-1.ps 'name="TSOIL"; level="L0.4-1";' -v 3
plot_data_plane
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
~/TSOIL_1-2.ps 'name="TSOIL"; level="L1-2";' -v 3
And the log messages confirm that we matched the right records:
DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
'TSOIL/L0.1-0' in GRIB2 record 491
DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
'TSOIL/L0.4-0.1' in GRIB2 record 493
DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
'TSOIL/L1-0.4' in GRIB2 record 495
DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
'TSOIL/L2-1' in GRIB2 record 497
On Tue, Jul 2, 2019 at 7:13 AM Mallory Row - NOAA Affiliate via RT <
met_help at ucar.edu> wrote:
>
> Tue Jul 02 07:13:20 2019: Request 90873 was acted upon.
> Transaction: Ticket created by mallory.row at noaa.gov
> Queue: met_help
> Subject: Defining decimal height levels
> Owner: Nobody
> Requestors: mallory.row at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=90873 >
>
>
> Hello,
>
> For the GFS we verify the moisture and temperature in the top soil
layer
> (0-10 cm). For the grib1 files in MET the variable information would
be
> defined, name="TSOIL"; level="Z0-10";...for example,
>
>
> *regrid_data_plane
>
/gpfs/hps3/emc/global/noscrub/Fanglin.Yang/stat/gfs/pgbf24.gfs.2019060100
> G002 test.nc <http://test.nc> -field 'name="TSOIL"; level="Z0-
10";'*DEBUG
> 1: Reading data file:
>
/gpfs/hps3/emc/global/noscrub/Fanglin.Yang/stat/gfs/pgbf24.gfs.2019060100
> DEBUG 2: Input grid: Projection: Lat/Lon Nx: 360 Ny: 181 lat_ll:
-90.000
> lon_ll: -0.000 delta_lat: 1.000 delta_lon: 1.000
> DEBUG 2: Output grid: Projection: Lat/Lon Nx: 144 Ny: 73 lat_ll:
-90.000
> lon_ll: -0.000 delta_lat: 2.500 delta_lon: 2.500
> DEBUG 2: Interpolation options: method = NEAREST, width = 1, shape =
> SQUARE, vld_thresh = 0.5
> DEBUG 2: Range of input data (name="TSOIL"; level="Z0-10";) is
205.948 to
> 317.748.
> DEBUG 2: Range of regridded data (name="TSOIL"; level="Z0-10";) is
206.148
> to 316.748.
> DEBUG 1: Writing output file: test.nc
>
> For the next version of the GFS the output is being saved as GRIB2
files.
> In the GRIB2 files the top soil layer is now defined in meters, so
0-0.1m,
> but I can't seem to get this variable to be recongized by MET. Can
the
> height layer be defined as decimals?
>
> *regrid_data_plane
>
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> G002 test.nc <http://test.nc> -field 'name="TSOIL"; level="Z0-
0.1";'*
> DEBUG 1: Reading data file:
>
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> DEBUG 2: Input grid: Projection: Lat/Lon Nx: 360 Ny: 181 lat_ll:
-90.000
> lon_ll: -0.000 delta_lat: 1.000 delta_lon: 1.000
> DEBUG 2: Output grid: Projection: Lat/Lon Nx: 144 Ny: 73 lat_ll:
-90.000
> lon_ll: -0.000 delta_lat: 2.500 delta_lon: 2.500
> DEBUG 2: Interpolation options: method = NEAREST, width = 1, shape =
> SQUARE, vld_thresh = 0.5
> WARNING:
> WARNING: MetGrib2DataFile::data_plane() -> No matching record found
for
> 'TSOIL/Z0.1-0'
> WARNING:
> ERROR :
> ERROR : process_data_file() -> trouble getting field "name="TSOIL";
> level="Z0-0.1";" from file
>
>
"/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2"
>
> Though the above examples are using the MET calls directly, I'm
actually
> using METplus to do the verification. I was just messing around with
> regrid_data_plane to try to see/figure out the variable level
information
> before setting in the METplus config file and running.
>
> Mallory
>
>
------------------------------------------------
Subject: Defining decimal height levels
From: Mallory Row - NOAA Affiliate
Time: Tue Jul 02 11:04:09 2019
Thanks! Yeah, I think that sounds good.
Mallory
On Tue, Jul 2, 2019 at 1:01 PM John Halley Gotway via RT
<met_help at ucar.edu>
wrote:
> Hi Mallory,
>
>
> The short answer is... just replace the 'Z' with an 'L' in your
level
> setting and it'll work. Listed below are examples of the
plot_data_plane
> commands I ran to test this.
>
>
> I see that the MET GRIB2 library code isn't interpreting this level
type
> value as being a "vertical level type" (in which case the 'Z' would
work).
>
>
> The level types are defined in this table:
>
> https://www.nco.ncep.noaa.gov/pmb/docs/grib2/grib2_doc/grib2_table4-
5.shtml
>
>
> In met-8.1, 'Z' matches level types 1, 101, 102, and 103.
>
> I propose that 'Z' should also match...
>
> *106: Depth Below Land Surface*
>
> *117: Mixed Layer Depth*
>
> *160: Depth Below Sea Level*
>
> *161: Depth Below Water Surface*
>
>
> Do you agree?
>
>
> Thanks,
> John
>
>
> Yes, met-8.1 should read those levels just fine. To confirm I
logged on to
> WCOSS and checked:
>
>
> [John.H.Gotway at t10a2 ~]$ module use
> /global/noscrub/Julie.Prestopnik/modulefiles
>
> [John.H.Gotway at t10a2 ~]$ module load met/8.1
>
> [John.H.Gotway at t10a2 ~]$ module load grib_util
>
> [John.H.Gotway at t10a2 ~]$ $WGRIB2
>
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> | grep TSOIL
>
> 491:24142266:d=2019062412:TSOIL:0-0.1 m below ground:42 hour fcst:
>
> 493:24196608:d=2019062412:TSOIL:0.1-0.4 m below ground:42 hour fcst:
>
> 495:24250540:d=2019062412:TSOIL:0.4-1 m below ground:42 hour fcst:
>
> 497:24302377:d=2019062412:TSOIL:1-2 m below ground:42 hour fcst:
>
>
> The TSOIL data we want is in records 491, 493, 495, and 497. You
can plot
> them like this:
>
>
> plot_data_plane
>
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> ~/TSOIL_0-0.1.ps 'name="TSOIL"; level="L0-0.1";' -v 3
>
>
> plot_data_plane
>
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> ~/TSOIL_0.1-0.4.ps 'name="TSOIL"; level="L0.1-0.4";' -v 3
>
>
> plot_data_plane
>
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> ~/TSOIL_0.4-1.ps 'name="TSOIL"; level="L0.4-1";' -v 3
>
>
> plot_data_plane
>
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> ~/TSOIL_1-2.ps 'name="TSOIL"; level="L1-2";' -v 3
>
>
> And the log messages confirm that we matched the right records:
>
> DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> 'TSOIL/L0.1-0' in GRIB2 record 491
>
> DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> 'TSOIL/L0.4-0.1' in GRIB2 record 493
>
> DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> 'TSOIL/L1-0.4' in GRIB2 record 495
>
> DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> 'TSOIL/L2-1' in GRIB2 record 497
>
> On Tue, Jul 2, 2019 at 7:13 AM Mallory Row - NOAA Affiliate via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Tue Jul 02 07:13:20 2019: Request 90873 was acted upon.
> > Transaction: Ticket created by mallory.row at noaa.gov
> > Queue: met_help
> > Subject: Defining decimal height levels
> > Owner: Nobody
> > Requestors: mallory.row at noaa.gov
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=90873 >
> >
> >
> > Hello,
> >
> > For the GFS we verify the moisture and temperature in the top soil
layer
> > (0-10 cm). For the grib1 files in MET the variable information
would be
> > defined, name="TSOIL"; level="Z0-10";...for example,
> >
> >
> > *regrid_data_plane
> >
/gpfs/hps3/emc/global/noscrub/Fanglin.Yang/stat/gfs/pgbf24.gfs.2019060100
> > G002 test.nc <http://test.nc> -field 'name="TSOIL";
> level="Z0-10";'*DEBUG
> > 1: Reading data file:
> >
/gpfs/hps3/emc/global/noscrub/Fanglin.Yang/stat/gfs/pgbf24.gfs.2019060100
> > DEBUG 2: Input grid: Projection: Lat/Lon Nx: 360 Ny: 181 lat_ll:
-90.000
> > lon_ll: -0.000 delta_lat: 1.000 delta_lon: 1.000
> > DEBUG 2: Output grid: Projection: Lat/Lon Nx: 144 Ny: 73 lat_ll:
-90.000
> > lon_ll: -0.000 delta_lat: 2.500 delta_lon: 2.500
> > DEBUG 2: Interpolation options: method = NEAREST, width = 1, shape
=
> > SQUARE, vld_thresh = 0.5
> > DEBUG 2: Range of input data (name="TSOIL"; level="Z0-10";) is
205.948 to
> > 317.748.
> > DEBUG 2: Range of regridded data (name="TSOIL"; level="Z0-10";) is
> 206.148
> > to 316.748.
> > DEBUG 1: Writing output file: test.nc
> >
> > For the next version of the GFS the output is being saved as GRIB2
files.
> > In the GRIB2 files the top soil layer is now defined in meters, so
> 0-0.1m,
> > but I can't seem to get this variable to be recongized by MET. Can
the
> > height layer be defined as decimals?
> >
> > *regrid_data_plane
> >
> >
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> > G002 test.nc <http://test.nc> -field 'name="TSOIL"; level="Z0-
0.1";'*
> > DEBUG 1: Reading data file:
> >
> >
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> > DEBUG 2: Input grid: Projection: Lat/Lon Nx: 360 Ny: 181 lat_ll:
-90.000
> > lon_ll: -0.000 delta_lat: 1.000 delta_lon: 1.000
> > DEBUG 2: Output grid: Projection: Lat/Lon Nx: 144 Ny: 73 lat_ll:
-90.000
> > lon_ll: -0.000 delta_lat: 2.500 delta_lon: 2.500
> > DEBUG 2: Interpolation options: method = NEAREST, width = 1, shape
=
> > SQUARE, vld_thresh = 0.5
> > WARNING:
> > WARNING: MetGrib2DataFile::data_plane() -> No matching record
found for
> > 'TSOIL/Z0.1-0'
> > WARNING:
> > ERROR :
> > ERROR : process_data_file() -> trouble getting field
"name="TSOIL";
> > level="Z0-0.1";" from file
> >
> >
>
"/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2"
> >
> > Though the above examples are using the MET calls directly, I'm
actually
> > using METplus to do the verification. I was just messing around
with
> > regrid_data_plane to try to see/figure out the variable level
information
> > before setting in the METplus config file and running.
> >
> > Mallory
> >
> >
>
>
------------------------------------------------
Subject: Defining decimal height levels
From: John Halley Gotway
Time: Tue Jul 02 12:04:29 2019
Mallory,
I created a GitHub issue for this (mostly just to track the change)
and
will include it in met-8.1.1:
https://github.com/NCAR/MET/issues/1159
Thanks,
John
On Tue, Jul 2, 2019 at 11:04 AM Mallory Row - NOAA Affiliate via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=90873 >
>
> Thanks! Yeah, I think that sounds good.
>
> Mallory
>
> On Tue, Jul 2, 2019 at 1:01 PM John Halley Gotway via RT <
> met_help at ucar.edu>
> wrote:
>
> > Hi Mallory,
> >
> >
> > The short answer is... just replace the 'Z' with an 'L' in your
level
> > setting and it'll work. Listed below are examples of the
plot_data_plane
> > commands I ran to test this.
> >
> >
> > I see that the MET GRIB2 library code isn't interpreting this
level type
> > value as being a "vertical level type" (in which case the 'Z'
would
> work).
> >
> >
> > The level types are defined in this table:
> >
> >
> https://www.nco.ncep.noaa.gov/pmb/docs/grib2/grib2_doc/grib2_table4-
5.shtml
> >
> >
> > In met-8.1, 'Z' matches level types 1, 101, 102, and 103.
> >
> > I propose that 'Z' should also match...
> >
> > *106: Depth Below Land Surface*
> >
> > *117: Mixed Layer Depth*
> >
> > *160: Depth Below Sea Level*
> >
> > *161: Depth Below Water Surface*
> >
> >
> > Do you agree?
> >
> >
> > Thanks,
> > John
> >
> >
> > Yes, met-8.1 should read those levels just fine. To confirm I
logged on
> to
> > WCOSS and checked:
> >
> >
> > [John.H.Gotway at t10a2 ~]$ module use
> > /global/noscrub/Julie.Prestopnik/modulefiles
> >
> > [John.H.Gotway at t10a2 ~]$ module load met/8.1
> >
> > [John.H.Gotway at t10a2 ~]$ module load grib_util
> >
> > [John.H.Gotway at t10a2 ~]$ $WGRIB2
> >
> >
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> > | grep TSOIL
> >
> > 491:24142266:d=2019062412:TSOIL:0-0.1 m below ground:42 hour fcst:
> >
> > 493:24196608:d=2019062412:TSOIL:0.1-0.4 m below ground:42 hour
fcst:
> >
> > 495:24250540:d=2019062412:TSOIL:0.4-1 m below ground:42 hour fcst:
> >
> > 497:24302377:d=2019062412:TSOIL:1-2 m below ground:42 hour fcst:
> >
> >
> > The TSOIL data we want is in records 491, 493, 495, and 497. You
can
> plot
> > them like this:
> >
> >
> > plot_data_plane
> >
> >
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> > ~/TSOIL_0-0.1.ps 'name="TSOIL"; level="L0-0.1";' -v 3
> >
> >
> > plot_data_plane
> >
> >
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> > ~/TSOIL_0.1-0.4.ps 'name="TSOIL"; level="L0.1-0.4";' -v 3
> >
> >
> > plot_data_plane
> >
> >
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> > ~/TSOIL_0.4-1.ps 'name="TSOIL"; level="L0.4-1";' -v 3
> >
> >
> > plot_data_plane
> >
> >
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> > ~/TSOIL_1-2.ps 'name="TSOIL"; level="L1-2";' -v 3
> >
> >
> > And the log messages confirm that we matched the right records:
> >
> > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> > 'TSOIL/L0.1-0' in GRIB2 record 491
> >
> > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> > 'TSOIL/L0.4-0.1' in GRIB2 record 493
> >
> > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> > 'TSOIL/L1-0.4' in GRIB2 record 495
> >
> > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> > 'TSOIL/L2-1' in GRIB2 record 497
> >
> > On Tue, Jul 2, 2019 at 7:13 AM Mallory Row - NOAA Affiliate via RT
<
> > met_help at ucar.edu> wrote:
> >
> > >
> > > Tue Jul 02 07:13:20 2019: Request 90873 was acted upon.
> > > Transaction: Ticket created by mallory.row at noaa.gov
> > > Queue: met_help
> > > Subject: Defining decimal height levels
> > > Owner: Nobody
> > > Requestors: mallory.row at noaa.gov
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=90873
> >
> > >
> > >
> > > Hello,
> > >
> > > For the GFS we verify the moisture and temperature in the top
soil
> layer
> > > (0-10 cm). For the grib1 files in MET the variable information
would be
> > > defined, name="TSOIL"; level="Z0-10";...for example,
> > >
> > >
> > > *regrid_data_plane
> > >
>
/gpfs/hps3/emc/global/noscrub/Fanglin.Yang/stat/gfs/pgbf24.gfs.2019060100
> > > G002 test.nc <http://test.nc> -field 'name="TSOIL";
> > level="Z0-10";'*DEBUG
> > > 1: Reading data file:
> > >
>
/gpfs/hps3/emc/global/noscrub/Fanglin.Yang/stat/gfs/pgbf24.gfs.2019060100
> > > DEBUG 2: Input grid: Projection: Lat/Lon Nx: 360 Ny: 181 lat_ll:
> -90.000
> > > lon_ll: -0.000 delta_lat: 1.000 delta_lon: 1.000
> > > DEBUG 2: Output grid: Projection: Lat/Lon Nx: 144 Ny: 73 lat_ll:
> -90.000
> > > lon_ll: -0.000 delta_lat: 2.500 delta_lon: 2.500
> > > DEBUG 2: Interpolation options: method = NEAREST, width = 1,
shape =
> > > SQUARE, vld_thresh = 0.5
> > > DEBUG 2: Range of input data (name="TSOIL"; level="Z0-10";) is
205.948
> to
> > > 317.748.
> > > DEBUG 2: Range of regridded data (name="TSOIL"; level="Z0-10";)
is
> > 206.148
> > > to 316.748.
> > > DEBUG 1: Writing output file: test.nc
> > >
> > > For the next version of the GFS the output is being saved as
GRIB2
> files.
> > > In the GRIB2 files the top soil layer is now defined in meters,
so
> > 0-0.1m,
> > > but I can't seem to get this variable to be recongized by MET.
Can the
> > > height layer be defined as decimals?
> > >
> > > *regrid_data_plane
> > >
> > >
> >
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> > > G002 test.nc <http://test.nc> -field 'name="TSOIL"; level="Z0-
0.1";'*
> > > DEBUG 1: Reading data file:
> > >
> > >
> >
>
/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2
> > > DEBUG 2: Input grid: Projection: Lat/Lon Nx: 360 Ny: 181 lat_ll:
> -90.000
> > > lon_ll: -0.000 delta_lat: 1.000 delta_lon: 1.000
> > > DEBUG 2: Output grid: Projection: Lat/Lon Nx: 144 Ny: 73 lat_ll:
> -90.000
> > > lon_ll: -0.000 delta_lat: 2.500 delta_lon: 2.500
> > > DEBUG 2: Interpolation options: method = NEAREST, width = 1,
shape =
> > > SQUARE, vld_thresh = 0.5
> > > WARNING:
> > > WARNING: MetGrib2DataFile::data_plane() -> No matching record
found for
> > > 'TSOIL/Z0.1-0'
> > > WARNING:
> > > ERROR :
> > > ERROR : process_data_file() -> trouble getting field
"name="TSOIL";
> > > level="Z0-0.1";" from file
> > >
> > >
> >
>
"/gpfs/dell2/emc/modeling/noscrub/emc.glopara/archive/v16rt0/pgbf42.gfs.2019062412.grib2"
> > >
> > > Though the above examples are using the MET calls directly, I'm
> actually
> > > using METplus to do the verification. I was just messing around
with
> > > regrid_data_plane to try to see/figure out the variable level
> information
> > > before setting in the METplus config file and running.
> > >
> > > Mallory
> > >
> > >
> >
> >
>
>
------------------------------------------------
More information about the Met_help
mailing list