[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