[Met_help] [rt.rap.ucar.edu #38811] History for Wind Thresholds
RAL HelpDesk {for Paul Oldenburg}
met_help at ucar.edu
Tue Jun 29 11:23:58 MDT 2010
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
I know you guys are probably sick of hearing from me but I have another
question. Now that I am finally running point stat with the updated version
of MET, I am not getting data for all the wind thresholds I want. Attached
is the config file I am using. When I had the threshold set to gt7.7 and
higher, I only got data for gt7.7. Then I added the NA (no threshold case
to the list) and now I only have data for NA case.
So, what am I doing wrong in the config file?
thanks
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
From: Paul Oldenburg
Time: Fri Jun 25 09:39:41 2010
Robert,
Can you please send your model output and observation file?
Thanks,
Paul
RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
> Fri Jun 25 09:37:58 2010: Request 38811 was acted upon.
> Transaction: Ticket created by craigr at offutt.af.mil
> Queue: met_help
> Subject: Wind Thresholds
> Owner: Nobody
> Requestors: craigr at offutt.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>
>
> I know you guys are probably sick of hearing from me but I have
another
> question. Now that I am finally running point stat with the updated
version
> of MET, I am not getting data for all the wind thresholds I want.
Attached
> is the config file I am using. When I had the threshold set to
gt7.7 and
> higher, I only got data for gt7.7. Then I added the NA (no
threshold case
> to the list) and now I only have data for NA case.
>
>
>
> So, what am I doing wrong in the config file?
>
>
>
> thanks
>
>
------------------------------------------------
Subject: Wind Thresholds
From: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Time: Fri Jun 25 10:14:18 2010
The attached WRF file is valid at 24 Jun 6Z. The output is listed
below.
/sphome/qcteam/METv2.0/bin/point_stat
/gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
/h/data/global/WXQC/da
ta/obs_prep/2010062406
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF -outdir
/h/data/global/WXQC/data/met/data -v 2
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=4257686604
Forecast File:
/gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
Climatology File: none
Configuration File:
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
Observation File: /h/data/global/WXQC/data/obs_prep/2010062406
----------------------------------------
Reading records for UGRD/Z10.
For UGRD/Z10 found 1 forecast levels and 0 climatology levels.
----------------------------------------
Reading records for VGRD/Z10.
For VGRD/Z10 found 1 forecast levels and 0 climatology levels.
----------------------------------------
----------------------------------------
Processing UGRD/Z10 versus UGRD/Z10, for observation type ANYSFC, over
region
FULL, for interpolation method MAX(4), using 2522
pairs.
Computing Continuous Statistics.
----------------------------------------
Processing VGRD/Z10 versus VGRD/Z10, for observation type ANYSFC, over
region
FULL, for interpolation method MAX(4), using 2522
pairs.
Computing Continuous Statistics.
Computing Vector Partial Sums, for forecast wind speed NA, for
observation
wind speed NA, using 2522 pairs.
----------------------------------------
Output file:
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V.stat
Output file:
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_ctc.txt
Output file:
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cts.txt
Output file:
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cnt.txt
Output file:
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_sl1l2.txt
Output file:
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_vl1l2.txt
Thanks
Bob
-----Original Message-----
From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
Sent: Friday, June 25, 2010 10:40 AM
To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
Robert,
Can you please send your model output and observation file?
Thanks,
Paul
RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
> Fri Jun 25 09:37:58 2010: Request 38811 was acted upon.
> Transaction: Ticket created by craigr at offutt.af.mil
> Queue: met_help
> Subject: Wind Thresholds
> Owner: Nobody
> Requestors: craigr at offutt.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>
>
> I know you guys are probably sick of hearing from me but I have
another
> question. Now that I am finally running point stat with the updated
version
> of MET, I am not getting data for all the wind thresholds I want.
Attached
> is the config file I am using. When I had the threshold set to
gt7.7 and
> higher, I only got data for gt7.7. Then I added the NA (no
threshold case
> to the list) and now I only have data for NA case.
>
>
>
> So, what am I doing wrong in the config file?
>
>
>
> thanks
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
From: Paul Oldenburg
Time: Fri Jun 25 10:23:19 2010
Bob,
I'm afraid the grib file you sent was corrupted, perhaps because it
was too
large for an email attachment. The obs file seems to be ok. Can you
please
post the grib file to our anonymous ftp site as described below:
ftp ftp.rap.ucar.edu
username = anonymous
password = craigr at offutt.af.mil
cd incoming/irap/met_help/craig_data
put qc_100622_us057g1011t02b060004800
bye
Thanks,
Paul
RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>
> The attached WRF file is valid at 24 Jun 6Z. The output is listed
below.
>
> /sphome/qcteam/METv2.0/bin/point_stat
> /gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
> /h/data/global/WXQC/da
> ta/obs_prep/2010062406
> /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF -outdir
> /h/data/global/WXQC/data/met/data -v 2
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=4257686604
> Forecast File:
> /gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
> Climatology File: none
> Configuration File:
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
> Observation File: /h/data/global/WXQC/data/obs_prep/2010062406
>
> ----------------------------------------
>
> Reading records for UGRD/Z10.
> For UGRD/Z10 found 1 forecast levels and 0 climatology levels.
>
> ----------------------------------------
>
> Reading records for VGRD/Z10.
> For VGRD/Z10 found 1 forecast levels and 0 climatology levels.
>
> ----------------------------------------
> ----------------------------------------
>
> Processing UGRD/Z10 versus UGRD/Z10, for observation type ANYSFC,
over region
> FULL, for interpolation method MAX(4), using 2522
> pairs.
> Computing Continuous Statistics.
>
> ----------------------------------------
>
> Processing VGRD/Z10 versus VGRD/Z10, for observation type ANYSFC,
over region
> FULL, for interpolation method MAX(4), using 2522
> pairs.
> Computing Continuous Statistics.
> Computing Vector Partial Sums, for forecast wind speed NA, for
observation
> wind speed NA, using 2522 pairs.
>
> ----------------------------------------
>
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V.stat
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_ctc.txt
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cts.txt
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cnt.txt
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_sl1l2.txt
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_vl1l2.txt
>
> Thanks
>
> Bob
>
> -----Original Message-----
> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent: Friday, June 25, 2010 10:40 AM
> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
>
> Robert,
>
> Can you please send your model output and observation file?
>
> Thanks,
>
> Paul
>
> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>> Fri Jun 25 09:37:58 2010: Request 38811 was acted upon.
>> Transaction: Ticket created by craigr at offutt.af.mil
>> Queue: met_help
>> Subject: Wind Thresholds
>> Owner: Nobody
>> Requestors: craigr at offutt.af.mil
>> Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>>
>>
>> I know you guys are probably sick of hearing from me but I have
another
>> question. Now that I am finally running point stat with the
updated version
>> of MET, I am not getting data for all the wind thresholds I want.
Attached
>> is the config file I am using. When I had the threshold set to
gt7.7 and
>> higher, I only got data for gt7.7. Then I added the NA (no
threshold case
>> to the list) and now I only have data for NA case.
>>
>>
>>
>> So, what am I doing wrong in the config file?
>>
>>
>>
>> thanks
>>
>>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #38811] Wind Thresholds
From: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Time: Fri Jun 25 11:25:43 2010
I ftped the file to your server.
Bob
-----Original Message-----
From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
Sent: Friday, June 25, 2010 11:23 AM
To: Craig, Robert J Civ USAF
AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
Bob,
I'm afraid the grib file you sent was corrupted, perhaps
because it was too
large for an email attachment. The obs file seems
to be ok. Can you please
post the grib file to our anonymous ftp
site as described below:
ftp ftp.rap.ucar.edu
username = anonymous
password = craigr at offutt.af.mil
cd incoming/irap/met_help/craig_data
put qc_100622_us057g1011t02b060004800
bye
Thanks,
Paul
RAL
HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>
> The
attached WRF file is valid at 24 Jun 6Z. The output is listed below.
>
> /sphome/qcteam/METv2.0/bin/point_stat
>
/gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
> /h/data/global/WXQC/da
> ta/obs_prep/2010062406
>
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF -outdir
>
/h/data/global/WXQC/data/met/data -v 2
> GSL_RNG_TYPE=mt19937
>
GSL_RNG_SEED=4257686604
> Forecast File:
>
/gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
>
Climatology File: none
> Configuration File:
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>
Observation File: /h/data/global/WXQC/data/obs_prep/2010062406
>
>
----------------------------------------
>
> Reading records for
UGRD/Z10.
> For UGRD/Z10 found 1 forecast levels and 0 climatology
levels.
>
> ----------------------------------------
>
> Reading
records for VGRD/Z10.
> For VGRD/Z10 found 1 forecast levels and 0
climatology levels.
>
> ----------------------------------------
>
----------------------------------------
>
> Processing UGRD/Z10
versus UGRD/Z10, for observation type ANYSFC, over region
> FULL,
for interpolation method MAX(4), using 2522
> pairs.
> Computing
Continuous Statistics.
>
> ----------------------------------------
>
> Processing VGRD/Z10 versus VGRD/Z10, for observation type
ANYSFC, over region
> FULL, for interpolation method MAX(4), using
2522
> pairs.
> Computing Continuous Statistics.
> Computing Vector
Partial Sums, for forecast wind speed NA, for observation
> wind
speed NA, using 2522 pairs.
>
>
----------------------------------------
>
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V.stat
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_ctc.txt
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cts.txt
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cnt.txt
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_sl1l2.txt
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_vl1l2.txt
>
> Thanks
>
> Bob
>
> -----Original Message-----
> From: RAL
HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent:
Friday, June 25, 2010 10:40 AM
> To: Craig, Robert J Civ USAF AFWA 16
WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
>
>
Robert,
>
> Can you please send your model output and observation
file?
>
> Thanks,
>
> Paul
>
> RAL HelpDesk {for Craig,
Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>> Fri Jun 25 09:37:58 2010:
Request 38811 was acted upon.
>> Transaction: Ticket created by
craigr at offutt.af.mil
>> Queue: met_help
>> Subject: Wind
Thresholds
>> Owner: Nobody
>> Requestors:
craigr at offutt.af.mil
>> Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>>
>>
>>
I know you guys are probably sick of hearing from me but I have
another
>> question. Now that I am finally running point stat with
the updated version
>> of MET, I am not getting data for all the wind
thresholds I want. Attached
>> is the config file I am using. When
I had the threshold set to gt7.7 and
>> higher, I only got data for
gt7.7. Then I added the NA (no threshold case
>> to the list) and
now I only have data for NA case.
>>
>>
>>
>> So, what am I doing
wrong in the config file?
>>
>>
>>
>> thanks
>>
>>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
From: Paul Oldenburg
Time: Fri Jun 25 11:42:51 2010
Bob,
The comment in the point_stat_config file shows the format of the
fcst_wind_thresh list should be as follows:
// e.g. fcst_wind_thresh[] = [ "NA", "ge1.0" ];
Your fcst_wind_thresh list is:
fcst_wind_thresh[] = [ "NA gt7.7 gt10.3 gt12.9 gt18" ];
I changed this to the following format, and point_stat was able to
find pairs at
the higher thresholds.
fcst_wind_thresh[] = [ "gt7.7", "gt10.3", "gt12.9", "gt18" ];
Paul
RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>
> I ftped the file to your server.
>
> Bob
>
> -----Original Message-----
> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent: Friday, June 25, 2010 11:23 AM
> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
>
> Bob,
>
> I'm afraid the grib file you sent was corrupted, perhaps because it
was too
> large for an email attachment. The obs file seems to be ok. Can
you please
> post the grib file to our anonymous ftp site as described below:
>
> ftp ftp.rap.ucar.edu
> username = anonymous
> password = craigr at offutt.af.mil
> cd incoming/irap/met_help/craig_data
> put qc_100622_us057g1011t02b060004800
> bye
>
> Thanks,
>
> Paul
>
> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>>
>> The attached WRF file is valid at 24 Jun 6Z. The output is listed
below.
>>
>> /sphome/qcteam/METv2.0/bin/point_stat
>> /gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
>> /h/data/global/WXQC/da
>> ta/obs_prep/2010062406
>> /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF -outdir
>> /h/data/global/WXQC/data/met/data -v 2
>> GSL_RNG_TYPE=mt19937
>> GSL_RNG_SEED=4257686604
>> Forecast File:
>> /gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
>> Climatology File: none
>> Configuration File:
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>> Observation File: /h/data/global/WXQC/data/obs_prep/2010062406
>>
>> ----------------------------------------
>>
>> Reading records for UGRD/Z10.
>> For UGRD/Z10 found 1 forecast levels and 0 climatology levels.
>>
>> ----------------------------------------
>>
>> Reading records for VGRD/Z10.
>> For VGRD/Z10 found 1 forecast levels and 0 climatology levels.
>>
>> ----------------------------------------
>> ----------------------------------------
>>
>> Processing UGRD/Z10 versus UGRD/Z10, for observation type ANYSFC,
over region
>> FULL, for interpolation method MAX(4), using 2522
>> pairs.
>> Computing Continuous Statistics.
>>
>> ----------------------------------------
>>
>> Processing VGRD/Z10 versus VGRD/Z10, for observation type ANYSFC,
over region
>> FULL, for interpolation method MAX(4), using 2522
>> pairs.
>> Computing Continuous Statistics.
>> Computing Vector Partial Sums, for forecast wind speed NA, for
observation
>> wind speed NA, using 2522 pairs.
>>
>> ----------------------------------------
>>
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V.stat
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_ctc.txt
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cts.txt
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cnt.txt
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_sl1l2.txt
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_vl1l2.txt
>>
>> Thanks
>>
>> Bob
>>
>> -----Original Message-----
>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent: Friday, June 25, 2010 10:40 AM
>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
>>
>> Robert,
>>
>> Can you please send your model output and observation file?
>>
>> Thanks,
>>
>> Paul
>>
>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>> Fri Jun 25 09:37:58 2010: Request 38811 was acted upon.
>>> Transaction: Ticket created by craigr at offutt.af.mil
>>> Queue: met_help
>>> Subject: Wind Thresholds
>>> Owner: Nobody
>>> Requestors: craigr at offutt.af.mil
>>> Status: new
>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>>>
>>>
>>> I know you guys are probably sick of hearing from me but I have
another
>>> question. Now that I am finally running point stat with the
updated version
>>> of MET, I am not getting data for all the wind thresholds I want.
Attached
>>> is the config file I am using. When I had the threshold set to
gt7.7 and
>>> higher, I only got data for gt7.7. Then I added the NA (no
threshold case
>>> to the list) and now I only have data for NA case.
>>>
>>>
>>>
>>> So, what am I doing wrong in the config file?
>>>
>>>
>>>
>>> thanks
>>>
>>>
>>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #38811] Wind Thresholds
From: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Time: Fri Jun 25 11:51:25 2010
I don't know how I missed that. By the way, where is the
"stat_job_aggr_wdir" stat analysis job?
Thanks
Bob
-----Original Message-----
From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
Sent: Friday, June 25, 2010 12:43 PM
To:
Craig, Robert J Civ USAF AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu
#38811] Wind Thresholds
Bob,
The comment in the
point_stat_config file shows the format of the
fcst_wind_thresh list
should be as follows:
// e.g. fcst_wind_thresh[] = [ "NA", "ge1.0"
];
Your fcst_wind_thresh list is:
fcst_wind_thresh[] = [ "NA
gt7.7 gt10.3 gt12.9 gt18" ];
I changed this to the following
format, and point_stat was able to find pairs at
the higher
thresholds.
fcst_wind_thresh[] = [ "gt7.7", "gt10.3", "gt12.9",
"gt18" ];
Paul
RAL HelpDesk {for Craig, Robert J Civ USAF AFWA
16 WS/WXN} wrote:
> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>
> I
ftped the file to your server.
>
> Bob
>
> -----Original
Message-----
> From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
> Sent: Friday, June 25, 2010 11:23 AM
>
To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
> Subject: Re:
[rt.rap.ucar.edu #38811] Wind Thresholds
>
> Bob,
>
> I'm afraid
the grib file you sent was corrupted, perhaps because it was too
>
large for an email attachment. The obs file seems to be ok. Can you
please
> post the grib file to our anonymous ftp site as described
below:
>
> ftp ftp.rap.ucar.edu
> username = anonymous
> password
= craigr at offutt.af.mil
> cd incoming/irap/met_help/craig_data
> put
qc_100622_us057g1011t02b060004800
> bye
>
> Thanks,
>
> Paul
>
> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>>
>> The attached WRF file is valid at 24 Jun 6Z. The output is listed
below.
>>
>> /sphome/qcteam/METv2.0/bin/point_stat
>>
/gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
>> /h/data/global/WXQC/da
>> ta/obs_prep/2010062406
>>
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF -outdir
>>
/h/data/global/WXQC/data/met/data -v 2
>> GSL_RNG_TYPE=mt19937
>>
GSL_RNG_SEED=4257686604
>> Forecast File:
>>
/gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
>> Climatology File: none
>> Configuration File:
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>
Observation File: /h/data/global/WXQC/data/obs_prep/2010062406
>>
>>
----------------------------------------
>>
>> Reading records for
UGRD/Z10.
>> For UGRD/Z10 found 1 forecast levels and 0 climatology
levels.
>>
>> ----------------------------------------
>>
>>
Reading records for VGRD/Z10.
>> For VGRD/Z10 found 1 forecast levels
and 0 climatology levels.
>>
>>
----------------------------------------
>>
----------------------------------------
>>
>> Processing UGRD/Z10
versus UGRD/Z10, for observation type ANYSFC, over region
>> FULL,
for interpolation method MAX(4), using 2522
>> pairs.
>> Computing
Continuous Statistics.
>>
>>
----------------------------------------
>>
>> Processing VGRD/Z10
versus VGRD/Z10, for observation type ANYSFC, over region
>> FULL,
for interpolation method MAX(4), using 2522
>> pairs.
>> Computing
Continuous Statistics.
>> Computing Vector Partial Sums, for forecast
wind speed NA, for observation
>> wind speed NA, using 2522 pairs.
>>
>> ----------------------------------------
>>
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V.stat
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_ctc.txt
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cts.txt
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cnt.txt
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_sl1l2.txt
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_vl1l2.txt
>>
>> Thanks
>>
>> Bob
>>
>> -----Original Message-----
>> From:
RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent:
Friday, June 25, 2010 10:40 AM
>> To: Craig, Robert J Civ USAF AFWA
16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
>>
>> Robert,
>>
>> Can you please send your model output and
observation file?
>>
>> Thanks,
>>
>> Paul
>>
>> RAL HelpDesk
{for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>> Fri Jun 25
09:37:58 2010: Request 38811 was acted upon.
>>> Transaction: Ticket
created by craigr at offutt.af.mil
>>> Queue: met_help
>>>
Subject: Wind Thresholds
>>> Owner: Nobody
>>> Requestors:
craigr at offutt.af.mil
>>> Status: new
>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>>>
>>>
>>> I know you guys are probably sick of hearing from me but I have
another
>>> question. Now that I am finally running point stat with
the updated version
>>> of MET, I am not getting data for all the
wind thresholds I want. Attached
>>> is the config file I am using.
When I had the threshold set to gt7.7 and
>>> higher, I only got data
for gt7.7. Then I added the NA (no threshold case
>>> to the list)
and now I only have data for NA case.
>>>
>>>
>>>
>>> So, what am
I doing wrong in the config file?
>>>
>>>
>>>
>>> thanks
>>>
>>>
>>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
From: Paul Oldenburg
Time: Fri Jun 25 11:57:11 2010
Bob,
Please see the following example from the MET Tutorial Practical
Session:
http://www.dtcenter.org/met/users/support/OnLinePractical/OnLinePractical_210/point_stat/wind_dir.php
Paul
RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>
> I don't know how I missed that. By the way, where is the
"stat_job_aggr_wdir" stat analysis job?
>
> Thanks
>
> Bob
>
> -----Original Message-----
> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent: Friday, June 25, 2010 12:43 PM
> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
>
> Bob,
>
> The comment in the point_stat_config file shows the format of the
> fcst_wind_thresh list should be as follows:
>
> // e.g. fcst_wind_thresh[] = [ "NA", "ge1.0" ];
>
> Your fcst_wind_thresh list is:
>
> fcst_wind_thresh[] = [ "NA gt7.7 gt10.3 gt12.9 gt18" ];
>
> I changed this to the following format, and point_stat was able to
find pairs at
> the higher thresholds.
>
> fcst_wind_thresh[] = [ "gt7.7", "gt10.3", "gt12.9", "gt18" ];
>
> Paul
>
>
> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>>
>> I ftped the file to your server.
>>
>> Bob
>>
>> -----Original Message-----
>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent: Friday, June 25, 2010 11:23 AM
>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
>>
>> Bob,
>>
>> I'm afraid the grib file you sent was corrupted, perhaps because it
was too
>> large for an email attachment. The obs file seems to be ok. Can
you please
>> post the grib file to our anonymous ftp site as described below:
>>
>> ftp ftp.rap.ucar.edu
>> username = anonymous
>> password = craigr at offutt.af.mil
>> cd incoming/irap/met_help/craig_data
>> put qc_100622_us057g1011t02b060004800
>> bye
>>
>> Thanks,
>>
>> Paul
>>
>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>>>
>>> The attached WRF file is valid at 24 Jun 6Z. The output is listed
below.
>>>
>>> /sphome/qcteam/METv2.0/bin/point_stat
>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
>>> /h/data/global/WXQC/da
>>> ta/obs_prep/2010062406
>>> /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF -outdir
>>> /h/data/global/WXQC/data/met/data -v 2
>>> GSL_RNG_TYPE=mt19937
>>> GSL_RNG_SEED=4257686604
>>> Forecast File:
>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
>>> Climatology File: none
>>> Configuration File:
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>> Observation File: /h/data/global/WXQC/data/obs_prep/2010062406
>>>
>>> ----------------------------------------
>>>
>>> Reading records for UGRD/Z10.
>>> For UGRD/Z10 found 1 forecast levels and 0 climatology levels.
>>>
>>> ----------------------------------------
>>>
>>> Reading records for VGRD/Z10.
>>> For VGRD/Z10 found 1 forecast levels and 0 climatology levels.
>>>
>>> ----------------------------------------
>>> ----------------------------------------
>>>
>>> Processing UGRD/Z10 versus UGRD/Z10, for observation type ANYSFC,
over region
>>> FULL, for interpolation method MAX(4), using 2522
>>> pairs.
>>> Computing Continuous Statistics.
>>>
>>> ----------------------------------------
>>>
>>> Processing VGRD/Z10 versus VGRD/Z10, for observation type ANYSFC,
over region
>>> FULL, for interpolation method MAX(4), using 2522
>>> pairs.
>>> Computing Continuous Statistics.
>>> Computing Vector Partial Sums, for forecast wind speed NA, for
observation
>>> wind speed NA, using 2522 pairs.
>>>
>>> ----------------------------------------
>>>
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V.stat
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_ctc.txt
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cts.txt
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cnt.txt
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_sl1l2.txt
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_vl1l2.txt
>>>
>>> Thanks
>>>
>>> Bob
>>>
>>> -----Original Message-----
>>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>>> Sent: Friday, June 25, 2010 10:40 AM
>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>> Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
>>>
>>> Robert,
>>>
>>> Can you please send your model output and observation file?
>>>
>>> Thanks,
>>>
>>> Paul
>>>
>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>>> Fri Jun 25 09:37:58 2010: Request 38811 was acted upon.
>>>> Transaction: Ticket created by craigr at offutt.af.mil
>>>> Queue: met_help
>>>> Subject: Wind Thresholds
>>>> Owner: Nobody
>>>> Requestors: craigr at offutt.af.mil
>>>> Status: new
>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>>>>
>>>>
>>>> I know you guys are probably sick of hearing from me but I have
another
>>>> question. Now that I am finally running point stat with the
updated version
>>>> of MET, I am not getting data for all the wind thresholds I want.
Attached
>>>> is the config file I am using. When I had the threshold set to
gt7.7 and
>>>> higher, I only got data for gt7.7. Then I added the NA (no
threshold case
>>>> to the list) and now I only have data for NA case.
>>>>
>>>>
>>>>
>>>> So, what am I doing wrong in the config file?
>>>>
>>>>
>>>>
>>>> thanks
>>>>
>>>>
>>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #38811] Wind Thresholds
From: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Time: Fri Jun 25 12:26:37 2010
Paul, I now remember now that I tried this awhile back. When I do I
get the following error:
ERROR: PointStatConfInfo::process_config()
-> The number fcst_thresh levels provided must match the number of
fiel
ds provided in fcst_field.
Do I need to have a U/V pair in
the Fcst Field for each threshold? How did you get my config file to
work with out getting this error?
Thanks
Bob
-----Original
Message-----
From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
Sent: Friday, June 25, 2010 12:43 PM
To:
Craig, Robert J Civ USAF AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu
#38811] Wind Thresholds
Bob,
The comment in the
point_stat_config file shows the format of the
fcst_wind_thresh list
should be as follows:
// e.g. fcst_wind_thresh[] = [ "NA", "ge1.0"
];
Your fcst_wind_thresh list is:
fcst_wind_thresh[] = [ "NA
gt7.7 gt10.3 gt12.9 gt18" ];
I changed this to the following
format, and point_stat was able to find pairs at
the higher
thresholds.
fcst_wind_thresh[] = [ "gt7.7", "gt10.3", "gt12.9",
"gt18" ];
Paul
RAL HelpDesk {for Craig, Robert J Civ USAF AFWA
16 WS/WXN} wrote:
> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>
> I
ftped the file to your server.
>
> Bob
>
> -----Original
Message-----
> From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
> Sent: Friday, June 25, 2010 11:23 AM
>
To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
> Subject: Re:
[rt.rap.ucar.edu #38811] Wind Thresholds
>
> Bob,
>
> I'm afraid
the grib file you sent was corrupted, perhaps because it was too
>
large for an email attachment. The obs file seems to be ok. Can you
please
> post the grib file to our anonymous ftp site as described
below:
>
> ftp ftp.rap.ucar.edu
> username = anonymous
> password
= craigr at offutt.af.mil
> cd incoming/irap/met_help/craig_data
> put
qc_100622_us057g1011t02b060004800
> bye
>
> Thanks,
>
> Paul
>
> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>>
>> The attached WRF file is valid at 24 Jun 6Z. The output is listed
below.
>>
>> /sphome/qcteam/METv2.0/bin/point_stat
>>
/gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
>> /h/data/global/WXQC/da
>> ta/obs_prep/2010062406
>>
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF -outdir
>>
/h/data/global/WXQC/data/met/data -v 2
>> GSL_RNG_TYPE=mt19937
>>
GSL_RNG_SEED=4257686604
>> Forecast File:
>>
/gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
>> Climatology File: none
>> Configuration File:
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>
Observation File: /h/data/global/WXQC/data/obs_prep/2010062406
>>
>>
----------------------------------------
>>
>> Reading records for
UGRD/Z10.
>> For UGRD/Z10 found 1 forecast levels and 0 climatology
levels.
>>
>> ----------------------------------------
>>
>>
Reading records for VGRD/Z10.
>> For VGRD/Z10 found 1 forecast levels
and 0 climatology levels.
>>
>>
----------------------------------------
>>
----------------------------------------
>>
>> Processing UGRD/Z10
versus UGRD/Z10, for observation type ANYSFC, over region
>> FULL,
for interpolation method MAX(4), using 2522
>> pairs.
>> Computing
Continuous Statistics.
>>
>>
----------------------------------------
>>
>> Processing VGRD/Z10
versus VGRD/Z10, for observation type ANYSFC, over region
>> FULL,
for interpolation method MAX(4), using 2522
>> pairs.
>> Computing
Continuous Statistics.
>> Computing Vector Partial Sums, for forecast
wind speed NA, for observation
>> wind speed NA, using 2522 pairs.
>>
>> ----------------------------------------
>>
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V.stat
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_ctc.txt
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cts.txt
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cnt.txt
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_sl1l2.txt
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_vl1l2.txt
>>
>> Thanks
>>
>> Bob
>>
>> -----Original Message-----
>> From:
RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent:
Friday, June 25, 2010 10:40 AM
>> To: Craig, Robert J Civ USAF AFWA
16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
>>
>> Robert,
>>
>> Can you please send your model output and
observation file?
>>
>> Thanks,
>>
>> Paul
>>
>> RAL HelpDesk
{for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>> Fri Jun 25
09:37:58 2010: Request 38811 was acted upon.
>>> Transaction: Ticket
created by craigr at offutt.af.mil
>>> Queue: met_help
>>>
Subject: Wind Thresholds
>>> Owner: Nobody
>>> Requestors:
craigr at offutt.af.mil
>>> Status: new
>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>>>
>>>
>>> I know you guys are probably sick of hearing from me but I have
another
>>> question. Now that I am finally running point stat with
the updated version
>>> of MET, I am not getting data for all the
wind thresholds I want. Attached
>>> is the config file I am using.
When I had the threshold set to gt7.7 and
>>> higher, I only got data
for gt7.7. Then I added the NA (no threshold case
>>> to the list)
and now I only have data for NA case.
>>>
>>>
>>>
>>> So, what am
I doing wrong in the config file?
>>>
>>>
>>>
>>> thanks
>>>
>>>
>>
>
>
------------------------------------------------
Subject: Wind Thresholds
From: Paul Oldenburg
Time: Fri Jun 25 12:36:39 2010
Bob,
I used the attached point_stat_config file with the model and obs data
you sent
me (slightly renamed) to produce the following output. The only
change was the
one that I mentioned to fcst_wind_thresh.
Paul
$ point_stat qc_100622_us057g1011t02b060004800.grb obs_2010062406.nc
point_stat_config_WRF
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=3153470540
Forecast File: qc_100622_us057g1011t02b060004800.grb
Climatology File: none
Configuration File: point_stat_config_WRF
Observation File: obs_2010062406.nc
----------------------------------------
Reading records for UGRD/Z10.
For UGRD/Z10 found 1 forecast levels and 0 climatology levels.
----------------------------------------
Reading records for VGRD/Z10.
For VGRD/Z10 found 1 forecast levels and 0 climatology levels.
----------------------------------------
Searching 227714 observations from 62279 PrepBufr messages.
----------------------------------------
Processing UGRD/Z10 versus UGRD/Z10, for observation type ANYSFC, over
region
FULL, for interpolation method MAX(4), using 2522 pairs.
Computing Continuous Statistics.
----------------------------------------
Processing VGRD/Z10 versus VGRD/Z10, for observation type ANYSFC, over
region
FULL, for interpolation method MAX(4), using 2522 pairs.
Computing Continuous Statistics.
Computing Vector Partial Sums, for forecast wind speed >7.700, for
observation
wind speed >7.700, using 32 pairs.
Computing Vector Partial Sums, for forecast wind speed >10.300, for
observation
wind speed >10.300, using 4 pairs.
Computing Vector Partial Sums, for forecast wind speed >12.900, for
observation
wind speed >12.900, using 0 pairs.
Computing Vector Partial Sums, for forecast wind speed >18.000, for
observation
wind speed >18.000, using 0 pairs.
----------------------------------------
Output file:
/d1/pgoldenb/opt/MET_builds/METv2.0_gnu4/out/point_stat/point_stat_WRF_480000L_20100624_060000V.stat
Output file:
/d1/pgoldenb/opt/MET_builds/METv2.0_gnu4/out/point_stat/point_stat_WRF_480000L_20100624_060000V_ctc.txt
Output file:
/d1/pgoldenb/opt/MET_builds/METv2.0_gnu4/out/point_stat/point_stat_WRF_480000L_20100624_060000V_cts.txt
Output file:
/d1/pgoldenb/opt/MET_builds/METv2.0_gnu4/out/point_stat/point_stat_WRF_480000L_20100624_060000V_cnt.txt
Output file:
/d1/pgoldenb/opt/MET_builds/METv2.0_gnu4/out/point_stat/point_stat_WRF_480000L_20100624_060000V_sl1l2.txt
Output file:
/d1/pgoldenb/opt/MET_builds/METv2.0_gnu4/out/point_stat/point_stat_WRF_480000L_20100624_060000V_vl1l2.txt
RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>
> Paul, I now remember now that I tried this awhile back. When I do I
get the following error:
>
> ERROR: PointStatConfInfo::process_config() -> The number fcst_thresh
levels provided must match the number of fiel
> ds provided in fcst_field.
>
> Do I need to have a U/V pair in the Fcst Field for each threshold?
How did you get my config file to work with out getting this error?
>
> Thanks
>
> Bob
>
>
> -----Original Message-----
> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent: Friday, June 25, 2010 12:43 PM
> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
>
> Bob,
>
> The comment in the point_stat_config file shows the format of the
> fcst_wind_thresh list should be as follows:
>
> // e.g. fcst_wind_thresh[] = [ "NA", "ge1.0" ];
>
> Your fcst_wind_thresh list is:
>
> fcst_wind_thresh[] = [ "NA gt7.7 gt10.3 gt12.9 gt18" ];
>
> I changed this to the following format, and point_stat was able to
find pairs at
> the higher thresholds.
>
> fcst_wind_thresh[] = [ "gt7.7", "gt10.3", "gt12.9", "gt18" ];
>
> Paul
>
>
> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>>
>> I ftped the file to your server.
>>
>> Bob
>>
>> -----Original Message-----
>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent: Friday, June 25, 2010 11:23 AM
>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
>>
>> Bob,
>>
>> I'm afraid the grib file you sent was corrupted, perhaps because it
was too
>> large for an email attachment. The obs file seems to be ok. Can
you please
>> post the grib file to our anonymous ftp site as described below:
>>
>> ftp ftp.rap.ucar.edu
>> username = anonymous
>> password = craigr at offutt.af.mil
>> cd incoming/irap/met_help/craig_data
>> put qc_100622_us057g1011t02b060004800
>> bye
>>
>> Thanks,
>>
>> Paul
>>
>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>>>
>>> The attached WRF file is valid at 24 Jun 6Z. The output is listed
below.
>>>
>>> /sphome/qcteam/METv2.0/bin/point_stat
>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
>>> /h/data/global/WXQC/da
>>> ta/obs_prep/2010062406
>>> /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF -outdir
>>> /h/data/global/WXQC/data/met/data -v 2
>>> GSL_RNG_TYPE=mt19937
>>> GSL_RNG_SEED=4257686604
>>> Forecast File:
>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100622_us057g1011t02b060004800
>>> Climatology File: none
>>> Configuration File:
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>> Observation File: /h/data/global/WXQC/data/obs_prep/2010062406
>>>
>>> ----------------------------------------
>>>
>>> Reading records for UGRD/Z10.
>>> For UGRD/Z10 found 1 forecast levels and 0 climatology levels.
>>>
>>> ----------------------------------------
>>>
>>> Reading records for VGRD/Z10.
>>> For VGRD/Z10 found 1 forecast levels and 0 climatology levels.
>>>
>>> ----------------------------------------
>>> ----------------------------------------
>>>
>>> Processing UGRD/Z10 versus UGRD/Z10, for observation type ANYSFC,
over region
>>> FULL, for interpolation method MAX(4), using 2522
>>> pairs.
>>> Computing Continuous Statistics.
>>>
>>> ----------------------------------------
>>>
>>> Processing VGRD/Z10 versus VGRD/Z10, for observation type ANYSFC,
over region
>>> FULL, for interpolation method MAX(4), using 2522
>>> pairs.
>>> Computing Continuous Statistics.
>>> Computing Vector Partial Sums, for forecast wind speed NA, for
observation
>>> wind speed NA, using 2522 pairs.
>>>
>>> ----------------------------------------
>>>
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V.stat
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_ctc.txt
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cts.txt
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_cnt.txt
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_sl1l2.txt
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_480000L_20100624_060000V_vl1l2.txt
>>>
>>> Thanks
>>>
>>> Bob
>>>
>>> -----Original Message-----
>>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>>> Sent: Friday, June 25, 2010 10:40 AM
>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>> Subject: Re: [rt.rap.ucar.edu #38811] Wind Thresholds
>>>
>>> Robert,
>>>
>>> Can you please send your model output and observation file?
>>>
>>> Thanks,
>>>
>>> Paul
>>>
>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>>> Fri Jun 25 09:37:58 2010: Request 38811 was acted upon.
>>>> Transaction: Ticket created by craigr at offutt.af.mil
>>>> Queue: met_help
>>>> Subject: Wind Thresholds
>>>> Owner: Nobody
>>>> Requestors: craigr at offutt.af.mil
>>>> Status: new
>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38811 >
>>>>
>>>>
>>>> I know you guys are probably sick of hearing from me but I have
another
>>>> question. Now that I am finally running point stat with the
updated version
>>>> of MET, I am not getting data for all the wind thresholds I want.
Attached
>>>> is the config file I am using. When I had the threshold set to
gt7.7 and
>>>> higher, I only got data for gt7.7. Then I added the NA (no
threshold case
>>>> to the list) and now I only have data for NA case.
>>>>
>>>>
>>>>
>>>> So, what am I doing wrong in the config file?
>>>>
>>>>
>>>>
>>>> thanks
>>>>
>>>>
>>
>
>
------------------------------------------------
Subject: Wind Thresholds
From: Paul Oldenburg
Time: Fri Jun 25 12:36:39 2010
////////////////////////////////////////////////////////////////////////////////
//
// Default point_stat configuration file
//
////////////////////////////////////////////////////////////////////////////////
//
// Specify a name to designate the model being verified. This name
will be
// written to the second column of the ASCII output generated.
//
model = "WRF";
//
// Beginning and ending time offset values in seconds for observations
// to be used. These time offsets are defined in reference to the
// forecast valid time, v. Observations with a valid time falling in
the
// window [v+beg_ds, v+end_ds] will be used.
// These selections are overridden by the command line arguments
// -valid_beg and -valid_end.
//
beg_ds = -600;
end_ds = 600;
//
// Specify a comma-separated list of fields to be verified. The
forecast and
// observation fields may be specified separately. If the obs_field
parameter
// is left blank, it will default to the contents of fcst_field.
//
// Each field is specified as a grib code or corresponding grib code
// abbreviation followed by an accumulation or vertical level
indicator.
//
// Each verification field is specified as one of the following:
// GC/ANNN for accumulation interval NNN
// GC/ZNNN for vertical level NNN
// GC/PNNN for pressure level NNN in hPa
// GC/PNNN-NNN for a range of pressure levels in hPa
// GC/LNNN for a generic level type
// GC/RNNN for a specific GRIB record number
// Where GC is the number of or abbreviation for the grib code
// to be verified.
// http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
//
// NOTE: To verify winds as vectors rather than scalars,
// specify UGRD (or 33) followd by VGRD (or 34) with the
// same level values.
//
// NOTE: To process a probability field, add "/PROB", such as
"POP/Z0/PROB".
//
// e.g. fcst_field[] = [ "SPFH/P500", "TMP/P500" ];
//
fcst_field[] = [ "33/Z10", "34/Z10"];
obs_field[] = [];
//
// Specify a comma-separated list of groups of thresholds to be
applied to the
// fields listed above. Thresholds for the forecast and observation
fields
// may be specified separately. If the obs_thresh parameter is left
blank,
// it will default to the contents of fcst_thresh.
//
// At least one threshold must be provided for each field listed
above. The
// lengths of the "fcst_field" and "fcst_thresh" arrays must match, as
must
// lengths of the "obs_field" and "obs_thresh" arrays. To apply
multiple
// thresholds to a field, separate the threshold values with a space.
//
// Each threshold must be preceded by a two letter indicator for the
type of
// thresholding to be performed:
// 'lt' for less than 'le' for less than or equal to
// 'eq' for equal to 'ne' for not equal to
// 'gt' for greater than 'ge' for greater than or equal to
//
// NOTE: Thresholds for probabilities must be preceeded by "ge".
//
// e.g. fcst_thresh[] = [ "gt80", "gt273" ];
//
fcst_thresh[] = [ "",""];
obs_thresh[] = [];
//
// Specify a comma-separated list of thresholds to be used when
computing
// VL1L2 and VAL1L2 partial sums for winds. The thresholds are
applied to the
// wind speed values derived from each U/V pair. Only those U/V pairs
which meet
// the wind speed threshold criteria are retained. If the
obs_wind_thresh
// parameter is left blank, it will default to the contents of
fcst_wind_thresh.
//
// To apply multiple wind speed thresholds, separate the threshold
values with a
// space. Use "NA" to indicate that no wind speed threshold should be
applied.
//
// Each threshold must be preceded by a two letter indicator for the
type of
// thresholding to be performed:
// 'lt' for less than 'le' for less than or equal to
// 'eq' for equal to 'ne' for not equal to
// 'gt' for greater than 'ge' for greater than or equal to
// 'NA' for no threshold
//
// e.g. fcst_wind_thresh[] = [ "NA", "ge1.0" ];
//
//fcst_wind_thresh[] = [ "NA gt7.7 gt10.3 gt12.9 gt18" ];
fcst_wind_thresh[] = [ "gt7.7", "gt10.3", "gt12.9", "gt18" ];
obs_wind_thresh[] = [];
//
// Specify a comma-separated list of PrepBufr message types with which
// to perform the verification. Statistics will be computed
separately
// for each message type specified. At least one PrepBufr message
type
// must be provided.
// List of valid message types:
// ADPUPA AIRCAR AIRCFT ADPSFC ERS1DA GOESND GPSIPW
// MSONET PROFLR QKSWND RASSDA SATEMP SATWND SFCBOG
// SFCSHP SPSSMI SYNDAT VADWND
// ANYAIR (= AIRCAR, AIRCFT)
// ANYSFC (= ADPSFC, SFCSHP, ADPUPA, PROFLR)
// ONLYSF (= ADPSFC, SFCSHP)
//
http://www.emc.ncep.noaa.gov/mmb/data_processing/prepbufr.doc/table_1.htm
//
// e.g. message_type[] = [ "ADPUPA", "AIRCAR" ];
//
message_type[] = [ "ANYSFC" ];
//
// Specify a comma-separated list of grids to be used in masking the
data over
// which to perform scoring. An empty list indicates that no masking
grid
// should be performed. The standard NCEP grids are named "GNNN"
where NNN
// indicates the three digit grid number. Enter "FULL" to score over
the
// entire domain.
// http://www.nco.ncep.noaa.gov/pmb/docs/on388/tableb.html
//
// e.g. mask_grid[] = [ "FULL" ];
//
mask_grid[] = [ "FULL" ];
//
// Specify a comma-separated list of masking regions to be applied.
// An empty list indicates that no additional masks should be used.
// The masking regions may be defined in one of 4 ways:
//
// (1) An ASCII file containing a lat/lon polygon.
// Latitude in degrees north and longitude in degrees east.
// By default, the first and last polygon points are connected.
// e.g. "MET_BASE/data/poly/EAST.poly" which consists of n points:
// "poly_name lat1 lon1 lat2 lon2... latn lonn"
//
// (2) The NetCDF output of the gen_poly_mask tool.
//
// (3) A NetCDF data file, followed by the name of the NetCDF variable
// to be used, and optionally, a threshold to be applied to the
field.
// e.g. "sample.nc var_name gt0.00"
//
// (4) A GRIB data file, followed by a description of the field
// to be used, and optionally, a threshold to be applied to the
field.
// e.g. "sample.grb APCP/A3 gt0.00"
//
// Any NetCDF or GRIB file used must have the same grid dimensions as
the
// data being verified.
//
// MET_BASE may be used in the path for the files above.
//
// e.g. mask_poly[] = [ "MET_BASE/data/poly/EAST.poly",
// "poly_mask.ncf",
// "sample.nc APCP",
// "sample.grb HGT/Z0 gt100.0" ];
//
mask_poly[]=[];
//
// Specify the name of an ASCII file containing a space-separated list
of
// station ID's at which to perform verification. Each station ID
specified
// is treated as an individual masking region.
//
// An empty list file name indicates that no station ID masks should
be used.
//
// MET_BASE may be used in the path for the station ID mask file name.
//
// e.g. mask_sid = "MET_BASE/data/stations/CONUS.stations";
//
mask_sid = "";
//
// Specify a comma-separated list of values for alpha to be used when
computing
// confidence intervals. Values of alpha must be between 0 and 1.
//
// e.g. ci_alpha[] = [ 0.05, 0.10 ];
//
ci_alpha[] = [ 0.05 ];
//
// Specify the method to be used for computing bootstrap confidence
intervals.
// The value for this is interpreted as follows:
// (0) Use the BCa interval method (computationally intensive)
// (1) Use the percentile interval method
//
boot_interval = 1;
//
// Specify a proportion between 0 and 1 to define the replicate sample
size
// to be used when computing percentile intervals. The replicate
sample
// size is set to boot_rep_prop * n, where n is the number of raw data
points.
//
// e.g boot_rep_prop = 0.80;
//
boot_rep_prop = 1.0;
//
// Specify the number of times each set of matched pair data should be
// resampled when computing bootstrap confidence intervals. A value
of
// zero disables the computation of bootstrap condifence intervals.
//
// e.g. n_boot_rep = 1000;
//
n_boot_rep = 1000;
//
// Specify the name of the random number generator to be used. See
the MET
// Users Guide for a list of possible random number generators.
//
boot_rng = "mt19937";
//
// Specify the seed value to be used when computing bootstrap
confidence
// intervals. If left unspecified, the seed will change for each run
and
// the computed bootstrap confidence intervals will not be
reproducable.
//
boot_seed = "";
//
// Specify a comma-separated list of interpolation method(s) to be
used
// for comparing the forecast grid to the observation points. String
values
// are interpreted as follows:
// MIN = Minimum in the neighborhood
// MAX = Maximum in the neighborhood
// MEDIAN = Median in the neighborhood
// UW_MEAN = Unweighted mean in the neighborhood
// DW_MEAN = Distance-weighted mean in the neighborhood
// LS_FIT = Least-squares fit in the neighborhood
//
// In all cases, vertical interpolation is performed in the natural
log
// of pressure of the levels above and below the observation.
//
// e.g. interp_method[] = [ "UW_MEAN", "MEDIAN" ];
//
interp_method[] = [ "MAX" ];
//
// Specify a comma-separated list of box widths to be used by the
// interpolation techniques listed above. A value of 1 indicates that
// the nearest neighbor approach should be used. For a value of n
// greater than 1, the n*n grid points closest to the observation
define
// the neighborhood.
//
// e.g. interp_width = [ 1, 3, 5 ];
//
interp_width[] = [ 2 ];
//
// When interpolating, compute a ratio of the number of valid data
points
// to the total number of points in the neighborhood. If that ratio
is
// less than this threshold, do not include the observation. This
// threshold must be between 0 and 1. Setting this threshold to 1
will
// require that each observation be surrounded by n*n valid forecast
// points.
//
// e.g. interp_thresh = 1.0;
//
interp_thresh = 1.0;
//
// Specify flags to indicate the type of data to be output:
// (1) STAT and FHO Text Files, Forecast, Hit, Observation Rates:
// Total (TOTAL),
// Forecast Rate (F_RATE),
// Hit Rate (H_RATE),
// Observation Rate (O_RATE)
//
// (2) STAT and CTC Text Files, Contingency Table Counts:
// Total (TOTAL),
// Forecast Yes and Observation Yes Count (FY_OY),
// Forecast Yes and Observation No Count (FY_ON),
// Forecast No and Observation Yes Count (FN_OY),
// Forecast No and Observation No Count (FN_ON)
//
// (3) STAT and CTS Text Files, Contingency Table Scores:
// Total (TOTAL),
// Base Rate (BASER), BASER_CL, BASER_CU,
// Forecast Mean (FMEAN), FMEAN_CL, FMEAN_CU,
// Accuracy (ACC), ACC_CL, ACC_CU,
// Frequency Bias (FBIAS),
// Probability of Detecting Yes (PODY), PODY_CL, PODY_CU,
// Probability of Detecting No (PODN), PODN_CL, PODN_CU,
// Probability of False Detection (POFD), POFD_CL, POFD_CU,
// False Alarm Ratio (FAR), FAR_CL, FAR_CU,
// Critical Success Index (CSI), CSI_CL, CSI_CU,
// Gilbert Skill Score (GSS),
// Hanssen and Kuipers Discriminant (HK), HK_CL, HK_CU,
// Heidke Skill Score (HSS),
// Odds Ratio (ODDS), ODDS_CL, ODDS_CU
//
// (4) STAT and CNT Text Files, Statistics of Continuous Variables:
// Total (TOTAL),
// Forecast Mean (FBAR), FBAR_CL, FBAR_CU,
// Forecast Standard Deviation (FSTDEV), FSTDEV_CL,
FSTDEV_CU
// Observation Mean (OBAR), OBAR_CL, OBAR_CU,
// Observation Standard Deviation (OSTDEV), OSTDEV_CL,
OSTDEV_CU,
// Pearson's Correlation Coefficient (PR_CORR), PR_CORR_CL,
PR_CORR_CU,
// Spearman's Rank Correlation Coefficient (SP_CORR),
// Kendall Tau Rank Correlation Coefficient (KT_CORR),
// Number of ranks compared (RANKS),
// Number of tied ranks in the forecast field (FRANK_TIES),
// Number of tied ranks in the observation field
(ORANK_TIES),
// Mean Error (ME), ME_CL, ME_CU,
// Standard Deviation of the Error (ESTDEV), ESTDEV_CL,
ESTDEV_CU,
// Bias (BIAS = FBAR - OBAR),
// Mean Absolute Error (MAE),
// Mean Squared Error (MSE),
// Bias-Corrected Mean Squared Error (BCMSE),
// Root Mean Squared Error (RMSE),
// Percentiles of the Error (E10, E25, E50, E75, E90)
//
// NOTE: CL and CU values define lower and upper
// confidence interval limits.
//
// (5) STAT and SL1L2 Text Files, Scalar Partial Sums:
// Total (TOTAL),
// Forecast Mean (FBAR),
// = mean(f)
// Observation Mean (OBAR),
// = mean(o)
// Forecast*Observation Product Mean (FOBAR),
// = mean(f*o)
// Forecast Squared Mean (FFBAR),
// = mean(f^2)
// Observation Squared Mean (OOBAR)
// = mean(o^2)
//
// (6) STAT and SAL1L2 Text Files, Scalar Anomaly Partial Sums:
// Total (TOTAL),
// Forecast Anomaly Mean (FABAR),
// = mean(f-c)
// Observation Anomaly Mean (OABAR),
// = mean(o-c)
// Product of Forecast and Observation Anomalies Mean
(FOABAR),
// = mean((f-c)*(o-c))
// Forecast Anomaly Squared Mean (FFABAR),
// = mean((f-c)^2)
// Observation Anomaly Squared Mean (OOABAR)
// = mean((o-c)^2)
//
// (7) STAT and VL1L2 Text Files, Vector Partial Sums:
// Total (TOTAL),
// U-Forecast Mean (UFBAR),
// = mean(uf)
// V-Forecast Mean (VFBAR),
// = mean(vf)
// U-Observation Mean (UOBAR),
// = mean(uo)
// V-Observation Mean (VOBAR),
// = mean(vo)
// U-Product Plus V-Product (UVFOBAR),
// = mean(uf*uo+vf*vo)
// U-Forecast Squared Plus V-Forecast Squared (UVFFBAR),
// = mean(uf^2+vf^2)
// U-Observation Squared Plus V-Observation Squared
(UVOOBAR)
// = mean(uo^2+vo^2)
//
// (8) STAT and VAL1L2 Text Files, Vector Anomaly Partial Sums:
// U-Forecast Anomaly Mean (UFABAR),
// = mean(uf-uc)
// V-Forecast Anomaly Mean (VFABAR),
// = mean(vf-vc)
// U-Observation Anomaly Mean (UOABAR),
// = mean(uo-uc)
// V-Observation Anomaly Mean (VOABAR),
// = mean(vo-vc)
// U-Anomaly Product Plus V-Anomaly Product (UVFOABAR),
// = mean((uf-uc)*(uo-uc)+(vf-vc)*(vo-vc))
// U-Forecast Anomaly Squared Plus V-Forecast Anomaly
Squared (UVFFABAR),
// = mean((uf-uc)^2+(vf-vc)^2)
// U-Observation Anomaly Squared Plus V-Observation Anomaly
Squared (UVOOABAR)
// = mean((uo-uc)^2+(vo-vc)^2)
//
// (9) STAT and PCT Text Files, Nx2 Probability Contingency Table
Counts:
// Total (TOTAL),
// Number of Forecast Probability Thresholds (N_THRESH),
// Probability Threshold Value (THRESH_i),
// Row Observation Yes Count (OY_i),
// Row Observation No Count (ON_i),
// NOTE: Previous 3 columns repeated for each row in the
table
// Last Probability Threshold Value (THRESH_n)
//
// (10) STAT and PSTD Text Files, Nx2 Probability Contingency Table
Scores:
// Total (TOTAL),
// Number of Forecast Probability Thresholds (N_THRESH),
// Reliability (RELIABILITY),
// Resolution (RESOLUTION),
// Uncertainty (UNCERTAINTY),
// Area Under the ROC Curve (ROC_AUC),
// Brier Score (BRIER), BRIER_NCL, BRIER_NCU,
// Probability Threshold Value (THRESH_i)
// NOTE: Previous column repeated for each probability
threshold
//
// (11) STAT and PJC Text Files, Joint/Continuous Statistics of
// Probabilistic Variables:
// Total (TOTAL),
// Number of Forecast Probability Thresholds (N_THRESH),
// Probability Threshold Value (THRESH_i),
// Observation Yes Count Divided by Total (OY_TP_i),
// Observation No Count Divided by Total (ON_TP_i),
// Calibration (CALIBRATION_i),
// Refinement (REFINEMENT_i),
// Likelikhood (LIKELIHOOD_i),
// Base Rate (BASER_i),
// NOTE: Previous 7 columns repeated for each row in the
table
// Last Probability Threshold Value (THRESH_n)
//
// (12) STAT and PRC Text Files, ROC Curve Points for
// Probabilistic Variables:
// Total (TOTAL),
// Number of Forecast Probability Thresholds (N_THRESH),
// Probability Threshold Value (THRESH_i),
// Probability of Detecting Yes (PODY_i),
// Probability of False Detection (POFD_i),
// NOTE: Previous 3 columns repeated for each row in the
table
// Last Probability Threshold Value (THRESH_n)
//
// (13) STAT and MPR Text Files, Matched Pair Data:
// Total (TOTAL),
// Index (INDEX),
// Latitude (LAT),
// Longitude (LON),
// Level (LEVEL),
// Forecast Value (FCST),
// Observation Value (OBS),
// Climatological Value (CLIMO),
// Interpolation Methold (INTERP_MTHD),
// Interpolation Points (INTERP_PNTS)
//
// In the expressions above, f are forecast values, o are observed
values,
// and c are climatological values.
//
// Values for these flags are interpreted as follows:
// (0) Do not generate output of this type
// (1) Write output to a STAT file
// (2) Write output to a STAT file and a text file
//
output_flag[] = [ 0, 2, 2, 2, 2, 0, 2, 0, 0, 0, 0, 0, 0 ];
//
// Flag to indicate whether Kendall's Tau and Spearman's Rank
Correlation
// Coefficients should be computed. Computing them over large
datasets is
// computationally intensive and slows down the runtime execution
significantly.
// (0) Do not compute these correlation coefficients
// (1) Compute these correlation coefficients
//
rank_corr_flag = 0;
//
// Specify the GRIB Table 2 parameter table version number to be used
// for interpreting GRIB codes.
// http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
//
grib_ptv = 2;
//
// Directory where temporary files should be written.
//
tmp_dir = "/tmp";
//
// Prefix to be used for the output file names.
//
output_prefix="WRF";
//
// Indicate a version number for the contents of this configuration
file.
// The value should generally not be modified.
//
version = "V 2.0";
------------------------------------------------
More information about the Met_help
mailing list