[Met_help] [rt.rap.ucar.edu #93905] History for Topo_mask for multiple nests in one file
John Halley Gotway via RT
met_help at ucar.edu
Wed Apr 1 11:01:59 MDT 2020
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hey John,
How do I specify the topo mask using python embedding if I have multiple
different nests represented in my field/name entry. Say I have
field = [
{name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_4a0331x0331_${INIT_TIME}_${FCST_HR
}_fcstfld ${WIND_FLAG}";},
{name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_3a0211x0211_${INIT_TIME}_${FCST_HR
}_fcstfld ${WIND_FLAG}";},
{name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_2a0271x0271_${INIT_TIME}_${FCST_HR
}_fcstfld ${WIND_FLAG}";},
{name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_1a0127x0127_${INIT_TIME}_${FCST_HR
}_fcstfld ${WIND_FLAG}";},
...
]
Would my topo_mask then look something lke this?
topo_mask = {
flag = TRUE;
FILE_NAME = ["PYTHON_NUMPY"];
field = [
{ name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
${DATA_DIR}/terrht_sfc_000000_000000_4a0331x0331_${INIT_TIME}_${FCST_HR}_fcs
tfld";},
{ name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
${DATA_DIR}/terrht_sfc_000000_000000_3a0211x0211_${INIT_TIME}_${FCST_HR}_fcs
tfld";},
{ name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
${DATA_DIR}/terrht_sfc_000000_000000_2a0271x0271_${INIT_TIME}_${FCST_HR}_fcs
tfld";},
{ name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
${DATA_DIR}/terrht_sfc_000000_000000_1a0127x0127_${INIT_TIME}_${FCST_HR}_fcs
tfld";},
.
]
regrid = { method = BILIN; width=2;}
use_obs_thresh = ge-500&&le500;
interp_fcst_thresh = ge-50&&le50;
}
Justin Tsu
Marine Meteorology Division
Data Assimilation/Mesoscale Modeling
Building 704 Room 212
Naval Research Laboratory, Code 7531
7 Grace Hopper Avenue
Monterey, CA 93943-5502
Ph. (831) 656-4111
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Topo_mask for multiple nests in one file
From: John Halley Gotway
Time: Fri Jan 24 09:21:52 2020
Justin,
I don't think this will work. When you run Point-Stat, you have to
pick a
single verification domain. By that I mean you pick a single grid on
which
the verification is performed. The "regrid" config file section
allows you
to be flexible on how you define that domain. But ultimately, you
pick a
single one. Things like the land/sea mask, topography, and masking
regions
must all be defined on that grid.
So if you want to evaluate different verification domains, you'll need
to
run Point-Stat separately for each one.
Hope that helps clarify.
Thanks,
John
On Thu, Jan 23, 2020 at 5:05 PM Tsu, Mr. Justin via RT
<met_help at ucar.edu>
wrote:
>
> Thu Jan 23 17:05:36 2020: Request 93905 was acted upon.
> Transaction: Ticket created by justin.tsu at nrlmry.navy.mil
> Queue: met_help
> Subject: Topo_mask for multiple nests in one file
> Owner: Nobody
> Requestors: justin.tsu at nrlmry.navy.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=93905 >
>
>
> Hey John,
>
>
>
> How do I specify the topo mask using python embedding if I have
multiple
> different nests represented in my field/name entry. Say I have
>
>
>
> field = [
>
> {name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_4a0331x0331_${INIT_TIME}_${FCST_HR
> }_fcstfld ${WIND_FLAG}";},
>
> {name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_3a0211x0211_${INIT_TIME}_${FCST_HR
> }_fcstfld ${WIND_FLAG}";},
>
> {name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_2a0271x0271_${INIT_TIME}_${FCST_HR
> }_fcstfld ${WIND_FLAG}";},
>
> {name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_1a0127x0127_${INIT_TIME}_${FCST_HR
> }_fcstfld ${WIND_FLAG}";},
>
> ...
>
> ]
>
>
>
> Would my topo_mask then look something lke this?
>
>
>
> topo_mask = {
>
> flag = TRUE;
>
> FILE_NAME = ["PYTHON_NUMPY"];
>
> field = [
>
> { name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/terrht_sfc_000000_000000_4a0331x0331_${INIT_TIME}_${FCST_HR}_fcs
> tfld";},
>
> { name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/terrht_sfc_000000_000000_3a0211x0211_${INIT_TIME}_${FCST_HR}_fcs
> tfld";},
>
> { name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/terrht_sfc_000000_000000_2a0271x0271_${INIT_TIME}_${FCST_HR}_fcs
> tfld";},
>
> { name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/terrht_sfc_000000_000000_1a0127x0127_${INIT_TIME}_${FCST_HR}_fcs
> tfld";},
>
> .
>
> ]
>
> regrid = { method = BILIN; width=2;}
>
> use_obs_thresh = ge-500&&le500;
>
> interp_fcst_thresh = ge-50&&le50;
>
> }
>
>
>
> Justin Tsu
>
> Marine Meteorology Division
>
> Data Assimilation/Mesoscale Modeling
>
> Building 704 Room 212
>
> Naval Research Laboratory, Code 7531
>
> 7 Grace Hopper Avenue
>
> Monterey, CA 93943-5502
>
>
>
> Ph. (831) 656-4111
>
>
>
>
>
------------------------------------------------
Subject: Topo_mask for multiple nests in one file
From: Tsu, Mr. Justin
Time: Fri Jan 24 12:14:31 2020
No problem. Thanks for that, I've repurposed my script so that it
will run
each nest individually with its own usage of point_stat. On a side
note I am
running point_stat with verbosity set to 0 and am still getting this
kind of
output:
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=18446744071784442718
Do you know how to suppress this?
Justin
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Friday, January 24, 2020 8:22 AM
To: Tsu, Mr. Justin
Subject: Re: [rt.rap.ucar.edu #93905] Topo_mask for multiple nests in
one file
Justin,
I don't think this will work. When you run Point-Stat, you have to
pick a
single verification domain. By that I mean you pick a single grid on
which
the verification is performed. The "regrid" config file section
allows you
to be flexible on how you define that domain. But ultimately, you
pick a
single one. Things like the land/sea mask, topography, and masking
regions
must all be defined on that grid.
So if you want to evaluate different verification domains, you'll need
to
run Point-Stat separately for each one.
Hope that helps clarify.
Thanks,
John
On Thu, Jan 23, 2020 at 5:05 PM Tsu, Mr. Justin via RT
<met_help at ucar.edu>
wrote:
>
> Thu Jan 23 17:05:36 2020: Request 93905 was acted upon.
> Transaction: Ticket created by justin.tsu at nrlmry.navy.mil
> Queue: met_help
> Subject: Topo_mask for multiple nests in one file
> Owner: Nobody
> Requestors: justin.tsu at nrlmry.navy.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=93905 >
>
>
> Hey John,
>
>
>
> How do I specify the topo mask using python embedding if I have
multiple
> different nests represented in my field/name entry. Say I have
>
>
>
> field = [
>
> {name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_4a0331x0331_${INIT_TIME}_${FCST_HR
> }_fcstfld ${WIND_FLAG}";},
>
> {name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_3a0211x0211_${INIT_TIME}_${FCST_HR
> }_fcstfld ${WIND_FLAG}";},
>
> {name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_2a0271x0271_${INIT_TIME}_${FCST_HR
> }_fcstfld ${WIND_FLAG}";},
>
> {name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_1a0127x0127_${INIT_TIME}_${FCST_HR
> }_fcstfld ${WIND_FLAG}";},
>
> ...
>
> ]
>
>
>
> Would my topo_mask then look something lke this?
>
>
>
> topo_mask = {
>
> flag = TRUE;
>
> FILE_NAME = ["PYTHON_NUMPY"];
>
> field = [
>
> { name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/terrht_sfc_000000_000000_4a0331x0331_${INIT_TIME}_${FCST_HR}_fcs
> tfld";},
>
> { name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/terrht_sfc_000000_000000_3a0211x0211_${INIT_TIME}_${FCST_HR}_fcs
> tfld";},
>
> { name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/terrht_sfc_000000_000000_2a0271x0271_${INIT_TIME}_${FCST_HR}_fcs
> tfld";},
>
> { name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
>
>
${DATA_DIR}/terrht_sfc_000000_000000_1a0127x0127_${INIT_TIME}_${FCST_HR}_fcs
> tfld";},
>
> .
>
> ]
>
> regrid = { method = BILIN; width=2;}
>
> use_obs_thresh = ge-500&&le500;
>
> interp_fcst_thresh = ge-50&&le50;
>
> }
>
>
>
> Justin Tsu
>
> Marine Meteorology Division
>
> Data Assimilation/Mesoscale Modeling
>
> Building 704 Room 212
>
> Naval Research Laboratory, Code 7531
>
> 7 Grace Hopper Avenue
>
> Monterey, CA 93943-5502
>
>
>
> Ph. (831) 656-4111
>
>
>
>
>
------------------------------------------------
Subject: Topo_mask for multiple nests in one file
From: John Halley Gotway
Time: Fri Jan 24 13:13:43 2020
Hah, yes, I agree that's very annoying. That's being written by the
GNU
Scientific Library when MET initializes a random number generator.
And
unfortunately, I don't know how to make it go away.
You could redirect all of the log output to /dev/null... but then you
wound't see MET's error messages.
John
On Fri, Jan 24, 2020 at 12:14 PM Tsu, Mr. Justin via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=93905 >
>
> No problem. Thanks for that, I've repurposed my script so that it
will
> run
> each nest individually with its own usage of point_stat. On a side
note I
> am
> running point_stat with verbosity set to 0 and am still getting this
kind
> of
> output:
>
>
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=18446744071784442718
>
> Do you know how to suppress this?
>
> Justin
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Friday, January 24, 2020 8:22 AM
> To: Tsu, Mr. Justin
> Subject: Re: [rt.rap.ucar.edu #93905] Topo_mask for multiple nests
in one
> file
>
> Justin,
>
> I don't think this will work. When you run Point-Stat, you have to
pick a
> single verification domain. By that I mean you pick a single grid
on which
> the verification is performed. The "regrid" config file section
allows you
> to be flexible on how you define that domain. But ultimately, you
pick a
> single one. Things like the land/sea mask, topography, and masking
regions
> must all be defined on that grid.
>
> So if you want to evaluate different verification domains, you'll
need to
> run Point-Stat separately for each one.
>
> Hope that helps clarify.
>
> Thanks,
> John
>
> On Thu, Jan 23, 2020 at 5:05 PM Tsu, Mr. Justin via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > Thu Jan 23 17:05:36 2020: Request 93905 was acted upon.
> > Transaction: Ticket created by justin.tsu at nrlmry.navy.mil
> > Queue: met_help
> > Subject: Topo_mask for multiple nests in one file
> > Owner: Nobody
> > Requestors: justin.tsu at nrlmry.navy.mil
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=93905 >
> >
> >
> > Hey John,
> >
> >
> >
> > How do I specify the topo mask using python embedding if I have
multiple
> > different nests represented in my field/name entry. Say I have
> >
> >
> >
> > field = [
> >
> > {name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
> >
> >
>
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_4a0331x0331_${INIT_TIME}_${FCST_HR
> > }_fcstfld ${WIND_FLAG}";},
> >
> > {name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
> >
> >
>
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_3a0211x0211_${INIT_TIME}_${FCST_HR
> > }_fcstfld ${WIND_FLAG}";},
> >
> > {name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
> >
> >
>
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_2a0271x0271_${INIT_TIME}_${FCST_HR
> > }_fcstfld ${WIND_FLAG}";},
> >
> > {name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
> >
> >
>
${DATA_DIR}/${VAR_NAME}_pre_000010_000000_1a0127x0127_${INIT_TIME}_${FCST_HR
> > }_fcstfld ${WIND_FLAG}";},
> >
> > ...
> >
> > ]
> >
> >
> >
> > Would my topo_mask then look something lke this?
> >
> >
> >
> > topo_mask = {
> >
> > flag = TRUE;
> >
> > FILE_NAME = ["PYTHON_NUMPY"];
> >
> > field = [
> >
> > { name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
> >
> >
>
${DATA_DIR}/terrht_sfc_000000_000000_4a0331x0331_${INIT_TIME}_${FCST_HR}_fcs
> > tfld";},
> >
> > { name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
> >
> >
>
${DATA_DIR}/terrht_sfc_000000_000000_3a0211x0211_${INIT_TIME}_${FCST_HR}_fcs
> > tfld";},
> >
> > { name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
> >
> >
>
${DATA_DIR}/terrht_sfc_000000_000000_2a0271x0271_${INIT_TIME}_${FCST_HR}_fcs
> > tfld";},
> >
> > { name = "/users/tsu/MET/work/PY_SRC/read_NRL_binary.py
> >
> >
>
${DATA_DIR}/terrht_sfc_000000_000000_1a0127x0127_${INIT_TIME}_${FCST_HR}_fcs
> > tfld";},
> >
> > .
> >
> > ]
> >
> > regrid = { method = BILIN; width=2;}
> >
> > use_obs_thresh = ge-500&&le500;
> >
> > interp_fcst_thresh = ge-50&&le50;
> >
> > }
> >
> >
> >
> > Justin Tsu
> >
> > Marine Meteorology Division
> >
> > Data Assimilation/Mesoscale Modeling
> >
> > Building 704 Room 212
> >
> > Naval Research Laboratory, Code 7531
> >
> > 7 Grace Hopper Avenue
> >
> > Monterey, CA 93943-5502
> >
> >
> >
> > Ph. (831) 656-4111
> >
> >
> >
> >
> >
>
>
>
------------------------------------------------
More information about the Met_help
mailing list