[Met_help] [rt.rap.ucar.edu #88041] History for UGRD and VGRD not in prebufr
John Halley Gotway via RT
met_help at ucar.edu
Tue Jul 9 12:07:14 MDT 2019
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
I am comparing model against observations retrieved from rda.ucar.edu <http://rda.ucar.edu/> I was able to extract TMP/L1, RH/L2, SPFH/L2 by point_stat program but I was not able to do the same with surface UGRD and VGRD. I tried to use level: L0, L1, L2, Z10, P1000, P1014 and P1013.
The UGRD was found in the forecast file but not in the observation file. The variables are
How can I know the leves in the prebufr file assigned to V and U variables?
Agustín García-Reynoso http://www.atmosfera.unam.mx/fqa
Centro de Ciencias de la Atmosfera http://www.atmosfera.unam.mx
52-55-56 22 43 96
Universidad Nacional Autonoma de Mexico 52-55-56 16 07 89 (fax)
Circuito Exterior s/n Mexico DF 04510
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: UGRD and VGRD not in prebufr
From: John Halley Gotway
Time: Wed Dec 05 14:40:55 2018
I see that you have a question about retrieving observations of type
UGRD
and VGRD from PREPBUFR files. For the surface, U and V winds are
typically
reported at 10 meters (i.e. level = "Z10";) and are stored in the
ADPSFC
message type.
However, there are some wrinkles depending on the PREPBUFR data
source.
For example, in GDAS (Global Data Assimilation), surface variables are
often assigned a quality marker value of 9... only to prevent them
from
being included in the data assimilation system. This is mentioned as
a
NOTE on the bottom of this page:
https://dtcenter.org/met/users/downloads/observation_data.php
If you're using GDAS, perhaps you're inadvertently filtering the
surface
obs out depending on how you have the quality_mark_thresh set in the
PB2NC
config file? You could try increasing that threshold to 9 to see if
you
get 10-meter winds:
quality_mark_thresh = 9;
If you have a more specific question, please either send me an example
data
file or point me to the exact one you're using online. Also please
send me
the MET configuration file and the command line you're running. If I
can
replicate the behavior you described, I can probably figure out why
it's
happening. Here's how to send us data:
https://dtcenter.org/met/users/support/met_help.php#ftp
Also, you may find it useful to run the PB2NC tool a single time using
the
"-index" command line option. It's very, very slow... but it'll dump
out
an inventory of the observation variables that are present in that
file
with a list of the message types in which they appear. You really
only
need to run this once for each data source to see what it contains and
what
variables to request in the PB2NC config file.
Hope that helps.
Thanks,
John Halley Gotway
------------------------------------------------
Subject: UGRD and VGRD not in prebufr
From: Jose Agustín García Reynoso
Time: Wed Dec 05 17:46:28 2018
John
I used the Dataset : ds337.0 and dowloaded the files in netcdf format
from
the site. Then I am not able to use the pb2nc in order to observe the
headers.
The quality_mark_thresh = 9; was included in the config file and the
same
result.
The met version is 8.0
I uploaded the files in the ftp site *ftp.rap.ucar.edu
<http://ftp.rap.ucar.edu> * in *incoming/irap/met_help/data_garcia *
directory.
El mié., 5 de dic. de 2018 a la(s) 15:40, John Halley Gotway via RT (
met_help at ucar.edu) escribió:
> I see that you have a question about retrieving observations of type
UGRD
> and VGRD from PREPBUFR files. For the surface, U and V winds are
typically
> reported at 10 meters (i.e. level = "Z10";) and are stored in the
ADPSFC
> message type.
>
> However, there are some wrinkles depending on the PREPBUFR data
source.
> For example, in GDAS (Global Data Assimilation), surface variables
are
> often assigned a quality marker value of 9... only to prevent them
from
> being included in the data assimilation system. This is mentioned
as a
> NOTE on the bottom of this page:
> https://dtcenter.org/met/users/downloads/observation_data.php
>
> If you're using GDAS, perhaps you're inadvertently filtering the
surface
> obs out depending on how you have the quality_mark_thresh set in the
PB2NC
> config file? You could try increasing that threshold to 9 to see if
you
> get 10-meter winds:
> quality_mark_thresh = 9;
>
> If you have a more specific question, please either send me an
example data
> file or point me to the exact one you're using online. Also please
send me
> the MET configuration file and the command line you're running. If
I can
> replicate the behavior you described, I can probably figure out why
it's
> happening. Here's how to send us data:
> https://dtcenter.org/met/users/support/met_help.php#ftp
>
> Also, you may find it useful to run the PB2NC tool a single time
using the
> "-index" command line option. It's very, very slow... but it'll
dump out
> an inventory of the observation variables that are present in that
file
> with a list of the message types in which they appear. You really
only
> need to run this once for each data source to see what it contains
and what
> variables to request in the PB2NC config file.
>
> Hope that helps.
>
> Thanks,
> John Halley Gotway
>
>
--
Dr. Jose Agustin Garcia Reynoso
Atmosfera <http://www.atmosfera.unam.mx/fqa>
*Centro de Ciencias de la Atmosfera, *UNAM
Tel (52) 55 *56 22 43 96* Fax (52) 55 56160789
Circuito de la Investigacion Cientifica s/n Mexico DF 04510
------------------------------------------------
Subject: UGRD and VGRD not in prebufr
From: Jose Agustín García Reynoso
Time: Wed Dec 05 18:09:34 2018
John
I re-dowloaded and I have U and V data by using:
Date Limits : 2018-10-20 06:00 2018-10-30 06:00
Latitude Limits : 18.0 N 32.0 N
Longitude Limits : 100.0 W 85.0 W
Parameters : SPFH TMP HGT UGRD VGRD DPT WIND RH MIXR
PRMSL
PREPBUFR Message Types : SFCSHP ADPSFC
PREPBUFR Report Types : 111 183 227 250 280 120 182 229 251 281
126 187 230 252 282 130 210 232
253 284 132 220 242 254 287 153 221 243 257 289 180 223 245 258 290
181 224 246 259
Input Report Types : 11 65 75 531 574 12 66 76 532 575 13 68 77
540 576 22 69 511 551 577 23
71 512 561 581 31 72 514 562 582 41 73 522 571 583 61 74 523 573 584
63
Quality Mark Threshold : 9
File Compression : gz
Data Output Format : NETCDF
El mié., 5 de dic. de 2018 a la(s) 15:40, John Halley Gotway via RT (
met_help at ucar.edu) escribió:
> I see that you have a question about retrieving observations of type
UGRD
> and VGRD from PREPBUFR files. For the surface, U and V winds are
typically
> reported at 10 meters (i.e. level = "Z10";) and are stored in the
ADPSFC
> message type.
>
> However, there are some wrinkles depending on the PREPBUFR data
source.
> For example, in GDAS (Global Data Assimilation), surface variables
are
> often assigned a quality marker value of 9... only to prevent them
from
> being included in the data assimilation system. This is mentioned
as a
> NOTE on the bottom of this page:
> https://dtcenter.org/met/users/downloads/observation_data.php
>
> If you're using GDAS, perhaps you're inadvertently filtering the
surface
> obs out depending on how you have the quality_mark_thresh set in the
PB2NC
> config file? You could try increasing that threshold to 9 to see if
you
> get 10-meter winds:
> quality_mark_thresh = 9;
>
> If you have a more specific question, please either send me an
example data
> file or point me to the exact one you're using online. Also please
send me
> the MET configuration file and the command line you're running. If
I can
> replicate the behavior you described, I can probably figure out why
it's
> happening. Here's how to send us data:
> https://dtcenter.org/met/users/support/met_help.php#ftp
>
> Also, you may find it useful to run the PB2NC tool a single time
using the
> "-index" command line option. It's very, very slow... but it'll
dump out
> an inventory of the observation variables that are present in that
file
> with a list of the message types in which they appear. You really
only
> need to run this once for each data source to see what it contains
and what
> variables to request in the PB2NC config file.
>
> Hope that helps.
>
> Thanks,
> John Halley Gotway
>
>
--
Dr. Jose Agustin Garcia Reynoso
Atmosfera <http://www.atmosfera.unam.mx/fqa>
*Centro de Ciencias de la Atmosfera, *UNAM
Tel (52) 55 *56 22 43 96* Fax (52) 55 56160789
Circuito de la Investigacion Cientifica s/n Mexico DF 04510
------------------------------------------------
Subject: UGRD and VGRD not in prebufr
From: John Halley Gotway
Time: Thu Dec 06 09:19:34 2018
Jose,
OK, so are you all set, or do you still have any questions about
verifying
10-meter winds?
Thanks,
John
On Wed, Dec 5, 2018 at 6:10 PM Jose Agustín García Reynoso via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88041 >
>
> John
>
> I re-dowloaded and I have U and V data by using:
>
> Date Limits : 2018-10-20 06:00 2018-10-30 06:00
> Latitude Limits : 18.0 N 32.0 N
> Longitude Limits : 100.0 W 85.0 W
> Parameters : SPFH TMP HGT UGRD VGRD DPT WIND RH MIXR
PRMSL
> PREPBUFR Message Types : SFCSHP ADPSFC
> PREPBUFR Report Types : 111 183 227 250 280 120 182 229 251 281
> 126 187 230 252 282 130 210 232
> 253 284 132 220 242 254 287 153 221 243 257 289 180 223 245 258 290
> 181 224 246 259
> Input Report Types : 11 65 75 531 574 12 66 76 532 575 13 68
77
> 540 576 22 69 511 551 577 23
> 71 512 561 581 31 72 514 562 582 41 73 522 571 583 61 74 523 573 584
63
> Quality Mark Threshold : 9
> File Compression : gz
> Data Output Format : NETCDF
>
>
> El mié., 5 de dic. de 2018 a la(s) 15:40, John Halley Gotway via RT
(
> met_help at ucar.edu) escribió:
>
> > I see that you have a question about retrieving observations of
type UGRD
> > and VGRD from PREPBUFR files. For the surface, U and V winds are
> typically
> > reported at 10 meters (i.e. level = "Z10";) and are stored in the
ADPSFC
> > message type.
> >
> > However, there are some wrinkles depending on the PREPBUFR data
source.
> > For example, in GDAS (Global Data Assimilation), surface variables
are
> > often assigned a quality marker value of 9... only to prevent them
from
> > being included in the data assimilation system. This is mentioned
as a
> > NOTE on the bottom of this page:
> > https://dtcenter.org/met/users/downloads/observation_data.php
> >
> > If you're using GDAS, perhaps you're inadvertently filtering the
surface
> > obs out depending on how you have the quality_mark_thresh set in
the
> PB2NC
> > config file? You could try increasing that threshold to 9 to see
if you
> > get 10-meter winds:
> > quality_mark_thresh = 9;
> >
> > If you have a more specific question, please either send me an
example
> data
> > file or point me to the exact one you're using online. Also
please send
> me
> > the MET configuration file and the command line you're running.
If I can
> > replicate the behavior you described, I can probably figure out
why it's
> > happening. Here's how to send us data:
> > https://dtcenter.org/met/users/support/met_help.php#ftp
> >
> > Also, you may find it useful to run the PB2NC tool a single time
using
> the
> > "-index" command line option. It's very, very slow... but it'll
dump out
> > an inventory of the observation variables that are present in that
file
> > with a list of the message types in which they appear. You really
only
> > need to run this once for each data source to see what it contains
and
> what
> > variables to request in the PB2NC config file.
> >
> > Hope that helps.
> >
> > Thanks,
> > John Halley Gotway
> >
> >
>
> --
> Dr. Jose Agustin Garcia Reynoso
> Atmosfera <http://www.atmosfera.unam.mx/fqa>
> *Centro de Ciencias de la Atmosfera, *UNAM
> Tel (52) 55 *56 22 43 96* Fax (52) 55 56160789
> Circuito de la Investigacion Cientifica s/n Mexico DF 04510
>
>
------------------------------------------------
Subject: UGRD and VGRD not in prebufr
From: Jose Agustín García Reynoso
Time: Thu Dec 06 09:41:07 2018
John
I used UGRD/L1, VGRD/L1 and WIND/L1 in the configuration file
> El 06/12/2018, a las 10:19, John Halley Gotway via RT
<met_help at ucar.edu> escribió:
>
> Jose,
>
> OK, so are you all set, or do you still have any questions about
verifying
> 10-meter winds?
>
> Thanks,
> John
>
> On Wed, Dec 5, 2018 at 6:10 PM Jose Agustín García Reynoso via RT <
> met_help at ucar.edu <mailto:met_help at ucar.edu>> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88041
<https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88041> >
>>
>> John
>>
>> I re-dowloaded and I have U and V data by using:
>>
>> Date Limits : 2018-10-20 06:00 2018-10-30 06:00
>> Latitude Limits : 18.0 N 32.0 N
>> Longitude Limits : 100.0 W 85.0 W
>> Parameters : SPFH TMP HGT UGRD VGRD DPT WIND RH MIXR
PRMSL
>> PREPBUFR Message Types : SFCSHP ADPSFC
>> PREPBUFR Report Types : 111 183 227 250 280 120 182 229 251 281
>> 126 187 230 252 282 130 210 232
>> 253 284 132 220 242 254 287 153 221 243 257 289 180 223 245 258 290
>> 181 224 246 259
>> Input Report Types : 11 65 75 531 574 12 66 76 532 575 13 68
77
>> 540 576 22 69 511 551 577 23
>> 71 512 561 581 31 72 514 562 582 41 73 522 571 583 61 74 523 573
584 63
>> Quality Mark Threshold : 9
>> File Compression : gz
>> Data Output Format : NETCDF
>>
>>
>> El mié., 5 de dic. de 2018 a la(s) 15:40, John Halley Gotway via RT
(
>> met_help at ucar.edu) escribió:
>>
>>> I see that you have a question about retrieving observations of
type UGRD
>>> and VGRD from PREPBUFR files. For the surface, U and V winds are
>> typically
>>> reported at 10 meters (i.e. level = "Z10";) and are stored in the
ADPSFC
>>> message type.
>>>
>>> However, there are some wrinkles depending on the PREPBUFR data
source.
>>> For example, in GDAS (Global Data Assimilation), surface variables
are
>>> often assigned a quality marker value of 9... only to prevent them
from
>>> being included in the data assimilation system. This is mentioned
as a
>>> NOTE on the bottom of this page:
>>> https://dtcenter.org/met/users/downloads/observation_data.php
>>>
>>> If you're using GDAS, perhaps you're inadvertently filtering the
surface
>>> obs out depending on how you have the quality_mark_thresh set in
the
>> PB2NC
>>> config file? You could try increasing that threshold to 9 to see
if you
>>> get 10-meter winds:
>>> quality_mark_thresh = 9;
>>>
>>> If you have a more specific question, please either send me an
example
>> data
>>> file or point me to the exact one you're using online. Also
please send
>> me
>>> the MET configuration file and the command line you're running.
If I can
>>> replicate the behavior you described, I can probably figure out
why it's
>>> happening. Here's how to send us data:
>>> https://dtcenter.org/met/users/support/met_help.php#ftp
>>>
>>> Also, you may find it useful to run the PB2NC tool a single time
using
>> the
>>> "-index" command line option. It's very, very slow... but it'll
dump out
>>> an inventory of the observation variables that are present in that
file
>>> with a list of the message types in which they appear. You really
only
>>> need to run this once for each data source to see what it contains
and
>> what
>>> variables to request in the PB2NC config file.
>>>
>>> Hope that helps.
>>>
>>> Thanks,
>>> John Halley Gotway
>>>
>>>
>>
>> --
>> Dr. Jose Agustin Garcia Reynoso
>> Atmosfera <http://www.atmosfera.unam.mx/fqa
<http://www.atmosfera.unam.mx/fqa>>
>> *Centro de Ciencias de la Atmosfera, *UNAM
>> Tel (52) 55 *56 22 43 96* Fax (52) 55 56160789
>> Circuito de la Investigacion Cientifica s/n Mexico DF 04510
------------------------------------------------
Subject: UGRD and VGRD not in prebufr
From: John Halley Gotway
Time: Thu Dec 06 10:08:27 2018
Jose,
If you're running Point-Stat to compare GRIB forecast files to point
observations, I would recommend increasing the logging level to "-v
3".
That'll include in the log output the GRIB records that are being
matched
by your level setting of "L1". And then use wgrib to inventory the
contents of the GRIB file and make sure Point-Stat is matching the
GRIB
records you expect.
That way you can be confident that it's evaluating the data you
expect.
Thanks,
John
On Thu, Dec 6, 2018 at 9:41 AM Jose Agustín García Reynoso via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88041 >
>
> John
>
> I used UGRD/L1, VGRD/L1 and WIND/L1 in the configuration file
>
>
> > El 06/12/2018, a las 10:19, John Halley Gotway via RT
<met_help at ucar.edu>
> escribió:
> >
> > Jose,
> >
> > OK, so are you all set, or do you still have any questions about
> verifying
> > 10-meter winds?
> >
> > Thanks,
> > John
> >
> > On Wed, Dec 5, 2018 at 6:10 PM Jose Agustín García Reynoso via RT
<
> > met_help at ucar.edu <mailto:met_help at ucar.edu>> wrote:
> >
> >>
> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88041 <
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88041> >
> >>
> >> John
> >>
> >> I re-dowloaded and I have U and V data by using:
> >>
> >> Date Limits : 2018-10-20 06:00 2018-10-30 06:00
> >> Latitude Limits : 18.0 N 32.0 N
> >> Longitude Limits : 100.0 W 85.0 W
> >> Parameters : SPFH TMP HGT UGRD VGRD DPT WIND RH
MIXR
> PRMSL
> >> PREPBUFR Message Types : SFCSHP ADPSFC
> >> PREPBUFR Report Types : 111 183 227 250 280 120 182 229 251
281
> >> 126 187 230 252 282 130 210 232
> >> 253 284 132 220 242 254 287 153 221 243 257 289 180 223 245 258
290
> >> 181 224 246 259
> >> Input Report Types : 11 65 75 531 574 12 66 76 532 575 13
68 77
> >> 540 576 22 69 511 551 577 23
> >> 71 512 561 581 31 72 514 562 582 41 73 522 571 583 61 74 523 573
584 63
> >> Quality Mark Threshold : 9
> >> File Compression : gz
> >> Data Output Format : NETCDF
> >>
> >>
> >> El mié., 5 de dic. de 2018 a la(s) 15:40, John Halley Gotway via
RT (
> >> met_help at ucar.edu) escribió:
> >>
> >>> I see that you have a question about retrieving observations of
type
> UGRD
> >>> and VGRD from PREPBUFR files. For the surface, U and V winds
are
> >> typically
> >>> reported at 10 meters (i.e. level = "Z10";) and are stored in
the
> ADPSFC
> >>> message type.
> >>>
> >>> However, there are some wrinkles depending on the PREPBUFR data
source.
> >>> For example, in GDAS (Global Data Assimilation), surface
variables are
> >>> often assigned a quality marker value of 9... only to prevent
them from
> >>> being included in the data assimilation system. This is
mentioned as a
> >>> NOTE on the bottom of this page:
> >>> https://dtcenter.org/met/users/downloads/observation_data.php
> >>>
> >>> If you're using GDAS, perhaps you're inadvertently filtering the
> surface
> >>> obs out depending on how you have the quality_mark_thresh set in
the
> >> PB2NC
> >>> config file? You could try increasing that threshold to 9 to
see if
> you
> >>> get 10-meter winds:
> >>> quality_mark_thresh = 9;
> >>>
> >>> If you have a more specific question, please either send me an
example
> >> data
> >>> file or point me to the exact one you're using online. Also
please
> send
> >> me
> >>> the MET configuration file and the command line you're running.
If I
> can
> >>> replicate the behavior you described, I can probably figure out
why
> it's
> >>> happening. Here's how to send us data:
> >>> https://dtcenter.org/met/users/support/met_help.php#ftp
> >>>
> >>> Also, you may find it useful to run the PB2NC tool a single time
using
> >> the
> >>> "-index" command line option. It's very, very slow... but it'll
dump
> out
> >>> an inventory of the observation variables that are present in
that file
> >>> with a list of the message types in which they appear. You
really only
> >>> need to run this once for each data source to see what it
contains and
> >> what
> >>> variables to request in the PB2NC config file.
> >>>
> >>> Hope that helps.
> >>>
> >>> Thanks,
> >>> John Halley Gotway
> >>>
> >>>
> >>
> >> --
> >> Dr. Jose Agustin Garcia Reynoso
> >> Atmosfera <http://www.atmosfera.unam.mx/fqa <
> http://www.atmosfera.unam.mx/fqa>>
> >> *Centro de Ciencias de la Atmosfera, *UNAM
> >> Tel (52) 55 *56 22 43 96* Fax (52) 55 56160789
> >> Circuito de la Investigacion Cientifica s/n Mexico DF 04510
>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #88041] UGRD and VGRD not in prebufr
From: Jose Agustín García Reynoso
Time: Thu Dec 06 10:15:40 2018
John
Thanks I did it up to -v 5 and that let me know that the observation
file has no useful data for UGRD and VGRN because all the rejected
lines were zero.
now:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Processing VGRD/L1 versus VGRD/L1, for observation type
SFCSHP, over region GLF, for interpolation method BILIN(4), using 160
pairs.
DEBUG 3: Number of matched pairs = 160
DEBUG 3: Observations processed = 43411
DEBUG 3: Rejected: SID exclusion = 0
DEBUG 3: Rejected: obs type = 39589
DEBUG 3: Rejected: valid time = 2437
DEBUG 3: Rejected: bad obs value = 0
DEBUG 3: Rejected: off the grid = 327
DEBUG 3: Rejected: level mismatch = 0
DEBUG 3: Rejected: quality marker = 0
DEBUG 3: Rejected: message type = 547
DEBUG 3: Rejected: masking region = 351
DEBUG 3: Rejected: bad fcst value = 0
DEBUG 3: Rejected: duplicates = 0
DEBUG 2: Computing Categorical Statistics.
DEBUG 2: Computing Continuous Statistics.
DEBUG 2: Computing Scalar Partial Sums.
there are some used and others rejected.
Thanks
> El 06/12/2018, a las 11:08, John Halley Gotway via RT
<met_help at ucar.edu> escribió:
>
> Jose,
>
> If you're running Point-Stat to compare GRIB forecast files to point
> observations, I would recommend increasing the logging level to "-v
3".
> That'll include in the log output the GRIB records that are being
matched
> by your level setting of "L1". And then use wgrib to inventory the
> contents of the GRIB file and make sure Point-Stat is matching the
GRIB
> records you expect.
>
> That way you can be confident that it's evaluating the data you
expect.
>
> Thanks,
> John
>
> On Thu, Dec 6, 2018 at 9:41 AM Jose Agustín García Reynoso via RT <
> met_help at ucar.edu> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88041 >
>>
>> John
>>
>> I used UGRD/L1, VGRD/L1 and WIND/L1 in the configuration file
>>
>>
>>> El 06/12/2018, a las 10:19, John Halley Gotway via RT
<met_help at ucar.edu>
>> escribió:
>>>
>>> Jose,
>>>
>>> OK, so are you all set, or do you still have any questions about
>> verifying
>>> 10-meter winds?
>>>
>>> Thanks,
>>> John
>>>
>>> On Wed, Dec 5, 2018 at 6:10 PM Jose Agustín García Reynoso via RT
<
>>> met_help at ucar.edu <mailto:met_help at ucar.edu>> wrote:
>>>
>>>>
>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88041 <
>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88041> >
>>>>
>>>> John
>>>>
>>>> I re-dowloaded and I have U and V data by using:
>>>>
>>>> Date Limits : 2018-10-20 06:00 2018-10-30 06:00
>>>> Latitude Limits : 18.0 N 32.0 N
>>>> Longitude Limits : 100.0 W 85.0 W
>>>> Parameters : SPFH TMP HGT UGRD VGRD DPT WIND RH
MIXR
>> PRMSL
>>>> PREPBUFR Message Types : SFCSHP ADPSFC
>>>> PREPBUFR Report Types : 111 183 227 250 280 120 182 229 251
281
>>>> 126 187 230 252 282 130 210 232
>>>> 253 284 132 220 242 254 287 153 221 243 257 289 180 223 245 258
290
>>>> 181 224 246 259
>>>> Input Report Types : 11 65 75 531 574 12 66 76 532 575 13
68 77
>>>> 540 576 22 69 511 551 577 23
>>>> 71 512 561 581 31 72 514 562 582 41 73 522 571 583 61 74 523 573
584 63
>>>> Quality Mark Threshold : 9
>>>> File Compression : gz
>>>> Data Output Format : NETCDF
>>>>
>>>>
>>>> El mié., 5 de dic. de 2018 a la(s) 15:40, John Halley Gotway via
RT (
>>>> met_help at ucar.edu) escribió:
>>>>
>>>>> I see that you have a question about retrieving observations of
type
>> UGRD
>>>>> and VGRD from PREPBUFR files. For the surface, U and V winds
are
>>>> typically
>>>>> reported at 10 meters (i.e. level = "Z10";) and are stored in
the
>> ADPSFC
>>>>> message type.
>>>>>
>>>>> However, there are some wrinkles depending on the PREPBUFR data
source.
>>>>> For example, in GDAS (Global Data Assimilation), surface
variables are
>>>>> often assigned a quality marker value of 9... only to prevent
them from
>>>>> being included in the data assimilation system. This is
mentioned as a
>>>>> NOTE on the bottom of this page:
>>>>> https://dtcenter.org/met/users/downloads/observation_data.php
>>>>>
>>>>> If you're using GDAS, perhaps you're inadvertently filtering the
>> surface
>>>>> obs out depending on how you have the quality_mark_thresh set in
the
>>>> PB2NC
>>>>> config file? You could try increasing that threshold to 9 to
see if
>> you
>>>>> get 10-meter winds:
>>>>> quality_mark_thresh = 9;
>>>>>
>>>>> If you have a more specific question, please either send me an
example
>>>> data
>>>>> file or point me to the exact one you're using online. Also
please
>> send
>>>> me
>>>>> the MET configuration file and the command line you're running.
If I
>> can
>>>>> replicate the behavior you described, I can probably figure out
why
>> it's
>>>>> happening. Here's how to send us data:
>>>>> https://dtcenter.org/met/users/support/met_help.php#ftp
>>>>>
>>>>> Also, you may find it useful to run the PB2NC tool a single time
using
>>>> the
>>>>> "-index" command line option. It's very, very slow... but it'll
dump
>> out
>>>>> an inventory of the observation variables that are present in
that file
>>>>> with a list of the message types in which they appear. You
really only
>>>>> need to run this once for each data source to see what it
contains and
>>>> what
>>>>> variables to request in the PB2NC config file.
>>>>>
>>>>> Hope that helps.
>>>>>
>>>>> Thanks,
>>>>> John Halley Gotway
>>>>>
>>>>>
>>>>
>>>> --
>>>> Dr. Jose Agustin Garcia Reynoso
>>>> Atmosfera <http://www.atmosfera.unam.mx/fqa <
>> http://www.atmosfera.unam.mx/fqa>>
>>>> *Centro de Ciencias de la Atmosfera, *UNAM
>>>> Tel (52) 55 *56 22 43 96* Fax (52) 55 56160789
>>>> Circuito de la Investigacion Cientifica s/n Mexico DF 04510
>>
>>
>>
>
------------------------------------------------
More information about the Met_help
mailing list