[Met_help] [rt.rap.ucar.edu #100423] History for Running Point-stat for precipitation verification
Julie Prestopnik via RT
met_help at ucar.edu
Tue Jul 6 11:19:54 MDT 2021
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Greetings,
My goal is to verify a GFS forecast of accumulated precipitation against point observations for an event. This is my first time using MET, so I followed the Point-stat tutorial successfully (http://dtcenter.org/met-online-tutorial-metv8-0/point-stat). However, when I try to run the same example with my own data I get no result (e.g. output stat files containing headers only). In addition, I have not being able to run MET with the mask produced by gen_vx_mask, I get an error saying that the file does not exit (while it does).
I have attached the files I am using to run my custom case. The directories have the same structure as the example in the tutorial.
By the way, I am using the docker version of MET.
Best regards,
Raúl Valenzuela
Assistant Professor
Institute of Engineering Sciences
University of O’Higgins
Rancagua, Chile
Tel +56 2332 86019
raul.valenzuela at uoh.cl
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Running Point-stat for precipitation verification
From: Julie Prestopnik
Time: Fri Jul 02 14:43:54 2021
Hi Raúl.
I see you are having trouble running Point-Stat with your data and
trouble
with running MET with the mask produced by gen_vx_mask. Thank you
for
sending your data. Can you please send me the command line calls that
you
are using so that I can be sure to run the same thing you are? Also,
can
you please tell me what version of MET you are running?
Thanks!
Julie
On Fri, Jul 2, 2021 at 1:57 PM Raúl Valenzuela Ruiz via RT <
met_help at ucar.edu> wrote:
>
> Fri Jul 02 13:57:01 2021: Request 100423 was acted upon.
> Transaction: Ticket created by raul.valenzuela at uoh.cl
> Queue: met_help
> Subject: Running Point-stat for precipitation verification
> Owner: Nobody
> Requestors: raul.valenzuela at uoh.cl
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423 >
>
>
> Greetings,
>
> My goal is to verify a GFS forecast of accumulated precipitation
against
> point observations for an event. This is my first time using MET, so
I
> followed the Point-stat tutorial successfully (
> http://dtcenter.org/met-online-tutorial-metv8-0/point-stat).
However,
> when I try to run the same example with my own data I get no result
(e.g.
> output stat files containing headers only). In addition, I have not
being
> able to run MET with the mask produced by gen_vx_mask, I get an
error
> saying that the file does not exit (while it does).
>
> I have attached the files I am using to run my custom case. The
> directories have the same structure as the example in the tutorial.
>
> By the way, I am using the docker version of MET.
>
> Best regards,
>
> Raúl Valenzuela
> Assistant Professor
> Institute of Engineering Sciences
> University of O’Higgins
> Rancagua, Chile
> Tel +56 2332 86019
> raul.valenzuela at uoh.cl
>
>
--
Julie Prestopnik (she/her)
Software Engineer
National Center for Atmospheric Research
Research Applications Laboratory
Email: jpresto at ucar.edu
My working day may not be your working day. Please do not feel
obliged to
reply to this email outside of your normal working hours.
------------------------------------------------
Subject: Running Point-stat for precipitation verification
From: Raúl Valenzuela Ruiz
Time: Fri Jul 02 15:04:39 2021
Thanks Julie!
I am running the following command:
point_stat \
$MET_RUN/input/fcst/gfs.0p25.2021012912.f036.grib2 \
$MET_RUN/output/ascii2nc/Acumulado_Precipitacion_2021-02-
01_23_00_Int_6.nc \
$MET_RUN/config/PointStatConfig_run01 \
-outdir $MET_RUN/output/point_stat -v 2
I thought that the problema was related to the masking area, so
replaced the mask grid with the following configuration:
mask = {
grid = [];
poly = [];
sid = [];
llpnt = [ { name = "BOX";
lat_thresh = >=-36&&<=-32;
lon_thresh = >=108&&<=111;
}
];
}
I realized that the GFS GRIB2 file had 0-360 longitude coordinates. I
tried to change them to -180 to 180 with CDO but I got some weird
results (longitudes 180 to 560). Then I read that GRIB2 only accepts
positive coordinates, so I kept the 0-360 longitude in the GFS GRIB2
file and shifted the longitude coordinates in my observations. No
success so far.
Best,
Raúl
Raúl Valenzuela
Assistant Professor
Institute of Engineering Sciences
University of O’Higgins
Rancagua, Chile
Tel +56 2332 86019
raul.valenzuela at uoh.cl
On Jul 2, 2021, 4:43 PM -0400, Julie Prestopnik via RT
<met_help at ucar.edu>, wrote:
> Hi Raúl.
>
> I see you are having trouble running Point-Stat with your data and
trouble
> with running MET with the mask produced by gen_vx_mask. Thank you
for
> sending your data. Can you please send me the command line calls
that you
> are using so that I can be sure to run the same thing you are? Also,
can
> you please tell me what version of MET you are running?
>
> Thanks!
>
> Julie
>
> On Fri, Jul 2, 2021 at 1:57 PM Raúl Valenzuela Ruiz via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Fri Jul 02 13:57:01 2021: Request 100423 was acted upon.
> > Transaction: Ticket created by raul.valenzuela at uoh.cl
> > Queue: met_help
> > Subject: Running Point-stat for precipitation verification
> > Owner: Nobody
> > Requestors: raul.valenzuela at uoh.cl
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423 >
> >
> >
> > Greetings,
> >
> > My goal is to verify a GFS forecast of accumulated precipitation
against
> > point observations for an event. This is my first time using MET,
so I
> > followed the Point-stat tutorial successfully (
> > http://dtcenter.org/met-online-tutorial-metv8-0/point-stat).
However,
> > when I try to run the same example with my own data I get no
result (e.g.
> > output stat files containing headers only). In addition, I have
not being
> > able to run MET with the mask produced by gen_vx_mask, I get an
error
> > saying that the file does not exit (while it does).
> >
> > I have attached the files I am using to run my custom case. The
> > directories have the same structure as the example in the
tutorial.
> >
> > By the way, I am using the docker version of MET.
> >
> > Best regards,
> >
> > Raúl Valenzuela
> > Assistant Professor
> > Institute of Engineering Sciences
> > University of O’Higgins
> > Rancagua, Chile
> > Tel +56 2332 86019
> > raul.valenzuela at uoh.cl
> >
> >
>
> --
> Julie Prestopnik (she/her)
> Software Engineer
> National Center for Atmospheric Research
> Research Applications Laboratory
> Email: jpresto at ucar.edu
>
> My working day may not be your working day. Please do not feel
obliged to
> reply to this email outside of your normal working hours.
>
------------------------------------------------
Subject: Running Point-stat for precipitation verification
From: Raúl Valenzuela Ruiz
Time: Fri Jul 02 15:08:55 2021
Sorry, I forgot the MET version, how should I verify that?
Raúl
Raúl Valenzuela
Assistant Professor
Institute of Engineering Sciences
University of O’Higgins
Rancagua, Chile
Tel +56 2332 86019
raul.valenzuela at uoh.cl
On Jul 2, 2021, 5:04 PM -0400, Raúl Valenzuela Ruiz
<raul.valenzuela at uoh.cl>, wrote:
> Thanks Julie!
>
> I am running the following command:
>
> point_stat \
> $MET_RUN/input/fcst/gfs.0p25.2021012912.f036.grib2 \
> $MET_RUN/output/ascii2nc/Acumulado_Precipitacion_2021-02-
01_23_00_Int_6.nc \
> $MET_RUN/config/PointStatConfig_run01 \
> -outdir $MET_RUN/output/point_stat -v 2
>
> I thought that the problema was related to the masking area, so
replaced the mask grid with the following configuration:
>
> mask = {
> grid = [];
> poly = [];
> sid = [];
> llpnt = [ { name = "BOX";
> lat_thresh = >=-36&&<=-32;
> lon_thresh = >=108&&<=111;
> }
> ];
> }
>
> I realized that the GFS GRIB2 file had 0-360 longitude coordinates.
I tried to change them to -180 to 180 with CDO but I got some weird
results (longitudes 180 to 560). Then I read that GRIB2 only accepts
positive coordinates, so I kept the 0-360 longitude in the GFS GRIB2
file and shifted the longitude coordinates in my observations. No
success so far.
>
> Best,
>
> Raúl
>
> Raúl Valenzuela
> Assistant Professor
> Institute of Engineering Sciences
> University of O’Higgins
> Rancagua, Chile
> Tel +56 2332 86019
> raul.valenzuela at uoh.cl
> On Jul 2, 2021, 4:43 PM -0400, Julie Prestopnik via RT
<met_help at ucar.edu>, wrote:
> > Hi Raúl.
> >
> > I see you are having trouble running Point-Stat with your data and
trouble
> > with running MET with the mask produced by gen_vx_mask. Thank you
for
> > sending your data. Can you please send me the command line calls
that you
> > are using so that I can be sure to run the same thing you are?
Also, can
> > you please tell me what version of MET you are running?
> >
> > Thanks!
> >
> > Julie
> >
> > On Fri, Jul 2, 2021 at 1:57 PM Raúl Valenzuela Ruiz via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > Fri Jul 02 13:57:01 2021: Request 100423 was acted upon.
> > > Transaction: Ticket created by raul.valenzuela at uoh.cl
> > > Queue: met_help
> > > Subject: Running Point-stat for precipitation verification
> > > Owner: Nobody
> > > Requestors: raul.valenzuela at uoh.cl
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423 >
> > >
> > >
> > > Greetings,
> > >
> > > My goal is to verify a GFS forecast of accumulated precipitation
against
> > > point observations for an event. This is my first time using
MET, so I
> > > followed the Point-stat tutorial successfully (
> > > http://dtcenter.org/met-online-tutorial-metv8-0/point-stat).
However,
> > > when I try to run the same example with my own data I get no
result (e.g.
> > > output stat files containing headers only). In addition, I have
not being
> > > able to run MET with the mask produced by gen_vx_mask, I get an
error
> > > saying that the file does not exit (while it does).
> > >
> > > I have attached the files I am using to run my custom case. The
> > > directories have the same structure as the example in the
tutorial.
> > >
> > > By the way, I am using the docker version of MET.
> > >
> > > Best regards,
> > >
> > > Raúl Valenzuela
> > > Assistant Professor
> > > Institute of Engineering Sciences
> > > University of O’Higgins
> > > Rancagua, Chile
> > > Tel +56 2332 86019
> > > raul.valenzuela at uoh.cl
> > >
> > >
> >
> > --
> > Julie Prestopnik (she/her)
> > Software Engineer
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Email: jpresto at ucar.edu
> >
> > My working day may not be your working day. Please do not feel
obliged to
> > reply to this email outside of your normal working hours.
> >
------------------------------------------------
Subject: Running Point-stat for precipitation verification
From: Julie Prestopnik
Time: Fri Jul 02 15:42:33 2021
Thank you for the information. You can run "--version" with any of
the MET
tools to see what version you are running. For example:
point_stat --version
Julie
On Fri, Jul 2, 2021 at 3:08 PM Raúl Valenzuela Ruiz via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423 >
>
> Sorry, I forgot the MET version, how should I verify that?
>
> Raúl
>
> Raúl Valenzuela
> Assistant Professor
> Institute of Engineering Sciences
> University of O’Higgins
> Rancagua, Chile
> Tel +56 2332 86019
> raul.valenzuela at uoh.cl
> On Jul 2, 2021, 5:04 PM -0400, Raúl Valenzuela Ruiz <
> raul.valenzuela at uoh.cl>, wrote:
> > Thanks Julie!
> >
> > I am running the following command:
> >
> > point_stat \
> > $MET_RUN/input/fcst/gfs.0p25.2021012912.f036.grib2 \
> >
> $MET_RUN/output/ascii2nc/Acumulado_Precipitacion_2021-02-
01_23_00_Int_6.nc \
> > $MET_RUN/config/PointStatConfig_run01 \
> > -outdir $MET_RUN/output/point_stat -v 2
> >
> > I thought that the problema was related to the masking area, so
replaced
> the mask grid with the following configuration:
> >
> > mask = {
> > grid = [];
> > poly = [];
> > sid = [];
> > llpnt = [ { name = "BOX";
> > lat_thresh = >=-36&&<=-32;
> > lon_thresh = >=108&&<=111;
> > }
> > ];
> > }
> >
> > I realized that the GFS GRIB2 file had 0-360 longitude
coordinates. I
> tried to change them to -180 to 180 with CDO but I got some weird
results
> (longitudes 180 to 560). Then I read that GRIB2 only accepts
positive
> coordinates, so I kept the 0-360 longitude in the GFS GRIB2 file and
> shifted the longitude coordinates in my observations. No success so
far.
> >
> > Best,
> >
> > Raúl
> >
> > Raúl Valenzuela
> > Assistant Professor
> > Institute of Engineering Sciences
> > University of O’Higgins
> > Rancagua, Chile
> > Tel +56 2332 86019
> > raul.valenzuela at uoh.cl
> > On Jul 2, 2021, 4:43 PM -0400, Julie Prestopnik via RT <
> met_help at ucar.edu>, wrote:
> > > Hi Raúl.
> > >
> > > I see you are having trouble running Point-Stat with your data
and
> trouble
> > > with running MET with the mask produced by gen_vx_mask. Thank
you for
> > > sending your data. Can you please send me the command line calls
that
> you
> > > are using so that I can be sure to run the same thing you are?
Also,
> can
> > > you please tell me what version of MET you are running?
> > >
> > > Thanks!
> > >
> > > Julie
> > >
> > > On Fri, Jul 2, 2021 at 1:57 PM Raúl Valenzuela Ruiz via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > Fri Jul 02 13:57:01 2021: Request 100423 was acted upon.
> > > > Transaction: Ticket created by raul.valenzuela at uoh.cl
> > > > Queue: met_help
> > > > Subject: Running Point-stat for precipitation verification
> > > > Owner: Nobody
> > > > Requestors: raul.valenzuela at uoh.cl
> > > > Status: new
> > > > Ticket <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423 >
> > > >
> > > >
> > > > Greetings,
> > > >
> > > > My goal is to verify a GFS forecast of accumulated
precipitation
> against
> > > > point observations for an event. This is my first time using
MET, so
> I
> > > > followed the Point-stat tutorial successfully (
> > > > http://dtcenter.org/met-online-tutorial-metv8-0/point-stat).
> However,
> > > > when I try to run the same example with my own data I get no
result
> (e.g.
> > > > output stat files containing headers only). In addition, I
have not
> being
> > > > able to run MET with the mask produced by gen_vx_mask, I get
an error
> > > > saying that the file does not exit (while it does).
> > > >
> > > > I have attached the files I am using to run my custom case.
The
> > > > directories have the same structure as the example in the
tutorial.
> > > >
> > > > By the way, I am using the docker version of MET.
> > > >
> > > > Best regards,
> > > >
> > > > Raúl Valenzuela
> > > > Assistant Professor
> > > > Institute of Engineering Sciences
> > > > University of O’Higgins
> > > > Rancagua, Chile
> > > > Tel +56 2332 86019
> > > > raul.valenzuela at uoh.cl
> > > >
> > > >
> > >
> > > --
> > > Julie Prestopnik (she/her)
> > > Software Engineer
> > > National Center for Atmospheric Research
> > > Research Applications Laboratory
> > > Email: jpresto at ucar.edu
> > >
> > > My working day may not be your working day. Please do not feel
obliged
> to
> > > reply to this email outside of your normal working hours.
> > >
>
>
--
Julie Prestopnik (she/her)
Software Engineer
National Center for Atmospheric Research
Research Applications Laboratory
Email: jpresto at ucar.edu
My working day may not be your working day. Please do not feel
obliged to
reply to this email outside of your normal working hours.
------------------------------------------------
Subject: Running Point-stat for precipitation verification
From: Raúl Valenzuela Ruiz
Time: Fri Jul 02 15:47:29 2021
Here is the output:
[root at 7c3def7bb5f7 met]# point_stat --version
MET Version: V10.0.0
Repository: https://github.com/dtcenter/MET
Revision: development
Change Date: development
Raúl
Raúl Valenzuela
Assistant Professor
Institute of Engineering Sciences
University of O’Higgins
Rancagua, Chile
Tel +56 2332 86019
raul.valenzuela at uoh.cl
On Jul 2, 2021, 5:42 PM -0400, Julie Prestopnik via RT
<met_help at ucar.edu>, wrote:
> Thank you for the information. You can run "--version" with any of
the MET
> tools to see what version you are running. For example:
>
> point_stat --version
>
> Julie
>
> On Fri, Jul 2, 2021 at 3:08 PM Raúl Valenzuela Ruiz via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423 >
> >
> > Sorry, I forgot the MET version, how should I verify that?
> >
> > Raúl
> >
> > Raúl Valenzuela
> > Assistant Professor
> > Institute of Engineering Sciences
> > University of O’Higgins
> > Rancagua, Chile
> > Tel +56 2332 86019
> > raul.valenzuela at uoh.cl
> > On Jul 2, 2021, 5:04 PM -0400, Raúl Valenzuela Ruiz <
> > raul.valenzuela at uoh.cl>, wrote:
> > > Thanks Julie!
> > >
> > > I am running the following command:
> > >
> > > point_stat \
> > > $MET_RUN/input/fcst/gfs.0p25.2021012912.f036.grib2 \
> > >
> > $MET_RUN/output/ascii2nc/Acumulado_Precipitacion_2021-02-
01_23_00_Int_6.nc \
> > > $MET_RUN/config/PointStatConfig_run01 \
> > > -outdir $MET_RUN/output/point_stat -v 2
> > >
> > > I thought that the problema was related to the masking area, so
replaced
> > the mask grid with the following configuration:
> > >
> > > mask = {
> > > grid = [];
> > > poly = [];
> > > sid = [];
> > > llpnt = [ { name = "BOX";
> > > lat_thresh = >=-36&&<=-32;
> > > lon_thresh = >=108&&<=111;
> > > }
> > > ];
> > > }
> > >
> > > I realized that the GFS GRIB2 file had 0-360 longitude
coordinates. I
> > tried to change them to -180 to 180 with CDO but I got some weird
results
> > (longitudes 180 to 560). Then I read that GRIB2 only accepts
positive
> > coordinates, so I kept the 0-360 longitude in the GFS GRIB2 file
and
> > shifted the longitude coordinates in my observations. No success
so far.
> > >
> > > Best,
> > >
> > > Raúl
> > >
> > > Raúl Valenzuela
> > > Assistant Professor
> > > Institute of Engineering Sciences
> > > University of O’Higgins
> > > Rancagua, Chile
> > > Tel +56 2332 86019
> > > raul.valenzuela at uoh.cl
> > > On Jul 2, 2021, 4:43 PM -0400, Julie Prestopnik via RT <
> > met_help at ucar.edu>, wrote:
> > > > Hi Raúl.
> > > >
> > > > I see you are having trouble running Point-Stat with your data
and
> > trouble
> > > > with running MET with the mask produced by gen_vx_mask. Thank
you for
> > > > sending your data. Can you please send me the command line
calls that
> > you
> > > > are using so that I can be sure to run the same thing you are?
Also,
> > can
> > > > you please tell me what version of MET you are running?
> > > >
> > > > Thanks!
> > > >
> > > > Julie
> > > >
> > > > On Fri, Jul 2, 2021 at 1:57 PM Raúl Valenzuela Ruiz via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > Fri Jul 02 13:57:01 2021: Request 100423 was acted upon.
> > > > > Transaction: Ticket created by raul.valenzuela at uoh.cl
> > > > > Queue: met_help
> > > > > Subject: Running Point-stat for precipitation verification
> > > > > Owner: Nobody
> > > > > Requestors: raul.valenzuela at uoh.cl
> > > > > Status: new
> > > > > Ticket <URL:
> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423 >
> > > > >
> > > > >
> > > > > Greetings,
> > > > >
> > > > > My goal is to verify a GFS forecast of accumulated
precipitation
> > against
> > > > > point observations for an event. This is my first time using
MET, so
> > I
> > > > > followed the Point-stat tutorial successfully (
> > > > > http://dtcenter.org/met-online-tutorial-metv8-0/point-stat).
> > However,
> > > > > when I try to run the same example with my own data I get no
result
> > (e.g.
> > > > > output stat files containing headers only). In addition, I
have not
> > being
> > > > > able to run MET with the mask produced by gen_vx_mask, I get
an error
> > > > > saying that the file does not exit (while it does).
> > > > >
> > > > > I have attached the files I am using to run my custom case.
The
> > > > > directories have the same structure as the example in the
tutorial.
> > > > >
> > > > > By the way, I am using the docker version of MET.
> > > > >
> > > > > Best regards,
> > > > >
> > > > > Raúl Valenzuela
> > > > > Assistant Professor
> > > > > Institute of Engineering Sciences
> > > > > University of O’Higgins
> > > > > Rancagua, Chile
> > > > > Tel +56 2332 86019
> > > > > raul.valenzuela at uoh.cl
> > > > >
> > > > >
> > > >
> > > > --
> > > > Julie Prestopnik (she/her)
> > > > Software Engineer
> > > > National Center for Atmospheric Research
> > > > Research Applications Laboratory
> > > > Email: jpresto at ucar.edu
> > > >
> > > > My working day may not be your working day. Please do not feel
obliged
> > to
> > > > reply to this email outside of your normal working hours.
> > > >
> >
> >
>
> --
> Julie Prestopnik (she/her)
> Software Engineer
> National Center for Atmospheric Research
> Research Applications Laboratory
> Email: jpresto at ucar.edu
>
> My working day may not be your working day. Please do not feel
obliged to
> reply to this email outside of your normal working hours.
>
------------------------------------------------
Subject: Running Point-stat for precipitation verification
From: Julie Prestopnik
Time: Fri Jul 02 15:54:15 2021
Running the files with a verbosity level of 2, I see the following
output,
which indicates that 1040 observations were processed, but all 1040
observations were rejected due to the valid time.
DEBUG 2: Processing APCP/A6 versus APCP/A6, for observation type
ADPSFC,
> over region BOX, for interpolation method NEAREST(1), using 0
matched pairs.
> DEBUG 2: Number of matched pairs = 0
> DEBUG 2: Observations processed = 1040
> DEBUG 2: Rejected: station id = 0
> DEBUG 2: Rejected: obs var name = 0
> DEBUG 2: Rejected: valid time = 1040
> DEBUG 2: Rejected: bad obs value = 0
> DEBUG 2: Rejected: off the grid = 0
> DEBUG 2: Rejected: topography = 0
> DEBUG 2: Rejected: level mismatch = 0
> DEBUG 2: Rejected: quality marker = 0
> DEBUG 2: Rejected: message type = 0
> DEBUG 2: Rejected: masking region = 0
> DEBUG 2: Rejected: bad fcst value = 0
> DEBUG 2: Rejected: bad climo mean = 0
> DEBUG 2: Rejected: bad climo stdev = 0
> DEBUG 2: Rejected: mpr filter = 0
> DEBUG 2: Rejected: duplicates = 0
In your fcst GRIB2 file I see the following valid time:
> > wgrib2 -vt gfs.0p25.2021012912.f036.grib2
> 1:0:vt=2021013100
In your obs NetCDF file I see the following valid time:
> hdr_vld_table =
> "20210201_230000" ;
I hope that helps!
Julie
On Fri, Jul 2, 2021 at 3:42 PM Julie Prestopnik <jpresto at ucar.edu>
wrote:
> Thank you for the information. You can run "--version" with any of
the
> MET tools to see what version you are running. For example:
>
> point_stat --version
>
> Julie
>
> On Fri, Jul 2, 2021 at 3:08 PM Raúl Valenzuela Ruiz via RT <
> met_help at ucar.edu> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423 >
>>
>> Sorry, I forgot the MET version, how should I verify that?
>>
>> Raúl
>>
>> Raúl Valenzuela
>> Assistant Professor
>> Institute of Engineering Sciences
>> University of O’Higgins
>> Rancagua, Chile
>> Tel +56 2332 86019
>> raul.valenzuela at uoh.cl
>> On Jul 2, 2021, 5:04 PM -0400, Raúl Valenzuela Ruiz <
>> raul.valenzuela at uoh.cl>, wrote:
>> > Thanks Julie!
>> >
>> > I am running the following command:
>> >
>> > point_stat \
>> > $MET_RUN/input/fcst/gfs.0p25.2021012912.f036.grib2 \
>> >
>> $MET_RUN/output/ascii2nc/Acumulado_Precipitacion_2021-02-
01_23_00_Int_6.nc \
>> > $MET_RUN/config/PointStatConfig_run01 \
>> > -outdir $MET_RUN/output/point_stat -v 2
>> >
>> > I thought that the problema was related to the masking area, so
>> replaced the mask grid with the following configuration:
>> >
>> > mask = {
>> > grid = [];
>> > poly = [];
>> > sid = [];
>> > llpnt = [ { name = "BOX";
>> > lat_thresh = >=-36&&<=-32;
>> > lon_thresh = >=108&&<=111;
>> > }
>> > ];
>> > }
>> >
>> > I realized that the GFS GRIB2 file had 0-360 longitude
coordinates. I
>> tried to change them to -180 to 180 with CDO but I got some weird
results
>> (longitudes 180 to 560). Then I read that GRIB2 only accepts
positive
>> coordinates, so I kept the 0-360 longitude in the GFS GRIB2 file
and
>> shifted the longitude coordinates in my observations. No success so
far.
>> >
>> > Best,
>> >
>> > Raúl
>> >
>> > Raúl Valenzuela
>> > Assistant Professor
>> > Institute of Engineering Sciences
>> > University of O’Higgins
>> > Rancagua, Chile
>> > Tel +56 2332 86019
>> > raul.valenzuela at uoh.cl
>> > On Jul 2, 2021, 4:43 PM -0400, Julie Prestopnik via RT <
>> met_help at ucar.edu>, wrote:
>> > > Hi Raúl.
>> > >
>> > > I see you are having trouble running Point-Stat with your data
and
>> trouble
>> > > with running MET with the mask produced by gen_vx_mask. Thank
you for
>> > > sending your data. Can you please send me the command line
calls that
>> you
>> > > are using so that I can be sure to run the same thing you are?
Also,
>> can
>> > > you please tell me what version of MET you are running?
>> > >
>> > > Thanks!
>> > >
>> > > Julie
>> > >
>> > > On Fri, Jul 2, 2021 at 1:57 PM Raúl Valenzuela Ruiz via RT <
>> > > met_help at ucar.edu> wrote:
>> > >
>> > > >
>> > > > Fri Jul 02 13:57:01 2021: Request 100423 was acted upon.
>> > > > Transaction: Ticket created by raul.valenzuela at uoh.cl
>> > > > Queue: met_help
>> > > > Subject: Running Point-stat for precipitation verification
>> > > > Owner: Nobody
>> > > > Requestors: raul.valenzuela at uoh.cl
>> > > > Status: new
>> > > > Ticket <URL:
>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423 >
>> > > >
>> > > >
>> > > > Greetings,
>> > > >
>> > > > My goal is to verify a GFS forecast of accumulated
precipitation
>> against
>> > > > point observations for an event. This is my first time using
MET,
>> so I
>> > > > followed the Point-stat tutorial successfully (
>> > > > http://dtcenter.org/met-online-tutorial-metv8-0/point-stat).
>> However,
>> > > > when I try to run the same example with my own data I get no
result
>> (e.g.
>> > > > output stat files containing headers only). In addition, I
have not
>> being
>> > > > able to run MET with the mask produced by gen_vx_mask, I get
an
>> error
>> > > > saying that the file does not exit (while it does).
>> > > >
>> > > > I have attached the files I am using to run my custom case.
The
>> > > > directories have the same structure as the example in the
tutorial.
>> > > >
>> > > > By the way, I am using the docker version of MET.
>> > > >
>> > > > Best regards,
>> > > >
>> > > > Raúl Valenzuela
>> > > > Assistant Professor
>> > > > Institute of Engineering Sciences
>> > > > University of O’Higgins
>> > > > Rancagua, Chile
>> > > > Tel +56 2332 86019
>> > > > raul.valenzuela at uoh.cl
>> > > >
>> > > >
>> > >
>> > > --
>> > > Julie Prestopnik (she/her)
>> > > Software Engineer
>> > > National Center for Atmospheric Research
>> > > Research Applications Laboratory
>> > > Email: jpresto at ucar.edu
>> > >
>> > > My working day may not be your working day. Please do not feel
>> obliged to
>> > > reply to this email outside of your normal working hours.
>> > >
>>
>>
>
> --
> Julie Prestopnik (she/her)
> Software Engineer
> National Center for Atmospheric Research
> Research Applications Laboratory
> Email: jpresto at ucar.edu
>
> My working day may not be your working day. Please do not feel
obliged to
> reply to this email outside of your normal working hours.
>
--
Julie Prestopnik (she/her)
Software Engineer
National Center for Atmospheric Research
Research Applications Laboratory
Email: jpresto at ucar.edu
My working day may not be your working day. Please do not feel
obliged to
reply to this email outside of your normal working hours.
------------------------------------------------
Subject: Running Point-stat for precipitation verification
From: Raúl Valenzuela Ruiz
Time: Fri Jul 02 16:28:53 2021
I knew it was a silly error. Thanks!
Raúl
Raúl Valenzuela
Assistant Professor
Institute of Engineering Sciences
University of O’Higgins
Rancagua, Chile
Tel +56 2332 86019
raul.valenzuela at uoh.cl
On Jul 2, 2021, 5:54 PM -0400, Julie Prestopnik via RT
<met_help at ucar.edu>, wrote:
> Running the files with a verbosity level of 2, I see the following
output,
> which indicates that 1040 observations were processed, but all 1040
> observations were rejected due to the valid time.
>
> DEBUG 2: Processing APCP/A6 versus APCP/A6, for observation type
ADPSFC,
> > over region BOX, for interpolation method NEAREST(1), using 0
matched pairs.
> > DEBUG 2: Number of matched pairs = 0
> > DEBUG 2: Observations processed = 1040
> > DEBUG 2: Rejected: station id = 0
> > DEBUG 2: Rejected: obs var name = 0
> > DEBUG 2: Rejected: valid time = 1040
> > DEBUG 2: Rejected: bad obs value = 0
> > DEBUG 2: Rejected: off the grid = 0
> > DEBUG 2: Rejected: topography = 0
> > DEBUG 2: Rejected: level mismatch = 0
> > DEBUG 2: Rejected: quality marker = 0
> > DEBUG 2: Rejected: message type = 0
> > DEBUG 2: Rejected: masking region = 0
> > DEBUG 2: Rejected: bad fcst value = 0
> > DEBUG 2: Rejected: bad climo mean = 0
> > DEBUG 2: Rejected: bad climo stdev = 0
> > DEBUG 2: Rejected: mpr filter = 0
> > DEBUG 2: Rejected: duplicates = 0
>
>
> In your fcst GRIB2 file I see the following valid time:
>
> > > wgrib2 -vt gfs.0p25.2021012912.f036.grib2
> > 1:0:vt=2021013100
>
>
> In your obs NetCDF file I see the following valid time:
>
> > hdr_vld_table =
> > "20210201_230000" ;
>
>
> I hope that helps!
>
> Julie
>
>
> On Fri, Jul 2, 2021 at 3:42 PM Julie Prestopnik <jpresto at ucar.edu>
wrote:
>
> > Thank you for the information. You can run "--version" with any of
the
> > MET tools to see what version you are running. For example:
> >
> > point_stat --version
> >
> > Julie
> >
> > On Fri, Jul 2, 2021 at 3:08 PM Raúl Valenzuela Ruiz via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423 >
> > >
> > > Sorry, I forgot the MET version, how should I verify that?
> > >
> > > Raúl
> > >
> > > Raúl Valenzuela
> > > Assistant Professor
> > > Institute of Engineering Sciences
> > > University of O’Higgins
> > > Rancagua, Chile
> > > Tel +56 2332 86019
> > > raul.valenzuela at uoh.cl
> > > On Jul 2, 2021, 5:04 PM -0400, Raúl Valenzuela Ruiz <
> > > raul.valenzuela at uoh.cl>, wrote:
> > > > Thanks Julie!
> > > >
> > > > I am running the following command:
> > > >
> > > > point_stat \
> > > > $MET_RUN/input/fcst/gfs.0p25.2021012912.f036.grib2 \
> > > >
> > > $MET_RUN/output/ascii2nc/Acumulado_Precipitacion_2021-02-
01_23_00_Int_6.nc \
> > > > $MET_RUN/config/PointStatConfig_run01 \
> > > > -outdir $MET_RUN/output/point_stat -v 2
> > > >
> > > > I thought that the problema was related to the masking area,
so
> > > replaced the mask grid with the following configuration:
> > > >
> > > > mask = {
> > > > grid = [];
> > > > poly = [];
> > > > sid = [];
> > > > llpnt = [ { name = "BOX";
> > > > lat_thresh = >=-36&&<=-32;
> > > > lon_thresh = >=108&&<=111;
> > > > }
> > > > ];
> > > > }
> > > >
> > > > I realized that the GFS GRIB2 file had 0-360 longitude
coordinates. I
> > > tried to change them to -180 to 180 with CDO but I got some
weird results
> > > (longitudes 180 to 560). Then I read that GRIB2 only accepts
positive
> > > coordinates, so I kept the 0-360 longitude in the GFS GRIB2 file
and
> > > shifted the longitude coordinates in my observations. No success
so far.
> > > >
> > > > Best,
> > > >
> > > > Raúl
> > > >
> > > > Raúl Valenzuela
> > > > Assistant Professor
> > > > Institute of Engineering Sciences
> > > > University of O’Higgins
> > > > Rancagua, Chile
> > > > Tel +56 2332 86019
> > > > raul.valenzuela at uoh.cl
> > > > On Jul 2, 2021, 4:43 PM -0400, Julie Prestopnik via RT <
> > > met_help at ucar.edu>, wrote:
> > > > > Hi Raúl.
> > > > >
> > > > > I see you are having trouble running Point-Stat with your
data and
> > > trouble
> > > > > with running MET with the mask produced by gen_vx_mask.
Thank you for
> > > > > sending your data. Can you please send me the command line
calls that
> > > you
> > > > > are using so that I can be sure to run the same thing you
are? Also,
> > > can
> > > > > you please tell me what version of MET you are running?
> > > > >
> > > > > Thanks!
> > > > >
> > > > > Julie
> > > > >
> > > > > On Fri, Jul 2, 2021 at 1:57 PM Raúl Valenzuela Ruiz via RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > >
> > > > > > Fri Jul 02 13:57:01 2021: Request 100423 was acted upon.
> > > > > > Transaction: Ticket created by raul.valenzuela at uoh.cl
> > > > > > Queue: met_help
> > > > > > Subject: Running Point-stat for precipitation verification
> > > > > > Owner: Nobody
> > > > > > Requestors: raul.valenzuela at uoh.cl
> > > > > > Status: new
> > > > > > Ticket <URL:
> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423 >
> > > > > >
> > > > > >
> > > > > > Greetings,
> > > > > >
> > > > > > My goal is to verify a GFS forecast of accumulated
precipitation
> > > against
> > > > > > point observations for an event. This is my first time
using MET,
> > > so I
> > > > > > followed the Point-stat tutorial successfully (
> > > > > > http://dtcenter.org/met-online-tutorial-metv8-0/point-
stat).
> > > However,
> > > > > > when I try to run the same example with my own data I get
no result
> > > (e.g.
> > > > > > output stat files containing headers only). In addition, I
have not
> > > being
> > > > > > able to run MET with the mask produced by gen_vx_mask, I
get an
> > > error
> > > > > > saying that the file does not exit (while it does).
> > > > > >
> > > > > > I have attached the files I am using to run my custom
case. The
> > > > > > directories have the same structure as the example in the
tutorial.
> > > > > >
> > > > > > By the way, I am using the docker version of MET.
> > > > > >
> > > > > > Best regards,
> > > > > >
> > > > > > Raúl Valenzuela
> > > > > > Assistant Professor
> > > > > > Institute of Engineering Sciences
> > > > > > University of O’Higgins
> > > > > > Rancagua, Chile
> > > > > > Tel +56 2332 86019
> > > > > > raul.valenzuela at uoh.cl
> > > > > >
> > > > > >
> > > > >
> > > > > --
> > > > > Julie Prestopnik (she/her)
> > > > > Software Engineer
> > > > > National Center for Atmospheric Research
> > > > > Research Applications Laboratory
> > > > > Email: jpresto at ucar.edu
> > > > >
> > > > > My working day may not be your working day. Please do not
feel
> > > obliged to
> > > > > reply to this email outside of your normal working hours.
> > > > >
> > >
> > >
> >
> > --
> > Julie Prestopnik (she/her)
> > Software Engineer
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Email: jpresto at ucar.edu
> >
> > My working day may not be your working day. Please do not feel
obliged to
> > reply to this email outside of your normal working hours.
> >
>
>
> --
> Julie Prestopnik (she/her)
> Software Engineer
> National Center for Atmospheric Research
> Research Applications Laboratory
> Email: jpresto at ucar.edu
>
> My working day may not be your working day. Please do not feel
obliged to
> reply to this email outside of your normal working hours.
>
------------------------------------------------
Subject: Running Point-stat for precipitation verification
From: Raúl Valenzuela Ruiz
Time: Tue Jul 06 08:30:23 2021
Hi Julia,
Should I open a new ticket for a new question or could I follow this
one?
Just in case I can follow this one:
I would like to compute the average precipitation of all observation
points inside a given forecast grid cell (e.g. verify average rain
gauge precipitation inside grid and grid forecast value). I went
through all the PDFs in the tutorial but I am not sure how to achieve
this. I though it was with Series-Analysis but the PDF only shows a
grid-grid example.
Best regards,
Raúl
Raúl Valenzuela
Assistant Professor
Institute of Engineering Sciences
University of O’Higgins
Rancagua, Chile
Tel +56 2332 86019
raul.valenzuela at uoh.cl
On 2 Jul 2021 18:28 -0400, Raúl Valenzuela Ruiz
<raul.valenzuela at uoh.cl>, wrote:
> I knew it was a silly error. Thanks!
>
> Raúl
>
> Raúl Valenzuela
> Assistant Professor
> Institute of Engineering Sciences
> University of O’Higgins
> Rancagua, Chile
> Tel +56 2332 86019
> raul.valenzuela at uoh.cl
> On Jul 2, 2021, 5:54 PM -0400, Julie Prestopnik via RT
<met_help at ucar.edu>, wrote:
> > Running the files with a verbosity level of 2, I see the following
output,
> > which indicates that 1040 observations were processed, but all
1040
> > observations were rejected due to the valid time.
> >
> > DEBUG 2: Processing APCP/A6 versus APCP/A6, for observation type
ADPSFC,
> > > over region BOX, for interpolation method NEAREST(1), using 0
matched pairs.
> > > DEBUG 2: Number of matched pairs = 0
> > > DEBUG 2: Observations processed = 1040
> > > DEBUG 2: Rejected: station id = 0
> > > DEBUG 2: Rejected: obs var name = 0
> > > DEBUG 2: Rejected: valid time = 1040
> > > DEBUG 2: Rejected: bad obs value = 0
> > > DEBUG 2: Rejected: off the grid = 0
> > > DEBUG 2: Rejected: topography = 0
> > > DEBUG 2: Rejected: level mismatch = 0
> > > DEBUG 2: Rejected: quality marker = 0
> > > DEBUG 2: Rejected: message type = 0
> > > DEBUG 2: Rejected: masking region = 0
> > > DEBUG 2: Rejected: bad fcst value = 0
> > > DEBUG 2: Rejected: bad climo mean = 0
> > > DEBUG 2: Rejected: bad climo stdev = 0
> > > DEBUG 2: Rejected: mpr filter = 0
> > > DEBUG 2: Rejected: duplicates = 0
> >
> >
> > In your fcst GRIB2 file I see the following valid time:
> >
> > > > wgrib2 -vt gfs.0p25.2021012912.f036.grib2
> > > 1:0:vt=2021013100
> >
> >
> > In your obs NetCDF file I see the following valid time:
> >
> > > hdr_vld_table =
> > > "20210201_230000" ;
> >
> >
> > I hope that helps!
> >
> > Julie
> >
> >
> > On Fri, Jul 2, 2021 at 3:42 PM Julie Prestopnik <jpresto at ucar.edu>
wrote:
> >
> > > Thank you for the information. You can run "--version" with any
of the
> > > MET tools to see what version you are running. For example:
> > >
> > > point_stat --version
> > >
> > > Julie
> > >
> > > On Fri, Jul 2, 2021 at 3:08 PM Raúl Valenzuela Ruiz via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423
>
> > > >
> > > > Sorry, I forgot the MET version, how should I verify that?
> > > >
> > > > Raúl
> > > >
> > > > Raúl Valenzuela
> > > > Assistant Professor
> > > > Institute of Engineering Sciences
> > > > University of O’Higgins
> > > > Rancagua, Chile
> > > > Tel +56 2332 86019
> > > > raul.valenzuela at uoh.cl
> > > > On Jul 2, 2021, 5:04 PM -0400, Raúl Valenzuela Ruiz <
> > > > raul.valenzuela at uoh.cl>, wrote:
> > > > > Thanks Julie!
> > > > >
> > > > > I am running the following command:
> > > > >
> > > > > point_stat \
> > > > > $MET_RUN/input/fcst/gfs.0p25.2021012912.f036.grib2 \
> > > > >
> > > > $MET_RUN/output/ascii2nc/Acumulado_Precipitacion_2021-02-
01_23_00_Int_6.nc \
> > > > > $MET_RUN/config/PointStatConfig_run01 \
> > > > > -outdir $MET_RUN/output/point_stat -v 2
> > > > >
> > > > > I thought that the problema was related to the masking area,
so
> > > > replaced the mask grid with the following configuration:
> > > > >
> > > > > mask = {
> > > > > grid = [];
> > > > > poly = [];
> > > > > sid = [];
> > > > > llpnt = [ { name = "BOX";
> > > > > lat_thresh = >=-36&&<=-32;
> > > > > lon_thresh = >=108&&<=111;
> > > > > }
> > > > > ];
> > > > > }
> > > > >
> > > > > I realized that the GFS GRIB2 file had 0-360 longitude
coordinates. I
> > > > tried to change them to -180 to 180 with CDO but I got some
weird results
> > > > (longitudes 180 to 560). Then I read that GRIB2 only accepts
positive
> > > > coordinates, so I kept the 0-360 longitude in the GFS GRIB2
file and
> > > > shifted the longitude coordinates in my observations. No
success so far.
> > > > >
> > > > > Best,
> > > > >
> > > > > Raúl
> > > > >
> > > > > Raúl Valenzuela
> > > > > Assistant Professor
> > > > > Institute of Engineering Sciences
> > > > > University of O’Higgins
> > > > > Rancagua, Chile
> > > > > Tel +56 2332 86019
> > > > > raul.valenzuela at uoh.cl
> > > > > On Jul 2, 2021, 4:43 PM -0400, Julie Prestopnik via RT <
> > > > met_help at ucar.edu>, wrote:
> > > > > > Hi Raúl.
> > > > > >
> > > > > > I see you are having trouble running Point-Stat with your
data and
> > > > trouble
> > > > > > with running MET with the mask produced by gen_vx_mask.
Thank you for
> > > > > > sending your data. Can you please send me the command line
calls that
> > > > you
> > > > > > are using so that I can be sure to run the same thing you
are? Also,
> > > > can
> > > > > > you please tell me what version of MET you are running?
> > > > > >
> > > > > > Thanks!
> > > > > >
> > > > > > Julie
> > > > > >
> > > > > > On Fri, Jul 2, 2021 at 1:57 PM Raúl Valenzuela Ruiz via RT
<
> > > > > > met_help at ucar.edu> wrote:
> > > > > >
> > > > > > >
> > > > > > > Fri Jul 02 13:57:01 2021: Request 100423 was acted upon.
> > > > > > > Transaction: Ticket created by raul.valenzuela at uoh.cl
> > > > > > > Queue: met_help
> > > > > > > Subject: Running Point-stat for precipitation
verification
> > > > > > > Owner: Nobody
> > > > > > > Requestors: raul.valenzuela at uoh.cl
> > > > > > > Status: new
> > > > > > > Ticket <URL:
> > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423 >
> > > > > > >
> > > > > > >
> > > > > > > Greetings,
> > > > > > >
> > > > > > > My goal is to verify a GFS forecast of accumulated
precipitation
> > > > against
> > > > > > > point observations for an event. This is my first time
using MET,
> > > > so I
> > > > > > > followed the Point-stat tutorial successfully (
> > > > > > > http://dtcenter.org/met-online-tutorial-metv8-0/point-
stat).
> > > > However,
> > > > > > > when I try to run the same example with my own data I
get no result
> > > > (e.g.
> > > > > > > output stat files containing headers only). In addition,
I have not
> > > > being
> > > > > > > able to run MET with the mask produced by gen_vx_mask, I
get an
> > > > error
> > > > > > > saying that the file does not exit (while it does).
> > > > > > >
> > > > > > > I have attached the files I am using to run my custom
case. The
> > > > > > > directories have the same structure as the example in
the tutorial.
> > > > > > >
> > > > > > > By the way, I am using the docker version of MET.
> > > > > > >
> > > > > > > Best regards,
> > > > > > >
> > > > > > > Raúl Valenzuela
> > > > > > > Assistant Professor
> > > > > > > Institute of Engineering Sciences
> > > > > > > University of O’Higgins
> > > > > > > Rancagua, Chile
> > > > > > > Tel +56 2332 86019
> > > > > > > raul.valenzuela at uoh.cl
> > > > > > >
> > > > > > >
> > > > > >
> > > > > > --
> > > > > > Julie Prestopnik (she/her)
> > > > > > Software Engineer
> > > > > > National Center for Atmospheric Research
> > > > > > Research Applications Laboratory
> > > > > > Email: jpresto at ucar.edu
> > > > > >
> > > > > > My working day may not be your working day. Please do not
feel
> > > > obliged to
> > > > > > reply to this email outside of your normal working hours.
> > > > > >
> > > >
> > > >
> > >
> > > --
> > > Julie Prestopnik (she/her)
> > > Software Engineer
> > > National Center for Atmospheric Research
> > > Research Applications Laboratory
> > > Email: jpresto at ucar.edu
> > >
> > > My working day may not be your working day. Please do not feel
obliged to
> > > reply to this email outside of your normal working hours.
> > >
> >
> >
> > --
> > Julie Prestopnik (she/her)
> > Software Engineer
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Email: jpresto at ucar.edu
> >
> > My working day may not be your working day. Please do not feel
obliged to
> > reply to this email outside of your normal working hours.
> >
------------------------------------------------
Subject: Running Point-stat for precipitation verification
From: Julie Prestopnik
Time: Tue Jul 06 10:02:45 2021
Normally, I would say, please open a new ticket. However, we are
transitioning to using GitHub Discussions as of today. If you were to
submit a new ticket to "met_help at ucar.edu" you would receive an auto-
reply
with the message below. Given that, if you could please post this to
the
Discussions Forum, we would appreciate it. Please let me know if you
have
any questions. Also, I should mention that to encourage community
support
of the Discussions, the MET team will not post a reply until the
following
day to allow others to answer.
Julie
"***** ALERT: THIS E-MAIL ADDRESS IS NO LONGER IN USE FOR SUPPORT OF
THE
METPLUS VERIFICATION SYSTEM. WE ARE NO LONGER SUPPORTING NEW
QUESTIONS OR
MONITORING THIS EMAIL.
PLEASE CREATE A FREE GITHUB ACCOUNT AND POST YOUR QUESTIONS TO THE
METPLUS
COMPONENTS DISCUSSION FORUM AT
https://github.com/dtcenter/METplus/discussions.
*****
Thank you,
met_help at ucar.edu"
On Tue, Jul 6, 2021 at 8:30 AM Raúl Valenzuela Ruiz via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423 >
>
> Hi Julia,
>
> Should I open a new ticket for a new question or could I follow this
one?
>
> Just in case I can follow this one:
>
> I would like to compute the average precipitation of all observation
> points inside a given forecast grid cell (e.g. verify average rain
gauge
> precipitation inside grid and grid forecast value). I went through
all the
> PDFs in the tutorial but I am not sure how to achieve this. I though
it was
> with Series-Analysis but the PDF only shows a grid-grid example.
>
> Best regards,
>
> Raúl
>
> Raúl Valenzuela
> Assistant Professor
> Institute of Engineering Sciences
> University of O’Higgins
> Rancagua, Chile
> Tel +56 2332 86019
> raul.valenzuela at uoh.cl
> On 2 Jul 2021 18:28 -0400, Raúl Valenzuela Ruiz
<raul.valenzuela at uoh.cl>,
> wrote:
> > I knew it was a silly error. Thanks!
> >
> > Raúl
> >
> > Raúl Valenzuela
> > Assistant Professor
> > Institute of Engineering Sciences
> > University of O’Higgins
> > Rancagua, Chile
> > Tel +56 2332 86019
> > raul.valenzuela at uoh.cl
> > On Jul 2, 2021, 5:54 PM -0400, Julie Prestopnik via RT <
> met_help at ucar.edu>, wrote:
> > > Running the files with a verbosity level of 2, I see the
following
> output,
> > > which indicates that 1040 observations were processed, but all
1040
> > > observations were rejected due to the valid time.
> > >
> > > DEBUG 2: Processing APCP/A6 versus APCP/A6, for observation type
> ADPSFC,
> > > > over region BOX, for interpolation method NEAREST(1), using 0
> matched pairs.
> > > > DEBUG 2: Number of matched pairs = 0
> > > > DEBUG 2: Observations processed = 1040
> > > > DEBUG 2: Rejected: station id = 0
> > > > DEBUG 2: Rejected: obs var name = 0
> > > > DEBUG 2: Rejected: valid time = 1040
> > > > DEBUG 2: Rejected: bad obs value = 0
> > > > DEBUG 2: Rejected: off the grid = 0
> > > > DEBUG 2: Rejected: topography = 0
> > > > DEBUG 2: Rejected: level mismatch = 0
> > > > DEBUG 2: Rejected: quality marker = 0
> > > > DEBUG 2: Rejected: message type = 0
> > > > DEBUG 2: Rejected: masking region = 0
> > > > DEBUG 2: Rejected: bad fcst value = 0
> > > > DEBUG 2: Rejected: bad climo mean = 0
> > > > DEBUG 2: Rejected: bad climo stdev = 0
> > > > DEBUG 2: Rejected: mpr filter = 0
> > > > DEBUG 2: Rejected: duplicates = 0
> > >
> > >
> > > In your fcst GRIB2 file I see the following valid time:
> > >
> > > > > wgrib2 -vt gfs.0p25.2021012912.f036.grib2
> > > > 1:0:vt=2021013100
> > >
> > >
> > > In your obs NetCDF file I see the following valid time:
> > >
> > > > hdr_vld_table =
> > > > "20210201_230000" ;
> > >
> > >
> > > I hope that helps!
> > >
> > > Julie
> > >
> > >
> > > On Fri, Jul 2, 2021 at 3:42 PM Julie Prestopnik
<jpresto at ucar.edu>
> wrote:
> > >
> > > > Thank you for the information. You can run "--version" with
any of
> the
> > > > MET tools to see what version you are running. For example:
> > > >
> > > > point_stat --version
> > > >
> > > > Julie
> > > >
> > > > On Fri, Jul 2, 2021 at 3:08 PM Raúl Valenzuela Ruiz via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423 >
> > > > >
> > > > > Sorry, I forgot the MET version, how should I verify that?
> > > > >
> > > > > Raúl
> > > > >
> > > > > Raúl Valenzuela
> > > > > Assistant Professor
> > > > > Institute of Engineering Sciences
> > > > > University of O’Higgins
> > > > > Rancagua, Chile
> > > > > Tel +56 2332 86019
> > > > > raul.valenzuela at uoh.cl
> > > > > On Jul 2, 2021, 5:04 PM -0400, Raúl Valenzuela Ruiz <
> > > > > raul.valenzuela at uoh.cl>, wrote:
> > > > > > Thanks Julie!
> > > > > >
> > > > > > I am running the following command:
> > > > > >
> > > > > > point_stat \
> > > > > > $MET_RUN/input/fcst/gfs.0p25.2021012912.f036.grib2 \
> > > > > >
> > > > >
> $MET_RUN/output/ascii2nc/Acumulado_Precipitacion_2021-02-
01_23_00_Int_6.nc \
> > > > > > $MET_RUN/config/PointStatConfig_run01 \
> > > > > > -outdir $MET_RUN/output/point_stat -v 2
> > > > > >
> > > > > > I thought that the problema was related to the masking
area, so
> > > > > replaced the mask grid with the following configuration:
> > > > > >
> > > > > > mask = {
> > > > > > grid = [];
> > > > > > poly = [];
> > > > > > sid = [];
> > > > > > llpnt = [ { name = "BOX";
> > > > > > lat_thresh = >=-36&&<=-32;
> > > > > > lon_thresh = >=108&&<=111;
> > > > > > }
> > > > > > ];
> > > > > > }
> > > > > >
> > > > > > I realized that the GFS GRIB2 file had 0-360 longitude
> coordinates. I
> > > > > tried to change them to -180 to 180 with CDO but I got some
weird
> results
> > > > > (longitudes 180 to 560). Then I read that GRIB2 only accepts
> positive
> > > > > coordinates, so I kept the 0-360 longitude in the GFS GRIB2
file
> and
> > > > > shifted the longitude coordinates in my observations. No
success
> so far.
> > > > > >
> > > > > > Best,
> > > > > >
> > > > > > Raúl
> > > > > >
> > > > > > Raúl Valenzuela
> > > > > > Assistant Professor
> > > > > > Institute of Engineering Sciences
> > > > > > University of O’Higgins
> > > > > > Rancagua, Chile
> > > > > > Tel +56 2332 86019
> > > > > > raul.valenzuela at uoh.cl
> > > > > > On Jul 2, 2021, 4:43 PM -0400, Julie Prestopnik via RT <
> > > > > met_help at ucar.edu>, wrote:
> > > > > > > Hi Raúl.
> > > > > > >
> > > > > > > I see you are having trouble running Point-Stat with
your data
> and
> > > > > trouble
> > > > > > > with running MET with the mask produced by gen_vx_mask.
Thank
> you for
> > > > > > > sending your data. Can you please send me the command
line
> calls that
> > > > > you
> > > > > > > are using so that I can be sure to run the same thing
you are?
> Also,
> > > > > can
> > > > > > > you please tell me what version of MET you are running?
> > > > > > >
> > > > > > > Thanks!
> > > > > > >
> > > > > > > Julie
> > > > > > >
> > > > > > > On Fri, Jul 2, 2021 at 1:57 PM Raúl Valenzuela Ruiz via
RT <
> > > > > > > met_help at ucar.edu> wrote:
> > > > > > >
> > > > > > > >
> > > > > > > > Fri Jul 02 13:57:01 2021: Request 100423 was acted
upon.
> > > > > > > > Transaction: Ticket created by raul.valenzuela at uoh.cl
> > > > > > > > Queue: met_help
> > > > > > > > Subject: Running Point-stat for precipitation
verification
> > > > > > > > Owner: Nobody
> > > > > > > > Requestors: raul.valenzuela at uoh.cl
> > > > > > > > Status: new
> > > > > > > > Ticket <URL:
> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100423 >
> > > > > > > >
> > > > > > > >
> > > > > > > > Greetings,
> > > > > > > >
> > > > > > > > My goal is to verify a GFS forecast of accumulated
> precipitation
> > > > > against
> > > > > > > > point observations for an event. This is my first time
using
> MET,
> > > > > so I
> > > > > > > > followed the Point-stat tutorial successfully (
> > > > > > > > http://dtcenter.org/met-online-tutorial-metv8-0/point-
stat).
> > > > > However,
> > > > > > > > when I try to run the same example with my own data I
get no
> result
> > > > > (e.g.
> > > > > > > > output stat files containing headers only). In
addition, I
> have not
> > > > > being
> > > > > > > > able to run MET with the mask produced by gen_vx_mask,
I get
> an
> > > > > error
> > > > > > > > saying that the file does not exit (while it does).
> > > > > > > >
> > > > > > > > I have attached the files I am using to run my custom
case.
> The
> > > > > > > > directories have the same structure as the example in
the
> tutorial.
> > > > > > > >
> > > > > > > > By the way, I am using the docker version of MET.
> > > > > > > >
> > > > > > > > Best regards,
> > > > > > > >
> > > > > > > > Raúl Valenzuela
> > > > > > > > Assistant Professor
> > > > > > > > Institute of Engineering Sciences
> > > > > > > > University of O’Higgins
> > > > > > > > Rancagua, Chile
> > > > > > > > Tel +56 2332 86019
> > > > > > > > raul.valenzuela at uoh.cl
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > > --
> > > > > > > Julie Prestopnik (she/her)
> > > > > > > Software Engineer
> > > > > > > National Center for Atmospheric Research
> > > > > > > Research Applications Laboratory
> > > > > > > Email: jpresto at ucar.edu
> > > > > > >
> > > > > > > My working day may not be your working day. Please do
not feel
> > > > > obliged to
> > > > > > > reply to this email outside of your normal working
hours.
> > > > > > >
> > > > >
> > > > >
> > > >
> > > > --
> > > > Julie Prestopnik (she/her)
> > > > Software Engineer
> > > > National Center for Atmospheric Research
> > > > Research Applications Laboratory
> > > > Email: jpresto at ucar.edu
> > > >
> > > > My working day may not be your working day. Please do not feel
> obliged to
> > > > reply to this email outside of your normal working hours.
> > > >
> > >
> > >
> > > --
> > > Julie Prestopnik (she/her)
> > > Software Engineer
> > > National Center for Atmospheric Research
> > > Research Applications Laboratory
> > > Email: jpresto at ucar.edu
> > >
> > > My working day may not be your working day. Please do not feel
obliged
> to
> > > reply to this email outside of your normal working hours.
> > >
>
>
--
Julie Prestopnik (she/her)
Software Engineer
National Center for Atmospheric Research
Research Applications Laboratory
Email: jpresto at ucar.edu
My working day may not be your working day. Please do not feel
obliged to
reply to this email outside of your normal working hours.
------------------------------------------------
More information about the Met_help
mailing list