[Met_help] [rt.rap.ucar.edu #83377] History for visibility in MET
John Halley Gotway via RT
met_help at ucar.edu
Thu Jan 4 15:21:18 MST 2018
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hi, there,
I am trying to verify visibility in MET. In my PB2NC config file, I have
the following:
//
// Observation values to retain or derive
//
obs_bufr_var = [ "POB", "TOB", "ZOB", "QOB", "UOB", "VOB", "HOVI" ];
//
// Mapping of BUFR variable name to GRIB name
//
obs_bufr_map = [
{ key = "HOVI"; val = "VIS";
key = "POB"; val = "PRMSL";
key = "ZOB"; val = "HGT";
key = "TOB"; val = "TMP";
key = "QOB"; val = "SPFH";
key = "UOB"; val = "UGRD";
key = "VOB"; val = "VGRD"; }
];
Then I have this in my point_stat config file:
fcst = {
field = [
{
name = "VIS";
level = "Z0";
cat_thresh = [ >=800, >=1600, >=8046, >=16090 ];
},
.
.
.
];
}
obs = fcst;
I still can't get VIS verification in my .stat file. Is there anything I
am missing here, that you can see? Do you need more information?
Thanks!
Perry
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: visibility in MET
From: Minna Win
Time: Thu Dec 21 11:55:46 2017
Hi Perry,
A lot of us are out of the office for the holidays, so you might
experience
a delay in your response until everyone returns on January 2. I don't
have
any expertise in point stat, but I think you need to have your level
enclosed in brackets
ie level = ["Z0"];
John or anyone else with more experience might be able to help you
out.
Regards,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax: 302-497-8401
On Thu, Dec 21, 2017 at 6:38 PM, perry.shafran at noaa.gov via RT <
met_help at ucar.edu> wrote:
>
> Thu Dec 21 11:38:09 2017: Request 83377 was acted upon.
> Transaction: Ticket created by perry.shafran at noaa.gov
> Queue: met_help
> Subject: visibility in MET
> Owner: Nobody
> Requestors: perry.shafran at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
>
>
> Hi, there,
>
> I am trying to verify visibility in MET. In my PB2NC config file, I
have
> the following:
>
> //
> // Observation values to retain or derive
> //
> obs_bufr_var = [ "POB", "TOB", "ZOB", "QOB", "UOB", "VOB", "HOVI" ];
>
> //
> // Mapping of BUFR variable name to GRIB name
> //
> obs_bufr_map = [
> { key = "HOVI"; val = "VIS";
> key = "POB"; val = "PRMSL";
> key = "ZOB"; val = "HGT";
> key = "TOB"; val = "TMP";
> key = "QOB"; val = "SPFH";
> key = "UOB"; val = "UGRD";
> key = "VOB"; val = "VGRD"; }
> ];
>
> Then I have this in my point_stat config file:
>
> fcst = {
> field = [
> {
> name = "VIS";
> level = "Z0";
> cat_thresh = [ >=800, >=1600, >=8046, >=16090 ];
> },
> .
> .
> .
> ];
>
> }
>
> obs = fcst;
>
> I still can't get VIS verification in my .stat file. Is there
anything I
> am missing here, that you can see? Do you need more information?
>
> Thanks!
>
> Perry
>
>
------------------------------------------------
Subject: visibility in MET
From: perry.shafran at noaa.gov
Time: Thu Dec 21 11:57:02 2017
Hi, Minna,
Thanks!
Perry
On Thu, Dec 21, 2017 at 1:55 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
> Hi Perry,
>
> A lot of us are out of the office for the holidays, so you might
experience
> a delay in your response until everyone returns on January 2. I
don't have
> any expertise in point stat, but I think you need to have your level
> enclosed in brackets
>
> ie level = ["Z0"];
>
>
> John or anyone else with more experience might be able to help you
out.
>
> Regards,
> Minna
>
> ---------------
> Minna Win
> NCAR
> Research Applications Lab
> Phone: 303-497-8423
> Fax: 302-497-8401
>
>
> On Thu, Dec 21, 2017 at 6:38 PM, perry.shafran at noaa.gov via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Thu Dec 21 11:38:09 2017: Request 83377 was acted upon.
> > Transaction: Ticket created by perry.shafran at noaa.gov
> > Queue: met_help
> > Subject: visibility in MET
> > Owner: Nobody
> > Requestors: perry.shafran at noaa.gov
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> >
> >
> > Hi, there,
> >
> > I am trying to verify visibility in MET. In my PB2NC config file,
I have
> > the following:
> >
> > //
> > // Observation values to retain or derive
> > //
> > obs_bufr_var = [ "POB", "TOB", "ZOB", "QOB", "UOB", "VOB", "HOVI"
];
> >
> > //
> > // Mapping of BUFR variable name to GRIB name
> > //
> > obs_bufr_map = [
> > { key = "HOVI"; val = "VIS";
> > key = "POB"; val = "PRMSL";
> > key = "ZOB"; val = "HGT";
> > key = "TOB"; val = "TMP";
> > key = "QOB"; val = "SPFH";
> > key = "UOB"; val = "UGRD";
> > key = "VOB"; val = "VGRD"; }
> > ];
> >
> > Then I have this in my point_stat config file:
> >
> > fcst = {
> > field = [
> > {
> > name = "VIS";
> > level = "Z0";
> > cat_thresh = [ >=800, >=1600, >=8046, >=16090 ];
> > },
> > .
> > .
> > .
> > ];
> >
> > }
> >
> > obs = fcst;
> >
> > I still can't get VIS verification in my .stat file. Is there
anything I
> > am missing here, that you can see? Do you need more information?
> >
> > Thanks!
> >
> > Perry
> >
> >
>
>
------------------------------------------------
Subject: visibility in MET
From: perry.shafran at noaa.gov
Time: Thu Dec 21 12:09:15 2017
Hi, Minna,
I'm not sure that really is an issue because I have other variables
(i.e.
TMP) with items not enclosed in brackets, and those appear in the
.stat
file. Part of the issue here I think has to do with the fact that the
variable has one mnemonic for obs (HOVI) and another for the model
(VIS).
I think I did it right the mapping of the BUFR name to the GRIB name,
but
I'm still not getting anything. I think VIS is considered Z0 but
maybe
it's something else? It didn't work when I set the level to L0
either.
Waiting until after the new year for assistance is fine.
Thanks!
Perry
On Thu, Dec 21, 2017 at 1:55 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
> Hi Perry,
>
> A lot of us are out of the office for the holidays, so you might
experience
> a delay in your response until everyone returns on January 2. I
don't have
> any expertise in point stat, but I think you need to have your level
> enclosed in brackets
>
> ie level = ["Z0"];
>
>
> John or anyone else with more experience might be able to help you
out.
>
> Regards,
> Minna
>
> ---------------
> Minna Win
> NCAR
> Research Applications Lab
> Phone: 303-497-8423
> Fax: 302-497-8401
>
>
> On Thu, Dec 21, 2017 at 6:38 PM, perry.shafran at noaa.gov via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Thu Dec 21 11:38:09 2017: Request 83377 was acted upon.
> > Transaction: Ticket created by perry.shafran at noaa.gov
> > Queue: met_help
> > Subject: visibility in MET
> > Owner: Nobody
> > Requestors: perry.shafran at noaa.gov
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> >
> >
> > Hi, there,
> >
> > I am trying to verify visibility in MET. In my PB2NC config file,
I have
> > the following:
> >
> > //
> > // Observation values to retain or derive
> > //
> > obs_bufr_var = [ "POB", "TOB", "ZOB", "QOB", "UOB", "VOB", "HOVI"
];
> >
> > //
> > // Mapping of BUFR variable name to GRIB name
> > //
> > obs_bufr_map = [
> > { key = "HOVI"; val = "VIS";
> > key = "POB"; val = "PRMSL";
> > key = "ZOB"; val = "HGT";
> > key = "TOB"; val = "TMP";
> > key = "QOB"; val = "SPFH";
> > key = "UOB"; val = "UGRD";
> > key = "VOB"; val = "VGRD"; }
> > ];
> >
> > Then I have this in my point_stat config file:
> >
> > fcst = {
> > field = [
> > {
> > name = "VIS";
> > level = "Z0";
> > cat_thresh = [ >=800, >=1600, >=8046, >=16090 ];
> > },
> > .
> > .
> > .
> > ];
> >
> > }
> >
> > obs = fcst;
> >
> > I still can't get VIS verification in my .stat file. Is there
anything I
> > am missing here, that you can see? Do you need more information?
> >
> > Thanks!
> >
> > Perry
> >
> >
>
>
------------------------------------------------
Subject: visibility in MET
From: Minna Win
Time: Thu Dec 21 12:29:45 2017
Hi Perry,
To check if VIS is Z0, you could use plot_data_plane to see if it is:
plot_data_plane \
<full path to input file> \
<full path and name of temporary ps file> \
'name="VIS"; level="Z0";' -v 4
Setting the verbosity to 4 provides enough information to indicate
whether
this variable is found, or if the level does/doesn't exist. This has
been a
helpful "trick" we
use on another project which has variable names and levels that are
not
well-documented.
It might be worth trying- it's up to you.
Best of luck!
Thanks,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax: 302-497-8401
On Thu, Dec 21, 2017 at 7:09 PM, perry.shafran at noaa.gov via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
>
> Hi, Minna,
>
> I'm not sure that really is an issue because I have other variables
(i.e.
> TMP) with items not enclosed in brackets, and those appear in the
.stat
> file. Part of the issue here I think has to do with the fact that
the
> variable has one mnemonic for obs (HOVI) and another for the model
(VIS).
> I think I did it right the mapping of the BUFR name to the GRIB
name, but
> I'm still not getting anything. I think VIS is considered Z0 but
maybe
> it's something else? It didn't work when I set the level to L0
either.
>
> Waiting until after the new year for assistance is fine.
>
> Thanks!
>
> Perry
>
> On Thu, Dec 21, 2017 at 1:55 PM, Minna Win via RT
<met_help at ucar.edu>
> wrote:
>
> > Hi Perry,
> >
> > A lot of us are out of the office for the holidays, so you might
> experience
> > a delay in your response until everyone returns on January 2. I
don't
> have
> > any expertise in point stat, but I think you need to have your
level
> > enclosed in brackets
> >
> > ie level = ["Z0"];
> >
> >
> > John or anyone else with more experience might be able to help you
out.
> >
> > Regards,
> > Minna
> >
> > ---------------
> > Minna Win
> > NCAR
> > Research Applications Lab
> > Phone: 303-497-8423
> > Fax: 302-497-8401
> >
> >
> > On Thu, Dec 21, 2017 at 6:38 PM, perry.shafran at noaa.gov via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > Thu Dec 21 11:38:09 2017: Request 83377 was acted upon.
> > > Transaction: Ticket created by perry.shafran at noaa.gov
> > > Queue: met_help
> > > Subject: visibility in MET
> > > Owner: Nobody
> > > Requestors: perry.shafran at noaa.gov
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377
> >
> > >
> > >
> > > Hi, there,
> > >
> > > I am trying to verify visibility in MET. In my PB2NC config
file, I
> have
> > > the following:
> > >
> > > //
> > > // Observation values to retain or derive
> > > //
> > > obs_bufr_var = [ "POB", "TOB", "ZOB", "QOB", "UOB", "VOB",
"HOVI" ];
> > >
> > > //
> > > // Mapping of BUFR variable name to GRIB name
> > > //
> > > obs_bufr_map = [
> > > { key = "HOVI"; val = "VIS";
> > > key = "POB"; val = "PRMSL";
> > > key = "ZOB"; val = "HGT";
> > > key = "TOB"; val = "TMP";
> > > key = "QOB"; val = "SPFH";
> > > key = "UOB"; val = "UGRD";
> > > key = "VOB"; val = "VGRD"; }
> > > ];
> > >
> > > Then I have this in my point_stat config file:
> > >
> > > fcst = {
> > > field = [
> > > {
> > > name = "VIS";
> > > level = "Z0";
> > > cat_thresh = [ >=800, >=1600, >=8046, >=16090 ];
> > > },
> > > .
> > > .
> > > .
> > > ];
> > >
> > > }
> > >
> > > obs = fcst;
> > >
> > > I still can't get VIS verification in my .stat file. Is there
> anything I
> > > am missing here, that you can see? Do you need more
information?
> > >
> > > Thanks!
> > >
> > > Perry
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: visibility in MET
From: John Halley Gotway
Time: Thu Dec 28 13:52:23 2017
Hi Perry,
Not sure if you're in the office today. If so, do you still have
questions
about this issue?
Thanks,
John
On Thu, Dec 21, 2017 at 12:29 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
>
> Hi Perry,
>
> To check if VIS is Z0, you could use plot_data_plane to see if it
is:
>
> plot_data_plane \
> <full path to input file> \
> <full path and name of temporary ps file> \
> 'name="VIS"; level="Z0";' -v 4
>
> Setting the verbosity to 4 provides enough information to indicate
whether
> this variable is found, or if the level does/doesn't exist. This has
been a
> helpful "trick" we
> use on another project which has variable names and levels that are
not
> well-documented.
>
> It might be worth trying- it's up to you.
>
> Best of luck!
>
> Thanks,
> Minna
>
>
>
>
> ---------------
> Minna Win
> NCAR
> Research Applications Lab
> Phone: 303-497-8423
> Fax: 302-497-8401
>
>
> On Thu, Dec 21, 2017 at 7:09 PM, perry.shafran at noaa.gov via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> >
> > Hi, Minna,
> >
> > I'm not sure that really is an issue because I have other
variables (i.e.
> > TMP) with items not enclosed in brackets, and those appear in the
.stat
> > file. Part of the issue here I think has to do with the fact that
the
> > variable has one mnemonic for obs (HOVI) and another for the model
(VIS).
> > I think I did it right the mapping of the BUFR name to the GRIB
name, but
> > I'm still not getting anything. I think VIS is considered Z0 but
maybe
> > it's something else? It didn't work when I set the level to L0
either.
> >
> > Waiting until after the new year for assistance is fine.
> >
> > Thanks!
> >
> > Perry
> >
> > On Thu, Dec 21, 2017 at 1:55 PM, Minna Win via RT
<met_help at ucar.edu>
> > wrote:
> >
> > > Hi Perry,
> > >
> > > A lot of us are out of the office for the holidays, so you might
> > experience
> > > a delay in your response until everyone returns on January 2. I
don't
> > have
> > > any expertise in point stat, but I think you need to have your
level
> > > enclosed in brackets
> > >
> > > ie level = ["Z0"];
> > >
> > >
> > > John or anyone else with more experience might be able to help
you out.
> > >
> > > Regards,
> > > Minna
> > >
> > > ---------------
> > > Minna Win
> > > NCAR
> > > Research Applications Lab
> > > Phone: 303-497-8423
> > > Fax: 302-497-8401
> > >
> > >
> > > On Thu, Dec 21, 2017 at 6:38 PM, perry.shafran at noaa.gov via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > Thu Dec 21 11:38:09 2017: Request 83377 was acted upon.
> > > > Transaction: Ticket created by perry.shafran at noaa.gov
> > > > Queue: met_help
> > > > Subject: visibility in MET
> > > > Owner: Nobody
> > > > Requestors: perry.shafran at noaa.gov
> > > > Status: new
> > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=83377
> > >
> > > >
> > > >
> > > > Hi, there,
> > > >
> > > > I am trying to verify visibility in MET. In my PB2NC config
file, I
> > have
> > > > the following:
> > > >
> > > > //
> > > > // Observation values to retain or derive
> > > > //
> > > > obs_bufr_var = [ "POB", "TOB", "ZOB", "QOB", "UOB", "VOB",
"HOVI" ];
> > > >
> > > > //
> > > > // Mapping of BUFR variable name to GRIB name
> > > > //
> > > > obs_bufr_map = [
> > > > { key = "HOVI"; val = "VIS";
> > > > key = "POB"; val = "PRMSL";
> > > > key = "ZOB"; val = "HGT";
> > > > key = "TOB"; val = "TMP";
> > > > key = "QOB"; val = "SPFH";
> > > > key = "UOB"; val = "UGRD";
> > > > key = "VOB"; val = "VGRD"; }
> > > > ];
> > > >
> > > > Then I have this in my point_stat config file:
> > > >
> > > > fcst = {
> > > > field = [
> > > > {
> > > > name = "VIS";
> > > > level = "Z0";
> > > > cat_thresh = [ >=800, >=1600, >=8046, >=16090 ];
> > > > },
> > > > .
> > > > .
> > > > .
> > > > ];
> > > >
> > > > }
> > > >
> > > > obs = fcst;
> > > >
> > > > I still can't get VIS verification in my .stat file. Is there
> > anything I
> > > > am missing here, that you can see? Do you need more
information?
> > > >
> > > > Thanks!
> > > >
> > > > Perry
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: visibility in MET
From: perry.shafran at noaa.gov
Time: Tue Jan 02 06:14:00 2018
Hi, John,
I was away all week so I might still have some questions. If you have
an
answer to the original question about the level for visibility, that
would
be great.
Thanks!
Perry
On Thu, Dec 28, 2017 at 3:52 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Hi Perry,
>
> Not sure if you're in the office today. If so, do you still have
questions
> about this issue?
>
> Thanks,
> John
>
> On Thu, Dec 21, 2017 at 12:29 PM, Minna Win via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> >
> > Hi Perry,
> >
> > To check if VIS is Z0, you could use plot_data_plane to see if it
is:
> >
> > plot_data_plane \
> > <full path to input file> \
> > <full path and name of temporary ps file> \
> > 'name="VIS"; level="Z0";' -v 4
> >
> > Setting the verbosity to 4 provides enough information to indicate
> whether
> > this variable is found, or if the level does/doesn't exist. This
has
> been a
> > helpful "trick" we
> > use on another project which has variable names and levels that
are not
> > well-documented.
> >
> > It might be worth trying- it's up to you.
> >
> > Best of luck!
> >
> > Thanks,
> > Minna
> >
> >
> >
> >
> > ---------------
> > Minna Win
> > NCAR
> > Research Applications Lab
> > Phone: 303-497-8423
> > Fax: 302-497-8401
> >
> >
> > On Thu, Dec 21, 2017 at 7:09 PM, perry.shafran at noaa.gov via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> > >
> > > Hi, Minna,
> > >
> > > I'm not sure that really is an issue because I have other
variables
> (i.e.
> > > TMP) with items not enclosed in brackets, and those appear in
the .stat
> > > file. Part of the issue here I think has to do with the fact
that the
> > > variable has one mnemonic for obs (HOVI) and another for the
model
> (VIS).
> > > I think I did it right the mapping of the BUFR name to the GRIB
name,
> but
> > > I'm still not getting anything. I think VIS is considered Z0
but maybe
> > > it's something else? It didn't work when I set the level to L0
either.
> > >
> > > Waiting until after the new year for assistance is fine.
> > >
> > > Thanks!
> > >
> > > Perry
> > >
> > > On Thu, Dec 21, 2017 at 1:55 PM, Minna Win via RT
<met_help at ucar.edu>
> > > wrote:
> > >
> > > > Hi Perry,
> > > >
> > > > A lot of us are out of the office for the holidays, so you
might
> > > experience
> > > > a delay in your response until everyone returns on January 2.
I
> don't
> > > have
> > > > any expertise in point stat, but I think you need to have your
level
> > > > enclosed in brackets
> > > >
> > > > ie level = ["Z0"];
> > > >
> > > >
> > > > John or anyone else with more experience might be able to help
you
> out.
> > > >
> > > > Regards,
> > > > Minna
> > > >
> > > > ---------------
> > > > Minna Win
> > > > NCAR
> > > > Research Applications Lab
> > > > Phone: 303-497-8423
> > > > Fax: 302-497-8401
> > > >
> > > >
> > > > On Thu, Dec 21, 2017 at 6:38 PM, perry.shafran at noaa.gov via RT
<
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > Thu Dec 21 11:38:09 2017: Request 83377 was acted upon.
> > > > > Transaction: Ticket created by perry.shafran at noaa.gov
> > > > > Queue: met_help
> > > > > Subject: visibility in MET
> > > > > Owner: Nobody
> > > > > Requestors: perry.shafran at noaa.gov
> > > > > Status: new
> > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > Ticket/Display.html?id=83377
> > > >
> > > > >
> > > > >
> > > > > Hi, there,
> > > > >
> > > > > I am trying to verify visibility in MET. In my PB2NC config
file,
> I
> > > have
> > > > > the following:
> > > > >
> > > > > //
> > > > > // Observation values to retain or derive
> > > > > //
> > > > > obs_bufr_var = [ "POB", "TOB", "ZOB", "QOB", "UOB", "VOB",
"HOVI"
> ];
> > > > >
> > > > > //
> > > > > // Mapping of BUFR variable name to GRIB name
> > > > > //
> > > > > obs_bufr_map = [
> > > > > { key = "HOVI"; val = "VIS";
> > > > > key = "POB"; val = "PRMSL";
> > > > > key = "ZOB"; val = "HGT";
> > > > > key = "TOB"; val = "TMP";
> > > > > key = "QOB"; val = "SPFH";
> > > > > key = "UOB"; val = "UGRD";
> > > > > key = "VOB"; val = "VGRD"; }
> > > > > ];
> > > > >
> > > > > Then I have this in my point_stat config file:
> > > > >
> > > > > fcst = {
> > > > > field = [
> > > > > {
> > > > > name = "VIS";
> > > > > level = "Z0";
> > > > > cat_thresh = [ >=800, >=1600, >=8046, >=16090 ];
> > > > > },
> > > > > .
> > > > > .
> > > > > .
> > > > > ];
> > > > >
> > > > > }
> > > > >
> > > > > obs = fcst;
> > > > >
> > > > > I still can't get VIS verification in my .stat file. Is
there
> > > anything I
> > > > > am missing here, that you can see? Do you need more
information?
> > > > >
> > > > > Thanks!
> > > > >
> > > > > Perry
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: visibility in MET
From: John Halley Gotway
Time: Tue Jan 02 09:26:45 2018
Hi Perry,
OK, I can take a closer look at this today. Since there are a lot of
potential reasons for 0 VIS matched pairs, I think it'd be best for me
to
test with the exact PB2NC and Point-Stat config files that you're
using.
Would you be able to either copy some data to theia for me to test or
post
it to our anonymous ftp site (
https://dtcenter.org/met/users/support/met_help.php#ftp)?
I'll need your PB2NC config file, Point-Stat config file, the PREPBUFR
input file for PB2NC, and the forecast GRIB file that you're
verifying.
I'll replicate the commands your running and then figure out why
you're
getting 0 matched pairs for visibility.
Thanks,
John
On Tue, Jan 2, 2018 at 6:14 AM, perry.shafran at noaa.gov via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
>
> Hi, John,
>
> I was away all week so I might still have some questions. If you
have an
> answer to the original question about the level for visibility, that
would
> be great.
>
> Thanks!
>
> Perry
>
> On Thu, Dec 28, 2017 at 3:52 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > Hi Perry,
> >
> > Not sure if you're in the office today. If so, do you still have
> questions
> > about this issue?
> >
> > Thanks,
> > John
> >
> > On Thu, Dec 21, 2017 at 12:29 PM, Minna Win via RT
<met_help at ucar.edu>
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> > >
> > > Hi Perry,
> > >
> > > To check if VIS is Z0, you could use plot_data_plane to see if
it is:
> > >
> > > plot_data_plane \
> > > <full path to input file> \
> > > <full path and name of temporary ps file> \
> > > 'name="VIS"; level="Z0";' -v 4
> > >
> > > Setting the verbosity to 4 provides enough information to
indicate
> > whether
> > > this variable is found, or if the level does/doesn't exist. This
has
> > been a
> > > helpful "trick" we
> > > use on another project which has variable names and levels that
are not
> > > well-documented.
> > >
> > > It might be worth trying- it's up to you.
> > >
> > > Best of luck!
> > >
> > > Thanks,
> > > Minna
> > >
> > >
> > >
> > >
> > > ---------------
> > > Minna Win
> > > NCAR
> > > Research Applications Lab
> > > Phone: 303-497-8423
> > > Fax: 302-497-8401
> > >
> > >
> > > On Thu, Dec 21, 2017 at 7:09 PM, perry.shafran at noaa.gov via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377
>
> > > >
> > > > Hi, Minna,
> > > >
> > > > I'm not sure that really is an issue because I have other
variables
> > (i.e.
> > > > TMP) with items not enclosed in brackets, and those appear in
the
> .stat
> > > > file. Part of the issue here I think has to do with the fact
that
> the
> > > > variable has one mnemonic for obs (HOVI) and another for the
model
> > (VIS).
> > > > I think I did it right the mapping of the BUFR name to the
GRIB name,
> > but
> > > > I'm still not getting anything. I think VIS is considered Z0
but
> maybe
> > > > it's something else? It didn't work when I set the level to
L0
> either.
> > > >
> > > > Waiting until after the new year for assistance is fine.
> > > >
> > > > Thanks!
> > > >
> > > > Perry
> > > >
> > > > On Thu, Dec 21, 2017 at 1:55 PM, Minna Win via RT
<met_help at ucar.edu
> >
> > > > wrote:
> > > >
> > > > > Hi Perry,
> > > > >
> > > > > A lot of us are out of the office for the holidays, so you
might
> > > > experience
> > > > > a delay in your response until everyone returns on January
2. I
> > don't
> > > > have
> > > > > any expertise in point stat, but I think you need to have
your
> level
> > > > > enclosed in brackets
> > > > >
> > > > > ie level = ["Z0"];
> > > > >
> > > > >
> > > > > John or anyone else with more experience might be able to
help you
> > out.
> > > > >
> > > > > Regards,
> > > > > Minna
> > > > >
> > > > > ---------------
> > > > > Minna Win
> > > > > NCAR
> > > > > Research Applications Lab
> > > > > Phone: 303-497-8423
> > > > > Fax: 302-497-8401
> > > > >
> > > > >
> > > > > On Thu, Dec 21, 2017 at 6:38 PM, perry.shafran at noaa.gov via
RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > >
> > > > > > Thu Dec 21 11:38:09 2017: Request 83377 was acted upon.
> > > > > > Transaction: Ticket created by perry.shafran at noaa.gov
> > > > > > Queue: met_help
> > > > > > Subject: visibility in MET
> > > > > > Owner: Nobody
> > > > > > Requestors: perry.shafran at noaa.gov
> > > > > > Status: new
> > > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > Ticket/Display.html?id=83377
> > > > >
> > > > > >
> > > > > >
> > > > > > Hi, there,
> > > > > >
> > > > > > I am trying to verify visibility in MET. In my PB2NC
config
> file,
> > I
> > > > have
> > > > > > the following:
> > > > > >
> > > > > > //
> > > > > > // Observation values to retain or derive
> > > > > > //
> > > > > > obs_bufr_var = [ "POB", "TOB", "ZOB", "QOB", "UOB", "VOB",
"HOVI"
> > ];
> > > > > >
> > > > > > //
> > > > > > // Mapping of BUFR variable name to GRIB name
> > > > > > //
> > > > > > obs_bufr_map = [
> > > > > > { key = "HOVI"; val = "VIS";
> > > > > > key = "POB"; val = "PRMSL";
> > > > > > key = "ZOB"; val = "HGT";
> > > > > > key = "TOB"; val = "TMP";
> > > > > > key = "QOB"; val = "SPFH";
> > > > > > key = "UOB"; val = "UGRD";
> > > > > > key = "VOB"; val = "VGRD"; }
> > > > > > ];
> > > > > >
> > > > > > Then I have this in my point_stat config file:
> > > > > >
> > > > > > fcst = {
> > > > > > field = [
> > > > > > {
> > > > > > name = "VIS";
> > > > > > level = "Z0";
> > > > > > cat_thresh = [ >=800, >=1600, >=8046, >=16090 ];
> > > > > > },
> > > > > > .
> > > > > > .
> > > > > > .
> > > > > > ];
> > > > > >
> > > > > > }
> > > > > >
> > > > > > obs = fcst;
> > > > > >
> > > > > > I still can't get VIS verification in my .stat file. Is
there
> > > > anything I
> > > > > > am missing here, that you can see? Do you need more
information?
> > > > > >
> > > > > > Thanks!
> > > > > >
> > > > > > Perry
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: visibility in MET
From: perry.shafran at noaa.gov
Time: Tue Jan 02 09:32:18 2018
Hi, John,
OK, great. I'll try to put something together for you to look at this
afternoon.
Thanks!
Perry
On Tue, Jan 2, 2018 at 11:26 AM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Hi Perry,
>
> OK, I can take a closer look at this today. Since there are a lot
of
> potential reasons for 0 VIS matched pairs, I think it'd be best for
me to
> test with the exact PB2NC and Point-Stat config files that you're
using.
>
> Would you be able to either copy some data to theia for me to test
or post
> it to our anonymous ftp site (
> https://dtcenter.org/met/users/support/met_help.php#ftp)?
>
> I'll need your PB2NC config file, Point-Stat config file, the
PREPBUFR
> input file for PB2NC, and the forecast GRIB file that you're
verifying.
>
> I'll replicate the commands your running and then figure out why
you're
> getting 0 matched pairs for visibility.
>
> Thanks,
> John
>
>
>
> On Tue, Jan 2, 2018 at 6:14 AM, perry.shafran at noaa.gov via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> >
> > Hi, John,
> >
> > I was away all week so I might still have some questions. If you
have an
> > answer to the original question about the level for visibility,
that
> would
> > be great.
> >
> > Thanks!
> >
> > Perry
> >
> > On Thu, Dec 28, 2017 at 3:52 PM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> > > Hi Perry,
> > >
> > > Not sure if you're in the office today. If so, do you still
have
> > questions
> > > about this issue?
> > >
> > > Thanks,
> > > John
> > >
> > > On Thu, Dec 21, 2017 at 12:29 PM, Minna Win via RT
<met_help at ucar.edu>
> > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377
>
> > > >
> > > > Hi Perry,
> > > >
> > > > To check if VIS is Z0, you could use plot_data_plane to see if
it is:
> > > >
> > > > plot_data_plane \
> > > > <full path to input file> \
> > > > <full path and name of temporary ps file> \
> > > > 'name="VIS"; level="Z0";' -v 4
> > > >
> > > > Setting the verbosity to 4 provides enough information to
indicate
> > > whether
> > > > this variable is found, or if the level does/doesn't exist.
This has
> > > been a
> > > > helpful "trick" we
> > > > use on another project which has variable names and levels
that are
> not
> > > > well-documented.
> > > >
> > > > It might be worth trying- it's up to you.
> > > >
> > > > Best of luck!
> > > >
> > > > Thanks,
> > > > Minna
> > > >
> > > >
> > > >
> > > >
> > > > ---------------
> > > > Minna Win
> > > > NCAR
> > > > Research Applications Lab
> > > > Phone: 303-497-8423
> > > > Fax: 302-497-8401
> > > >
> > > >
> > > > On Thu, Dec 21, 2017 at 7:09 PM, perry.shafran at noaa.gov via RT
<
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> > > > >
> > > > > Hi, Minna,
> > > > >
> > > > > I'm not sure that really is an issue because I have other
variables
> > > (i.e.
> > > > > TMP) with items not enclosed in brackets, and those appear
in the
> > .stat
> > > > > file. Part of the issue here I think has to do with the
fact that
> > the
> > > > > variable has one mnemonic for obs (HOVI) and another for the
model
> > > (VIS).
> > > > > I think I did it right the mapping of the BUFR name to the
GRIB
> name,
> > > but
> > > > > I'm still not getting anything. I think VIS is considered
Z0 but
> > maybe
> > > > > it's something else? It didn't work when I set the level to
L0
> > either.
> > > > >
> > > > > Waiting until after the new year for assistance is fine.
> > > > >
> > > > > Thanks!
> > > > >
> > > > > Perry
> > > > >
> > > > > On Thu, Dec 21, 2017 at 1:55 PM, Minna Win via RT <
> met_help at ucar.edu
> > >
> > > > > wrote:
> > > > >
> > > > > > Hi Perry,
> > > > > >
> > > > > > A lot of us are out of the office for the holidays, so you
might
> > > > > experience
> > > > > > a delay in your response until everyone returns on January
2. I
> > > don't
> > > > > have
> > > > > > any expertise in point stat, but I think you need to have
your
> > level
> > > > > > enclosed in brackets
> > > > > >
> > > > > > ie level = ["Z0"];
> > > > > >
> > > > > >
> > > > > > John or anyone else with more experience might be able to
help
> you
> > > out.
> > > > > >
> > > > > > Regards,
> > > > > > Minna
> > > > > >
> > > > > > ---------------
> > > > > > Minna Win
> > > > > > NCAR
> > > > > > Research Applications Lab
> > > > > > Phone: 303-497-8423
> > > > > > Fax: 302-497-8401
> > > > > >
> > > > > >
> > > > > > On Thu, Dec 21, 2017 at 6:38 PM, perry.shafran at noaa.gov
via RT <
> > > > > > met_help at ucar.edu> wrote:
> > > > > >
> > > > > > >
> > > > > > > Thu Dec 21 11:38:09 2017: Request 83377 was acted upon.
> > > > > > > Transaction: Ticket created by perry.shafran at noaa.gov
> > > > > > > Queue: met_help
> > > > > > > Subject: visibility in MET
> > > > > > > Owner: Nobody
> > > > > > > Requestors: perry.shafran at noaa.gov
> > > > > > > Status: new
> > > > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > Ticket/Display.html?id=83377
> > > > > >
> > > > > > >
> > > > > > >
> > > > > > > Hi, there,
> > > > > > >
> > > > > > > I am trying to verify visibility in MET. In my PB2NC
config
> > file,
> > > I
> > > > > have
> > > > > > > the following:
> > > > > > >
> > > > > > > //
> > > > > > > // Observation values to retain or derive
> > > > > > > //
> > > > > > > obs_bufr_var = [ "POB", "TOB", "ZOB", "QOB", "UOB",
"VOB",
> "HOVI"
> > > ];
> > > > > > >
> > > > > > > //
> > > > > > > // Mapping of BUFR variable name to GRIB name
> > > > > > > //
> > > > > > > obs_bufr_map = [
> > > > > > > { key = "HOVI"; val = "VIS";
> > > > > > > key = "POB"; val = "PRMSL";
> > > > > > > key = "ZOB"; val = "HGT";
> > > > > > > key = "TOB"; val = "TMP";
> > > > > > > key = "QOB"; val = "SPFH";
> > > > > > > key = "UOB"; val = "UGRD";
> > > > > > > key = "VOB"; val = "VGRD"; }
> > > > > > > ];
> > > > > > >
> > > > > > > Then I have this in my point_stat config file:
> > > > > > >
> > > > > > > fcst = {
> > > > > > > field = [
> > > > > > > {
> > > > > > > name = "VIS";
> > > > > > > level = "Z0";
> > > > > > > cat_thresh = [ >=800, >=1600, >=8046, >=16090 ];
> > > > > > > },
> > > > > > > .
> > > > > > > .
> > > > > > > .
> > > > > > > ];
> > > > > > >
> > > > > > > }
> > > > > > >
> > > > > > > obs = fcst;
> > > > > > >
> > > > > > > I still can't get VIS verification in my .stat file. Is
there
> > > > > anything I
> > > > > > > am missing here, that you can see? Do you need more
> information?
> > > > > > >
> > > > > > > Thanks!
> > > > > > >
> > > > > > > Perry
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: visibility in MET
From: perry.shafran at noaa.gov
Time: Tue Jan 02 11:38:30 2018
Hi, John,
Should have everything you need on theia in the directory
/scratch4/NCEPDEV/meso/noscrub/Perry.Shafran/for_john.
Thanks!
Perry
On Tue, Jan 2, 2018 at 11:26 AM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Hi Perry,
>
> OK, I can take a closer look at this today. Since there are a lot
of
> potential reasons for 0 VIS matched pairs, I think it'd be best for
me to
> test with the exact PB2NC and Point-Stat config files that you're
using.
>
> Would you be able to either copy some data to theia for me to test
or post
> it to our anonymous ftp site (
> https://dtcenter.org/met/users/support/met_help.php#ftp)?
>
> I'll need your PB2NC config file, Point-Stat config file, the
PREPBUFR
> input file for PB2NC, and the forecast GRIB file that you're
verifying.
>
> I'll replicate the commands your running and then figure out why
you're
> getting 0 matched pairs for visibility.
>
> Thanks,
> John
>
>
>
> On Tue, Jan 2, 2018 at 6:14 AM, perry.shafran at noaa.gov via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> >
> > Hi, John,
> >
> > I was away all week so I might still have some questions. If you
have an
> > answer to the original question about the level for visibility,
that
> would
> > be great.
> >
> > Thanks!
> >
> > Perry
> >
> > On Thu, Dec 28, 2017 at 3:52 PM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> > > Hi Perry,
> > >
> > > Not sure if you're in the office today. If so, do you still
have
> > questions
> > > about this issue?
> > >
> > > Thanks,
> > > John
> > >
> > > On Thu, Dec 21, 2017 at 12:29 PM, Minna Win via RT
<met_help at ucar.edu>
> > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377
>
> > > >
> > > > Hi Perry,
> > > >
> > > > To check if VIS is Z0, you could use plot_data_plane to see if
it is:
> > > >
> > > > plot_data_plane \
> > > > <full path to input file> \
> > > > <full path and name of temporary ps file> \
> > > > 'name="VIS"; level="Z0";' -v 4
> > > >
> > > > Setting the verbosity to 4 provides enough information to
indicate
> > > whether
> > > > this variable is found, or if the level does/doesn't exist.
This has
> > > been a
> > > > helpful "trick" we
> > > > use on another project which has variable names and levels
that are
> not
> > > > well-documented.
> > > >
> > > > It might be worth trying- it's up to you.
> > > >
> > > > Best of luck!
> > > >
> > > > Thanks,
> > > > Minna
> > > >
> > > >
> > > >
> > > >
> > > > ---------------
> > > > Minna Win
> > > > NCAR
> > > > Research Applications Lab
> > > > Phone: 303-497-8423
> > > > Fax: 302-497-8401
> > > >
> > > >
> > > > On Thu, Dec 21, 2017 at 7:09 PM, perry.shafran at noaa.gov via RT
<
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> > > > >
> > > > > Hi, Minna,
> > > > >
> > > > > I'm not sure that really is an issue because I have other
variables
> > > (i.e.
> > > > > TMP) with items not enclosed in brackets, and those appear
in the
> > .stat
> > > > > file. Part of the issue here I think has to do with the
fact that
> > the
> > > > > variable has one mnemonic for obs (HOVI) and another for the
model
> > > (VIS).
> > > > > I think I did it right the mapping of the BUFR name to the
GRIB
> name,
> > > but
> > > > > I'm still not getting anything. I think VIS is considered
Z0 but
> > maybe
> > > > > it's something else? It didn't work when I set the level to
L0
> > either.
> > > > >
> > > > > Waiting until after the new year for assistance is fine.
> > > > >
> > > > > Thanks!
> > > > >
> > > > > Perry
> > > > >
> > > > > On Thu, Dec 21, 2017 at 1:55 PM, Minna Win via RT <
> met_help at ucar.edu
> > >
> > > > > wrote:
> > > > >
> > > > > > Hi Perry,
> > > > > >
> > > > > > A lot of us are out of the office for the holidays, so you
might
> > > > > experience
> > > > > > a delay in your response until everyone returns on January
2. I
> > > don't
> > > > > have
> > > > > > any expertise in point stat, but I think you need to have
your
> > level
> > > > > > enclosed in brackets
> > > > > >
> > > > > > ie level = ["Z0"];
> > > > > >
> > > > > >
> > > > > > John or anyone else with more experience might be able to
help
> you
> > > out.
> > > > > >
> > > > > > Regards,
> > > > > > Minna
> > > > > >
> > > > > > ---------------
> > > > > > Minna Win
> > > > > > NCAR
> > > > > > Research Applications Lab
> > > > > > Phone: 303-497-8423
> > > > > > Fax: 302-497-8401
> > > > > >
> > > > > >
> > > > > > On Thu, Dec 21, 2017 at 6:38 PM, perry.shafran at noaa.gov
via RT <
> > > > > > met_help at ucar.edu> wrote:
> > > > > >
> > > > > > >
> > > > > > > Thu Dec 21 11:38:09 2017: Request 83377 was acted upon.
> > > > > > > Transaction: Ticket created by perry.shafran at noaa.gov
> > > > > > > Queue: met_help
> > > > > > > Subject: visibility in MET
> > > > > > > Owner: Nobody
> > > > > > > Requestors: perry.shafran at noaa.gov
> > > > > > > Status: new
> > > > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > Ticket/Display.html?id=83377
> > > > > >
> > > > > > >
> > > > > > >
> > > > > > > Hi, there,
> > > > > > >
> > > > > > > I am trying to verify visibility in MET. In my PB2NC
config
> > file,
> > > I
> > > > > have
> > > > > > > the following:
> > > > > > >
> > > > > > > //
> > > > > > > // Observation values to retain or derive
> > > > > > > //
> > > > > > > obs_bufr_var = [ "POB", "TOB", "ZOB", "QOB", "UOB",
"VOB",
> "HOVI"
> > > ];
> > > > > > >
> > > > > > > //
> > > > > > > // Mapping of BUFR variable name to GRIB name
> > > > > > > //
> > > > > > > obs_bufr_map = [
> > > > > > > { key = "HOVI"; val = "VIS";
> > > > > > > key = "POB"; val = "PRMSL";
> > > > > > > key = "ZOB"; val = "HGT";
> > > > > > > key = "TOB"; val = "TMP";
> > > > > > > key = "QOB"; val = "SPFH";
> > > > > > > key = "UOB"; val = "UGRD";
> > > > > > > key = "VOB"; val = "VGRD"; }
> > > > > > > ];
> > > > > > >
> > > > > > > Then I have this in my point_stat config file:
> > > > > > >
> > > > > > > fcst = {
> > > > > > > field = [
> > > > > > > {
> > > > > > > name = "VIS";
> > > > > > > level = "Z0";
> > > > > > > cat_thresh = [ >=800, >=1600, >=8046, >=16090 ];
> > > > > > > },
> > > > > > > .
> > > > > > > .
> > > > > > > .
> > > > > > > ];
> > > > > > >
> > > > > > > }
> > > > > > >
> > > > > > > obs = fcst;
> > > > > > >
> > > > > > > I still can't get VIS verification in my .stat file. Is
there
> > > > > anything I
> > > > > > > am missing here, that you can see? Do you need more
> information?
> > > > > > >
> > > > > > > Thanks!
> > > > > > >
> > > > > > > Perry
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: visibility in MET
From: John Halley Gotway
Time: Tue Jan 02 14:14:13 2018
Perry,
Great, thanks for sending that to me. I took a look at the PB2NC
config
file and immediately saw a problem. The "obs_bufr_map" is supposed to
be
an array of dictionaries. Each dictionary is surrounded by curly
braces
and defines a key and a value. So basically, the "obs_bufr_map" entry
isn't formatted correctly. Instead, it should look like this:
obs_bufr_map = [
{ key = "HOVI"; val = "VIS"; },
{ key = "POB"; val = "PRMSL"; },
{ key = "ZOB"; val = "HGT"; },
{ key = "TOB"; val = "TMP"; },
{ key = "QOB"; val = "SPFH"; },
{ key = "UOB"; val = "UGRD"; },
{ key = "VOB"; val = "VGRD"; }
];
To make sure this did the trick, I ran the following commands:
# Load MET module
module use /contrib/modulefiles
module load met/6.1
# Run PB2NC using your original settings
pb2nc prepda.2018010100 prepda.2018010100_ORIG.nc PB2NCConfig
# Run plot_point_obs to plot VIS observations (I expect there to be 0
of
them)
plot_point_obs prepda.2018010100_ORIG.nc prepda.2018010100_ORIG.ps
-obs_var
VIS
# I log messages from plot_point_obs confirm that there are no
observations
named VIS in the output:
# DEBUG 2: Finished plotting 0 locations.
# Edit the config file to define obs_bufr_map as listed above and then
rerun PB2NC
pb2nc prepda.2018010100 prepda.2018010100_FIX.nc PB2NCConfig_FIX
# Rerun plot_point_obs and make sure that there really are VIS
observations
present
plot_point_obs prepda.2018010100_FIX.nc prepda.2018010100_FIX.ps
-obs_var
VIS
# And the log message now confirm that we do have observations named
VIS:
# DEBUG 2: Finished plotting 3657 locations.
So that change to the PB2NC config file enabled PB2NC to output
observations named VIS as 3657 locations. Hopefully this change will
fix
your issue.
Thanks,
John
Thanks,
John
On Tue, Jan 2, 2018 at 11:38 AM, perry.shafran at noaa.gov via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
>
> Hi, John,
>
> Should have everything you need on theia in the directory
> /scratch4/NCEPDEV/meso/noscrub/Perry.Shafran/for_john.
>
> Thanks!
>
> Perry
>
> On Tue, Jan 2, 2018 at 11:26 AM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > Hi Perry,
> >
> > OK, I can take a closer look at this today. Since there are a lot
of
> > potential reasons for 0 VIS matched pairs, I think it'd be best
for me to
> > test with the exact PB2NC and Point-Stat config files that you're
using.
> >
> > Would you be able to either copy some data to theia for me to test
or
> post
> > it to our anonymous ftp site (
> > https://dtcenter.org/met/users/support/met_help.php#ftp)?
> >
> > I'll need your PB2NC config file, Point-Stat config file, the
PREPBUFR
> > input file for PB2NC, and the forecast GRIB file that you're
verifying.
> >
> > I'll replicate the commands your running and then figure out why
you're
> > getting 0 matched pairs for visibility.
> >
> > Thanks,
> > John
> >
> >
> >
> > On Tue, Jan 2, 2018 at 6:14 AM, perry.shafran at noaa.gov via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> > >
> > > Hi, John,
> > >
> > > I was away all week so I might still have some questions. If
you have
> an
> > > answer to the original question about the level for visibility,
that
> > would
> > > be great.
> > >
> > > Thanks!
> > >
> > > Perry
> > >
> > > On Thu, Dec 28, 2017 at 3:52 PM, John Halley Gotway via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > > Hi Perry,
> > > >
> > > > Not sure if you're in the office today. If so, do you still
have
> > > questions
> > > > about this issue?
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > > On Thu, Dec 21, 2017 at 12:29 PM, Minna Win via RT <
> met_help at ucar.edu>
> > > > wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> > > > >
> > > > > Hi Perry,
> > > > >
> > > > > To check if VIS is Z0, you could use plot_data_plane to see
if it
> is:
> > > > >
> > > > > plot_data_plane \
> > > > > <full path to input file> \
> > > > > <full path and name of temporary ps file> \
> > > > > 'name="VIS"; level="Z0";' -v 4
> > > > >
> > > > > Setting the verbosity to 4 provides enough information to
indicate
> > > > whether
> > > > > this variable is found, or if the level does/doesn't exist.
This
> has
> > > > been a
> > > > > helpful "trick" we
> > > > > use on another project which has variable names and levels
that are
> > not
> > > > > well-documented.
> > > > >
> > > > > It might be worth trying- it's up to you.
> > > > >
> > > > > Best of luck!
> > > > >
> > > > > Thanks,
> > > > > Minna
> > > > >
> > > > >
> > > > >
> > > > >
> > > > > ---------------
> > > > > Minna Win
> > > > > NCAR
> > > > > Research Applications Lab
> > > > > Phone: 303-497-8423
> > > > > Fax: 302-497-8401
> > > > >
> > > > >
> > > > > On Thu, Dec 21, 2017 at 7:09 PM, perry.shafran at noaa.gov via
RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> > > > > >
> > > > > > Hi, Minna,
> > > > > >
> > > > > > I'm not sure that really is an issue because I have other
> variables
> > > > (i.e.
> > > > > > TMP) with items not enclosed in brackets, and those appear
in the
> > > .stat
> > > > > > file. Part of the issue here I think has to do with the
fact
> that
> > > the
> > > > > > variable has one mnemonic for obs (HOVI) and another for
the
> model
> > > > (VIS).
> > > > > > I think I did it right the mapping of the BUFR name to the
GRIB
> > name,
> > > > but
> > > > > > I'm still not getting anything. I think VIS is considered
Z0 but
> > > maybe
> > > > > > it's something else? It didn't work when I set the level
to L0
> > > either.
> > > > > >
> > > > > > Waiting until after the new year for assistance is fine.
> > > > > >
> > > > > > Thanks!
> > > > > >
> > > > > > Perry
> > > > > >
> > > > > > On Thu, Dec 21, 2017 at 1:55 PM, Minna Win via RT <
> > met_help at ucar.edu
> > > >
> > > > > > wrote:
> > > > > >
> > > > > > > Hi Perry,
> > > > > > >
> > > > > > > A lot of us are out of the office for the holidays, so
you
> might
> > > > > > experience
> > > > > > > a delay in your response until everyone returns on
January 2.
> I
> > > > don't
> > > > > > have
> > > > > > > any expertise in point stat, but I think you need to
have your
> > > level
> > > > > > > enclosed in brackets
> > > > > > >
> > > > > > > ie level = ["Z0"];
> > > > > > >
> > > > > > >
> > > > > > > John or anyone else with more experience might be able
to help
> > you
> > > > out.
> > > > > > >
> > > > > > > Regards,
> > > > > > > Minna
> > > > > > >
> > > > > > > ---------------
> > > > > > > Minna Win
> > > > > > > NCAR
> > > > > > > Research Applications Lab
> > > > > > > Phone: 303-497-8423
> > > > > > > Fax: 302-497-8401
> > > > > > >
> > > > > > >
> > > > > > > On Thu, Dec 21, 2017 at 6:38 PM, perry.shafran at noaa.gov
via
> RT <
> > > > > > > met_help at ucar.edu> wrote:
> > > > > > >
> > > > > > > >
> > > > > > > > Thu Dec 21 11:38:09 2017: Request 83377 was acted
upon.
> > > > > > > > Transaction: Ticket created by perry.shafran at noaa.gov
> > > > > > > > Queue: met_help
> > > > > > > > Subject: visibility in MET
> > > > > > > > Owner: Nobody
> > > > > > > > Requestors: perry.shafran at noaa.gov
> > > > > > > > Status: new
> > > > > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > > Ticket/Display.html?id=83377
> > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > Hi, there,
> > > > > > > >
> > > > > > > > I am trying to verify visibility in MET. In my PB2NC
config
> > > file,
> > > > I
> > > > > > have
> > > > > > > > the following:
> > > > > > > >
> > > > > > > > //
> > > > > > > > // Observation values to retain or derive
> > > > > > > > //
> > > > > > > > obs_bufr_var = [ "POB", "TOB", "ZOB", "QOB", "UOB",
"VOB",
> > "HOVI"
> > > > ];
> > > > > > > >
> > > > > > > > //
> > > > > > > > // Mapping of BUFR variable name to GRIB name
> > > > > > > > //
> > > > > > > > obs_bufr_map = [
> > > > > > > > { key = "HOVI"; val = "VIS";
> > > > > > > > key = "POB"; val = "PRMSL";
> > > > > > > > key = "ZOB"; val = "HGT";
> > > > > > > > key = "TOB"; val = "TMP";
> > > > > > > > key = "QOB"; val = "SPFH";
> > > > > > > > key = "UOB"; val = "UGRD";
> > > > > > > > key = "VOB"; val = "VGRD"; }
> > > > > > > > ];
> > > > > > > >
> > > > > > > > Then I have this in my point_stat config file:
> > > > > > > >
> > > > > > > > fcst = {
> > > > > > > > field = [
> > > > > > > > {
> > > > > > > > name = "VIS";
> > > > > > > > level = "Z0";
> > > > > > > > cat_thresh = [ >=800, >=1600, >=8046, >=16090
];
> > > > > > > > },
> > > > > > > > .
> > > > > > > > .
> > > > > > > > .
> > > > > > > > ];
> > > > > > > >
> > > > > > > > }
> > > > > > > >
> > > > > > > > obs = fcst;
> > > > > > > >
> > > > > > > > I still can't get VIS verification in my .stat file.
Is
> there
> > > > > > anything I
> > > > > > > > am missing here, that you can see? Do you need more
> > information?
> > > > > > > >
> > > > > > > > Thanks!
> > > > > > > >
> > > > > > > > Perry
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: visibility in MET
From: perry.shafran at noaa.gov
Time: Tue Jan 02 14:23:45 2018
Hi, John,
Ah, I need to use the curly braces. But I did get verification points
for
some of the others though (such as TMP - TOB). Interesting. But the
curly
braces should help.
I am not sure whether, when I run point_stat, whether the level should
be
Z0 or L0. What level is it?
Thanks!
Perry
On Tue, Jan 2, 2018 at 4:14 PM, John Halley Gotway via RT
<met_help at ucar.edu
> wrote:
> Perry,
>
> Great, thanks for sending that to me. I took a look at the PB2NC
config
> file and immediately saw a problem. The "obs_bufr_map" is supposed
to be
> an array of dictionaries. Each dictionary is surrounded by curly
braces
> and defines a key and a value. So basically, the "obs_bufr_map"
entry
> isn't formatted correctly. Instead, it should look like this:
>
> obs_bufr_map = [
> { key = "HOVI"; val = "VIS"; },
> { key = "POB"; val = "PRMSL"; },
> { key = "ZOB"; val = "HGT"; },
> { key = "TOB"; val = "TMP"; },
> { key = "QOB"; val = "SPFH"; },
> { key = "UOB"; val = "UGRD"; },
> { key = "VOB"; val = "VGRD"; }
> ];
>
> To make sure this did the trick, I ran the following commands:
>
> # Load MET module
> module use /contrib/modulefiles
> module load met/6.1
>
> # Run PB2NC using your original settings
> pb2nc prepda.2018010100 prepda.2018010100_ORIG.nc PB2NCConfig
>
> # Run plot_point_obs to plot VIS observations (I expect there to be
0 of
> them)
> plot_point_obs prepda.2018010100_ORIG.nc prepda.2018010100_ORIG.ps
> -obs_var
> VIS
>
> # I log messages from plot_point_obs confirm that there are no
observations
> named VIS in the output:
> # DEBUG 2: Finished plotting 0 locations.
>
> # Edit the config file to define obs_bufr_map as listed above and
then
> rerun PB2NC
> pb2nc prepda.2018010100 prepda.2018010100_FIX.nc PB2NCConfig_FIX
>
> # Rerun plot_point_obs and make sure that there really are VIS
observations
> present
> plot_point_obs prepda.2018010100_FIX.nc prepda.2018010100_FIX.ps
-obs_var
> VIS
>
> # And the log message now confirm that we do have observations named
VIS:
> # DEBUG 2: Finished plotting 3657 locations.
>
> So that change to the PB2NC config file enabled PB2NC to output
> observations named VIS as 3657 locations. Hopefully this change
will fix
> your issue.
>
> Thanks,
> John
>
> Thanks,
> John
>
>
> On Tue, Jan 2, 2018 at 11:38 AM, perry.shafran at noaa.gov via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> >
> > Hi, John,
> >
> > Should have everything you need on theia in the directory
> > /scratch4/NCEPDEV/meso/noscrub/Perry.Shafran/for_john.
> >
> > Thanks!
> >
> > Perry
> >
> > On Tue, Jan 2, 2018 at 11:26 AM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> > > Hi Perry,
> > >
> > > OK, I can take a closer look at this today. Since there are a
lot of
> > > potential reasons for 0 VIS matched pairs, I think it'd be best
for me
> to
> > > test with the exact PB2NC and Point-Stat config files that
you're
> using.
> > >
> > > Would you be able to either copy some data to theia for me to
test or
> > post
> > > it to our anonymous ftp site (
> > > https://dtcenter.org/met/users/support/met_help.php#ftp)?
> > >
> > > I'll need your PB2NC config file, Point-Stat config file, the
PREPBUFR
> > > input file for PB2NC, and the forecast GRIB file that you're
verifying.
> > >
> > > I'll replicate the commands your running and then figure out why
you're
> > > getting 0 matched pairs for visibility.
> > >
> > > Thanks,
> > > John
> > >
> > >
> > >
> > > On Tue, Jan 2, 2018 at 6:14 AM, perry.shafran at noaa.gov via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377
>
> > > >
> > > > Hi, John,
> > > >
> > > > I was away all week so I might still have some questions. If
you
> have
> > an
> > > > answer to the original question about the level for
visibility, that
> > > would
> > > > be great.
> > > >
> > > > Thanks!
> > > >
> > > > Perry
> > > >
> > > > On Thu, Dec 28, 2017 at 3:52 PM, John Halley Gotway via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > > Hi Perry,
> > > > >
> > > > > Not sure if you're in the office today. If so, do you still
have
> > > > questions
> > > > > about this issue?
> > > > >
> > > > > Thanks,
> > > > > John
> > > > >
> > > > > On Thu, Dec 21, 2017 at 12:29 PM, Minna Win via RT <
> > met_help at ucar.edu>
> > > > > wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> > > > > >
> > > > > > Hi Perry,
> > > > > >
> > > > > > To check if VIS is Z0, you could use plot_data_plane to
see if it
> > is:
> > > > > >
> > > > > > plot_data_plane \
> > > > > > <full path to input file> \
> > > > > > <full path and name of temporary ps file> \
> > > > > > 'name="VIS"; level="Z0";' -v 4
> > > > > >
> > > > > > Setting the verbosity to 4 provides enough information to
> indicate
> > > > > whether
> > > > > > this variable is found, or if the level does/doesn't
exist. This
> > has
> > > > > been a
> > > > > > helpful "trick" we
> > > > > > use on another project which has variable names and levels
that
> are
> > > not
> > > > > > well-documented.
> > > > > >
> > > > > > It might be worth trying- it's up to you.
> > > > > >
> > > > > > Best of luck!
> > > > > >
> > > > > > Thanks,
> > > > > > Minna
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > > ---------------
> > > > > > Minna Win
> > > > > > NCAR
> > > > > > Research Applications Lab
> > > > > > Phone: 303-497-8423
> > > > > > Fax: 302-497-8401
> > > > > >
> > > > > >
> > > > > > On Thu, Dec 21, 2017 at 7:09 PM, perry.shafran at noaa.gov
via RT <
> > > > > > met_help at ucar.edu> wrote:
> > > > > >
> > > > > > >
> > > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377
> >
> > > > > > >
> > > > > > > Hi, Minna,
> > > > > > >
> > > > > > > I'm not sure that really is an issue because I have
other
> > variables
> > > > > (i.e.
> > > > > > > TMP) with items not enclosed in brackets, and those
appear in
> the
> > > > .stat
> > > > > > > file. Part of the issue here I think has to do with the
fact
> > that
> > > > the
> > > > > > > variable has one mnemonic for obs (HOVI) and another for
the
> > model
> > > > > (VIS).
> > > > > > > I think I did it right the mapping of the BUFR name to
the GRIB
> > > name,
> > > > > but
> > > > > > > I'm still not getting anything. I think VIS is
considered Z0
> but
> > > > maybe
> > > > > > > it's something else? It didn't work when I set the
level to L0
> > > > either.
> > > > > > >
> > > > > > > Waiting until after the new year for assistance is fine.
> > > > > > >
> > > > > > > Thanks!
> > > > > > >
> > > > > > > Perry
> > > > > > >
> > > > > > > On Thu, Dec 21, 2017 at 1:55 PM, Minna Win via RT <
> > > met_help at ucar.edu
> > > > >
> > > > > > > wrote:
> > > > > > >
> > > > > > > > Hi Perry,
> > > > > > > >
> > > > > > > > A lot of us are out of the office for the holidays, so
you
> > might
> > > > > > > experience
> > > > > > > > a delay in your response until everyone returns on
January 2.
> > I
> > > > > don't
> > > > > > > have
> > > > > > > > any expertise in point stat, but I think you need to
have
> your
> > > > level
> > > > > > > > enclosed in brackets
> > > > > > > >
> > > > > > > > ie level = ["Z0"];
> > > > > > > >
> > > > > > > >
> > > > > > > > John or anyone else with more experience might be able
to
> help
> > > you
> > > > > out.
> > > > > > > >
> > > > > > > > Regards,
> > > > > > > > Minna
> > > > > > > >
> > > > > > > > ---------------
> > > > > > > > Minna Win
> > > > > > > > NCAR
> > > > > > > > Research Applications Lab
> > > > > > > > Phone: 303-497-8423
> > > > > > > > Fax: 302-497-8401
> > > > > > > >
> > > > > > > >
> > > > > > > > On Thu, Dec 21, 2017 at 6:38 PM,
perry.shafran at noaa.gov via
> > RT <
> > > > > > > > met_help at ucar.edu> wrote:
> > > > > > > >
> > > > > > > > >
> > > > > > > > > Thu Dec 21 11:38:09 2017: Request 83377 was acted
upon.
> > > > > > > > > Transaction: Ticket created by
perry.shafran at noaa.gov
> > > > > > > > > Queue: met_help
> > > > > > > > > Subject: visibility in MET
> > > > > > > > > Owner: Nobody
> > > > > > > > > Requestors: perry.shafran at noaa.gov
> > > > > > > > > Status: new
> > > > > > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > > > Ticket/Display.html?id=83377
> > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > > Hi, there,
> > > > > > > > >
> > > > > > > > > I am trying to verify visibility in MET. In my
PB2NC
> config
> > > > file,
> > > > > I
> > > > > > > have
> > > > > > > > > the following:
> > > > > > > > >
> > > > > > > > > //
> > > > > > > > > // Observation values to retain or derive
> > > > > > > > > //
> > > > > > > > > obs_bufr_var = [ "POB", "TOB", "ZOB", "QOB", "UOB",
"VOB",
> > > "HOVI"
> > > > > ];
> > > > > > > > >
> > > > > > > > > //
> > > > > > > > > // Mapping of BUFR variable name to GRIB name
> > > > > > > > > //
> > > > > > > > > obs_bufr_map = [
> > > > > > > > > { key = "HOVI"; val = "VIS";
> > > > > > > > > key = "POB"; val = "PRMSL";
> > > > > > > > > key = "ZOB"; val = "HGT";
> > > > > > > > > key = "TOB"; val = "TMP";
> > > > > > > > > key = "QOB"; val = "SPFH";
> > > > > > > > > key = "UOB"; val = "UGRD";
> > > > > > > > > key = "VOB"; val = "VGRD"; }
> > > > > > > > > ];
> > > > > > > > >
> > > > > > > > > Then I have this in my point_stat config file:
> > > > > > > > >
> > > > > > > > > fcst = {
> > > > > > > > > field = [
> > > > > > > > > {
> > > > > > > > > name = "VIS";
> > > > > > > > > level = "Z0";
> > > > > > > > > cat_thresh = [ >=800, >=1600, >=8046, >=16090
];
> > > > > > > > > },
> > > > > > > > > .
> > > > > > > > > .
> > > > > > > > > .
> > > > > > > > > ];
> > > > > > > > >
> > > > > > > > > }
> > > > > > > > >
> > > > > > > > > obs = fcst;
> > > > > > > > >
> > > > > > > > > I still can't get VIS verification in my .stat file.
Is
> > there
> > > > > > > anything I
> > > > > > > > > am missing here, that you can see? Do you need more
> > > information?
> > > > > > > > >
> > > > > > > > > Thanks!
> > > > > > > > >
> > > > > > > > > Perry
> > > > > > > > >
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: visibility in MET
From: John Halley Gotway
Time: Tue Jan 02 14:43:29 2018
Perry,
When you run pb2nc, it first reads the *default* config file from
here:
/contrib/met/6.1/met-6.1/share/met/config/PB2NCConfig_default
That default config file contains the "obs_prepbufr_map" entry which
maps
TOB to TMP. So the only *new* thing you added was the mapping of HOVI
to
VIS. Since it wasn't formatted correctly, that's why your setting had
no
effect. But the other mappings were still defined from the default
config
file.
I ran point_stat using the config file you sent me and it produced
3518
matched pairs for VIS:
point_stat AWIP3D00.tm00 prepda.2018010100_FIX.nc
PointStatConfig_cv
-outdir out
DEBUG 2: Processing VIS/Z0 versus VIS/Z0, for observation type ADPSFC,
over
region G236, for interpolation method NEAREST(1), using 3518 pairs.
So it looks like 'Z0' works just fine.
FYI, I'm working in this directory:
/scratch4/BMC/dtc/John.H.Gotway/MET/MET_Help/shafran_data_20180102
Thanks,
John
On Tue, Jan 2, 2018 at 2:23 PM, perry.shafran at noaa.gov via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
>
> Hi, John,
>
> Ah, I need to use the curly braces. But I did get verification
points for
> some of the others though (such as TMP - TOB). Interesting. But
the curly
> braces should help.
>
> I am not sure whether, when I run point_stat, whether the level
should be
> Z0 or L0. What level is it?
>
> Thanks!
>
> Perry
>
> On Tue, Jan 2, 2018 at 4:14 PM, John Halley Gotway via RT <
> met_help at ucar.edu
> > wrote:
>
> > Perry,
> >
> > Great, thanks for sending that to me. I took a look at the PB2NC
config
> > file and immediately saw a problem. The "obs_bufr_map" is
supposed to be
> > an array of dictionaries. Each dictionary is surrounded by curly
braces
> > and defines a key and a value. So basically, the "obs_bufr_map"
entry
> > isn't formatted correctly. Instead, it should look like this:
> >
> > obs_bufr_map = [
> > { key = "HOVI"; val = "VIS"; },
> > { key = "POB"; val = "PRMSL"; },
> > { key = "ZOB"; val = "HGT"; },
> > { key = "TOB"; val = "TMP"; },
> > { key = "QOB"; val = "SPFH"; },
> > { key = "UOB"; val = "UGRD"; },
> > { key = "VOB"; val = "VGRD"; }
> > ];
> >
> > To make sure this did the trick, I ran the following commands:
> >
> > # Load MET module
> > module use /contrib/modulefiles
> > module load met/6.1
> >
> > # Run PB2NC using your original settings
> > pb2nc prepda.2018010100 prepda.2018010100_ORIG.nc PB2NCConfig
> >
> > # Run plot_point_obs to plot VIS observations (I expect there to
be 0 of
> > them)
> > plot_point_obs prepda.2018010100_ORIG.nc prepda.2018010100_ORIG.ps
> > -obs_var
> > VIS
> >
> > # I log messages from plot_point_obs confirm that there are no
> observations
> > named VIS in the output:
> > # DEBUG 2: Finished plotting 0 locations.
> >
> > # Edit the config file to define obs_bufr_map as listed above and
then
> > rerun PB2NC
> > pb2nc prepda.2018010100 prepda.2018010100_FIX.nc PB2NCConfig_FIX
> >
> > # Rerun plot_point_obs and make sure that there really are VIS
> observations
> > present
> > plot_point_obs prepda.2018010100_FIX.nc prepda.2018010100_FIX.ps
> -obs_var
> > VIS
> >
> > # And the log message now confirm that we do have observations
named VIS:
> > # DEBUG 2: Finished plotting 3657 locations.
> >
> > So that change to the PB2NC config file enabled PB2NC to output
> > observations named VIS as 3657 locations. Hopefully this change
will fix
> > your issue.
> >
> > Thanks,
> > John
> >
> > Thanks,
> > John
> >
> >
> > On Tue, Jan 2, 2018 at 11:38 AM, perry.shafran at noaa.gov via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> > >
> > > Hi, John,
> > >
> > > Should have everything you need on theia in the directory
> > > /scratch4/NCEPDEV/meso/noscrub/Perry.Shafran/for_john.
> > >
> > > Thanks!
> > >
> > > Perry
> > >
> > > On Tue, Jan 2, 2018 at 11:26 AM, John Halley Gotway via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > > Hi Perry,
> > > >
> > > > OK, I can take a closer look at this today. Since there are a
lot of
> > > > potential reasons for 0 VIS matched pairs, I think it'd be
best for
> me
> > to
> > > > test with the exact PB2NC and Point-Stat config files that
you're
> > using.
> > > >
> > > > Would you be able to either copy some data to theia for me to
test or
> > > post
> > > > it to our anonymous ftp site (
> > > > https://dtcenter.org/met/users/support/met_help.php#ftp)?
> > > >
> > > > I'll need your PB2NC config file, Point-Stat config file, the
> PREPBUFR
> > > > input file for PB2NC, and the forecast GRIB file that you're
> verifying.
> > > >
> > > > I'll replicate the commands your running and then figure out
why
> you're
> > > > getting 0 matched pairs for visibility.
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > >
> > > >
> > > > On Tue, Jan 2, 2018 at 6:14 AM, perry.shafran at noaa.gov via RT
<
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> > > > >
> > > > > Hi, John,
> > > > >
> > > > > I was away all week so I might still have some questions.
If you
> > have
> > > an
> > > > > answer to the original question about the level for
visibility,
> that
> > > > would
> > > > > be great.
> > > > >
> > > > > Thanks!
> > > > >
> > > > > Perry
> > > > >
> > > > > On Thu, Dec 28, 2017 at 3:52 PM, John Halley Gotway via RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > > Hi Perry,
> > > > > >
> > > > > > Not sure if you're in the office today. If so, do you
still have
> > > > > questions
> > > > > > about this issue?
> > > > > >
> > > > > > Thanks,
> > > > > > John
> > > > > >
> > > > > > On Thu, Dec 21, 2017 at 12:29 PM, Minna Win via RT <
> > > met_help at ucar.edu>
> > > > > > wrote:
> > > > > >
> > > > > > >
> > > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377
> >
> > > > > > >
> > > > > > > Hi Perry,
> > > > > > >
> > > > > > > To check if VIS is Z0, you could use plot_data_plane to
see if
> it
> > > is:
> > > > > > >
> > > > > > > plot_data_plane \
> > > > > > > <full path to input file> \
> > > > > > > <full path and name of temporary ps file> \
> > > > > > > 'name="VIS"; level="Z0";' -v 4
> > > > > > >
> > > > > > > Setting the verbosity to 4 provides enough information
to
> > indicate
> > > > > > whether
> > > > > > > this variable is found, or if the level does/doesn't
exist.
> This
> > > has
> > > > > > been a
> > > > > > > helpful "trick" we
> > > > > > > use on another project which has variable names and
levels that
> > are
> > > > not
> > > > > > > well-documented.
> > > > > > >
> > > > > > > It might be worth trying- it's up to you.
> > > > > > >
> > > > > > > Best of luck!
> > > > > > >
> > > > > > > Thanks,
> > > > > > > Minna
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > ---------------
> > > > > > > Minna Win
> > > > > > > NCAR
> > > > > > > Research Applications Lab
> > > > > > > Phone: 303-497-8423
> > > > > > > Fax: 302-497-8401
> > > > > > >
> > > > > > >
> > > > > > > On Thu, Dec 21, 2017 at 7:09 PM, perry.shafran at noaa.gov
via
> RT <
> > > > > > > met_help at ucar.edu> wrote:
> > > > > > >
> > > > > > > >
> > > > > > > > <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=83377
> > >
> > > > > > > >
> > > > > > > > Hi, Minna,
> > > > > > > >
> > > > > > > > I'm not sure that really is an issue because I have
other
> > > variables
> > > > > > (i.e.
> > > > > > > > TMP) with items not enclosed in brackets, and those
appear in
> > the
> > > > > .stat
> > > > > > > > file. Part of the issue here I think has to do with
the fact
> > > that
> > > > > the
> > > > > > > > variable has one mnemonic for obs (HOVI) and another
for the
> > > model
> > > > > > (VIS).
> > > > > > > > I think I did it right the mapping of the BUFR name to
the
> GRIB
> > > > name,
> > > > > > but
> > > > > > > > I'm still not getting anything. I think VIS is
considered Z0
> > but
> > > > > maybe
> > > > > > > > it's something else? It didn't work when I set the
level to
> L0
> > > > > either.
> > > > > > > >
> > > > > > > > Waiting until after the new year for assistance is
fine.
> > > > > > > >
> > > > > > > > Thanks!
> > > > > > > >
> > > > > > > > Perry
> > > > > > > >
> > > > > > > > On Thu, Dec 21, 2017 at 1:55 PM, Minna Win via RT <
> > > > met_help at ucar.edu
> > > > > >
> > > > > > > > wrote:
> > > > > > > >
> > > > > > > > > Hi Perry,
> > > > > > > > >
> > > > > > > > > A lot of us are out of the office for the holidays,
so you
> > > might
> > > > > > > > experience
> > > > > > > > > a delay in your response until everyone returns on
January
> 2.
> > > I
> > > > > > don't
> > > > > > > > have
> > > > > > > > > any expertise in point stat, but I think you need to
have
> > your
> > > > > level
> > > > > > > > > enclosed in brackets
> > > > > > > > >
> > > > > > > > > ie level = ["Z0"];
> > > > > > > > >
> > > > > > > > >
> > > > > > > > > John or anyone else with more experience might be
able to
> > help
> > > > you
> > > > > > out.
> > > > > > > > >
> > > > > > > > > Regards,
> > > > > > > > > Minna
> > > > > > > > >
> > > > > > > > > ---------------
> > > > > > > > > Minna Win
> > > > > > > > > NCAR
> > > > > > > > > Research Applications Lab
> > > > > > > > > Phone: 303-497-8423
> > > > > > > > > Fax: 302-497-8401
> > > > > > > > >
> > > > > > > > >
> > > > > > > > > On Thu, Dec 21, 2017 at 6:38 PM,
perry.shafran at noaa.gov
> via
> > > RT <
> > > > > > > > > met_help at ucar.edu> wrote:
> > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > Thu Dec 21 11:38:09 2017: Request 83377 was acted
upon.
> > > > > > > > > > Transaction: Ticket created by
perry.shafran at noaa.gov
> > > > > > > > > > Queue: met_help
> > > > > > > > > > Subject: visibility in MET
> > > > > > > > > > Owner: Nobody
> > > > > > > > > > Requestors: perry.shafran at noaa.gov
> > > > > > > > > > Status: new
> > > > > > > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > > > > Ticket/Display.html?id=83377
> > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > Hi, there,
> > > > > > > > > >
> > > > > > > > > > I am trying to verify visibility in MET. In my
PB2NC
> > config
> > > > > file,
> > > > > > I
> > > > > > > > have
> > > > > > > > > > the following:
> > > > > > > > > >
> > > > > > > > > > //
> > > > > > > > > > // Observation values to retain or derive
> > > > > > > > > > //
> > > > > > > > > > obs_bufr_var = [ "POB", "TOB", "ZOB", "QOB",
"UOB",
> "VOB",
> > > > "HOVI"
> > > > > > ];
> > > > > > > > > >
> > > > > > > > > > //
> > > > > > > > > > // Mapping of BUFR variable name to GRIB name
> > > > > > > > > > //
> > > > > > > > > > obs_bufr_map = [
> > > > > > > > > > { key = "HOVI"; val = "VIS";
> > > > > > > > > > key = "POB"; val = "PRMSL";
> > > > > > > > > > key = "ZOB"; val = "HGT";
> > > > > > > > > > key = "TOB"; val = "TMP";
> > > > > > > > > > key = "QOB"; val = "SPFH";
> > > > > > > > > > key = "UOB"; val = "UGRD";
> > > > > > > > > > key = "VOB"; val = "VGRD"; }
> > > > > > > > > > ];
> > > > > > > > > >
> > > > > > > > > > Then I have this in my point_stat config file:
> > > > > > > > > >
> > > > > > > > > > fcst = {
> > > > > > > > > > field = [
> > > > > > > > > > {
> > > > > > > > > > name = "VIS";
> > > > > > > > > > level = "Z0";
> > > > > > > > > > cat_thresh = [ >=800, >=1600, >=8046,
>=16090 ];
> > > > > > > > > > },
> > > > > > > > > > .
> > > > > > > > > > .
> > > > > > > > > > .
> > > > > > > > > > ];
> > > > > > > > > >
> > > > > > > > > > }
> > > > > > > > > >
> > > > > > > > > > obs = fcst;
> > > > > > > > > >
> > > > > > > > > > I still can't get VIS verification in my .stat
file. Is
> > > there
> > > > > > > > anything I
> > > > > > > > > > am missing here, that you can see? Do you need
more
> > > > information?
> > > > > > > > > >
> > > > > > > > > > Thanks!
> > > > > > > > > >
> > > > > > > > > > Perry
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: visibility in MET
From: perry.shafran at noaa.gov
Time: Wed Jan 03 07:01:07 2018
Hi, John,
I'll add in the formatting of the braces properly in the control file.
Thanks!
Perry
On Tue, Jan 2, 2018 at 4:43 PM, John Halley Gotway via RT
<met_help at ucar.edu
> wrote:
> Perry,
>
> When you run pb2nc, it first reads the *default* config file from
here:
> /contrib/met/6.1/met-6.1/share/met/config/PB2NCConfig_default
>
> That default config file contains the "obs_prepbufr_map" entry which
maps
> TOB to TMP. So the only *new* thing you added was the mapping of
HOVI to
> VIS. Since it wasn't formatted correctly, that's why your setting
had no
> effect. But the other mappings were still defined from the default
config
> file.
>
> I ran point_stat using the config file you sent me and it produced
3518
> matched pairs for VIS:
> point_stat AWIP3D00.tm00 prepda.2018010100_FIX.nc
PointStatConfig_cv
> -outdir out
>
> DEBUG 2: Processing VIS/Z0 versus VIS/Z0, for observation type
ADPSFC, over
> region G236, for interpolation method NEAREST(1), using 3518 pairs.
>
> So it looks like 'Z0' works just fine.
>
> FYI, I'm working in this directory:
> /scratch4/BMC/dtc/John.H.Gotway/MET/MET_Help/shafran_data_20180102
>
> Thanks,
> John
>
>
>
> On Tue, Jan 2, 2018 at 2:23 PM, perry.shafran at noaa.gov via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> >
> > Hi, John,
> >
> > Ah, I need to use the curly braces. But I did get verification
points
> for
> > some of the others though (such as TMP - TOB). Interesting. But
the
> curly
> > braces should help.
> >
> > I am not sure whether, when I run point_stat, whether the level
should be
> > Z0 or L0. What level is it?
> >
> > Thanks!
> >
> > Perry
> >
> > On Tue, Jan 2, 2018 at 4:14 PM, John Halley Gotway via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > > Perry,
> > >
> > > Great, thanks for sending that to me. I took a look at the
PB2NC
> config
> > > file and immediately saw a problem. The "obs_bufr_map" is
supposed to
> be
> > > an array of dictionaries. Each dictionary is surrounded by
curly
> braces
> > > and defines a key and a value. So basically, the "obs_bufr_map"
entry
> > > isn't formatted correctly. Instead, it should look like this:
> > >
> > > obs_bufr_map = [
> > > { key = "HOVI"; val = "VIS"; },
> > > { key = "POB"; val = "PRMSL"; },
> > > { key = "ZOB"; val = "HGT"; },
> > > { key = "TOB"; val = "TMP"; },
> > > { key = "QOB"; val = "SPFH"; },
> > > { key = "UOB"; val = "UGRD"; },
> > > { key = "VOB"; val = "VGRD"; }
> > > ];
> > >
> > > To make sure this did the trick, I ran the following commands:
> > >
> > > # Load MET module
> > > module use /contrib/modulefiles
> > > module load met/6.1
> > >
> > > # Run PB2NC using your original settings
> > > pb2nc prepda.2018010100 prepda.2018010100_ORIG.nc PB2NCConfig
> > >
> > > # Run plot_point_obs to plot VIS observations (I expect there to
be 0
> of
> > > them)
> > > plot_point_obs prepda.2018010100_ORIG.nc
prepda.2018010100_ORIG.ps
> > > -obs_var
> > > VIS
> > >
> > > # I log messages from plot_point_obs confirm that there are no
> > observations
> > > named VIS in the output:
> > > # DEBUG 2: Finished plotting 0 locations.
> > >
> > > # Edit the config file to define obs_bufr_map as listed above
and then
> > > rerun PB2NC
> > > pb2nc prepda.2018010100 prepda.2018010100_FIX.nc PB2NCConfig_FIX
> > >
> > > # Rerun plot_point_obs and make sure that there really are VIS
> > observations
> > > present
> > > plot_point_obs prepda.2018010100_FIX.nc prepda.2018010100_FIX.ps
> > -obs_var
> > > VIS
> > >
> > > # And the log message now confirm that we do have observations
named
> VIS:
> > > # DEBUG 2: Finished plotting 3657 locations.
> > >
> > > So that change to the PB2NC config file enabled PB2NC to output
> > > observations named VIS as 3657 locations. Hopefully this change
will
> fix
> > > your issue.
> > >
> > > Thanks,
> > > John
> > >
> > > Thanks,
> > > John
> > >
> > >
> > > On Tue, Jan 2, 2018 at 11:38 AM, perry.shafran at noaa.gov via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377
>
> > > >
> > > > Hi, John,
> > > >
> > > > Should have everything you need on theia in the directory
> > > > /scratch4/NCEPDEV/meso/noscrub/Perry.Shafran/for_john.
> > > >
> > > > Thanks!
> > > >
> > > > Perry
> > > >
> > > > On Tue, Jan 2, 2018 at 11:26 AM, John Halley Gotway via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > > Hi Perry,
> > > > >
> > > > > OK, I can take a closer look at this today. Since there are
a lot
> of
> > > > > potential reasons for 0 VIS matched pairs, I think it'd be
best for
> > me
> > > to
> > > > > test with the exact PB2NC and Point-Stat config files that
you're
> > > using.
> > > > >
> > > > > Would you be able to either copy some data to theia for me
to test
> or
> > > > post
> > > > > it to our anonymous ftp site (
> > > > > https://dtcenter.org/met/users/support/met_help.php#ftp)?
> > > > >
> > > > > I'll need your PB2NC config file, Point-Stat config file,
the
> > PREPBUFR
> > > > > input file for PB2NC, and the forecast GRIB file that you're
> > verifying.
> > > > >
> > > > > I'll replicate the commands your running and then figure out
why
> > you're
> > > > > getting 0 matched pairs for visibility.
> > > > >
> > > > > Thanks,
> > > > > John
> > > > >
> > > > >
> > > > >
> > > > > On Tue, Jan 2, 2018 at 6:14 AM, perry.shafran at noaa.gov via
RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83377 >
> > > > > >
> > > > > > Hi, John,
> > > > > >
> > > > > > I was away all week so I might still have some questions.
If you
> > > have
> > > > an
> > > > > > answer to the original question about the level for
visibility,
> > that
> > > > > would
> > > > > > be great.
> > > > > >
> > > > > > Thanks!
> > > > > >
> > > > > > Perry
> > > > > >
> > > > > > On Thu, Dec 28, 2017 at 3:52 PM, John Halley Gotway via RT
<
> > > > > > met_help at ucar.edu> wrote:
> > > > > >
> > > > > > > Hi Perry,
> > > > > > >
> > > > > > > Not sure if you're in the office today. If so, do you
still
> have
> > > > > > questions
> > > > > > > about this issue?
> > > > > > >
> > > > > > > Thanks,
> > > > > > > John
> > > > > > >
> > > > > > > On Thu, Dec 21, 2017 at 12:29 PM, Minna Win via RT <
> > > > met_help at ucar.edu>
> > > > > > > wrote:
> > > > > > >
> > > > > > > >
> > > > > > > > <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=83377
> > >
> > > > > > > >
> > > > > > > > Hi Perry,
> > > > > > > >
> > > > > > > > To check if VIS is Z0, you could use plot_data_plane
to see
> if
> > it
> > > > is:
> > > > > > > >
> > > > > > > > plot_data_plane \
> > > > > > > > <full path to input file> \
> > > > > > > > <full path and name of temporary ps file> \
> > > > > > > > 'name="VIS"; level="Z0";' -v 4
> > > > > > > >
> > > > > > > > Setting the verbosity to 4 provides enough information
to
> > > indicate
> > > > > > > whether
> > > > > > > > this variable is found, or if the level does/doesn't
exist.
> > This
> > > > has
> > > > > > > been a
> > > > > > > > helpful "trick" we
> > > > > > > > use on another project which has variable names and
levels
> that
> > > are
> > > > > not
> > > > > > > > well-documented.
> > > > > > > >
> > > > > > > > It might be worth trying- it's up to you.
> > > > > > > >
> > > > > > > > Best of luck!
> > > > > > > >
> > > > > > > > Thanks,
> > > > > > > > Minna
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > ---------------
> > > > > > > > Minna Win
> > > > > > > > NCAR
> > > > > > > > Research Applications Lab
> > > > > > > > Phone: 303-497-8423
> > > > > > > > Fax: 302-497-8401
> > > > > > > >
> > > > > > > >
> > > > > > > > On Thu, Dec 21, 2017 at 7:09 PM,
perry.shafran at noaa.gov via
> > RT <
> > > > > > > > met_help at ucar.edu> wrote:
> > > > > > > >
> > > > > > > > >
> > > > > > > > > <URL: https://rt.rap.ucar.edu/rt/
> > Ticket/Display.html?id=83377
> > > >
> > > > > > > > >
> > > > > > > > > Hi, Minna,
> > > > > > > > >
> > > > > > > > > I'm not sure that really is an issue because I have
other
> > > > variables
> > > > > > > (i.e.
> > > > > > > > > TMP) with items not enclosed in brackets, and those
appear
> in
> > > the
> > > > > > .stat
> > > > > > > > > file. Part of the issue here I think has to do with
the
> fact
> > > > that
> > > > > > the
> > > > > > > > > variable has one mnemonic for obs (HOVI) and another
for
> the
> > > > model
> > > > > > > (VIS).
> > > > > > > > > I think I did it right the mapping of the BUFR name
to the
> > GRIB
> > > > > name,
> > > > > > > but
> > > > > > > > > I'm still not getting anything. I think VIS is
considered
> Z0
> > > but
> > > > > > maybe
> > > > > > > > > it's something else? It didn't work when I set the
level
> to
> > L0
> > > > > > either.
> > > > > > > > >
> > > > > > > > > Waiting until after the new year for assistance is
fine.
> > > > > > > > >
> > > > > > > > > Thanks!
> > > > > > > > >
> > > > > > > > > Perry
> > > > > > > > >
> > > > > > > > > On Thu, Dec 21, 2017 at 1:55 PM, Minna Win via RT <
> > > > > met_help at ucar.edu
> > > > > > >
> > > > > > > > > wrote:
> > > > > > > > >
> > > > > > > > > > Hi Perry,
> > > > > > > > > >
> > > > > > > > > > A lot of us are out of the office for the
holidays, so
> you
> > > > might
> > > > > > > > > experience
> > > > > > > > > > a delay in your response until everyone returns on
> January
> > 2.
> > > > I
> > > > > > > don't
> > > > > > > > > have
> > > > > > > > > > any expertise in point stat, but I think you need
to have
> > > your
> > > > > > level
> > > > > > > > > > enclosed in brackets
> > > > > > > > > >
> > > > > > > > > > ie level = ["Z0"];
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > John or anyone else with more experience might be
able to
> > > help
> > > > > you
> > > > > > > out.
> > > > > > > > > >
> > > > > > > > > > Regards,
> > > > > > > > > > Minna
> > > > > > > > > >
> > > > > > > > > > ---------------
> > > > > > > > > > Minna Win
> > > > > > > > > > NCAR
> > > > > > > > > > Research Applications Lab
> > > > > > > > > > Phone: 303-497-8423
> > > > > > > > > > Fax: 302-497-8401
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > On Thu, Dec 21, 2017 at 6:38 PM,
perry.shafran at noaa.gov
> > via
> > > > RT <
> > > > > > > > > > met_help at ucar.edu> wrote:
> > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > > Thu Dec 21 11:38:09 2017: Request 83377 was
acted upon.
> > > > > > > > > > > Transaction: Ticket created by
perry.shafran at noaa.gov
> > > > > > > > > > > Queue: met_help
> > > > > > > > > > > Subject: visibility in MET
> > > > > > > > > > > Owner: Nobody
> > > > > > > > > > > Requestors: perry.shafran at noaa.gov
> > > > > > > > > > > Status: new
> > > > > > > > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > > > > > > Ticket/Display.html?id=83377
> > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > > Hi, there,
> > > > > > > > > > >
> > > > > > > > > > > I am trying to verify visibility in MET. In my
PB2NC
> > > config
> > > > > > file,
> > > > > > > I
> > > > > > > > > have
> > > > > > > > > > > the following:
> > > > > > > > > > >
> > > > > > > > > > > //
> > > > > > > > > > > // Observation values to retain or derive
> > > > > > > > > > > //
> > > > > > > > > > > obs_bufr_var = [ "POB", "TOB", "ZOB", "QOB",
"UOB",
> > "VOB",
> > > > > "HOVI"
> > > > > > > ];
> > > > > > > > > > >
> > > > > > > > > > > //
> > > > > > > > > > > // Mapping of BUFR variable name to GRIB name
> > > > > > > > > > > //
> > > > > > > > > > > obs_bufr_map = [
> > > > > > > > > > > { key = "HOVI"; val = "VIS";
> > > > > > > > > > > key = "POB"; val = "PRMSL";
> > > > > > > > > > > key = "ZOB"; val = "HGT";
> > > > > > > > > > > key = "TOB"; val = "TMP";
> > > > > > > > > > > key = "QOB"; val = "SPFH";
> > > > > > > > > > > key = "UOB"; val = "UGRD";
> > > > > > > > > > > key = "VOB"; val = "VGRD"; }
> > > > > > > > > > > ];
> > > > > > > > > > >
> > > > > > > > > > > Then I have this in my point_stat config file:
> > > > > > > > > > >
> > > > > > > > > > > fcst = {
> > > > > > > > > > > field = [
> > > > > > > > > > > {
> > > > > > > > > > > name = "VIS";
> > > > > > > > > > > level = "Z0";
> > > > > > > > > > > cat_thresh = [ >=800, >=1600, >=8046,
>=16090 ];
> > > > > > > > > > > },
> > > > > > > > > > > .
> > > > > > > > > > > .
> > > > > > > > > > > .
> > > > > > > > > > > ];
> > > > > > > > > > >
> > > > > > > > > > > }
> > > > > > > > > > >
> > > > > > > > > > > obs = fcst;
> > > > > > > > > > >
> > > > > > > > > > > I still can't get VIS verification in my .stat
file.
> Is
> > > > there
> > > > > > > > > anything I
> > > > > > > > > > > am missing here, that you can see? Do you need
more
> > > > > information?
> > > > > > > > > > >
> > > > > > > > > > > Thanks!
> > > > > > > > > > >
> > > > > > > > > > > Perry
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
More information about the Met_help
mailing list