[Met_help] [rt.rap.ucar.edu #97265] History for Isolate one raob station per vertical level
Minna Win via RT
met_help at ucar.edu
Thu Mar 18 10:01:06 MDT 2021
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello,
Is it possible to force point_stat to only use one sounding to verify a particular vertical level. Essentially, I want the "TOTAL" column in _cnt.txt to be 1, and hence FBAR will equal FCST in _mpr.txt. If that isn't possible, is there a way to print out station (ob) and model (fcst) outputs for each of the soundings that went into a vertical level verification?
Thank you
Justin
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Isolate one raob station per vertical level
From: Minna Win
Time: Thu Oct 29 09:27:50 2020
Hi Justin,
It doesn't look like there is a way to indicate only one sounding per
vertical level. It looks like you'll need modify your input files if
you
want to do this. I don't know if the obs and fcst information which
you
are interested are printed. Perhaps if the debug level is set to a
high
value? I will need to ask one of the engineers who is familiar with
the
code to see if this information can be extracted from the log.
Regards,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
---------------
Pronouns: she/her
On Wed, Oct 28, 2020 at 6:17 PM Tsu, Mr. Justin via RT
<met_help at ucar.edu>
wrote:
>
> Wed Oct 28 18:17:10 2020: Request 97265 was acted upon.
> Transaction: Ticket created by justin.tsu at nrlmry.navy.mil
> Queue: met_help
> Subject: Isolate one raob station per vertical level
> Owner: Nobody
> Requestors: justin.tsu at nrlmry.navy.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=97265 >
>
>
> Hello,
>
> Is it possible to force point_stat to only use one sounding to
verify a
> particular vertical level. Essentially, I want the "TOTAL" column
in
> _cnt.txt to be 1, and hence FBAR will equal FCST in _mpr.txt. If
that
> isn't possible, is there a way to print out station (ob) and model
(fcst)
> outputs for each of the soundings that went into a vertical level
> verification?
>
> Thank you
> Justin
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #97265] Isolate one raob station per vertical level
From: Tsu, Mr. Justin
Time: Thu Oct 29 10:39:52 2020
I would think that at least the latter option (printing out station
information along side the matched pairs) would be possible, since
station IDs are part of the input obs netcdf file. Please let me know
what you find
Thanks!
Justin
-----Original Message-----
From: Minna Win via RT [mailto:met_help at ucar.edu]
Sent: Thursday, October 29, 2020 8:28 AM
To: Tsu, Mr. Justin
Subject: Re: [rt.rap.ucar.edu #97265] Isolate one raob station per
vertical level
Hi Justin,
It doesn't look like there is a way to indicate only one sounding per
vertical level. It looks like you'll need modify your input files if
you
want to do this. I don't know if the obs and fcst information which
you
are interested are printed. Perhaps if the debug level is set to a
high
value? I will need to ask one of the engineers who is familiar with
the
code to see if this information can be extracted from the log.
Regards,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
---------------
Pronouns: she/her
On Wed, Oct 28, 2020 at 6:17 PM Tsu, Mr. Justin via RT
<met_help at ucar.edu>
wrote:
>
> Wed Oct 28 18:17:10 2020: Request 97265 was acted upon.
> Transaction: Ticket created by justin.tsu at nrlmry.navy.mil
> Queue: met_help
> Subject: Isolate one raob station per vertical level
> Owner: Nobody
> Requestors: justin.tsu at nrlmry.navy.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=97265 >
>
>
> Hello,
>
> Is it possible to force point_stat to only use one sounding to
verify a
> particular vertical level. Essentially, I want the "TOTAL" column
in
> _cnt.txt to be 1, and hence FBAR will equal FCST in _mpr.txt. If
that
> isn't possible, is there a way to print out station (ob) and model
(fcst)
> outputs for each of the soundings that went into a vertical level
> verification?
>
> Thank you
> Justin
>
>
------------------------------------------------
Subject: Isolate one raob station per vertical level
From: John Halley Gotway
Time: Thu Oct 29 15:14:25 2020
Minna and Justin,
Here's an option to consider. Point-Stat includes a "mask.sid" option
for
Station ID masking.
Listed below is an explanation taken from:
https://dtcenter.github.io/MET/latest/Users_Guide/data_io.html
If you set:
mask = {
sid = [ "KDEN", "KCOS" ];
...
}
Then those stations will be processed individually and you'll get CNT
lines
in the output where:
VX_MASK = KDEN and VX_MASK = KCOS
So if you know the station names you want ahead of time, that could
work.
Another approach would be running Point-Stat to generate the MPR
output
lines. And then run Stat-Analysis to convert each one to a CNT line:
stat_analysis -lookin point_stat_output.stat -job aggregate_stat
-line_type MPR -out_line_type CNT -out_stat mpr_to_cnt.stat -by
FCST_VAR,FCST_LEV,OBS_SID
Hope that helps.
John
// - The "sid" entry is an array of strings which define groups of
// observation station ID's over which to compute statistics.
Each entry
// in the array is either a filename of a comma-separated list.
// - For a filename, the strings are whitespace-separated. The
first
// string is the mask "name" and the remaining strings are the
station
// ID's to be used.
// - For a comma-separated list, optionally use a colon to
specify a name.
// For "MY_LIST:SID1,SID2", name = MY_LIST and values = SID1
and SID2.
// - For a comma-separated list of length one with no name
specified, the
// mask "name" and value are both set to the single station ID
string.
// For "SID1", name = SID1 and value = SID1.
// - For a comma-separated list of length greater than one with
no name
// specified, the name is set to MASK_SID and the values are
the station
// ID's to be used.
// For "SID1,SID2", name = MASK_SID and values = SID1 and
SID2.
// - The "name" of the station ID mask is written to the VX_MASK
column
// of the MET output files.
On Thu, Oct 29, 2020 at 10:40 AM Tsu, Mr. Justin via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=97265 >
>
> I would think that at least the latter option (printing out station
> information along side the matched pairs) would be possible, since
station
> IDs are part of the input obs netcdf file. Please let me know what
you find
>
>
> Thanks!
> Justin
>
> -----Original Message-----
> From: Minna Win via RT [mailto:met_help at ucar.edu]
> Sent: Thursday, October 29, 2020 8:28 AM
> To: Tsu, Mr. Justin
> Subject: Re: [rt.rap.ucar.edu #97265] Isolate one raob station per
> vertical level
>
> Hi Justin,
>
> It doesn't look like there is a way to indicate only one sounding
per
> vertical level. It looks like you'll need modify your input files
if you
> want to do this. I don't know if the obs and fcst information which
you
> are interested are printed. Perhaps if the debug level is set to a
high
> value? I will need to ask one of the engineers who is familiar with
the
> code to see if this information can be extracted from the log.
>
> Regards,
> Minna
> ---------------
> Minna Win
> National Center for Atmospheric Research
> Developmental Testbed Center
> Phone: 303-497-8423
> Fax: 303-497-8401
> ---------------
> Pronouns: she/her
>
>
> On Wed, Oct 28, 2020 at 6:17 PM Tsu, Mr. Justin via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > Wed Oct 28 18:17:10 2020: Request 97265 was acted upon.
> > Transaction: Ticket created by justin.tsu at nrlmry.navy.mil
> > Queue: met_help
> > Subject: Isolate one raob station per vertical level
> > Owner: Nobody
> > Requestors: justin.tsu at nrlmry.navy.mil
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=97265 >
> >
> >
> > Hello,
> >
> > Is it possible to force point_stat to only use one sounding to
verify a
> > particular vertical level. Essentially, I want the "TOTAL" column
in
> > _cnt.txt to be 1, and hence FBAR will equal FCST in _mpr.txt. If
that
> > isn't possible, is there a way to print out station (ob) and model
(fcst)
> > outputs for each of the soundings that went into a vertical level
> > verification?
> >
> > Thank you
> > Justin
> >
> >
>
>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #97265] Isolate one raob station per vertical level
From: Tsu, Mr. Justin
Time: Thu Oct 29 18:15:32 2020
I ended up using the mask approach and got two lines for each variable
in my cnt file. Looking at the total columns, one line has some
number like 74 and the second line has a TOTAL of 1 ( my mask sid is
one element). Is this second line showing me the ob value for that
one specified station id with the forecast interpolated to that ob
location?
Justin
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Thursday, October 29, 2020 2:14 PM
To: Tsu, Mr. Justin
Subject: Re: [rt.rap.ucar.edu #97265] Isolate one raob station per
vertical level
Minna and Justin,
Here's an option to consider. Point-Stat includes a "mask.sid" option
for
Station ID masking.
Listed below is an explanation taken from:
https://dtcenter.github.io/MET/latest/Users_Guide/data_io.html
If you set:
mask = {
sid = [ "KDEN", "KCOS" ];
...
}
Then those stations will be processed individually and you'll get CNT
lines
in the output where:
VX_MASK = KDEN and VX_MASK = KCOS
So if you know the station names you want ahead of time, that could
work.
Another approach would be running Point-Stat to generate the MPR
output
lines. And then run Stat-Analysis to convert each one to a CNT line:
stat_analysis -lookin point_stat_output.stat -job aggregate_stat
-line_type MPR -out_line_type CNT -out_stat mpr_to_cnt.stat -by
FCST_VAR,FCST_LEV,OBS_SID
Hope that helps.
John
// - The "sid" entry is an array of strings which define groups of
// observation station ID's over which to compute statistics.
Each entry
// in the array is either a filename of a comma-separated list.
// - For a filename, the strings are whitespace-separated. The
first
// string is the mask "name" and the remaining strings are the
station
// ID's to be used.
// - For a comma-separated list, optionally use a colon to
specify a name.
// For "MY_LIST:SID1,SID2", name = MY_LIST and values = SID1
and SID2.
// - For a comma-separated list of length one with no name
specified, the
// mask "name" and value are both set to the single station ID
string.
// For "SID1", name = SID1 and value = SID1.
// - For a comma-separated list of length greater than one with
no name
// specified, the name is set to MASK_SID and the values are
the station
// ID's to be used.
// For "SID1,SID2", name = MASK_SID and values = SID1 and
SID2.
// - The "name" of the station ID mask is written to the VX_MASK
column
// of the MET output files.
On Thu, Oct 29, 2020 at 10:40 AM Tsu, Mr. Justin via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=97265 >
>
> I would think that at least the latter option (printing out station
> information along side the matched pairs) would be possible, since
station
> IDs are part of the input obs netcdf file. Please let me know what
you find
>
>
> Thanks!
> Justin
>
> -----Original Message-----
> From: Minna Win via RT [mailto:met_help at ucar.edu]
> Sent: Thursday, October 29, 2020 8:28 AM
> To: Tsu, Mr. Justin
> Subject: Re: [rt.rap.ucar.edu #97265] Isolate one raob station per
> vertical level
>
> Hi Justin,
>
> It doesn't look like there is a way to indicate only one sounding
per
> vertical level. It looks like you'll need modify your input files
if you
> want to do this. I don't know if the obs and fcst information which
you
> are interested are printed. Perhaps if the debug level is set to a
high
> value? I will need to ask one of the engineers who is familiar with
the
> code to see if this information can be extracted from the log.
>
> Regards,
> Minna
> ---------------
> Minna Win
> National Center for Atmospheric Research
> Developmental Testbed Center
> Phone: 303-497-8423
> Fax: 303-497-8401
> ---------------
> Pronouns: she/her
>
>
> On Wed, Oct 28, 2020 at 6:17 PM Tsu, Mr. Justin via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > Wed Oct 28 18:17:10 2020: Request 97265 was acted upon.
> > Transaction: Ticket created by justin.tsu at nrlmry.navy.mil
> > Queue: met_help
> > Subject: Isolate one raob station per vertical level
> > Owner: Nobody
> > Requestors: justin.tsu at nrlmry.navy.mil
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=97265 >
> >
> >
> > Hello,
> >
> > Is it possible to force point_stat to only use one sounding to
verify a
> > particular vertical level. Essentially, I want the "TOTAL" column
in
> > _cnt.txt to be 1, and hence FBAR will equal FCST in _mpr.txt. If
that
> > isn't possible, is there a way to print out station (ob) and model
(fcst)
> > outputs for each of the soundings that went into a vertical level
> > verification?
> >
> > Thank you
> > Justin
> >
> >
>
>
>
>
------------------------------------------------
Subject: Isolate one raob station per vertical level
From: John Halley Gotway
Time: Fri Oct 30 12:33:36 2020
Justin,
Without seeing your exact config file and output files, I can't answer
with
certainty. You describe seeing TOTAL counts in the output with values
of 74
and 1. My guess is that the lines where TOTAL = 1 corresponds to the
station id mask(s) you've defined. And the line where TOTAL = 74
corresponds to some other spatial aggregation, like VX_MASK = FULL or
CONUS
or something.
For each CNT line, look closely at the TOTAL and VX_MASK columns.
John
On Thu, Oct 29, 2020 at 6:16 PM Tsu, Mr. Justin via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=97265 >
>
> I ended up using the mask approach and got two lines for each
variable in
> my cnt file. Looking at the total columns, one line has some number
like
> 74 and the second line has a TOTAL of 1 ( my mask sid is one
element). Is
> this second line showing me the ob value for that one specified
station id
> with the forecast interpolated to that ob location?
>
> Justin
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Thursday, October 29, 2020 2:14 PM
> To: Tsu, Mr. Justin
> Subject: Re: [rt.rap.ucar.edu #97265] Isolate one raob station per
> vertical level
>
> Minna and Justin,
>
> Here's an option to consider. Point-Stat includes a "mask.sid"
option for
> Station ID masking.
> Listed below is an explanation taken from:
> https://dtcenter.github.io/MET/latest/Users_Guide/data_io.html
>
> If you set:
> mask = {
> sid = [ "KDEN", "KCOS" ];
> ...
> }
>
> Then those stations will be processed individually and you'll get
CNT lines
> in the output where:
> VX_MASK = KDEN and VX_MASK = KCOS
>
> So if you know the station names you want ahead of time, that could
work.
>
> Another approach would be running Point-Stat to generate the MPR
output
> lines. And then run Stat-Analysis to convert each one to a CNT line:
> stat_analysis -lookin point_stat_output.stat -job aggregate_stat
> -line_type MPR -out_line_type CNT -out_stat mpr_to_cnt.stat -by
> FCST_VAR,FCST_LEV,OBS_SID
>
> Hope that helps.
>
> John
>
> // - The "sid" entry is an array of strings which define groups
of
> // observation station ID's over which to compute statistics.
Each
> entry
> // in the array is either a filename of a comma-separated list.
> // - For a filename, the strings are whitespace-separated. The
first
> // string is the mask "name" and the remaining strings are
the
> station
> // ID's to be used.
> // - For a comma-separated list, optionally use a colon to
specify a
> name.
> // For "MY_LIST:SID1,SID2", name = MY_LIST and values = SID1
and
> SID2.
> // - For a comma-separated list of length one with no name
> specified, the
> // mask "name" and value are both set to the single station
ID
> string.
> // For "SID1", name = SID1 and value = SID1.
> // - For a comma-separated list of length greater than one
with no
> name
> // specified, the name is set to MASK_SID and the values are
the
> station
> // ID's to be used.
> // For "SID1,SID2", name = MASK_SID and values = SID1 and
SID2.
> // - The "name" of the station ID mask is written to the
VX_MASK
> column
> // of the MET output files.
>
>
> On Thu, Oct 29, 2020 at 10:40 AM Tsu, Mr. Justin via RT
<met_help at ucar.edu
> >
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=97265 >
> >
> > I would think that at least the latter option (printing out
station
> > information along side the matched pairs) would be possible, since
> station
> > IDs are part of the input obs netcdf file. Please let me know
what you
> find
> >
> >
> > Thanks!
> > Justin
> >
> > -----Original Message-----
> > From: Minna Win via RT [mailto:met_help at ucar.edu]
> > Sent: Thursday, October 29, 2020 8:28 AM
> > To: Tsu, Mr. Justin
> > Subject: Re: [rt.rap.ucar.edu #97265] Isolate one raob station per
> > vertical level
> >
> > Hi Justin,
> >
> > It doesn't look like there is a way to indicate only one sounding
per
> > vertical level. It looks like you'll need modify your input files
if you
> > want to do this. I don't know if the obs and fcst information
which you
> > are interested are printed. Perhaps if the debug level is set to
a high
> > value? I will need to ask one of the engineers who is familiar
with the
> > code to see if this information can be extracted from the log.
> >
> > Regards,
> > Minna
> > ---------------
> > Minna Win
> > National Center for Atmospheric Research
> > Developmental Testbed Center
> > Phone: 303-497-8423
> > Fax: 303-497-8401
> > ---------------
> > Pronouns: she/her
> >
> >
> > On Wed, Oct 28, 2020 at 6:17 PM Tsu, Mr. Justin via RT <
> met_help at ucar.edu>
> > wrote:
> >
> > >
> > > Wed Oct 28 18:17:10 2020: Request 97265 was acted upon.
> > > Transaction: Ticket created by justin.tsu at nrlmry.navy.mil
> > > Queue: met_help
> > > Subject: Isolate one raob station per vertical level
> > > Owner: Nobody
> > > Requestors: justin.tsu at nrlmry.navy.mil
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=97265
> >
> > >
> > >
> > > Hello,
> > >
> > > Is it possible to force point_stat to only use one sounding to
verify a
> > > particular vertical level. Essentially, I want the "TOTAL"
column in
> > > _cnt.txt to be 1, and hence FBAR will equal FCST in _mpr.txt.
If that
> > > isn't possible, is there a way to print out station (ob) and
model
> (fcst)
> > > outputs for each of the soundings that went into a vertical
level
> > > verification?
> > >
> > > Thank you
> > > Justin
> > >
> > >
> >
> >
> >
> >
>
>
>
>
------------------------------------------------
More information about the Met_help
mailing list