[Met_help] [rt.rap.ucar.edu #75487] History for Trouble with Point-stat tool for some variables

John Halley Gotway via RT met_help at ucar.edu
Mon Mar 14 15:16:26 MDT 2016


----------------------------------------------------------------
  Initial Request
----------------------------------------------------------------

I am having a problem using the point-stat tool with some different variables in the Prepbufr format.  For example the variable "vs" (forecast) for 10m meridional wind I believe should compare with VGRD level Z10 in both the Prepbufr point observations and the climatology.  Yet when I run the tool it doesn't process any matches, although it does find a forecast and climatology level.  I can run the grid-stat tool using the same forecast and climatology fields, but a gridded analysis instead of point observations, and it works fine.  I get the same problem with 200 mb level winds (P200), although I'm almost positive I have been able to get statistics such as RMSE with these variables in the point-stat tool a while ago.  When I run the tool I get this:

DEBUG 3: Number of matched pairs  = 0                                         
DEBUG 3: Observations processed   = 45684                                       
DEBUG 3: Rejected: GRIB code      = 45684                                                                                                       
I must have the wrong variable name.  I thought for zonal winds at 200 mb would have the variable name "UGRD" and the level is "P200" from this website http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html.  Do you know if this is incorrect or where I could get more info?


----------------------------------------------------------------
  Complete Ticket History
----------------------------------------------------------------

Subject: Trouble with Point-stat tool for some variables
From: Shannon Rees - NOAA Affiliate
Time: Thu Mar 10 12:21:03 2016

Here is the configuration file.  I am also uploading some obs and fcst
data
to the ftp site in the "REES" directory.

Thank you,
Shannon


On Thu, Mar 10, 2016 at 1:42 PM, met_help at ucar.edu via RT
<met_help at ucar.edu
> wrote:

> Greetings,
>
> This message has been automatically generated in response to the
> creation of a trouble ticket regarding:
>         "Trouble with Point-stat tool for some variables",
> a summary of which appears below.
>
> There is no need to reply to this message right now.  Your ticket
has been
> assigned an ID of [rt.rap.ucar.edu #75487].
>
> Please include the string:
>
>          [rt.rap.ucar.edu #75487]
>
> in the subject line of all future correspondence about this issue.
To do
> so,
> you may reply to this message.
>
>                         Thank you,
>                         met_help at ucar.edu
>
>
-------------------------------------------------------------------------
> I am having a problem using the point-stat tool with some different
> variables in the Prepbufr format.  For example the variable "vs"
(forecast)
> for 10m meridional wind I believe should compare with VGRD level Z10
in
> both the Prepbufr point observations and the climatology.  Yet when
I run
> the tool it doesn't process any matches, although it does find a
forecast
> and climatology level.  I can run the grid-stat tool using the same
> forecast and climatology fields, but a gridded analysis instead of
point
> observations, and it works fine.  I get the same problem with 200 mb
level
> winds (P200), although I'm almost positive I have been able to get
> statistics such as RMSE with these variables in the point-stat tool
a while
> ago.  When I run the tool I get this:
>
> DEBUG 3: Number of matched pairs  = 0
> DEBUG 3: Observations processed   = 45684
> DEBUG 3: Rejected: GRIB code      = 45684
> I must have the wrong variable name.  I thought for zonal winds at
200 mb
> would have the variable name "UGRD" and the level is "P200" from
this
> website http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html.  Do
you
> know if this is incorrect or where I could get more info?
>
>


--
Shannon Rees
Programming Scientist
Engility
Geophysical Fluid Dynamics Lab
201 Forrestal Rd Princeton, NJ
609-452-5384

------------------------------------------------
Subject: Trouble with Point-stat tool for some variables
From: John Halley Gotway
Time: Thu Mar 10 12:27:22 2016

Shannon,

Looking in your config file, I see the following:

grib_fields = [
  {
    name  = "${obs_var}";
    level = "${obs_lev}";
    cat_thresh = [  ];
  }
];

obs = {
   message_type = [ "${mess_type}" ];
   sid_exc      = [];
   field        = grib_fields;
}

Do you have a script to wrap these calls to Point-Stat?  I'm wondering
how
${obs_var} is getting set.

John


On Thu, Mar 10, 2016 at 12:21 PM, Shannon Rees - NOAA Affiliate via RT
<
met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75487 >
>
> Here is the configuration file.  I am also uploading some obs and
fcst data
> to the ftp site in the "REES" directory.
>
> Thank you,
> Shannon
>
>
> On Thu, Mar 10, 2016 at 1:42 PM, met_help at ucar.edu via RT <
> met_help at ucar.edu
> > wrote:
>
> > Greetings,
> >
> > This message has been automatically generated in response to the
> > creation of a trouble ticket regarding:
> >         "Trouble with Point-stat tool for some variables",
> > a summary of which appears below.
> >
> > There is no need to reply to this message right now.  Your ticket
has
> been
> > assigned an ID of [rt.rap.ucar.edu #75487].
> >
> > Please include the string:
> >
> >          [rt.rap.ucar.edu #75487]
> >
> > in the subject line of all future correspondence about this issue.
To do
> > so,
> > you may reply to this message.
> >
> >                         Thank you,
> >                         met_help at ucar.edu
> >
> >
-------------------------------------------------------------------------
> > I am having a problem using the point-stat tool with some
different
> > variables in the Prepbufr format.  For example the variable "vs"
> (forecast)
> > for 10m meridional wind I believe should compare with VGRD level
Z10 in
> > both the Prepbufr point observations and the climatology.  Yet
when I run
> > the tool it doesn't process any matches, although it does find a
forecast
> > and climatology level.  I can run the grid-stat tool using the
same
> > forecast and climatology fields, but a gridded analysis instead of
point
> > observations, and it works fine.  I get the same problem with 200
mb
> level
> > winds (P200), although I'm almost positive I have been able to get
> > statistics such as RMSE with these variables in the point-stat
tool a
> while
> > ago.  When I run the tool I get this:
> >
> > DEBUG 3: Number of matched pairs  = 0
> > DEBUG 3: Observations processed   = 45684
> > DEBUG 3: Rejected: GRIB code      = 45684
> > I must have the wrong variable name.  I thought for zonal winds at
200 mb
> > would have the variable name "UGRD" and the level is "P200" from
this
> > website http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html.
Do you
> > know if this is incorrect or where I could get more info?
> >
> >
>
>
> --
> Shannon Rees
> Programming Scientist
> Engility
> Geophysical Fluid Dynamics Lab
> 201 Forrestal Rd Princeton, NJ
> 609-452-5384
>
>

------------------------------------------------
Subject: Trouble with Point-stat tool for some variables
From: Shannon Rees - NOAA Affiliate
Time: Thu Mar 10 13:40:34 2016

Yes it is attached here.

On Thu, Mar 10, 2016 at 2:27 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:

> Shannon,
>
> Looking in your config file, I see the following:
>
> grib_fields = [
>   {
>     name  = "${obs_var}";
>     level = "${obs_lev}";
>     cat_thresh = [  ];
>   }
> ];
>
> obs = {
>    message_type = [ "${mess_type}" ];
>    sid_exc      = [];
>    field        = grib_fields;
> }
>
> Do you have a script to wrap these calls to Point-Stat?  I'm
wondering how
> ${obs_var} is getting set.
>
> John
>
>
> On Thu, Mar 10, 2016 at 12:21 PM, Shannon Rees - NOAA Affiliate via
RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75487 >
> >
> > Here is the configuration file.  I am also uploading some obs and
fcst
> data
> > to the ftp site in the "REES" directory.
> >
> > Thank you,
> > Shannon
> >
> >
> > On Thu, Mar 10, 2016 at 1:42 PM, met_help at ucar.edu via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > > Greetings,
> > >
> > > This message has been automatically generated in response to the
> > > creation of a trouble ticket regarding:
> > >         "Trouble with Point-stat tool for some variables",
> > > a summary of which appears below.
> > >
> > > There is no need to reply to this message right now.  Your
ticket has
> > been
> > > assigned an ID of [rt.rap.ucar.edu #75487].
> > >
> > > Please include the string:
> > >
> > >          [rt.rap.ucar.edu #75487]
> > >
> > > in the subject line of all future correspondence about this
issue. To
> do
> > > so,
> > > you may reply to this message.
> > >
> > >                         Thank you,
> > >                         met_help at ucar.edu
> > >
> > >
>
-------------------------------------------------------------------------
> > > I am having a problem using the point-stat tool with some
different
> > > variables in the Prepbufr format.  For example the variable "vs"
> > (forecast)
> > > for 10m meridional wind I believe should compare with VGRD level
Z10 in
> > > both the Prepbufr point observations and the climatology.  Yet
when I
> run
> > > the tool it doesn't process any matches, although it does find a
> forecast
> > > and climatology level.  I can run the grid-stat tool using the
same
> > > forecast and climatology fields, but a gridded analysis instead
of
> point
> > > observations, and it works fine.  I get the same problem with
200 mb
> > level
> > > winds (P200), although I'm almost positive I have been able to
get
> > > statistics such as RMSE with these variables in the point-stat
tool a
> > while
> > > ago.  When I run the tool I get this:
> > >
> > > DEBUG 3: Number of matched pairs  = 0
> > > DEBUG 3: Observations processed   = 45684
> > > DEBUG 3: Rejected: GRIB code      = 45684
> > > I must have the wrong variable name.  I thought for zonal winds
at 200
> mb
> > > would have the variable name "UGRD" and the level is "P200" from
this
> > > website http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html.
Do
> you
> > > know if this is incorrect or where I could get more info?
> > >
> > >
> >
> >
> > --
> > Shannon Rees
> > Programming Scientist
> > Engility
> > Geophysical Fluid Dynamics Lab
> > 201 Forrestal Rd Princeton, NJ
> > 609-452-5384
> >
> >
>
>


--
Shannon Rees
Programming Scientist
Engility
Geophysical Fluid Dynamics Lab
201 Forrestal Rd Princeton, NJ
609-452-5384

------------------------------------------------
Subject: Trouble with Point-stat tool for some variables
From: Shannon Rees - NOAA Affiliate
Time: Thu Mar 10 13:49:59 2016

And here is the description of how to use it:


On Thu, Mar 10, 2016 at 3:40 PM, Shannon Rees - NOAA Affiliate <
shannon.rees at noaa.gov> wrote:

> Yes it is attached here.
>
> On Thu, Mar 10, 2016 at 2:27 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
>> Shannon,
>>
>> Looking in your config file, I see the following:
>>
>> grib_fields = [
>>   {
>>     name  = "${obs_var}";
>>     level = "${obs_lev}";
>>     cat_thresh = [  ];
>>   }
>> ];
>>
>> obs = {
>>    message_type = [ "${mess_type}" ];
>>    sid_exc      = [];
>>    field        = grib_fields;
>> }
>>
>> Do you have a script to wrap these calls to Point-Stat?  I'm
wondering how
>> ${obs_var} is getting set.
>>
>> John
>>
>>
>> On Thu, Mar 10, 2016 at 12:21 PM, Shannon Rees - NOAA Affiliate via
RT <
>> met_help at ucar.edu> wrote:
>>
>> >
>> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75487 >
>> >
>> > Here is the configuration file.  I am also uploading some obs and
fcst
>> data
>> > to the ftp site in the "REES" directory.
>> >
>> > Thank you,
>> > Shannon
>> >
>> >
>> > On Thu, Mar 10, 2016 at 1:42 PM, met_help at ucar.edu via RT <
>> > met_help at ucar.edu
>> > > wrote:
>> >
>> > > Greetings,
>> > >
>> > > This message has been automatically generated in response to
the
>> > > creation of a trouble ticket regarding:
>> > >         "Trouble with Point-stat tool for some variables",
>> > > a summary of which appears below.
>> > >
>> > > There is no need to reply to this message right now.  Your
ticket has
>> > been
>> > > assigned an ID of [rt.rap.ucar.edu #75487].
>> > >
>> > > Please include the string:
>> > >
>> > >          [rt.rap.ucar.edu #75487]
>> > >
>> > > in the subject line of all future correspondence about this
issue. To
>> do
>> > > so,
>> > > you may reply to this message.
>> > >
>> > >                         Thank you,
>> > >                         met_help at ucar.edu
>> > >
>> > >
>>
-------------------------------------------------------------------------
>> > > I am having a problem using the point-stat tool with some
different
>> > > variables in the Prepbufr format.  For example the variable
"vs"
>> > (forecast)
>> > > for 10m meridional wind I believe should compare with VGRD
level Z10
>> in
>> > > both the Prepbufr point observations and the climatology.  Yet
when I
>> run
>> > > the tool it doesn't process any matches, although it does find
a
>> forecast
>> > > and climatology level.  I can run the grid-stat tool using the
same
>> > > forecast and climatology fields, but a gridded analysis instead
of
>> point
>> > > observations, and it works fine.  I get the same problem with
200 mb
>> > level
>> > > winds (P200), although I'm almost positive I have been able to
get
>> > > statistics such as RMSE with these variables in the point-stat
tool a
>> > while
>> > > ago.  When I run the tool I get this:
>> > >
>> > > DEBUG 3: Number of matched pairs  = 0
>> > > DEBUG 3: Observations processed   = 45684
>> > > DEBUG 3: Rejected: GRIB code      = 45684
>> > > I must have the wrong variable name.  I thought for zonal winds
at
>> 200 mb
>> > > would have the variable name "UGRD" and the level is "P200"
from this
>> > > website
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html.  Do
>> you
>> > > know if this is incorrect or where I could get more info?
>> > >
>> > >
>> >
>> >
>> > --
>> > Shannon Rees
>> > Programming Scientist
>> > Engility
>> > Geophysical Fluid Dynamics Lab
>> > 201 Forrestal Rd Princeton, NJ
>> > 609-452-5384
>> >
>> >
>>
>>
>
>
> --
> Shannon Rees
> Programming Scientist
> Engility
> Geophysical Fluid Dynamics Lab
> 201 Forrestal Rd Princeton, NJ
> 609-452-5384
>



--
Shannon Rees
Programming Scientist
Engility
Geophysical Fluid Dynamics Lab
201 Forrestal Rd Princeton, NJ
609-452-5384

------------------------------------------------
Subject: Trouble with Point-stat tool for some variables
From: John Halley Gotway
Time: Fri Mar 11 14:31:32 2016

Shannon,

I ran with your data and was able to reproduce the behavior you're
seeing:

DEBUG 3: Observations processed   = 46351
DEBUG 3: Rejected: SID exclusion  = 0
DEBUG 3: Rejected: GRIB code      = 46351

Next, I ran plot_point_obs to see if your data has any wind obs in it:

plot_point_obs prepbufr.gdas.2015080100.nc prepbufr.gdas.2015080100.ps
-gc
34
...
DEBUG 2: Finished plotting 0 locations.
DEBUG 2: Skipped 0 locations off the grid.

And it does not... thus no matches.  Since you're using GDAS point
observations, I likely has something to do with the quality marker
you're
using.  Please look at the NOTE on the bottom of this page about GDAS
quality markers for surface observations:
   http://www.dtcenter.org/met/users/downloads/observation_data.php

Thanks,
John

On Thu, Mar 10, 2016 at 1:49 PM, Shannon Rees - NOAA Affiliate via RT
<
met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75487 >
>
> And here is the description of how to use it:
>
>
> On Thu, Mar 10, 2016 at 3:40 PM, Shannon Rees - NOAA Affiliate <
> shannon.rees at noaa.gov> wrote:
>
> > Yes it is attached here.
> >
> > On Thu, Mar 10, 2016 at 2:27 PM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> >> Shannon,
> >>
> >> Looking in your config file, I see the following:
> >>
> >> grib_fields = [
> >>   {
> >>     name  = "${obs_var}";
> >>     level = "${obs_lev}";
> >>     cat_thresh = [  ];
> >>   }
> >> ];
> >>
> >> obs = {
> >>    message_type = [ "${mess_type}" ];
> >>    sid_exc      = [];
> >>    field        = grib_fields;
> >> }
> >>
> >> Do you have a script to wrap these calls to Point-Stat?  I'm
wondering
> how
> >> ${obs_var} is getting set.
> >>
> >> John
> >>
> >>
> >> On Thu, Mar 10, 2016 at 12:21 PM, Shannon Rees - NOAA Affiliate
via RT <
> >> met_help at ucar.edu> wrote:
> >>
> >> >
> >> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75487 >
> >> >
> >> > Here is the configuration file.  I am also uploading some obs
and fcst
> >> data
> >> > to the ftp site in the "REES" directory.
> >> >
> >> > Thank you,
> >> > Shannon
> >> >
> >> >
> >> > On Thu, Mar 10, 2016 at 1:42 PM, met_help at ucar.edu via RT <
> >> > met_help at ucar.edu
> >> > > wrote:
> >> >
> >> > > Greetings,
> >> > >
> >> > > This message has been automatically generated in response to
the
> >> > > creation of a trouble ticket regarding:
> >> > >         "Trouble with Point-stat tool for some variables",
> >> > > a summary of which appears below.
> >> > >
> >> > > There is no need to reply to this message right now.  Your
ticket
> has
> >> > been
> >> > > assigned an ID of [rt.rap.ucar.edu #75487].
> >> > >
> >> > > Please include the string:
> >> > >
> >> > >          [rt.rap.ucar.edu #75487]
> >> > >
> >> > > in the subject line of all future correspondence about this
issue.
> To
> >> do
> >> > > so,
> >> > > you may reply to this message.
> >> > >
> >> > >                         Thank you,
> >> > >                         met_help at ucar.edu
> >> > >
> >> > >
> >>
>
-------------------------------------------------------------------------
> >> > > I am having a problem using the point-stat tool with some
different
> >> > > variables in the Prepbufr format.  For example the variable
"vs"
> >> > (forecast)
> >> > > for 10m meridional wind I believe should compare with VGRD
level Z10
> >> in
> >> > > both the Prepbufr point observations and the climatology.
Yet when
> I
> >> run
> >> > > the tool it doesn't process any matches, although it does
find a
> >> forecast
> >> > > and climatology level.  I can run the grid-stat tool using
the same
> >> > > forecast and climatology fields, but a gridded analysis
instead of
> >> point
> >> > > observations, and it works fine.  I get the same problem with
200 mb
> >> > level
> >> > > winds (P200), although I'm almost positive I have been able
to get
> >> > > statistics such as RMSE with these variables in the point-
stat tool
> a
> >> > while
> >> > > ago.  When I run the tool I get this:
> >> > >
> >> > > DEBUG 3: Number of matched pairs  = 0
> >> > > DEBUG 3: Observations processed   = 45684
> >> > > DEBUG 3: Rejected: GRIB code      = 45684
> >> > > I must have the wrong variable name.  I thought for zonal
winds at
> >> 200 mb
> >> > > would have the variable name "UGRD" and the level is "P200"
from
> this
> >> > > website
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html.
> Do
> >> you
> >> > > know if this is incorrect or where I could get more info?
> >> > >
> >> > >
> >> >
> >> >
> >> > --
> >> > Shannon Rees
> >> > Programming Scientist
> >> > Engility
> >> > Geophysical Fluid Dynamics Lab
> >> > 201 Forrestal Rd Princeton, NJ
> >> > 609-452-5384
> >> >
> >> >
> >>
> >>
> >
> >
> > --
> > Shannon Rees
> > Programming Scientist
> > Engility
> > Geophysical Fluid Dynamics Lab
> > 201 Forrestal Rd Princeton, NJ
> > 609-452-5384
> >
>
>
>
> --
> Shannon Rees
> Programming Scientist
> Engility
> Geophysical Fluid Dynamics Lab
> 201 Forrestal Rd Princeton, NJ
> 609-452-5384
>
>

------------------------------------------------
Subject: Trouble with Point-stat tool for some variables
From: Shannon Rees - NOAA Affiliate
Time: Mon Mar 14 14:05:28 2016

Hi John,

I think I downloaded all of the GDAS data with quality marker = 9
which
should preserve the surface variables.  I am downloading some more
GDAS
data with quality marker = 15 to test it and I'll let you know how
that
goes.

Thanks,

On Fri, Mar 11, 2016 at 4:31 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:

> Shannon,
>
> I ran with your data and was able to reproduce the behavior you're
seeing:
>
> DEBUG 3: Observations processed   = 46351
> DEBUG 3: Rejected: SID exclusion  = 0
> DEBUG 3: Rejected: GRIB code      = 46351
>
> Next, I ran plot_point_obs to see if your data has any wind obs in
it:
>
> plot_point_obs prepbufr.gdas.2015080100.nc
prepbufr.gdas.2015080100.ps -gc
> 34
> ...
> DEBUG 2: Finished plotting 0 locations.
> DEBUG 2: Skipped 0 locations off the grid.
>
> And it does not... thus no matches.  Since you're using GDAS point
> observations, I likely has something to do with the quality marker
you're
> using.  Please look at the NOTE on the bottom of this page about
GDAS
> quality markers for surface observations:
>    http://www.dtcenter.org/met/users/downloads/observation_data.php
>
> Thanks,
> John
>
> On Thu, Mar 10, 2016 at 1:49 PM, Shannon Rees - NOAA Affiliate via
RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75487 >
> >
> > And here is the description of how to use it:
> >
> >
> > On Thu, Mar 10, 2016 at 3:40 PM, Shannon Rees - NOAA Affiliate <
> > shannon.rees at noaa.gov> wrote:
> >
> > > Yes it is attached here.
> > >
> > > On Thu, Mar 10, 2016 at 2:27 PM, John Halley Gotway via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > >> Shannon,
> > >>
> > >> Looking in your config file, I see the following:
> > >>
> > >> grib_fields = [
> > >>   {
> > >>     name  = "${obs_var}";
> > >>     level = "${obs_lev}";
> > >>     cat_thresh = [  ];
> > >>   }
> > >> ];
> > >>
> > >> obs = {
> > >>    message_type = [ "${mess_type}" ];
> > >>    sid_exc      = [];
> > >>    field        = grib_fields;
> > >> }
> > >>
> > >> Do you have a script to wrap these calls to Point-Stat?  I'm
wondering
> > how
> > >> ${obs_var} is getting set.
> > >>
> > >> John
> > >>
> > >>
> > >> On Thu, Mar 10, 2016 at 12:21 PM, Shannon Rees - NOAA Affiliate
via
> RT <
> > >> met_help at ucar.edu> wrote:
> > >>
> > >> >
> > >> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75487
>
> > >> >
> > >> > Here is the configuration file.  I am also uploading some obs
and
> fcst
> > >> data
> > >> > to the ftp site in the "REES" directory.
> > >> >
> > >> > Thank you,
> > >> > Shannon
> > >> >
> > >> >
> > >> > On Thu, Mar 10, 2016 at 1:42 PM, met_help at ucar.edu via RT <
> > >> > met_help at ucar.edu
> > >> > > wrote:
> > >> >
> > >> > > Greetings,
> > >> > >
> > >> > > This message has been automatically generated in response
to the
> > >> > > creation of a trouble ticket regarding:
> > >> > >         "Trouble with Point-stat tool for some variables",
> > >> > > a summary of which appears below.
> > >> > >
> > >> > > There is no need to reply to this message right now.  Your
ticket
> > has
> > >> > been
> > >> > > assigned an ID of [rt.rap.ucar.edu #75487].
> > >> > >
> > >> > > Please include the string:
> > >> > >
> > >> > >          [rt.rap.ucar.edu #75487]
> > >> > >
> > >> > > in the subject line of all future correspondence about this
issue.
> > To
> > >> do
> > >> > > so,
> > >> > > you may reply to this message.
> > >> > >
> > >> > >                         Thank you,
> > >> > >                         met_help at ucar.edu
> > >> > >
> > >> > >
> > >>
> >
-------------------------------------------------------------------------
> > >> > > I am having a problem using the point-stat tool with some
> different
> > >> > > variables in the Prepbufr format.  For example the variable
"vs"
> > >> > (forecast)
> > >> > > for 10m meridional wind I believe should compare with VGRD
level
> Z10
> > >> in
> > >> > > both the Prepbufr point observations and the climatology.
Yet
> when
> > I
> > >> run
> > >> > > the tool it doesn't process any matches, although it does
find a
> > >> forecast
> > >> > > and climatology level.  I can run the grid-stat tool using
the
> same
> > >> > > forecast and climatology fields, but a gridded analysis
instead of
> > >> point
> > >> > > observations, and it works fine.  I get the same problem
with 200
> mb
> > >> > level
> > >> > > winds (P200), although I'm almost positive I have been able
to get
> > >> > > statistics such as RMSE with these variables in the point-
stat
> tool
> > a
> > >> > while
> > >> > > ago.  When I run the tool I get this:
> > >> > >
> > >> > > DEBUG 3: Number of matched pairs  = 0
> > >> > > DEBUG 3: Observations processed   = 45684
> > >> > > DEBUG 3: Rejected: GRIB code      = 45684
> > >> > > I must have the wrong variable name.  I thought for zonal
winds at
> > >> 200 mb
> > >> > > would have the variable name "UGRD" and the level is "P200"
from
> > this
> > >> > > website
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html.
> > Do
> > >> you
> > >> > > know if this is incorrect or where I could get more info?
> > >> > >
> > >> > >
> > >> >
> > >> >
> > >> > --
> > >> > Shannon Rees
> > >> > Programming Scientist
> > >> > Engility
> > >> > Geophysical Fluid Dynamics Lab
> > >> > 201 Forrestal Rd Princeton, NJ
> > >> > 609-452-5384
> > >> >
> > >> >
> > >>
> > >>
> > >
> > >
> > > --
> > > Shannon Rees
> > > Programming Scientist
> > > Engility
> > > Geophysical Fluid Dynamics Lab
> > > 201 Forrestal Rd Princeton, NJ
> > > 609-452-5384
> > >
> >
> >
> >
> > --
> > Shannon Rees
> > Programming Scientist
> > Engility
> > Geophysical Fluid Dynamics Lab
> > 201 Forrestal Rd Princeton, NJ
> > 609-452-5384
> >
> >
>
>


--
Shannon Rees
Programming Scientist
Engility
Geophysical Fluid Dynamics Lab
201 Forrestal Rd Princeton, NJ
609-452-5384

------------------------------------------------
Subject: Trouble with Point-stat tool for some variables
From: John Halley Gotway
Time: Mon Mar 14 15:16:06 2016

Shannon,

Sounds good.  I'll go ahead and resolve this ticket now, but just let
us
know if more issues arise.

Thanks,
John

On Mon, Mar 14, 2016 at 2:05 PM, Shannon Rees - NOAA Affiliate via RT
<
met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75487 >
>
> Hi John,
>
> I think I downloaded all of the GDAS data with quality marker = 9
which
> should preserve the surface variables.  I am downloading some more
GDAS
> data with quality marker = 15 to test it and I'll let you know how
that
> goes.
>
> Thanks,
>
> On Fri, Mar 11, 2016 at 4:31 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > Shannon,
> >
> > I ran with your data and was able to reproduce the behavior you're
> seeing:
> >
> > DEBUG 3: Observations processed   = 46351
> > DEBUG 3: Rejected: SID exclusion  = 0
> > DEBUG 3: Rejected: GRIB code      = 46351
> >
> > Next, I ran plot_point_obs to see if your data has any wind obs in
it:
> >
> > plot_point_obs prepbufr.gdas.2015080100.nc
prepbufr.gdas.2015080100.ps
> -gc
> > 34
> > ...
> > DEBUG 2: Finished plotting 0 locations.
> > DEBUG 2: Skipped 0 locations off the grid.
> >
> > And it does not... thus no matches.  Since you're using GDAS point
> > observations, I likely has something to do with the quality marker
you're
> > using.  Please look at the NOTE on the bottom of this page about
GDAS
> > quality markers for surface observations:
> >
http://www.dtcenter.org/met/users/downloads/observation_data.php
> >
> > Thanks,
> > John
> >
> > On Thu, Mar 10, 2016 at 1:49 PM, Shannon Rees - NOAA Affiliate via
RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75487 >
> > >
> > > And here is the description of how to use it:
> > >
> > >
> > > On Thu, Mar 10, 2016 at 3:40 PM, Shannon Rees - NOAA Affiliate <
> > > shannon.rees at noaa.gov> wrote:
> > >
> > > > Yes it is attached here.
> > > >
> > > > On Thu, Mar 10, 2016 at 2:27 PM, John Halley Gotway via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > >> Shannon,
> > > >>
> > > >> Looking in your config file, I see the following:
> > > >>
> > > >> grib_fields = [
> > > >>   {
> > > >>     name  = "${obs_var}";
> > > >>     level = "${obs_lev}";
> > > >>     cat_thresh = [  ];
> > > >>   }
> > > >> ];
> > > >>
> > > >> obs = {
> > > >>    message_type = [ "${mess_type}" ];
> > > >>    sid_exc      = [];
> > > >>    field        = grib_fields;
> > > >> }
> > > >>
> > > >> Do you have a script to wrap these calls to Point-Stat?  I'm
> wondering
> > > how
> > > >> ${obs_var} is getting set.
> > > >>
> > > >> John
> > > >>
> > > >>
> > > >> On Thu, Mar 10, 2016 at 12:21 PM, Shannon Rees - NOAA
Affiliate via
> > RT <
> > > >> met_help at ucar.edu> wrote:
> > > >>
> > > >> >
> > > >> > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=75487 >
> > > >> >
> > > >> > Here is the configuration file.  I am also uploading some
obs and
> > fcst
> > > >> data
> > > >> > to the ftp site in the "REES" directory.
> > > >> >
> > > >> > Thank you,
> > > >> > Shannon
> > > >> >
> > > >> >
> > > >> > On Thu, Mar 10, 2016 at 1:42 PM, met_help at ucar.edu via RT <
> > > >> > met_help at ucar.edu
> > > >> > > wrote:
> > > >> >
> > > >> > > Greetings,
> > > >> > >
> > > >> > > This message has been automatically generated in response
to the
> > > >> > > creation of a trouble ticket regarding:
> > > >> > >         "Trouble with Point-stat tool for some
variables",
> > > >> > > a summary of which appears below.
> > > >> > >
> > > >> > > There is no need to reply to this message right now.
Your
> ticket
> > > has
> > > >> > been
> > > >> > > assigned an ID of [rt.rap.ucar.edu #75487].
> > > >> > >
> > > >> > > Please include the string:
> > > >> > >
> > > >> > >          [rt.rap.ucar.edu #75487]
> > > >> > >
> > > >> > > in the subject line of all future correspondence about
this
> issue.
> > > To
> > > >> do
> > > >> > > so,
> > > >> > > you may reply to this message.
> > > >> > >
> > > >> > >                         Thank you,
> > > >> > >                         met_help at ucar.edu
> > > >> > >
> > > >> > >
> > > >>
> > >
>
-------------------------------------------------------------------------
> > > >> > > I am having a problem using the point-stat tool with some
> > different
> > > >> > > variables in the Prepbufr format.  For example the
variable "vs"
> > > >> > (forecast)
> > > >> > > for 10m meridional wind I believe should compare with
VGRD level
> > Z10
> > > >> in
> > > >> > > both the Prepbufr point observations and the climatology.
Yet
> > when
> > > I
> > > >> run
> > > >> > > the tool it doesn't process any matches, although it does
find a
> > > >> forecast
> > > >> > > and climatology level.  I can run the grid-stat tool
using the
> > same
> > > >> > > forecast and climatology fields, but a gridded analysis
instead
> of
> > > >> point
> > > >> > > observations, and it works fine.  I get the same problem
with
> 200
> > mb
> > > >> > level
> > > >> > > winds (P200), although I'm almost positive I have been
able to
> get
> > > >> > > statistics such as RMSE with these variables in the
point-stat
> > tool
> > > a
> > > >> > while
> > > >> > > ago.  When I run the tool I get this:
> > > >> > >
> > > >> > > DEBUG 3: Number of matched pairs  = 0
> > > >> > > DEBUG 3: Observations processed   = 45684
> > > >> > > DEBUG 3: Rejected: GRIB code      = 45684
> > > >> > > I must have the wrong variable name.  I thought for zonal
winds
> at
> > > >> 200 mb
> > > >> > > would have the variable name "UGRD" and the level is
"P200" from
> > > this
> > > >> > > website
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
> .
> > > Do
> > > >> you
> > > >> > > know if this is incorrect or where I could get more info?
> > > >> > >
> > > >> > >
> > > >> >
> > > >> >
> > > >> > --
> > > >> > Shannon Rees
> > > >> > Programming Scientist
> > > >> > Engility
> > > >> > Geophysical Fluid Dynamics Lab
> > > >> > 201 Forrestal Rd Princeton, NJ
> > > >> > 609-452-5384
> > > >> >
> > > >> >
> > > >>
> > > >>
> > > >
> > > >
> > > > --
> > > > Shannon Rees
> > > > Programming Scientist
> > > > Engility
> > > > Geophysical Fluid Dynamics Lab
> > > > 201 Forrestal Rd Princeton, NJ
> > > > 609-452-5384
> > > >
> > >
> > >
> > >
> > > --
> > > Shannon Rees
> > > Programming Scientist
> > > Engility
> > > Geophysical Fluid Dynamics Lab
> > > 201 Forrestal Rd Princeton, NJ
> > > 609-452-5384
> > >
> > >
> >
> >
>
>
> --
> Shannon Rees
> Programming Scientist
> Engility
> Geophysical Fluid Dynamics Lab
> 201 Forrestal Rd Princeton, NJ
> 609-452-5384
>
>

------------------------------------------------


More information about the Met_help mailing list