[Met_help] [rt.rap.ucar.edu #84814] History for METV6.1: point_stat not producing SL1L2 stat lines for relative humidity (RH)
Julie Prestopnik via RT
met_help at ucar.edu
Tue Apr 17 14:54:02 MDT 2018
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello,
I am running METV6.1. I noticed that not SL1L2 lines were being produced
for relative humidity (RH). From the output, I see that MET is finding the
RH variables from the model forecast file, but no pairs are being found.
*DEBUG 2: Reading data for RH/P700.DEBUG 2: For RH/P700 found 1 forecast
levels, 0 climatology mean levels, and 0 climatology standard deviation
levels.DEBUG 2: Processing RH/P700 versus RH/P700, for observation type
ADPUPA, over region FULL, for interpolation method BILIN(4), using 0 pairs.*
I ran pb2nc prior to running point_stat with* { key = "D_RH"; val =
"RH"; }* in *obs_prefbufr_map*. In my point_stat config file, I have
*fcst = { field = [ { name = "RH"; level = [ "P1000", "P925",
"P850", "P700", "P500", "P400", "P300" ]; } ];}obs = fcst;*
When I do the same process with METV6.0, I get the SL1L2 lines for RH. I
know that pb2nc got a revamp in METV6.1 and this may be causing an issue?
Maybe the way I have my point_stat config file set up isn't correct with
the new way pb2nc saves information?
Thanks in advance for the help!
Mallory
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: METV6.1: point_stat not producing SL1L2 stat lines for relative humidity (RH)
From: Julie Prestopnik
Time: Tue Apr 17 10:51:22 2018
Hi Mallory.
Could you please tell me where your both your data and config files
are
located for both pb2nc and also for point_stat? Could you also please
send
me the command lines that you are using to run both pb2nc and
point_stat?
John is out this week. I will do my best to help out and try to
figure out
what is going on.
Thank you!
Julie
On Tue, Apr 17, 2018 at 9:31 AM, Mallory Row - NOAA Affiliate via RT <
met_help at ucar.edu> wrote:
>
> Tue Apr 17 09:31:49 2018: Request 84814 was acted upon.
> Transaction: Ticket created by mallory.row at noaa.gov
> Queue: met_help
> Subject: METV6.1: point_stat not producing SL1L2 stat lines for
> relative humidity (RH)
> Owner: Nobody
> Requestors: mallory.row at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84814 >
>
>
> Hello,
>
> I am running METV6.1. I noticed that not SL1L2 lines were being
produced
> for relative humidity (RH). From the output, I see that MET is
finding the
> RH variables from the model forecast file, but no pairs are being
found.
>
>
>
> *DEBUG 2: Reading data for RH/P700.DEBUG 2: For RH/P700 found 1
forecast
> levels, 0 climatology mean levels, and 0 climatology standard
deviation
> levels.DEBUG 2: Processing RH/P700 versus RH/P700, for observation
type
> ADPUPA, over region FULL, for interpolation method BILIN(4), using 0
> pairs.*
>
> I ran pb2nc prior to running point_stat with* { key = "D_RH"; val
=
> "RH"; }* in *obs_prefbufr_map*. In my point_stat config file, I
have
>
>
>
>
>
>
>
> *fcst = { field = [ { name = "RH"; level = [ "P1000",
"P925",
> "P850", "P700", "P500", "P400", "P300" ]; } ];}obs = fcst;*
>
> When I do the same process with METV6.0, I get the SL1L2 lines for
RH. I
> know that pb2nc got a revamp in METV6.1 and this may be causing an
issue?
> Maybe the way I have my point_stat config file set up isn't correct
with
> the new way pb2nc saves information?
>
> Thanks in advance for the help!
>
> Mallory
>
>
------------------------------------------------
Subject: METV6.1: point_stat not producing SL1L2 stat lines for relative humidity (RH)
From: Mallory Row - NOAA Affiliate
Time: Tue Apr 17 13:28:44 2018
Hi Julie,
I'm working on Theia in
/scratch4/NCEPDEV/global/save/Mallory.Row/VRFY/RH_test. Within that,
there
are two directories metv60 and metv61. The corresponding config files
are
in there along with the output I ran and commands (in commands.txt).
For METV6.0:
module use /contrib/modulefiles
module load met/6.0
pb2nc
/scratch4/NCEPDEV/global/noscrub/stat/prepbufr/gdas/prepbufr.gdas.2018041000
prepbufr.gdas.2018041000.nc PB2NCConfig_upper_air
point_stat
/scratch4/NCEPDEV/global/noscrub/stat/gfs/pgbf24.gfs.2018040900
prepbufr.gdas.2018041000.nc PointStatConfig_upper_air
For METV6.1:
module use /contrib/modulefiles
module load met/6.1
pb2nc
/scratch4/NCEPDEV/global/noscrub/stat/prepbufr/gdas/prepbufr.gdas.2018041000
prepbufr.gdas.2018041000.nc PB2NConfig
point_stat
/scratch4/NCEPDEV/global/noscrub/stat/gfs/pgbf24.gfs.2018040900
prepbufr.gdas.2018041000.nc PointStatConfig
I'll note that this was first noticed in the grid-to-obs use case from
MET+. To take MET+ out of the equation, I tested to MET commands
outside of
it.
Thanks!
Mallory
On Tue, Apr 17, 2018 at 12:51 PM, Julie Prestopnik via RT
<met_help at ucar.edu
> wrote:
> Hi Mallory.
>
> Could you please tell me where your both your data and config files
are
> located for both pb2nc and also for point_stat? Could you also
please send
> me the command lines that you are using to run both pb2nc and
point_stat?
> John is out this week. I will do my best to help out and try to
figure out
> what is going on.
>
> Thank you!
>
> Julie
>
> On Tue, Apr 17, 2018 at 9:31 AM, Mallory Row - NOAA Affiliate via RT
<
> met_help at ucar.edu> wrote:
>
> >
> > Tue Apr 17 09:31:49 2018: Request 84814 was acted upon.
> > Transaction: Ticket created by mallory.row at noaa.gov
> > Queue: met_help
> > Subject: METV6.1: point_stat not producing SL1L2 stat lines
for
> > relative humidity (RH)
> > Owner: Nobody
> > Requestors: mallory.row at noaa.gov
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84814 >
> >
> >
> > Hello,
> >
> > I am running METV6.1. I noticed that not SL1L2 lines were being
produced
> > for relative humidity (RH). From the output, I see that MET is
finding
> the
> > RH variables from the model forecast file, but no pairs are being
found.
> >
> >
> >
> > *DEBUG 2: Reading data for RH/P700.DEBUG 2: For RH/P700 found 1
forecast
> > levels, 0 climatology mean levels, and 0 climatology standard
deviation
> > levels.DEBUG 2: Processing RH/P700 versus RH/P700, for observation
type
> > ADPUPA, over region FULL, for interpolation method BILIN(4), using
0
> > pairs.*
> >
> > I ran pb2nc prior to running point_stat with* { key = "D_RH";
val =
> > "RH"; }* in *obs_prefbufr_map*. In my point_stat config file,
I have
> >
> >
> >
> >
> >
> >
> >
> > *fcst = { field = [ { name = "RH"; level = [ "P1000",
"P925",
> > "P850", "P700", "P500", "P400", "P300" ]; } ];}obs = fcst;*
> >
> > When I do the same process with METV6.0, I get the SL1L2 lines for
RH. I
> > know that pb2nc got a revamp in METV6.1 and this may be causing an
issue?
> > Maybe the way I have my point_stat config file set up isn't
correct with
> > the new way pb2nc saves information?
> >
> > Thanks in advance for the help!
> >
> > Mallory
> >
> >
>
>
------------------------------------------------
Subject: METV6.1: point_stat not producing SL1L2 stat lines for relative humidity (RH)
From: Julie Prestopnik
Time: Tue Apr 17 14:30:27 2018
Hi Mallory.
I copied over your config files and made some modifications.
I made your 6.1 config file for both pb2nc and for point_stat look
more
like your config file for 6.0. When I was doing that I noticed that
you
did not have RH in obs_bufr_var. obs_bufr_var is the list of BUFR
variable
names to retain or derive. I modified obs_bufr_var to include the
same
vars that were in your obs_grib_code for 6.0.
These are all of the things I changed to make your 6.1 pb2nc config
file
look more like your 6.0 file:
Changed *message_type = []; *to *message_type = ["ADPUPA"];*
Changed *grid = "";* to *grid = "G003";*
Changed* end = 511;* to *end = 255;*
Changed *obs_bufr_var = [ "QOB", "TOB", "ZOB", "UOB", "VOB" ]; * to
*obs_bufr_var = [ "TOB", "D_RH", "UOB", "VOB", "ZOB" ];*
I believe I am seeing the output that you wanted to see for RH, after
adding RH to obs_bufr_var. Please take a look and let me know if the
output in my directory looks more like what you would expect and if
that
resolves the problem you were having.
/scratch4/BMC/dtc/Julie.Prestopnik/met_help/row_20180417/6.1/point_stat_240000L_20180410_000000V.stat
I hope that helps!
Julie
On Tue, Apr 17, 2018 at 1:28 PM, Mallory Row - NOAA Affiliate via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84814 >
>
> Hi Julie,
>
> I'm working on Theia in
> /scratch4/NCEPDEV/global/save/Mallory.Row/VRFY/RH_test. Within that,
there
> are two directories metv60 and metv61. The corresponding config
files are
> in there along with the output I ran and commands (in commands.txt).
>
> For METV6.0:
> module use /contrib/modulefiles
> module load met/6.0
> pb2nc
> /scratch4/NCEPDEV/global/noscrub/stat/prepbufr/gdas/
> prepbufr.gdas.2018041000
> prepbufr.gdas.2018041000.nc PB2NCConfig_upper_air
> point_stat
/scratch4/NCEPDEV/global/noscrub/stat/gfs/pgbf24.gfs.2018040900
> prepbufr.gdas.2018041000.nc PointStatConfig_upper_air
>
> For METV6.1:
> module use /contrib/modulefiles
> module load met/6.1
> pb2nc
> /scratch4/NCEPDEV/global/noscrub/stat/prepbufr/gdas/
> prepbufr.gdas.2018041000
> prepbufr.gdas.2018041000.nc PB2NConfig
> point_stat
/scratch4/NCEPDEV/global/noscrub/stat/gfs/pgbf24.gfs.2018040900
> prepbufr.gdas.2018041000.nc PointStatConfig
>
> I'll note that this was first noticed in the grid-to-obs use case
from
> MET+. To take MET+ out of the equation, I tested to MET commands
outside of
> it.
>
> Thanks!
>
> Mallory
>
> On Tue, Apr 17, 2018 at 12:51 PM, Julie Prestopnik via RT <
> met_help at ucar.edu
> > wrote:
>
> > Hi Mallory.
> >
> > Could you please tell me where your both your data and config
files are
> > located for both pb2nc and also for point_stat? Could you also
please
> send
> > me the command lines that you are using to run both pb2nc and
point_stat?
> > John is out this week. I will do my best to help out and try to
figure
> out
> > what is going on.
> >
> > Thank you!
> >
> > Julie
> >
> > On Tue, Apr 17, 2018 at 9:31 AM, Mallory Row - NOAA Affiliate via
RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > Tue Apr 17 09:31:49 2018: Request 84814 was acted upon.
> > > Transaction: Ticket created by mallory.row at noaa.gov
> > > Queue: met_help
> > > Subject: METV6.1: point_stat not producing SL1L2 stat lines
for
> > > relative humidity (RH)
> > > Owner: Nobody
> > > Requestors: mallory.row at noaa.gov
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84814
> >
> > >
> > >
> > > Hello,
> > >
> > > I am running METV6.1. I noticed that not SL1L2 lines were being
> produced
> > > for relative humidity (RH). From the output, I see that MET is
finding
> > the
> > > RH variables from the model forecast file, but no pairs are
being
> found.
> > >
> > >
> > >
> > > *DEBUG 2: Reading data for RH/P700.DEBUG 2: For RH/P700 found 1
> forecast
> > > levels, 0 climatology mean levels, and 0 climatology standard
deviation
> > > levels.DEBUG 2: Processing RH/P700 versus RH/P700, for
observation type
> > > ADPUPA, over region FULL, for interpolation method BILIN(4),
using 0
> > > pairs.*
> > >
> > > I ran pb2nc prior to running point_stat with* { key = "D_RH";
val =
> > > "RH"; }* in *obs_prefbufr_map*. In my point_stat config
file, I
> have
> > >
> > >
> > >
> > >
> > >
> > >
> > >
> > > *fcst = { field = [ { name = "RH"; level = [ "P1000",
> "P925",
> > > "P850", "P700", "P500", "P400", "P300" ]; } ];}obs = fcst;*
> > >
> > > When I do the same process with METV6.0, I get the SL1L2 lines
for RH.
> I
> > > know that pb2nc got a revamp in METV6.1 and this may be causing
an
> issue?
> > > Maybe the way I have my point_stat config file set up isn't
correct
> with
> > > the new way pb2nc saves information?
> > >
> > > Thanks in advance for the help!
> > >
> > > Mallory
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: METV6.1: point_stat not producing SL1L2 stat lines for relative humidity (RH)
From: Mallory Row - NOAA Affiliate
Time: Tue Apr 17 14:32:56 2018
OHHHH, I see now! Okay, this makes sense! Thank you!
Mallory
On Tue, Apr 17, 2018 at 4:30 PM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Mallory.
>
> I copied over your config files and made some modifications.
>
> I made your 6.1 config file for both pb2nc and for point_stat look
more
> like your config file for 6.0. When I was doing that I noticed that
you
> did not have RH in obs_bufr_var. obs_bufr_var is the list of BUFR
variable
> names to retain or derive. I modified obs_bufr_var to include the
same
> vars that were in your obs_grib_code for 6.0.
>
> These are all of the things I changed to make your 6.1 pb2nc config
file
> look more like your 6.0 file:
>
> Changed *message_type = []; *to *message_type = ["ADPUPA"];*
> Changed *grid = "";* to *grid = "G003";*
> Changed* end = 511;* to *end = 255;*
> Changed *obs_bufr_var = [ "QOB", "TOB", "ZOB", "UOB", "VOB" ]; * to
> *obs_bufr_var = [ "TOB", "D_RH", "UOB", "VOB", "ZOB" ];*
> I believe I am seeing the output that you wanted to see for RH,
after
> adding RH to obs_bufr_var. Please take a look and let me know if
the
> output in my directory looks more like what you would expect and if
that
> resolves the problem you were having.
>
>
/scratch4/BMC/dtc/Julie.Prestopnik/met_help/row_20180417/6.1/point_stat_
> 240000L_20180410_000000V.stat
>
> I hope that helps!
>
> Julie
>
> On Tue, Apr 17, 2018 at 1:28 PM, Mallory Row - NOAA Affiliate via RT
<
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84814 >
> >
> > Hi Julie,
> >
> > I'm working on Theia in
> > /scratch4/NCEPDEV/global/save/Mallory.Row/VRFY/RH_test. Within
that,
> there
> > are two directories metv60 and metv61. The corresponding config
files are
> > in there along with the output I ran and commands (in
commands.txt).
> >
> > For METV6.0:
> > module use /contrib/modulefiles
> > module load met/6.0
> > pb2nc
> > /scratch4/NCEPDEV/global/noscrub/stat/prepbufr/gdas/
> > prepbufr.gdas.2018041000
> > prepbufr.gdas.2018041000.nc PB2NCConfig_upper_air
> > point_stat /scratch4/NCEPDEV/global/noscrub/stat/gfs/pgbf24.gfs.
> 2018040900
> > prepbufr.gdas.2018041000.nc PointStatConfig_upper_air
> >
> > For METV6.1:
> > module use /contrib/modulefiles
> > module load met/6.1
> > pb2nc
> > /scratch4/NCEPDEV/global/noscrub/stat/prepbufr/gdas/
> > prepbufr.gdas.2018041000
> > prepbufr.gdas.2018041000.nc PB2NConfig
> > point_stat /scratch4/NCEPDEV/global/noscrub/stat/gfs/pgbf24.gfs.
> 2018040900
> > prepbufr.gdas.2018041000.nc PointStatConfig
> >
> > I'll note that this was first noticed in the grid-to-obs use case
from
> > MET+. To take MET+ out of the equation, I tested to MET commands
outside
> of
> > it.
> >
> > Thanks!
> >
> > Mallory
> >
> > On Tue, Apr 17, 2018 at 12:51 PM, Julie Prestopnik via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > > Hi Mallory.
> > >
> > > Could you please tell me where your both your data and config
files are
> > > located for both pb2nc and also for point_stat? Could you also
please
> > send
> > > me the command lines that you are using to run both pb2nc and
> point_stat?
> > > John is out this week. I will do my best to help out and try to
figure
> > out
> > > what is going on.
> > >
> > > Thank you!
> > >
> > > Julie
> > >
> > > On Tue, Apr 17, 2018 at 9:31 AM, Mallory Row - NOAA Affiliate
via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > Tue Apr 17 09:31:49 2018: Request 84814 was acted upon.
> > > > Transaction: Ticket created by mallory.row at noaa.gov
> > > > Queue: met_help
> > > > Subject: METV6.1: point_stat not producing SL1L2 stat
lines for
> > > > relative humidity (RH)
> > > > Owner: Nobody
> > > > Requestors: mallory.row at noaa.gov
> > > > Status: new
> > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=84814
> > >
> > > >
> > > >
> > > > Hello,
> > > >
> > > > I am running METV6.1. I noticed that not SL1L2 lines were
being
> > produced
> > > > for relative humidity (RH). From the output, I see that MET is
> finding
> > > the
> > > > RH variables from the model forecast file, but no pairs are
being
> > found.
> > > >
> > > >
> > > >
> > > > *DEBUG 2: Reading data for RH/P700.DEBUG 2: For RH/P700 found
1
> > forecast
> > > > levels, 0 climatology mean levels, and 0 climatology standard
> deviation
> > > > levels.DEBUG 2: Processing RH/P700 versus RH/P700, for
observation
> type
> > > > ADPUPA, over region FULL, for interpolation method BILIN(4),
using 0
> > > > pairs.*
> > > >
> > > > I ran pb2nc prior to running point_stat with* { key = "D_RH";
val
> =
> > > > "RH"; }* in *obs_prefbufr_map*. In my point_stat config
file, I
> > have
> > > >
> > > >
> > > >
> > > >
> > > >
> > > >
> > > >
> > > > *fcst = { field = [ { name = "RH"; level = [
"P1000",
> > "P925",
> > > > "P850", "P700", "P500", "P400", "P300" ]; } ];}obs = fcst;*
> > > >
> > > > When I do the same process with METV6.0, I get the SL1L2 lines
for
> RH.
> > I
> > > > know that pb2nc got a revamp in METV6.1 and this may be
causing an
> > issue?
> > > > Maybe the way I have my point_stat config file set up isn't
correct
> > with
> > > > the new way pb2nc saves information?
> > > >
> > > > Thanks in advance for the help!
> > > >
> > > > Mallory
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: METV6.1: point_stat not producing SL1L2 stat lines for relative humidity (RH)
From: Minna Win
Time: Tue Apr 17 14:36:00 2018
Thanks Julie- I can't believe I missed it, too!
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax: 302-497-8401
On Tue, Apr 17, 2018 at 2:32 PM, Mallory Row - NOAA Affiliate via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84814 >
>
> OHHHH, I see now! Okay, this makes sense! Thank you!
>
> Mallory
>
> On Tue, Apr 17, 2018 at 4:30 PM, Julie Prestopnik via RT <
> met_help at ucar.edu>
> wrote:
>
> > Hi Mallory.
> >
> > I copied over your config files and made some modifications.
> >
> > I made your 6.1 config file for both pb2nc and for point_stat look
more
> > like your config file for 6.0. When I was doing that I noticed
that you
> > did not have RH in obs_bufr_var. obs_bufr_var is the list of BUFR
> variable
> > names to retain or derive. I modified obs_bufr_var to include the
same
> > vars that were in your obs_grib_code for 6.0.
> >
> > These are all of the things I changed to make your 6.1 pb2nc
config file
> > look more like your 6.0 file:
> >
> > Changed *message_type = []; *to *message_type = ["ADPUPA"];*
> > Changed *grid = "";* to *grid = "G003";*
> > Changed* end = 511;* to *end = 255;*
> > Changed *obs_bufr_var = [ "QOB", "TOB", "ZOB", "UOB", "VOB" ]; *
to
> > *obs_bufr_var = [ "TOB", "D_RH", "UOB", "VOB", "ZOB" ];*
> > I believe I am seeing the output that you wanted to see for RH,
after
> > adding RH to obs_bufr_var. Please take a look and let me know if
the
> > output in my directory looks more like what you would expect and
if that
> > resolves the problem you were having.
> >
> >
/scratch4/BMC/dtc/Julie.Prestopnik/met_help/row_20180417/6.1/point_stat_
> > 240000L_20180410_000000V.stat
> >
> > I hope that helps!
> >
> > Julie
> >
> > On Tue, Apr 17, 2018 at 1:28 PM, Mallory Row - NOAA Affiliate via
RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84814 >
> > >
> > > Hi Julie,
> > >
> > > I'm working on Theia in
> > > /scratch4/NCEPDEV/global/save/Mallory.Row/VRFY/RH_test. Within
that,
> > there
> > > are two directories metv60 and metv61. The corresponding config
files
> are
> > > in there along with the output I ran and commands (in
commands.txt).
> > >
> > > For METV6.0:
> > > module use /contrib/modulefiles
> > > module load met/6.0
> > > pb2nc
> > > /scratch4/NCEPDEV/global/noscrub/stat/prepbufr/gdas/
> > > prepbufr.gdas.2018041000
> > > prepbufr.gdas.2018041000.nc PB2NCConfig_upper_air
> > > point_stat /scratch4/NCEPDEV/global/noscrub/stat/gfs/pgbf24.gfs.
> > 2018040900
> > > prepbufr.gdas.2018041000.nc PointStatConfig_upper_air
> > >
> > > For METV6.1:
> > > module use /contrib/modulefiles
> > > module load met/6.1
> > > pb2nc
> > > /scratch4/NCEPDEV/global/noscrub/stat/prepbufr/gdas/
> > > prepbufr.gdas.2018041000
> > > prepbufr.gdas.2018041000.nc PB2NConfig
> > > point_stat /scratch4/NCEPDEV/global/noscrub/stat/gfs/pgbf24.gfs.
> > 2018040900
> > > prepbufr.gdas.2018041000.nc PointStatConfig
> > >
> > > I'll note that this was first noticed in the grid-to-obs use
case from
> > > MET+. To take MET+ out of the equation, I tested to MET commands
> outside
> > of
> > > it.
> > >
> > > Thanks!
> > >
> > > Mallory
> > >
> > > On Tue, Apr 17, 2018 at 12:51 PM, Julie Prestopnik via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > > Hi Mallory.
> > > >
> > > > Could you please tell me where your both your data and config
files
> are
> > > > located for both pb2nc and also for point_stat? Could you
also
> please
> > > send
> > > > me the command lines that you are using to run both pb2nc and
> > point_stat?
> > > > John is out this week. I will do my best to help out and try
to
> figure
> > > out
> > > > what is going on.
> > > >
> > > > Thank you!
> > > >
> > > > Julie
> > > >
> > > > On Tue, Apr 17, 2018 at 9:31 AM, Mallory Row - NOAA Affiliate
via RT
> <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > Tue Apr 17 09:31:49 2018: Request 84814 was acted upon.
> > > > > Transaction: Ticket created by mallory.row at noaa.gov
> > > > > Queue: met_help
> > > > > Subject: METV6.1: point_stat not producing SL1L2 stat
lines
> for
> > > > > relative humidity (RH)
> > > > > Owner: Nobody
> > > > > Requestors: mallory.row at noaa.gov
> > > > > Status: new
> > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > Ticket/Display.html?id=84814
> > > >
> > > > >
> > > > >
> > > > > Hello,
> > > > >
> > > > > I am running METV6.1. I noticed that not SL1L2 lines were
being
> > > produced
> > > > > for relative humidity (RH). From the output, I see that MET
is
> > finding
> > > > the
> > > > > RH variables from the model forecast file, but no pairs are
being
> > > found.
> > > > >
> > > > >
> > > > >
> > > > > *DEBUG 2: Reading data for RH/P700.DEBUG 2: For RH/P700
found 1
> > > forecast
> > > > > levels, 0 climatology mean levels, and 0 climatology
standard
> > deviation
> > > > > levels.DEBUG 2: Processing RH/P700 versus RH/P700, for
observation
> > type
> > > > > ADPUPA, over region FULL, for interpolation method BILIN(4),
using
> 0
> > > > > pairs.*
> > > > >
> > > > > I ran pb2nc prior to running point_stat with* { key =
"D_RH";
> val
> > =
> > > > > "RH"; }* in *obs_prefbufr_map*. In my point_stat config
file, I
> > > have
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > > *fcst = { field = [ { name = "RH"; level = [
"P1000",
> > > "P925",
> > > > > "P850", "P700", "P500", "P400", "P300" ]; } ];}obs =
fcst;*
> > > > >
> > > > > When I do the same process with METV6.0, I get the SL1L2
lines for
> > RH.
> > > I
> > > > > know that pb2nc got a revamp in METV6.1 and this may be
causing an
> > > issue?
> > > > > Maybe the way I have my point_stat config file set up isn't
correct
> > > with
> > > > > the new way pb2nc saves information?
> > > > >
> > > > > Thanks in advance for the help!
> > > > >
> > > > > Mallory
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: METV6.1: point_stat not producing SL1L2 stat lines for relative humidity (RH)
From: Julie Prestopnik
Time: Tue Apr 17 14:52:46 2018
Great! Thank you! I'll go ahead and resolve this ticket now. Please
let
us know if you encounter any other issues.
Thanks,
Julie
On Tue, Apr 17, 2018 at 2:32 PM, Mallory Row - NOAA Affiliate via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84814 >
>
> OHHHH, I see now! Okay, this makes sense! Thank you!
>
> Mallory
>
> On Tue, Apr 17, 2018 at 4:30 PM, Julie Prestopnik via RT <
> met_help at ucar.edu>
> wrote:
>
> > Hi Mallory.
> >
> > I copied over your config files and made some modifications.
> >
> > I made your 6.1 config file for both pb2nc and for point_stat look
more
> > like your config file for 6.0. When I was doing that I noticed
that you
> > did not have RH in obs_bufr_var. obs_bufr_var is the list of BUFR
> variable
> > names to retain or derive. I modified obs_bufr_var to include the
same
> > vars that were in your obs_grib_code for 6.0.
> >
> > These are all of the things I changed to make your 6.1 pb2nc
config file
> > look more like your 6.0 file:
> >
> > Changed *message_type = []; *to *message_type = ["ADPUPA"];*
> > Changed *grid = "";* to *grid = "G003";*
> > Changed* end = 511;* to *end = 255;*
> > Changed *obs_bufr_var = [ "QOB", "TOB", "ZOB", "UOB", "VOB" ]; *
to
> > *obs_bufr_var = [ "TOB", "D_RH", "UOB", "VOB", "ZOB" ];*
> > I believe I am seeing the output that you wanted to see for RH,
after
> > adding RH to obs_bufr_var. Please take a look and let me know if
the
> > output in my directory looks more like what you would expect and
if that
> > resolves the problem you were having.
> >
> >
/scratch4/BMC/dtc/Julie.Prestopnik/met_help/row_20180417/6.1/point_stat_
> > 240000L_20180410_000000V.stat
> >
> > I hope that helps!
> >
> > Julie
> >
> > On Tue, Apr 17, 2018 at 1:28 PM, Mallory Row - NOAA Affiliate via
RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84814 >
> > >
> > > Hi Julie,
> > >
> > > I'm working on Theia in
> > > /scratch4/NCEPDEV/global/save/Mallory.Row/VRFY/RH_test. Within
that,
> > there
> > > are two directories metv60 and metv61. The corresponding config
files
> are
> > > in there along with the output I ran and commands (in
commands.txt).
> > >
> > > For METV6.0:
> > > module use /contrib/modulefiles
> > > module load met/6.0
> > > pb2nc
> > > /scratch4/NCEPDEV/global/noscrub/stat/prepbufr/gdas/
> > > prepbufr.gdas.2018041000
> > > prepbufr.gdas.2018041000.nc PB2NCConfig_upper_air
> > > point_stat /scratch4/NCEPDEV/global/noscrub/stat/gfs/pgbf24.gfs.
> > 2018040900
> > > prepbufr.gdas.2018041000.nc PointStatConfig_upper_air
> > >
> > > For METV6.1:
> > > module use /contrib/modulefiles
> > > module load met/6.1
> > > pb2nc
> > > /scratch4/NCEPDEV/global/noscrub/stat/prepbufr/gdas/
> > > prepbufr.gdas.2018041000
> > > prepbufr.gdas.2018041000.nc PB2NConfig
> > > point_stat /scratch4/NCEPDEV/global/noscrub/stat/gfs/pgbf24.gfs.
> > 2018040900
> > > prepbufr.gdas.2018041000.nc PointStatConfig
> > >
> > > I'll note that this was first noticed in the grid-to-obs use
case from
> > > MET+. To take MET+ out of the equation, I tested to MET commands
> outside
> > of
> > > it.
> > >
> > > Thanks!
> > >
> > > Mallory
> > >
> > > On Tue, Apr 17, 2018 at 12:51 PM, Julie Prestopnik via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > > Hi Mallory.
> > > >
> > > > Could you please tell me where your both your data and config
files
> are
> > > > located for both pb2nc and also for point_stat? Could you
also
> please
> > > send
> > > > me the command lines that you are using to run both pb2nc and
> > point_stat?
> > > > John is out this week. I will do my best to help out and try
to
> figure
> > > out
> > > > what is going on.
> > > >
> > > > Thank you!
> > > >
> > > > Julie
> > > >
> > > > On Tue, Apr 17, 2018 at 9:31 AM, Mallory Row - NOAA Affiliate
via RT
> <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > Tue Apr 17 09:31:49 2018: Request 84814 was acted upon.
> > > > > Transaction: Ticket created by mallory.row at noaa.gov
> > > > > Queue: met_help
> > > > > Subject: METV6.1: point_stat not producing SL1L2 stat
lines
> for
> > > > > relative humidity (RH)
> > > > > Owner: Nobody
> > > > > Requestors: mallory.row at noaa.gov
> > > > > Status: new
> > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > Ticket/Display.html?id=84814
> > > >
> > > > >
> > > > >
> > > > > Hello,
> > > > >
> > > > > I am running METV6.1. I noticed that not SL1L2 lines were
being
> > > produced
> > > > > for relative humidity (RH). From the output, I see that MET
is
> > finding
> > > > the
> > > > > RH variables from the model forecast file, but no pairs are
being
> > > found.
> > > > >
> > > > >
> > > > >
> > > > > *DEBUG 2: Reading data for RH/P700.DEBUG 2: For RH/P700
found 1
> > > forecast
> > > > > levels, 0 climatology mean levels, and 0 climatology
standard
> > deviation
> > > > > levels.DEBUG 2: Processing RH/P700 versus RH/P700, for
observation
> > type
> > > > > ADPUPA, over region FULL, for interpolation method BILIN(4),
using
> 0
> > > > > pairs.*
> > > > >
> > > > > I ran pb2nc prior to running point_stat with* { key =
"D_RH";
> val
> > =
> > > > > "RH"; }* in *obs_prefbufr_map*. In my point_stat config
file, I
> > > have
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > > *fcst = { field = [ { name = "RH"; level = [
"P1000",
> > > "P925",
> > > > > "P850", "P700", "P500", "P400", "P300" ]; } ];}obs =
fcst;*
> > > > >
> > > > > When I do the same process with METV6.0, I get the SL1L2
lines for
> > RH.
> > > I
> > > > > know that pb2nc got a revamp in METV6.1 and this may be
causing an
> > > issue?
> > > > > Maybe the way I have my point_stat config file set up isn't
correct
> > > with
> > > > > the new way pb2nc saves information?
> > > > >
> > > > > Thanks in advance for the help!
> > > > >
> > > > > Mallory
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
More information about the Met_help
mailing list