[Met_help] [rt.rap.ucar.edu #38629] History for Point Stat Question
RAL HelpDesk {for Paul Oldenburg}
met_help at ucar.edu
Thu Jun 24 15:57:24 MDT 2010
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
I am running point stat on WRF winds but I am not getting any stats. My
code to run this is:
/sphome/qcteam/METv2.0/bin/point_stat
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
/h/data/global/WXQC/data/obs_pr
ep/2010061512 /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
-outdir /h/data/global/WXQC/data/met/data -v 2
Attached is the config file along with the ob file and model file.
The following is the output I get. Why the 0 pairs? Thanks Bob Craig
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=3818134092
Forecast File:
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
Climatology File: none
Configuration File: /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
----------------------------------------
Reading records for UGRD/L10.
For UGRD/L10 found 1 forecast levels and 0 climatology levels.
----------------------------------------
Reading records for VGRD/L10.
For VGRD/L10 found 1 forecast levels and 0 climatology levels.
----------------------------------------
Searching 10473 observations from 170 PrepBufr messages.
----------------------------------------
Processing UGRD/L10 versus UGRD/L10, for observation type ADPSFC, over
region FULL, for interpolation method MAX(4), using 0 pairs.
----------------------------------------
Processing VGRD/L10 versus VGRD/L10, for observation type ADPSFC, over
region FULL, for interpolation method MAX(4), using 0 pairs.
----------------------------------------
Output file:
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
at
Output file:
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
c.txt
Output file:
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
s.txt
Output file:
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
t.txt
Output file:
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
1l2.txt
Output file:
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
1l2.txt
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
From: Paul Oldenburg
Time: Fri Jun 18 10:56:21 2010
Robert,
Your observation file has a message type of PROFLR, but your config
file is
specifying a message type of ADPSFC. This is one potential problem.
Another
problem is that the model output appears to be valid at 20100615_12
and the
observations are valid 24 hours later. You could change your config
file to
open the valid time window to a larger value using beg_ds and end_ds.
Please
let us know if you have any other questions.
Thanks,
Paul
RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
> Thu Jun 17 15:45:21 2010: Request 38629 was acted upon.
> Transaction: Ticket created by craigr at offutt.af.mil
> Queue: met_help
> Subject: Point Stat Question
> Owner: Nobody
> Requestors: craigr at offutt.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>
>
> I am running point stat on WRF winds but I am not getting any stats.
My
> code to run this is:
>
>
>
> /sphome/qcteam/METv2.0/bin/point_stat
> /gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
> /h/data/global/WXQC/data/obs_pr
>
> ep/2010061512
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
> -outdir /h/data/global/WXQC/data/met/data -v 2
>
>
>
>
>
> Attached is the config file along with the ob file and model file.
>
>
>
>
>
> The following is the output I get. Why the 0 pairs? Thanks Bob
Craig
>
>
>
>
>
> GSL_RNG_TYPE=mt19937
>
> GSL_RNG_SEED=3818134092
>
> Forecast File:
> /gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>
> Climatology File: none
>
> Configuration File:
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>
> Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>
>
>
> ----------------------------------------
>
>
>
> Reading records for UGRD/L10.
>
> For UGRD/L10 found 1 forecast levels and 0 climatology levels.
>
>
>
> ----------------------------------------
>
>
>
> Reading records for VGRD/L10.
>
> For VGRD/L10 found 1 forecast levels and 0 climatology levels.
>
>
>
> ----------------------------------------
>
>
>
> Searching 10473 observations from 170 PrepBufr messages.
>
>
>
> ----------------------------------------
>
>
>
> Processing UGRD/L10 versus UGRD/L10, for observation type ADPSFC,
over
> region FULL, for interpolation method MAX(4), using 0 pairs.
>
>
>
> ----------------------------------------
>
>
>
> Processing VGRD/L10 versus VGRD/L10, for observation type ADPSFC,
over
> region FULL, for interpolation method MAX(4), using 0 pairs.
>
>
>
> ----------------------------------------
>
>
>
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
> at
>
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
> c.txt
>
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
> s.txt
>
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
> t.txt
>
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
> 1l2.txt
>
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
> 1l2.txt
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #38629] Point Stat Question
From: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Time: Tue Jun 22 07:56:31 2010
Thanks for looking into this. I inadvertently left off my config file
name from the command string, hence it defaulted to profiler data. My
config file has "ANYSFC" which the documentation says includes ADPUPA
(RAOB data).
For the observations, I am pulling the prepbufr files
from ftp.ncep.noaa.gov in the
/pub/data/nccf/com/gfs/prod/gdas.20100615 for the 15 Jun files for the
case below. An example of the filenames I am pulling are
gdas2.t12z.prepbufr.nr. I rename these files to 20100615_12 in this
case. Since this file came out of the 20100615 directory I would
expect it to be valid for 20100615 at 12Z. Is there a way I can dump
the header of the prepbufr file to see the date that's in the file?
Bob
-----Original Message-----
From: RAL HelpDesk {for Paul
Oldenburg} [mailto:met_help at ucar.edu]
Sent: Friday, June 18, 2010
11:56 AM
To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Subject: Re:
[rt.rap.ucar.edu #38629] Point Stat Question
Robert,
Your
observation file has a message type of PROFLR, but your config file is
specifying a message type of ADPSFC. This is one potential problem.
Another
problem is that the model output appears to be valid at
20100615_12 and the
observations are valid 24 hours later. You could
change your config file to
open the valid time window to a larger
value using beg_ds and end_ds. Please
let us know if you have any
other questions.
Thanks,
Paul
RAL HelpDesk {for Craig,
Robert J Civ USAF AFWA 16 WS/WXN} wrote:
> Thu Jun 17 15:45:21 2010:
Request 38629 was acted upon.
> Transaction: Ticket created by
craigr at offutt.af.mil
> Queue: met_help
> Subject: Point
Stat Question
> Owner: Nobody
> Requestors:
craigr at offutt.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>
>
> I
am running point stat on WRF winds but I am not getting any stats. My
> code to run this is:
>
>
>
>
/sphome/qcteam/METv2.0/bin/point_stat
>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>
/h/data/global/WXQC/data/obs_pr
>
> ep/2010061512
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
> -outdir
/h/data/global/WXQC/data/met/data -v 2
>
>
>
>
>
>
Attached is the config file along with the ob file and model file.
>
>
>
>
>
> The following is the output I get. Why the 0
pairs? Thanks Bob Craig
>
>
>
>
>
> GSL_RNG_TYPE=mt19937
>
> GSL_RNG_SEED=3818134092
>
> Forecast File:
>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>
> Climatology File: none
>
> Configuration File:
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>
>
Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>
>
>
> ----------------------------------------
>
>
>
> Reading
records for UGRD/L10.
>
> For UGRD/L10 found 1 forecast levels and
0 climatology levels.
>
>
>
>
----------------------------------------
>
>
>
> Reading
records for VGRD/L10.
>
> For VGRD/L10 found 1 forecast levels and
0 climatology levels.
>
>
>
>
----------------------------------------
>
>
>
> Searching
10473 observations from 170 PrepBufr messages.
>
>
>
>
----------------------------------------
>
>
>
> Processing
UGRD/L10 versus UGRD/L10, for observation type ADPSFC, over
> region
FULL, for interpolation method MAX(4), using 0 pairs.
>
>
>
>
----------------------------------------
>
>
>
> Processing
VGRD/L10 versus VGRD/L10, for observation type ADPSFC, over
> region
FULL, for interpolation method MAX(4), using 0 pairs.
>
>
>
>
----------------------------------------
>
>
>
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
> at
>
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
> c.txt
>
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
> s.txt
>
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
> t.txt
>
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
> 1l2.txt
>
> Output file:
>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
> 1l2.txt
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
From: Paul Oldenburg
Time: Tue Jun 22 14:30:07 2010
Bob,
We looked more closely at the observation data that you included in
your
original email, and we noticed that they don't appear to be surface
observations. When verifying 2m temp or 10m winds in point_stat, only
observations of type ADPSFC or SFCSHP can be used (or ONLYSF for
either).
Please change your point_stat config fcst_field setting to the
following value:
fcst_field[] = [ "33/Z10", "34/Z10"];
also, please change the message type:
message_type[] = [ "ADPSFC" ];
Please let us know if you have any more questions.
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=38629 >
>
> Thanks for looking into this. I inadvertently left off my config
file name from the command string, hence it defaulted to profiler
data. My config file has "ANYSFC" which the documentation says
includes ADPUPA (RAOB data).
>
> For the observations, I am pulling the prepbufr files from
ftp.ncep.noaa.gov in the /pub/data/nccf/com/gfs/prod/gdas.20100615 for
the 15 Jun files for the case below. An example of the filenames I am
pulling are gdas2.t12z.prepbufr.nr. I rename these files to
20100615_12 in this case. Since this file came out of the 20100615
directory I would expect it to be valid for 20100615 at 12Z. Is there
a way I can dump the header of the prepbufr file to see the date
that's in the file?
>
> Bob
>
>
> -----Original Message-----
> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent: Friday, June 18, 2010 11:56 AM
> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>
> Robert,
>
> Your observation file has a message type of PROFLR, but your config
file is
> specifying a message type of ADPSFC. This is one potential problem.
Another
> problem is that the model output appears to be valid at 20100615_12
and the
> observations are valid 24 hours later. You could change your config
file to
> open the valid time window to a larger value using beg_ds and
end_ds. Please
> let us know if you have any other questions.
>
> Thanks,
>
> Paul
>
>
>
> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>> Thu Jun 17 15:45:21 2010: Request 38629 was acted upon.
>> Transaction: Ticket created by craigr at offutt.af.mil
>> Queue: met_help
>> Subject: Point Stat Question
>> Owner: Nobody
>> Requestors: craigr at offutt.af.mil
>> Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>
>>
>> I am running point stat on WRF winds but I am not getting any
stats. My
>> code to run this is:
>>
>>
>>
>> /sphome/qcteam/METv2.0/bin/point_stat
>> /gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>> /h/data/global/WXQC/data/obs_pr
>>
>> ep/2010061512
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>> -outdir /h/data/global/WXQC/data/met/data -v 2
>>
>>
>>
>>
>>
>> Attached is the config file along with the ob file and model file.
>>
>>
>>
>>
>>
>> The following is the output I get. Why the 0 pairs? Thanks Bob
Craig
>>
>>
>>
>>
>>
>> GSL_RNG_TYPE=mt19937
>>
>> GSL_RNG_SEED=3818134092
>>
>> Forecast File:
>> /gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>
>> Climatology File: none
>>
>> Configuration File:
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>
>> Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>>
>>
>>
>> ----------------------------------------
>>
>>
>>
>> Reading records for UGRD/L10.
>>
>> For UGRD/L10 found 1 forecast levels and 0 climatology levels.
>>
>>
>>
>> ----------------------------------------
>>
>>
>>
>> Reading records for VGRD/L10.
>>
>> For VGRD/L10 found 1 forecast levels and 0 climatology levels.
>>
>>
>>
>> ----------------------------------------
>>
>>
>>
>> Searching 10473 observations from 170 PrepBufr messages.
>>
>>
>>
>> ----------------------------------------
>>
>>
>>
>> Processing UGRD/L10 versus UGRD/L10, for observation type ADPSFC,
over
>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>
>>
>>
>> ----------------------------------------
>>
>>
>>
>> Processing VGRD/L10 versus VGRD/L10, for observation type ADPSFC,
over
>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>
>>
>>
>> ----------------------------------------
>>
>>
>>
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
>> at
>>
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>> c.txt
>>
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>> s.txt
>>
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
>> t.txt
>>
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
>> 1l2.txt
>>
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
>> 1l2.txt
>>
>
>
>
------------------------------------------------
Subject: Point Stat Question
From: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Time: Tue Jun 22 15:16:40 2010
I am getting the same 0 pair error with the attached ob file. I built
this
file using the attached config file. Is there surface observations in
my
file? Is the valid time of the file correct? I am using the -dump
option but
I am not getting any file showing up in the directory I set for this.
Bob
-----Original Message-----
From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
Sent: Tuesday, June 22, 2010 3:30 PM
To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
Bob,
We looked more closely at the observation data that you included in
your
original email, and we noticed that they don't appear to be surface
observations. When verifying 2m temp or 10m winds in point_stat, only
observations of type ADPSFC or SFCSHP can be used (or ONLYSF for
either).
Please change your point_stat config fcst_field setting to the
following
value:
fcst_field[] = [ "33/Z10", "34/Z10"];
also, please change the message type:
message_type[] = [ "ADPSFC" ];
Please let us know if you have any more questions.
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=38629 >
>
> Thanks for looking into this. I inadvertently left off my config
file name
> from the command string, hence it defaulted to profiler data. My
config file
> has "ANYSFC" which the documentation says includes ADPUPA (RAOB
data).
>
> For the observations, I am pulling the prepbufr files from
ftp.ncep.noaa.gov
> in the /pub/data/nccf/com/gfs/prod/gdas.20100615 for the 15 Jun
files for
> the case below. An example of the filenames I am pulling are
> gdas2.t12z.prepbufr.nr. I rename these files to 20100615_12 in
this case.
> Since this file came out of the 20100615 directory I would expect it
to be
> valid for 20100615 at 12Z. Is there a way I can dump the header of
the
> prepbufr file to see the date that's in the file?
>
> Bob
>
>
> -----Original Message-----
> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent: Friday, June 18, 2010 11:56 AM
> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>
> Robert,
>
> Your observation file has a message type of PROFLR, but your config
file is
> specifying a message type of ADPSFC. This is one potential problem.
> Another
> problem is that the model output appears to be valid at 20100615_12
and the
> observations are valid 24 hours later. You could change your config
file to
> open the valid time window to a larger value using beg_ds and
end_ds.
> Please
> let us know if you have any other questions.
>
> Thanks,
>
> Paul
>
>
>
> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>> Thu Jun 17 15:45:21 2010: Request 38629 was acted upon.
>> Transaction: Ticket created by craigr at offutt.af.mil
>> Queue: met_help
>> Subject: Point Stat Question
>> Owner: Nobody
>> Requestors: craigr at offutt.af.mil
>> Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>
>>
>> I am running point stat on WRF winds but I am not getting any
stats. My
>> code to run this is:
>>
>>
>>
>> /sphome/qcteam/METv2.0/bin/point_stat
>> /gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>> /h/data/global/WXQC/data/obs_pr
>>
>> ep/2010061512
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>> -outdir /h/data/global/WXQC/data/met/data -v 2
>>
>>
>>
>>
>>
>> Attached is the config file along with the ob file and model file.
>>
>>
>>
>>
>>
>> The following is the output I get. Why the 0 pairs? Thanks Bob
Craig
>>
>>
>>
>>
>>
>> GSL_RNG_TYPE=mt19937
>>
>> GSL_RNG_SEED=3818134092
>>
>> Forecast File:
>> /gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>
>> Climatology File: none
>>
>> Configuration File:
>> /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>
>> Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>>
>>
>>
>> ----------------------------------------
>>
>>
>>
>> Reading records for UGRD/L10.
>>
>> For UGRD/L10 found 1 forecast levels and 0 climatology levels.
>>
>>
>>
>> ----------------------------------------
>>
>>
>>
>> Reading records for VGRD/L10.
>>
>> For VGRD/L10 found 1 forecast levels and 0 climatology levels.
>>
>>
>>
>> ----------------------------------------
>>
>>
>>
>> Searching 10473 observations from 170 PrepBufr messages.
>>
>>
>>
>> ----------------------------------------
>>
>>
>>
>> Processing UGRD/L10 versus UGRD/L10, for observation type ADPSFC,
over
>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>
>>
>>
>> ----------------------------------------
>>
>>
>>
>> Processing VGRD/L10 versus VGRD/L10, for observation type ADPSFC,
over
>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>
>>
>>
>> ----------------------------------------
>>
>>
>>
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
>> at
>>
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>> c.txt
>>
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>> s.txt
>>
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
>> t.txt
>>
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
>> 1l2.txt
>>
>> Output file:
>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
>> 1l2.txt
>>
>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
From: Paul Oldenburg
Time: Wed Jun 23 09:57:59 2010
John,
I'm trying to look into Bob's problems, but I'm struggling with
ncdump. Am I
doing something wrong below where I attempt to view the variable
values?
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -h
obs_2010062200.nc
netcdf obs_2010062200 {
dimensions:
mxstr = 16 ;
hdr_arr_len = 3 ;
obs_arr_len = 5 ;
nobs = UNLIMITED ; // (10004 currently)
nhdr = 162 ;
variables:
float obs_arr(nobs, obs_arr_len) ;
obs_arr:long_name = "array of observation values" ;
obs_arr:_fill_value = -9999.f ;
obs_arr:columns = "hdr_id gc lvl hgt ob" ;
obs_arr:hdr_id_long_name = "index of matching header
data" ;
obs_arr:gc_long_name = "grib code corresponding to the
observation type" ;
obs_arr:lvl_long_name = "pressure level (hPa) or
accumulation
interval (h)" ;
obs_arr:hgt_long_name = "height in meters above sea
level (msl)" ;
obs_arr:ob_long_name = "observation value" ;
char hdr_typ(nhdr, mxstr) ;
hdr_typ:long_name = "message type" ;
char hdr_sid(nhdr, mxstr) ;
hdr_sid:long_name = "station identification" ;
char hdr_vld(nhdr, mxstr) ;
hdr_vld:long_name = "valid time" ;
hdr_vld:units = "YYYYMMDD_HHMMSS" ;
float hdr_arr(nhdr, hdr_arr_len) ;
hdr_arr:long_name = "array of observation station
header values" ;
hdr_arr:_fill_value = -9999.f ;
hdr_arr:columns = "lat lon elv" ;
hdr_arr:lat_long_name = "latitude" ;
hdr_arr:lat_units = "degrees_north" ;
hdr_arr:lon_long_name = "longitude" ;
hdr_arr:lon_units = "degrees_east" ;
hdr_arr:elv_long_name = "elevation" ;
hdr_arr:elv_units = "meters above sea level (msl)" ;
// global attributes:
:FileOrigins = "File
/h/data/global/WXQC/data/obs_prep/2010062200 generated 20100622_184241
UTC on
host a31n09.afwa.af.mil by the PB2NC tool" ;
}
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_arr:long_name
obs_2010062200.nc
ncdump: NetCDF: Variable not found
netcdf obs_2010062200 {
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
obs_2010062200.nc
ncdump: NetCDF: Variable not found
netcdf obs_2010062200 {
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
obs_2010062200.nc
ncdump: NetCDF: Variable not found
netcdf obs_2010062200 {
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_typ:long_name
obs_2010062200.nc
ncdump: NetCDF: Variable not found
netcdf obs_2010062200 {
RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>
> I am getting the same 0 pair error with the attached ob file. I
built this
> file using the attached config file. Is there surface observations
in my
> file? Is the valid time of the file correct? I am using the -dump
option but
> I am not getting any file showing up in the directory I set for
this.
>
> Bob
>
> -----Original Message-----
> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent: Tuesday, June 22, 2010 3:30 PM
> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>
> Bob,
>
> We looked more closely at the observation data that you included in
your
> original email, and we noticed that they don't appear to be surface
> observations. When verifying 2m temp or 10m winds in point_stat,
only
> observations of type ADPSFC or SFCSHP can be used (or ONLYSF for
either).
>
> Please change your point_stat config fcst_field setting to the
following
> value:
>
> fcst_field[] = [ "33/Z10", "34/Z10"];
>
> also, please change the message type:
>
> message_type[] = [ "ADPSFC" ];
>
> Please let us know if you have any more questions.
>
> 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=38629 >
>>
>> Thanks for looking into this. I inadvertently left off my config
file name
>> from the command string, hence it defaulted to profiler data. My
config file
>> has "ANYSFC" which the documentation says includes ADPUPA (RAOB
data).
>>
>> For the observations, I am pulling the prepbufr files from
ftp.ncep.noaa.gov
>> in the /pub/data/nccf/com/gfs/prod/gdas.20100615 for the 15 Jun
files for
>> the case below. An example of the filenames I am pulling are
>> gdas2.t12z.prepbufr.nr. I rename these files to 20100615_12 in
this case.
>> Since this file came out of the 20100615 directory I would expect
it to be
>> valid for 20100615 at 12Z. Is there a way I can dump the header of
the
>> prepbufr file to see the date that's in the file?
>>
>> Bob
>>
>>
>> -----Original Message-----
>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent: Friday, June 18, 2010 11:56 AM
>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>
>> Robert,
>>
>> Your observation file has a message type of PROFLR, but your config
file is
>> specifying a message type of ADPSFC. This is one potential
problem.
>> Another
>> problem is that the model output appears to be valid at 20100615_12
and the
>> observations are valid 24 hours later. You could change your
config file to
>> open the valid time window to a larger value using beg_ds and
end_ds.
>> Please
>> let us know if you have any other questions.
>>
>> Thanks,
>>
>> Paul
>>
>>
>>
>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>> Thu Jun 17 15:45:21 2010: Request 38629 was acted upon.
>>> Transaction: Ticket created by craigr at offutt.af.mil
>>> Queue: met_help
>>> Subject: Point Stat Question
>>> Owner: Nobody
>>> Requestors: craigr at offutt.af.mil
>>> Status: new
>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>
>>>
>>> I am running point stat on WRF winds but I am not getting any
stats. My
>>> code to run this is:
>>>
>>>
>>>
>>> /sphome/qcteam/METv2.0/bin/point_stat
>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>> /h/data/global/WXQC/data/obs_pr
>>>
>>> ep/2010061512
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>> -outdir /h/data/global/WXQC/data/met/data -v 2
>>>
>>>
>>>
>>>
>>>
>>> Attached is the config file along with the ob file and model file.
>>>
>>>
>>>
>>>
>>>
>>> The following is the output I get. Why the 0 pairs? Thanks Bob
Craig
>>>
>>>
>>>
>>>
>>>
>>> GSL_RNG_TYPE=mt19937
>>>
>>> GSL_RNG_SEED=3818134092
>>>
>>> Forecast File:
>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>
>>> Climatology File: none
>>>
>>> Configuration File:
>>> /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>
>>> Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>>>
>>>
>>>
>>> ----------------------------------------
>>>
>>>
>>>
>>> Reading records for UGRD/L10.
>>>
>>> For UGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>
>>>
>>>
>>> ----------------------------------------
>>>
>>>
>>>
>>> Reading records for VGRD/L10.
>>>
>>> For VGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>
>>>
>>>
>>> ----------------------------------------
>>>
>>>
>>>
>>> Searching 10473 observations from 170 PrepBufr messages.
>>>
>>>
>>>
>>> ----------------------------------------
>>>
>>>
>>>
>>> Processing UGRD/L10 versus UGRD/L10, for observation type ADPSFC,
over
>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>
>>>
>>>
>>> ----------------------------------------
>>>
>>>
>>>
>>> Processing VGRD/L10 versus VGRD/L10, for observation type ADPSFC,
over
>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>
>>>
>>>
>>> ----------------------------------------
>>>
>>>
>>>
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
>>> at
>>>
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>> c.txt
>>>
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>> s.txt
>>>
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
>>> t.txt
>>>
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
>>> 1l2.txt
>>>
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
>>> 1l2.txt
>>>
>>
>>
>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #38629] Point Stat Question
From: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Time: Thu Jun 24 09:36:39 2010
John, you sent this to me - did you intend someone else?
Bob
-----Original Message-----
From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
Sent: Wednesday, June 23, 2010 10:58 AM
To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Subject: Re:
[rt.rap.ucar.edu #38629] Point Stat Question
John,
I'm trying to
look into Bob's problems, but I'm struggling with ncdump. Am I
doing
something wrong below where I attempt to view the variable values?
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -h
obs_2010062200.nc
netcdf obs_2010062200 {
dimensions:
mxstr
= 16 ;
hdr_arr_len = 3 ;
obs_arr_len = 5 ;
nobs = UNLIMITED ; // (10004 currently)
nhdr = 162 ;
variables:
float obs_arr(nobs, obs_arr_len) ;
obs_arr:long_name = "array of observation values" ;
obs_arr:_fill_value = -9999.f ;
obs_arr:columns =
"hdr_id gc lvl hgt ob" ;
obs_arr:hdr_id_long_name =
"index of matching header data" ;
obs_arr:gc_long_name = "grib code corresponding to the
observation
type" ;
obs_arr:lvl_long_name = "pressure level (hPa)
or accumulation
interval (h)" ;
obs_arr:hgt_long_name = "height in meters above sea level (msl)" ;
obs_arr:ob_long_name = "observation value" ;
char
hdr_typ(nhdr, mxstr) ;
hdr_typ:long_name = "message
type" ;
char hdr_sid(nhdr, mxstr) ;
hdr_sid:long_name = "station identification" ;
char
hdr_vld(nhdr, mxstr) ;
hdr_vld:long_name = "valid
time" ;
hdr_vld:units = "YYYYMMDD_HHMMSS" ;
float hdr_arr(nhdr, hdr_arr_len) ;
hdr_arr:long_name
= "array of observation station header values" ;
hdr_arr:_fill_value = -9999.f ;
hdr_arr:columns =
"lat lon elv" ;
hdr_arr:lat_long_name = "latitude" ;
hdr_arr:lat_units = "degrees_north" ;
hdr_arr:lon_long_name = "longitude" ;
hdr_arr:lon_units = "degrees_east" ;
hdr_arr:elv_long_name = "elevation" ;
hdr_arr:elv_units = "meters above sea level (msl)" ;
// global
attributes:
:FileOrigins = "File
/h/data/global/WXQC/data/obs_prep/2010062200 generated 20100622_184241
UTC on
host a31n09.afwa.af.mil by the PB2NC tool" ;
}
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_arr:long_name
obs_2010062200.nc
ncdump: NetCDF: Variable not
found
netcdf obs_2010062200 {
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
obs_2010062200.nc
ncdump: NetCDF: Variable not
found
netcdf obs_2010062200 {
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
obs_2010062200.nc
ncdump: NetCDF: Variable not
found
netcdf obs_2010062200 {
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_typ:long_name
obs_2010062200.nc
ncdump: NetCDF: Variable not
found
netcdf obs_2010062200 {
RAL HelpDesk {for Craig, Robert J
Civ USAF AFWA 16 WS/WXN} wrote:
> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>
> I am
getting the same 0 pair error with the attached ob file. I built this
> file using the attached config file. Is there surface observations
in my
> file? Is the valid time of the file correct? I am using
the -dump option but
> I am not getting any file showing up in the
directory I set for this.
>
> Bob
>
> -----Original Message-----
> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent: Tuesday, June 22, 2010 3:30 PM
> To: Craig, Robert J Civ USAF
AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat
Question
>
> Bob,
>
> We looked more closely at the observation
data that you included in your
> original email, and we noticed that
they don't appear to be surface
> observations. When verifying 2m
temp or 10m winds in point_stat, only
> observations of type ADPSFC
or SFCSHP can be used (or ONLYSF for either).
>
> Please change
your point_stat config fcst_field setting to the following
> value:
>
> fcst_field[] = [ "33/Z10", "34/Z10"];
>
> also, please change
the message type:
>
> message_type[] = [ "ADPSFC" ];
>
> Please
let us know if you have any more questions.
>
> 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=38629 >
>>
>>
Thanks for looking into this. I inadvertently left off my config file
name
>> from the command string, hence it defaulted to profiler
data. My config file
>> has "ANYSFC" which the documentation says
includes ADPUPA (RAOB data).
>>
>> For the observations, I am
pulling the prepbufr files from ftp.ncep.noaa.gov
>> in the
/pub/data/nccf/com/gfs/prod/gdas.20100615 for the 15 Jun files for
>> the case below. An example of the filenames I am pulling are
>>
gdas2.t12z.prepbufr.nr. I rename these files to 20100615_12 in this
case.
>> Since this file came out of the 20100615 directory I would
expect it to be
>> valid for 20100615 at 12Z. Is there a way I can
dump the header of the
>> prepbufr file to see the date that's in
the file?
>>
>> Bob
>>
>>
>> -----Original Message-----
>> From:
RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent:
Friday, June 18, 2010 11:56 AM
>> To: Craig, Robert J Civ USAF AFWA
16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat
Question
>>
>> Robert,
>>
>> Your observation file has a message
type of PROFLR, but your config file is
>> specifying a message type
of ADPSFC. This is one potential problem.
>> Another
>> problem is
that the model output appears to be valid at 20100615_12 and the
>>
observations are valid 24 hours later. You could change your config
file to
>> open the valid time window to a larger value using beg_ds
and end_ds.
>> Please
>> let us know if you have any other
questions.
>>
>> Thanks,
>>
>> Paul
>>
>>
>>
>> RAL HelpDesk
{for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>> Thu Jun 17
15:45:21 2010: Request 38629 was acted upon.
>>> Transaction: Ticket
created by craigr at offutt.af.mil
>>> Queue: met_help
>>>
Subject: Point Stat Question
>>> Owner: Nobody
>>>
Requestors: craigr at offutt.af.mil
>>> Status: new
>>> Ticket
<URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>
>>>
>>> I am running point stat on WRF winds but I am not getting any
stats. My
>>> code to run this is:
>>>
>>>
>>>
>>>
/sphome/qcteam/METv2.0/bin/point_stat
>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>> /h/data/global/WXQC/data/obs_pr
>>>
>>> ep/2010061512
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>> -outdir
/h/data/global/WXQC/data/met/data -v 2
>>>
>>>
>>>
>>>
>>>
>>>
Attached is the config file along with the ob file and model file.
>>>
>>>
>>>
>>>
>>>
>>> The following is the output I get. Why
the 0 pairs? Thanks Bob Craig
>>>
>>>
>>>
>>>
>>>
>>>
GSL_RNG_TYPE=mt19937
>>>
>>> GSL_RNG_SEED=3818134092
>>>
>>>
Forecast File:
>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>
>>> Climatology File: none
>>>
>>> Configuration File:
>>>
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>
>>>
Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>>>
>>>
>>>
>>> ----------------------------------------
>>>
>>>
>>>
>>> Reading records for UGRD/L10.
>>>
>>> For UGRD/L10 found 1
forecast levels and 0 climatology levels.
>>>
>>>
>>>
>>>
----------------------------------------
>>>
>>>
>>>
>>> Reading
records for VGRD/L10.
>>>
>>> For VGRD/L10 found 1 forecast levels
and 0 climatology levels.
>>>
>>>
>>>
>>>
----------------------------------------
>>>
>>>
>>>
>>> Searching
10473 observations from 170 PrepBufr messages.
>>>
>>>
>>>
>>>
----------------------------------------
>>>
>>>
>>>
>>>
Processing UGRD/L10 versus UGRD/L10, for observation type ADPSFC, over
>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>
>>>
>>>
>>> ----------------------------------------
>>>
>>>
>>>
>>> Processing VGRD/L10 versus VGRD/L10, for observation type ADPSFC,
over
>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>
>>>
>>>
>>> ----------------------------------------
>>>
>>>
>>>
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
>>> at
>>>
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>> c.txt
>>>
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>> s.txt
>>>
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
>>> t.txt
>>>
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
>>> 1l2.txt
>>>
>>> Output file:
>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
>>> 1l2.txt
>>>
>>
>>
>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
From: Paul Oldenburg
Time: Thu Jun 24 10:21:57 2010
Bob,
I sent that email to John through the met_help ticketing system. It
was not
intended for you, sorry.
We tested the PB2NCConfig file that you sent, and we found that it did
not
process the message types that are needed by point_stat. The reason
was that
the quality_mark_thresh setting was set to 1 instead of 2, and the
message types
ADPSFC and ADPUPA were not processed as a result. Please try the
following
setting in your PB2NCConfig file:
quality_mark_thresh = 2;
Also, it is important to make sure that you are running the latest MET
patches.
If you have not done so, please download and deploy the patches from
the MET
website at:
http://www.dtcenter.org/met/users/support/known_issues/METv2.0/index.php
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=38629 >
>
> John, you sent this to me - did you intend someone else?
>
> Bob
>
> -----Original Message-----
> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent: Wednesday, June 23, 2010 10:58 AM
> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>
> John,
>
> I'm trying to look into Bob's problems, but I'm struggling with
ncdump. Am I
> doing something wrong below where I attempt to view the variable
values?
>
> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -h
obs_2010062200.nc
> netcdf obs_2010062200 {
> dimensions:
> mxstr = 16 ;
> hdr_arr_len = 3 ;
> obs_arr_len = 5 ;
> nobs = UNLIMITED ; // (10004 currently)
> nhdr = 162 ;
> variables:
> float obs_arr(nobs, obs_arr_len) ;
> obs_arr:long_name = "array of observation values" ;
> obs_arr:_fill_value = -9999.f ;
> obs_arr:columns = "hdr_id gc lvl hgt ob" ;
> obs_arr:hdr_id_long_name = "index of matching header
data" ;
> obs_arr:gc_long_name = "grib code corresponding to
the
> observation type" ;
> obs_arr:lvl_long_name = "pressure level (hPa) or
accumulation
> interval (h)" ;
> obs_arr:hgt_long_name = "height in meters above sea
level (msl)" ;
> obs_arr:ob_long_name = "observation value" ;
> char hdr_typ(nhdr, mxstr) ;
> hdr_typ:long_name = "message type" ;
> char hdr_sid(nhdr, mxstr) ;
> hdr_sid:long_name = "station identification" ;
> char hdr_vld(nhdr, mxstr) ;
> hdr_vld:long_name = "valid time" ;
> hdr_vld:units = "YYYYMMDD_HHMMSS" ;
> float hdr_arr(nhdr, hdr_arr_len) ;
> hdr_arr:long_name = "array of observation station
header values" ;
> hdr_arr:_fill_value = -9999.f ;
> hdr_arr:columns = "lat lon elv" ;
> hdr_arr:lat_long_name = "latitude" ;
> hdr_arr:lat_units = "degrees_north" ;
> hdr_arr:lon_long_name = "longitude" ;
> hdr_arr:lon_units = "degrees_east" ;
> hdr_arr:elv_long_name = "elevation" ;
> hdr_arr:elv_units = "meters above sea level (msl)" ;
>
> // global attributes:
> :FileOrigins = "File
> /h/data/global/WXQC/data/obs_prep/2010062200 generated
20100622_184241 UTC on
> host a31n09.afwa.af.mil by the PB2NC tool" ;
> }
> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_arr:long_name
> obs_2010062200.nc
> ncdump: NetCDF: Variable not found
> netcdf obs_2010062200 {
> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
> obs_2010062200.nc
> ncdump: NetCDF: Variable not found
> netcdf obs_2010062200 {
> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
> obs_2010062200.nc
> ncdump: NetCDF: Variable not found
> netcdf obs_2010062200 {
> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_typ:long_name
> obs_2010062200.nc
> ncdump: NetCDF: Variable not found
> netcdf obs_2010062200 {
>
>
> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>
>> I am getting the same 0 pair error with the attached ob file. I
built this
>> file using the attached config file. Is there surface observations
in my
>> file? Is the valid time of the file correct? I am using the -dump
option but
>> I am not getting any file showing up in the directory I set for
this.
>>
>> Bob
>>
>> -----Original Message-----
>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent: Tuesday, June 22, 2010 3:30 PM
>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>
>> Bob,
>>
>> We looked more closely at the observation data that you included in
your
>> original email, and we noticed that they don't appear to be surface
>> observations. When verifying 2m temp or 10m winds in point_stat,
only
>> observations of type ADPSFC or SFCSHP can be used (or ONLYSF for
either).
>>
>> Please change your point_stat config fcst_field setting to the
following
>> value:
>>
>> fcst_field[] = [ "33/Z10", "34/Z10"];
>>
>> also, please change the message type:
>>
>> message_type[] = [ "ADPSFC" ];
>>
>> Please let us know if you have any more questions.
>>
>> 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=38629 >
>>>
>>> Thanks for looking into this. I inadvertently left off my config
file name
>>> from the command string, hence it defaulted to profiler data. My
config file
>>> has "ANYSFC" which the documentation says includes ADPUPA (RAOB
data).
>>>
>>> For the observations, I am pulling the prepbufr files from
ftp.ncep.noaa.gov
>>> in the /pub/data/nccf/com/gfs/prod/gdas.20100615 for the 15 Jun
files for
>>> the case below. An example of the filenames I am pulling are
>>> gdas2.t12z.prepbufr.nr. I rename these files to 20100615_12 in
this case.
>>> Since this file came out of the 20100615 directory I would expect
it to be
>>> valid for 20100615 at 12Z. Is there a way I can dump the header
of the
>>> prepbufr file to see the date that's in the file?
>>>
>>> Bob
>>>
>>>
>>> -----Original Message-----
>>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>>> Sent: Friday, June 18, 2010 11:56 AM
>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>
>>> Robert,
>>>
>>> Your observation file has a message type of PROFLR, but your
config file is
>>> specifying a message type of ADPSFC. This is one potential
problem.
>>> Another
>>> problem is that the model output appears to be valid at
20100615_12 and the
>>> observations are valid 24 hours later. You could change your
config file to
>>> open the valid time window to a larger value using beg_ds and
end_ds.
>>> Please
>>> let us know if you have any other questions.
>>>
>>> Thanks,
>>>
>>> Paul
>>>
>>>
>>>
>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>>> Thu Jun 17 15:45:21 2010: Request 38629 was acted upon.
>>>> Transaction: Ticket created by craigr at offutt.af.mil
>>>> Queue: met_help
>>>> Subject: Point Stat Question
>>>> Owner: Nobody
>>>> Requestors: craigr at offutt.af.mil
>>>> Status: new
>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>
>>>>
>>>> I am running point stat on WRF winds but I am not getting any
stats. My
>>>> code to run this is:
>>>>
>>>>
>>>>
>>>> /sphome/qcteam/METv2.0/bin/point_stat
>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>> /h/data/global/WXQC/data/obs_pr
>>>>
>>>> ep/2010061512
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>> -outdir /h/data/global/WXQC/data/met/data -v 2
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> Attached is the config file along with the ob file and model
file.
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> The following is the output I get. Why the 0 pairs? Thanks Bob
Craig
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> GSL_RNG_TYPE=mt19937
>>>>
>>>> GSL_RNG_SEED=3818134092
>>>>
>>>> Forecast File:
>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>
>>>> Climatology File: none
>>>>
>>>> Configuration File:
>>>> /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>
>>>> Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>>>>
>>>>
>>>>
>>>> ----------------------------------------
>>>>
>>>>
>>>>
>>>> Reading records for UGRD/L10.
>>>>
>>>> For UGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>
>>>>
>>>>
>>>> ----------------------------------------
>>>>
>>>>
>>>>
>>>> Reading records for VGRD/L10.
>>>>
>>>> For VGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>
>>>>
>>>>
>>>> ----------------------------------------
>>>>
>>>>
>>>>
>>>> Searching 10473 observations from 170 PrepBufr messages.
>>>>
>>>>
>>>>
>>>> ----------------------------------------
>>>>
>>>>
>>>>
>>>> Processing UGRD/L10 versus UGRD/L10, for observation type ADPSFC,
over
>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>
>>>>
>>>>
>>>> ----------------------------------------
>>>>
>>>>
>>>>
>>>> Processing VGRD/L10 versus VGRD/L10, for observation type ADPSFC,
over
>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>
>>>>
>>>>
>>>> ----------------------------------------
>>>>
>>>>
>>>>
>>>> Output file:
>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
>>>> at
>>>>
>>>> Output file:
>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>> c.txt
>>>>
>>>> Output file:
>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>> s.txt
>>>>
>>>> Output file:
>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
>>>> t.txt
>>>>
>>>> Output file:
>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
>>>> 1l2.txt
>>>>
>>>> Output file:
>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
>>>> 1l2.txt
>>>>
>>>
>>
>>
>
>
------------------------------------------------
Subject: Point Stat Question
From: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Time: Thu Jun 24 12:28:17 2010
I will reset the quality threshold.
I download the updates to MET and placed in the top level MET dir.
When I run
a Make clean I get the following errors:
[496 qcteam a31n09 /sphome/qcteam/METv2.0] make clean
"Makefile", line 191: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 193: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 194: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 196: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 197: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 199: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 200: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 202: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 203: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 205: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 206: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 208: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 209: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 211: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 212: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 214: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 215: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 217: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 218: make: 1254-055 Dependency line needs colon or
double
colon operator.
"Makefile", line 220: make: 1254-055 Dependency line needs colon or
double
colon operator.
make: 1254-058 Fatal errors encountered -- cannot continue.
The make file is attached. What is causing this error?
Thanks
Bob
-----Original Message-----
From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
Sent: Thursday, June 24, 2010 11:22 AM
To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
Bob,
I sent that email to John through the met_help ticketing system. It
was not
intended for you, sorry.
We tested the PB2NCConfig file that you sent, and we found that it did
not
process the message types that are needed by point_stat. The reason
was that
the quality_mark_thresh setting was set to 1 instead of 2, and the
message
types
ADPSFC and ADPUPA were not processed as a result. Please try the
following
setting in your PB2NCConfig file:
quality_mark_thresh = 2;
Also, it is important to make sure that you are running the latest MET
patches.
If you have not done so, please download and deploy the patches from
the MET
website at:
http://www.dtcenter.org/met/users/support/known_issues/METv2.0/index.php
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=38629 >
>
> John, you sent this to me - did you intend someone else?
>
> Bob
>
> -----Original Message-----
> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent: Wednesday, June 23, 2010 10:58 AM
> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>
> John,
>
> I'm trying to look into Bob's problems, but I'm struggling with
ncdump. Am
> I
> doing something wrong below where I attempt to view the variable
values?
>
> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -h
obs_2010062200.nc
> netcdf obs_2010062200 {
> dimensions:
> mxstr = 16 ;
> hdr_arr_len = 3 ;
> obs_arr_len = 5 ;
> nobs = UNLIMITED ; // (10004 currently)
> nhdr = 162 ;
> variables:
> float obs_arr(nobs, obs_arr_len) ;
> obs_arr:long_name = "array of observation values" ;
> obs_arr:_fill_value = -9999.f ;
> obs_arr:columns = "hdr_id gc lvl hgt ob" ;
> obs_arr:hdr_id_long_name = "index of matching header
data" ;
> obs_arr:gc_long_name = "grib code corresponding to
the
> observation type" ;
> obs_arr:lvl_long_name = "pressure level (hPa) or
> accumulation
> interval (h)" ;
> obs_arr:hgt_long_name = "height in meters above sea
level
> (msl)" ;
> obs_arr:ob_long_name = "observation value" ;
> char hdr_typ(nhdr, mxstr) ;
> hdr_typ:long_name = "message type" ;
> char hdr_sid(nhdr, mxstr) ;
> hdr_sid:long_name = "station identification" ;
> char hdr_vld(nhdr, mxstr) ;
> hdr_vld:long_name = "valid time" ;
> hdr_vld:units = "YYYYMMDD_HHMMSS" ;
> float hdr_arr(nhdr, hdr_arr_len) ;
> hdr_arr:long_name = "array of observation station
header
> values" ;
> hdr_arr:_fill_value = -9999.f ;
> hdr_arr:columns = "lat lon elv" ;
> hdr_arr:lat_long_name = "latitude" ;
> hdr_arr:lat_units = "degrees_north" ;
> hdr_arr:lon_long_name = "longitude" ;
> hdr_arr:lon_units = "degrees_east" ;
> hdr_arr:elv_long_name = "elevation" ;
> hdr_arr:elv_units = "meters above sea level (msl)" ;
>
> // global attributes:
> :FileOrigins = "File
> /h/data/global/WXQC/data/obs_prep/2010062200 generated
20100622_184241 UTC
> on
> host a31n09.afwa.af.mil by the PB2NC tool" ;
> }
> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_arr:long_name
> obs_2010062200.nc
> ncdump: NetCDF: Variable not found
> netcdf obs_2010062200 {
> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
> obs_2010062200.nc
> ncdump: NetCDF: Variable not found
> netcdf obs_2010062200 {
> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
> obs_2010062200.nc
> ncdump: NetCDF: Variable not found
> netcdf obs_2010062200 {
> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_typ:long_name
> obs_2010062200.nc
> ncdump: NetCDF: Variable not found
> netcdf obs_2010062200 {
>
>
> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>
>> I am getting the same 0 pair error with the attached ob file. I
built this
>> file using the attached config file. Is there surface observations
in my
>> file? Is the valid time of the file correct? I am using the -dump
option
>> but
>> I am not getting any file showing up in the directory I set for
this.
>>
>> Bob
>>
>> -----Original Message-----
>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent: Tuesday, June 22, 2010 3:30 PM
>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>
>> Bob,
>>
>> We looked more closely at the observation data that you included in
your
>> original email, and we noticed that they don't appear to be surface
>> observations. When verifying 2m temp or 10m winds in point_stat,
only
>> observations of type ADPSFC or SFCSHP can be used (or ONLYSF for
either).
>>
>> Please change your point_stat config fcst_field setting to the
following
>> value:
>>
>> fcst_field[] = [ "33/Z10", "34/Z10"];
>>
>> also, please change the message type:
>>
>> message_type[] = [ "ADPSFC" ];
>>
>> Please let us know if you have any more questions.
>>
>> 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=38629 >
>>>
>>> Thanks for looking into this. I inadvertently left off my config
file
>>> name
>>> from the command string, hence it defaulted to profiler data. My
config
>>> file
>>> has "ANYSFC" which the documentation says includes ADPUPA (RAOB
data).
>>>
>>> For the observations, I am pulling the prepbufr files from
>>> ftp.ncep.noaa.gov
>>> in the /pub/data/nccf/com/gfs/prod/gdas.20100615 for the 15 Jun
files for
>>> the case below. An example of the filenames I am pulling are
>>> gdas2.t12z.prepbufr.nr. I rename these files to 20100615_12 in
this
>>> case.
>>> Since this file came out of the 20100615 directory I would expect
it to be
>>> valid for 20100615 at 12Z. Is there a way I can dump the header
of the
>>> prepbufr file to see the date that's in the file?
>>>
>>> Bob
>>>
>>>
>>> -----Original Message-----
>>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>>> Sent: Friday, June 18, 2010 11:56 AM
>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>
>>> Robert,
>>>
>>> Your observation file has a message type of PROFLR, but your
config file
>>> is
>>> specifying a message type of ADPSFC. This is one potential
problem.
>>> Another
>>> problem is that the model output appears to be valid at
20100615_12 and
>>> the
>>> observations are valid 24 hours later. You could change your
config file
>>> to
>>> open the valid time window to a larger value using beg_ds and
end_ds.
>>> Please
>>> let us know if you have any other questions.
>>>
>>> Thanks,
>>>
>>> Paul
>>>
>>>
>>>
>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>>> Thu Jun 17 15:45:21 2010: Request 38629 was acted upon.
>>>> Transaction: Ticket created by craigr at offutt.af.mil
>>>> Queue: met_help
>>>> Subject: Point Stat Question
>>>> Owner: Nobody
>>>> Requestors: craigr at offutt.af.mil
>>>> Status: new
>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>
>>>>
>>>> I am running point stat on WRF winds but I am not getting any
stats. My
>>>> code to run this is:
>>>>
>>>>
>>>>
>>>> /sphome/qcteam/METv2.0/bin/point_stat
>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>> /h/data/global/WXQC/data/obs_pr
>>>>
>>>> ep/2010061512
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>> -outdir /h/data/global/WXQC/data/met/data -v 2
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> Attached is the config file along with the ob file and model
file.
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> The following is the output I get. Why the 0 pairs? Thanks Bob
Craig
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> GSL_RNG_TYPE=mt19937
>>>>
>>>> GSL_RNG_SEED=3818134092
>>>>
>>>> Forecast File:
>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>
>>>> Climatology File: none
>>>>
>>>> Configuration File:
>>>> /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>
>>>> Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>>>>
>>>>
>>>>
>>>> ----------------------------------------
>>>>
>>>>
>>>>
>>>> Reading records for UGRD/L10.
>>>>
>>>> For UGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>
>>>>
>>>>
>>>> ----------------------------------------
>>>>
>>>>
>>>>
>>>> Reading records for VGRD/L10.
>>>>
>>>> For VGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>
>>>>
>>>>
>>>> ----------------------------------------
>>>>
>>>>
>>>>
>>>> Searching 10473 observations from 170 PrepBufr messages.
>>>>
>>>>
>>>>
>>>> ----------------------------------------
>>>>
>>>>
>>>>
>>>> Processing UGRD/L10 versus UGRD/L10, for observation type ADPSFC,
over
>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>
>>>>
>>>>
>>>> ----------------------------------------
>>>>
>>>>
>>>>
>>>> Processing VGRD/L10 versus VGRD/L10, for observation type ADPSFC,
over
>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>
>>>>
>>>>
>>>> ----------------------------------------
>>>>
>>>>
>>>>
>>>> Output file:
>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
>>>> at
>>>>
>>>> Output file:
>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>> c.txt
>>>>
>>>> Output file:
>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>> s.txt
>>>>
>>>> Output file:
>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
>>>> t.txt
>>>>
>>>> Output file:
>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
>>>> 1l2.txt
>>>>
>>>> Output file:
>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
>>>> 1l2.txt
>>>>
>>>
>>
>>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
From: Paul Oldenburg
Time: Thu Jun 24 13:00:23 2010
Bob,
These errors look like errors that are output by the IBM version of
make. MET
needs to be built with GNU make, which is called gmake on our IBM
machine.
Please look on your IBM machine to see if you have gmake installed and
use it
instead.
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=38629 >
>
> I will reset the quality threshold.
>
> I download the updates to MET and placed in the top level MET dir.
When I run
> a Make clean I get the following errors:
>
> [496 qcteam a31n09 /sphome/qcteam/METv2.0] make clean
> "Makefile", line 191: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 193: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 194: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 196: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 197: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 199: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 200: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 202: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 203: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 205: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 206: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 208: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 209: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 211: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 212: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 214: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 215: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 217: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 218: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> "Makefile", line 220: make: 1254-055 Dependency line needs colon or
double
> colon operator.
> make: 1254-058 Fatal errors encountered -- cannot continue.
>
> The make file is attached. What is causing this error?
>
> Thanks
>
> Bob
>
>
>
> -----Original Message-----
> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent: Thursday, June 24, 2010 11:22 AM
> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>
> Bob,
>
> I sent that email to John through the met_help ticketing system. It
was not
> intended for you, sorry.
>
> We tested the PB2NCConfig file that you sent, and we found that it
did not
> process the message types that are needed by point_stat. The reason
was that
> the quality_mark_thresh setting was set to 1 instead of 2, and the
message
> types
> ADPSFC and ADPUPA were not processed as a result. Please try the
following
> setting in your PB2NCConfig file:
>
> quality_mark_thresh = 2;
>
> Also, it is important to make sure that you are running the latest
MET
> patches.
> If you have not done so, please download and deploy the patches
from the MET
> website at:
>
>
http://www.dtcenter.org/met/users/support/known_issues/METv2.0/index.php
>
> 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=38629 >
>>
>> John, you sent this to me - did you intend someone else?
>>
>> Bob
>>
>> -----Original Message-----
>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent: Wednesday, June 23, 2010 10:58 AM
>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>
>> John,
>>
>> I'm trying to look into Bob's problems, but I'm struggling with
ncdump. Am
>> I
>> doing something wrong below where I attempt to view the variable
values?
>>
>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -h
obs_2010062200.nc
>> netcdf obs_2010062200 {
>> dimensions:
>> mxstr = 16 ;
>> hdr_arr_len = 3 ;
>> obs_arr_len = 5 ;
>> nobs = UNLIMITED ; // (10004 currently)
>> nhdr = 162 ;
>> variables:
>> float obs_arr(nobs, obs_arr_len) ;
>> obs_arr:long_name = "array of observation values" ;
>> obs_arr:_fill_value = -9999.f ;
>> obs_arr:columns = "hdr_id gc lvl hgt ob" ;
>> obs_arr:hdr_id_long_name = "index of matching
header data" ;
>> obs_arr:gc_long_name = "grib code corresponding to
the
>> observation type" ;
>> obs_arr:lvl_long_name = "pressure level (hPa) or
>> accumulation
>> interval (h)" ;
>> obs_arr:hgt_long_name = "height in meters above sea
level
>> (msl)" ;
>> obs_arr:ob_long_name = "observation value" ;
>> char hdr_typ(nhdr, mxstr) ;
>> hdr_typ:long_name = "message type" ;
>> char hdr_sid(nhdr, mxstr) ;
>> hdr_sid:long_name = "station identification" ;
>> char hdr_vld(nhdr, mxstr) ;
>> hdr_vld:long_name = "valid time" ;
>> hdr_vld:units = "YYYYMMDD_HHMMSS" ;
>> float hdr_arr(nhdr, hdr_arr_len) ;
>> hdr_arr:long_name = "array of observation station
header
>> values" ;
>> hdr_arr:_fill_value = -9999.f ;
>> hdr_arr:columns = "lat lon elv" ;
>> hdr_arr:lat_long_name = "latitude" ;
>> hdr_arr:lat_units = "degrees_north" ;
>> hdr_arr:lon_long_name = "longitude" ;
>> hdr_arr:lon_units = "degrees_east" ;
>> hdr_arr:elv_long_name = "elevation" ;
>> hdr_arr:elv_units = "meters above sea level (msl)"
;
>>
>> // global attributes:
>> :FileOrigins = "File
>> /h/data/global/WXQC/data/obs_prep/2010062200 generated
20100622_184241 UTC
>> on
>> host a31n09.afwa.af.mil by the PB2NC tool" ;
>> }
>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_arr:long_name
>> obs_2010062200.nc
>> ncdump: NetCDF: Variable not found
>> netcdf obs_2010062200 {
>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>> obs_2010062200.nc
>> ncdump: NetCDF: Variable not found
>> netcdf obs_2010062200 {
>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>> obs_2010062200.nc
>> ncdump: NetCDF: Variable not found
>> netcdf obs_2010062200 {
>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_typ:long_name
>> obs_2010062200.nc
>> ncdump: NetCDF: Variable not found
>> netcdf obs_2010062200 {
>>
>>
>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>
>>> I am getting the same 0 pair error with the attached ob file. I
built this
>>> file using the attached config file. Is there surface
observations in my
>>> file? Is the valid time of the file correct? I am using the
-dump option
>>> but
>>> I am not getting any file showing up in the directory I set for
this.
>>>
>>> Bob
>>>
>>> -----Original Message-----
>>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>>> Sent: Tuesday, June 22, 2010 3:30 PM
>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>
>>> Bob,
>>>
>>> We looked more closely at the observation data that you included
in your
>>> original email, and we noticed that they don't appear to be
surface
>>> observations. When verifying 2m temp or 10m winds in point_stat,
only
>>> observations of type ADPSFC or SFCSHP can be used (or ONLYSF for
either).
>>>
>>> Please change your point_stat config fcst_field setting to the
following
>>> value:
>>>
>>> fcst_field[] = [ "33/Z10", "34/Z10"];
>>>
>>> also, please change the message type:
>>>
>>> message_type[] = [ "ADPSFC" ];
>>>
>>> Please let us know if you have any more questions.
>>>
>>> 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=38629 >
>>>>
>>>> Thanks for looking into this. I inadvertently left off my config
file
>>>> name
>>>> from the command string, hence it defaulted to profiler data. My
config
>>>> file
>>>> has "ANYSFC" which the documentation says includes ADPUPA (RAOB
data).
>>>>
>>>> For the observations, I am pulling the prepbufr files from
>>>> ftp.ncep.noaa.gov
>>>> in the /pub/data/nccf/com/gfs/prod/gdas.20100615 for the 15 Jun
files for
>>>> the case below. An example of the filenames I am pulling are
>>>> gdas2.t12z.prepbufr.nr. I rename these files to 20100615_12 in
this
>>>> case.
>>>> Since this file came out of the 20100615 directory I would expect
it to be
>>>> valid for 20100615 at 12Z. Is there a way I can dump the header
of the
>>>> prepbufr file to see the date that's in the file?
>>>>
>>>> Bob
>>>>
>>>>
>>>> -----Original Message-----
>>>> From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
>>>> Sent: Friday, June 18, 2010 11:56 AM
>>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>
>>>> Robert,
>>>>
>>>> Your observation file has a message type of PROFLR, but your
config file
>>>> is
>>>> specifying a message type of ADPSFC. This is one potential
problem.
>>>> Another
>>>> problem is that the model output appears to be valid at
20100615_12 and
>>>> the
>>>> observations are valid 24 hours later. You could change your
config file
>>>> to
>>>> open the valid time window to a larger value using beg_ds and
end_ds.
>>>> Please
>>>> let us know if you have any other questions.
>>>>
>>>> Thanks,
>>>>
>>>> Paul
>>>>
>>>>
>>>>
>>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>>>> Thu Jun 17 15:45:21 2010: Request 38629 was acted upon.
>>>>> Transaction: Ticket created by craigr at offutt.af.mil
>>>>> Queue: met_help
>>>>> Subject: Point Stat Question
>>>>> Owner: Nobody
>>>>> Requestors: craigr at offutt.af.mil
>>>>> Status: new
>>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>>
>>>>>
>>>>> I am running point stat on WRF winds but I am not getting any
stats. My
>>>>> code to run this is:
>>>>>
>>>>>
>>>>>
>>>>> /sphome/qcteam/METv2.0/bin/point_stat
>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>> /h/data/global/WXQC/data/obs_pr
>>>>>
>>>>> ep/2010061512
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>> -outdir /h/data/global/WXQC/data/met/data -v 2
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Attached is the config file along with the ob file and model
file.
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> The following is the output I get. Why the 0 pairs? Thanks
Bob Craig
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> GSL_RNG_TYPE=mt19937
>>>>>
>>>>> GSL_RNG_SEED=3818134092
>>>>>
>>>>> Forecast File:
>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>>
>>>>> Climatology File: none
>>>>>
>>>>> Configuration File:
>>>>> /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>
>>>>> Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>>>>>
>>>>>
>>>>>
>>>>> ----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>> Reading records for UGRD/L10.
>>>>>
>>>>> For UGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>
>>>>>
>>>>>
>>>>> ----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>> Reading records for VGRD/L10.
>>>>>
>>>>> For VGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>
>>>>>
>>>>>
>>>>> ----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>> Searching 10473 observations from 170 PrepBufr messages.
>>>>>
>>>>>
>>>>>
>>>>> ----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>> Processing UGRD/L10 versus UGRD/L10, for observation type
ADPSFC, over
>>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>>
>>>>>
>>>>>
>>>>> ----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>> Processing VGRD/L10 versus VGRD/L10, for observation type
ADPSFC, over
>>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>>
>>>>>
>>>>>
>>>>> ----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
>>>>> at
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>> c.txt
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>> s.txt
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
>>>>> t.txt
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
>>>>> 1l2.txt
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
>>>>> 1l2.txt
>>>>>
>>>
>>
>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #38629] Point Stat Question
From: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Time: Thu Jun 24 13:20:20 2010
Gmake solved that problem - compiled with the following errors (these
were at the bottom of the listing):
ld: 0711-224 WARNING: Duplicate
symbol: .strlenFPCc
ld: 0711-345 Use the -bloadmap or -bnoquiet
option to obtain more information.
ld: 0711-317 ERROR: Undefined
symbol: .openbf_
ld: 0711-317 ERROR: Undefined symbol: .datelen_
ld:
0711-317 ERROR: Undefined symbol: .ireadns_
ld: 0711-317 ERROR:
Undefined symbol: .ufbint_
ld: 0711-317 ERROR: Undefined symbol:
.ufbevn_
ld: 0711-317 ERROR: Undefined symbol: .ufbtab_
gmake[3]:
*** [pb2nc] Error 8
gmake[3]: Leaving directory
`/h/sphome/qcteam/METv2.0/src/pb2nc'
gmake[2]: *** [all] Error 2
gmake[2]: Leaving directory `/h/sphome/qcteam/METv2.0/src/pb2nc'
gmake[1]: *** [targets] Error 2
gmake[1]: Leaving directory
`/h/sphome/qcteam/METv2.0'
gmake: *** [all] Error 2
Should I worry
about these - what is the fix?
Thanks
-----Original Message-----
From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
Sent: Thursday, June 24, 2010 2:00 PM
To: Craig, Robert J Civ USAF
AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat
Question
Bob,
These errors look like errors that are output by
the IBM version of make. MET
needs to be built with GNU make, which
is called gmake on our IBM machine.
Please look on your IBM machine
to see if you have gmake installed and use it
instead.
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=38629 >
>
> I will
reset the quality threshold.
>
> I download the updates to MET and
placed in the top level MET dir. When I run
> a Make clean I get the
following errors:
>
> [496 qcteam a31n09 /sphome/qcteam/METv2.0]
make clean
> "Makefile", line 191: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 193:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 194: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 196:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 197: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 199:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 200: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 202:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 203: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 205:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 206: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 208:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 209: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 211:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 212: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 214:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 215: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 217:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 218: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 220:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> make: 1254-058 Fatal errors encountered -- cannot
continue.
>
> The make file is attached. What is causing this
error?
>
> Thanks
>
> Bob
>
>
>
> -----Original Message-----
> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent: Thursday, June 24, 2010 11:22 AM
> To: Craig, Robert J Civ
USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38629] Point
Stat Question
>
> Bob,
>
> I sent that email to John through the
met_help ticketing system. It was not
> intended for you, sorry.
>
> We tested the PB2NCConfig file that you sent, and we found that it
did not
> process the message types that are needed by point_stat.
The reason was that
> the quality_mark_thresh setting was set to 1
instead of 2, and the message
> types
> ADPSFC and ADPUPA were not
processed as a result. Please try the following
> setting in your
PB2NCConfig file:
>
> quality_mark_thresh = 2;
>
> Also, it is
important to make sure that you are running the latest MET
> patches.
> If you have not done so, please download and deploy the patches
from the MET
> website at:
>
>
http://www.dtcenter.org/met/users/support/known_issues/METv2.0/index.php
>
> 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=38629 >
>>
>>
John, you sent this to me - did you intend someone else?
>>
>> Bob
>>
>> -----Original Message-----
>> From: RAL HelpDesk {for Paul
Oldenburg} [mailto:met_help at ucar.edu]
>> Sent: Wednesday, June 23,
2010 10:58 AM
>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>
>>
John,
>>
>> I'm trying to look into Bob's problems, but I'm
struggling with ncdump. Am
>> I
>> doing something wrong below
where I attempt to view the variable values?
>>
>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -h
obs_2010062200.nc
>> netcdf obs_2010062200 {
>> dimensions:
>>
mxstr = 16 ;
>> hdr_arr_len = 3 ;
>> obs_arr_len = 5
;
>> nobs = UNLIMITED ; // (10004 currently)
>> nhdr
= 162 ;
>> variables:
>> float obs_arr(nobs, obs_arr_len) ;
>> obs_arr:long_name = "array of observation values" ;
>> obs_arr:_fill_value = -9999.f ;
>>
obs_arr:columns = "hdr_id gc lvl hgt ob" ;
>>
obs_arr:hdr_id_long_name = "index of matching header data" ;
>>
obs_arr:gc_long_name = "grib code corresponding to the
>> observation
type" ;
>> obs_arr:lvl_long_name = "pressure level
(hPa) or
>> accumulation
>> interval (h)" ;
>>
obs_arr:hgt_long_name = "height in meters above sea level
>> (msl)" ;
>> obs_arr:ob_long_name = "observation value" ;
>>
char hdr_typ(nhdr, mxstr) ;
>> hdr_typ:long_name =
"message type" ;
>> char hdr_sid(nhdr, mxstr) ;
>>
hdr_sid:long_name = "station identification" ;
>> char
hdr_vld(nhdr, mxstr) ;
>> hdr_vld:long_name = "valid
time" ;
>> hdr_vld:units = "YYYYMMDD_HHMMSS" ;
>>
float hdr_arr(nhdr, hdr_arr_len) ;
>>
hdr_arr:long_name = "array of observation station header
>> values" ;
>> hdr_arr:_fill_value = -9999.f ;
>>
hdr_arr:columns = "lat lon elv" ;
>>
hdr_arr:lat_long_name = "latitude" ;
>>
hdr_arr:lat_units = "degrees_north" ;
>>
hdr_arr:lon_long_name = "longitude" ;
>>
hdr_arr:lon_units = "degrees_east" ;
>>
hdr_arr:elv_long_name = "elevation" ;
>>
hdr_arr:elv_units = "meters above sea level (msl)" ;
>>
>> // global
attributes:
>> :FileOrigins = "File
>>
/h/data/global/WXQC/data/obs_prep/2010062200 generated 20100622_184241
UTC
>> on
>> host a31n09.afwa.af.mil by the PB2NC tool" ;
>> }
>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_arr:long_name
>> obs_2010062200.nc
>> ncdump: NetCDF: Variable
not found
>> netcdf obs_2010062200 {
>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>> obs_2010062200.nc
>> ncdump: NetCDF:
Variable not found
>> netcdf obs_2010062200 {
>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>> obs_2010062200.nc
>> ncdump: NetCDF:
Variable not found
>> netcdf obs_2010062200 {
>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_typ:long_name
>> obs_2010062200.nc
>> ncdump: NetCDF: Variable
not found
>> netcdf obs_2010062200 {
>>
>>
>> RAL HelpDesk {for
Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>
>>> I
am getting the same 0 pair error with the attached ob file. I built
this
>>> file using the attached config file. Is there surface
observations in my
>>> file? Is the valid time of the file correct?
I am using the -dump option
>>> but
>>> I am not getting any file
showing up in the directory I set for this.
>>>
>>> Bob
>>>
>>>
-----Original Message-----
>>> From: RAL HelpDesk {for Paul
Oldenburg} [mailto:met_help at ucar.edu]
>>> Sent: Tuesday, June 22,
2010 3:30 PM
>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>
>>>
Bob,
>>>
>>> We looked more closely at the observation data that you
included in your
>>> original email, and we noticed that they don't
appear to be surface
>>> observations. When verifying 2m temp or 10m
winds in point_stat, only
>>> observations of type ADPSFC or SFCSHP
can be used (or ONLYSF for either).
>>>
>>> Please change your
point_stat config fcst_field setting to the following
>>> value:
>>>
>>> fcst_field[] = [ "33/Z10", "34/Z10"];
>>>
>>> also, please
change the message type:
>>>
>>> message_type[] = [ "ADPSFC" ];
>>>
>>> Please let us know if you have any more questions.
>>>
>>>
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=38629 >
>>>>
>>>>
Thanks for looking into this. I inadvertently left off my config file
>>>> name
>>>> from the command string, hence it defaulted to
profiler data. My config
>>>> file
>>>> has "ANYSFC" which the
documentation says includes ADPUPA (RAOB data).
>>>>
>>>> For the
observations, I am pulling the prepbufr files from
>>>>
ftp.ncep.noaa.gov
>>>> in the
/pub/data/nccf/com/gfs/prod/gdas.20100615 for the 15 Jun files for
>>>> the case below. An example of the filenames I am pulling are
>>>> gdas2.t12z.prepbufr.nr. I rename these files to 20100615_12 in
this
>>>> case.
>>>> Since this file came out of the 20100615
directory I would expect it to be
>>>> valid for 20100615 at 12Z. Is
there a way I can dump the header of the
>>>> prepbufr file to see
the date that's in the file?
>>>>
>>>> Bob
>>>>
>>>>
>>>>
-----Original Message-----
>>>> From: RAL HelpDesk {for Paul
Oldenburg} [mailto:met_help at ucar.edu]
>>>> Sent: Friday, June 18,
2010 11:56 AM
>>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>>
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>
>>>>
Robert,
>>>>
>>>> Your observation file has a message type of
PROFLR, but your config file
>>>> is
>>>> specifying a message type
of ADPSFC. This is one potential problem.
>>>> Another
>>>> problem
is that the model output appears to be valid at 20100615_12 and
>>>>
the
>>>> observations are valid 24 hours later. You could change
your config file
>>>> to
>>>> open the valid time window to a larger
value using beg_ds and end_ds.
>>>> Please
>>>> let us know if you
have any other questions.
>>>>
>>>> Thanks,
>>>>
>>>> Paul
>>>>
>>>>
>>>>
>>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16
WS/WXN} wrote:
>>>>> Thu Jun 17 15:45:21 2010: Request 38629 was
acted upon.
>>>>> Transaction: Ticket created by craigr at offutt.af.mil
>>>>> Queue: met_help
>>>>> Subject: Point Stat Question
>>>>> Owner: Nobody
>>>>> Requestors: craigr at offutt.af.mil
>>>>> Status: new
>>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>>
>>>>>
>>>>> I am running point stat on WRF winds but I am not getting
any stats. My
>>>>> code to run this is:
>>>>>
>>>>>
>>>>>
>>>>>
/sphome/qcteam/METv2.0/bin/point_stat
>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>> /h/data/global/WXQC/data/obs_pr
>>>>>
>>>>> ep/2010061512
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>
-outdir /h/data/global/WXQC/data/met/data -v 2
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Attached is the config file along with the ob file
and model file.
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> The
following is the output I get. Why the 0 pairs? Thanks Bob Craig
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> GSL_RNG_TYPE=mt19937
>>>>>
>>>>> GSL_RNG_SEED=3818134092
>>>>>
>>>>> Forecast File:
>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>>
>>>>> Climatology File: none
>>>>>
>>>>> Configuration File:
>>>>> /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>
>>>>> Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>>>>>
>>>>>
>>>>>
>>>>> ----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>> Reading records for UGRD/L10.
>>>>>
>>>>>
For UGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>
>>>>>
>>>>>
>>>>> ----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>> Reading records for VGRD/L10.
>>>>>
>>>>> For
VGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>
>>>>>
>>>>>
>>>>> ----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>> Searching 10473 observations from 170 PrepBufr
messages.
>>>>>
>>>>>
>>>>>
>>>>>
----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>>
Processing UGRD/L10 versus UGRD/L10, for observation type ADPSFC, over
>>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>>
>>>>>
>>>>>
>>>>> ----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>> Processing VGRD/L10 versus VGRD/L10, for
observation type ADPSFC, over
>>>>> region FULL, for interpolation
method MAX(4), using 0 pairs.
>>>>>
>>>>>
>>>>>
>>>>>
----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>>
Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
>>>>> at
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>> c.txt
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>> s.txt
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
>>>>> t.txt
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
>>>>> 1l2.txt
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
>>>>> 1l2.txt
>>>>>
>>>
>>
>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
From: Paul Oldenburg
Time: Thu Jun 24 13:21:40 2010
Bob,
Can you please redirect the make output to a log file and send us the
entire
make output? That will help us diagnose.
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=38629 >
>
> Gmake solved that problem - compiled with the following errors
(these were at the bottom of the listing):
>
> ld: 0711-224 WARNING: Duplicate symbol: .strlenFPCc
> ld: 0711-345 Use the -bloadmap or -bnoquiet option to obtain more
information.
> ld: 0711-317 ERROR: Undefined symbol: .openbf_
> ld: 0711-317 ERROR: Undefined symbol: .datelen_
> ld: 0711-317 ERROR: Undefined symbol: .ireadns_
> ld: 0711-317 ERROR: Undefined symbol: .ufbint_
> ld: 0711-317 ERROR: Undefined symbol: .ufbevn_
> ld: 0711-317 ERROR: Undefined symbol: .ufbtab_
> gmake[3]: *** [pb2nc] Error 8
> gmake[3]: Leaving directory `/h/sphome/qcteam/METv2.0/src/pb2nc'
> gmake[2]: *** [all] Error 2
> gmake[2]: Leaving directory `/h/sphome/qcteam/METv2.0/src/pb2nc'
> gmake[1]: *** [targets] Error 2
> gmake[1]: Leaving directory `/h/sphome/qcteam/METv2.0'
> gmake: *** [all] Error 2
>
> Should I worry about these - what is the fix?
>
> Thanks
>
> -----Original Message-----
> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent: Thursday, June 24, 2010 2:00 PM
> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>
> Bob,
>
> These errors look like errors that are output by the IBM version of
make. MET
> needs to be built with GNU make, which is called gmake on our IBM
machine.
> Please look on your IBM machine to see if you have gmake installed
and use it
> instead.
>
> 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=38629 >
>>
>> I will reset the quality threshold.
>>
>> I download the updates to MET and placed in the top level MET dir.
When I run
>> a Make clean I get the following errors:
>>
>> [496 qcteam a31n09 /sphome/qcteam/METv2.0] make clean
>> "Makefile", line 191: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 193: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 194: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 196: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 197: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 199: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 200: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 202: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 203: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 205: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 206: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 208: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 209: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 211: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 212: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 214: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 215: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 217: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 218: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 220: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> make: 1254-058 Fatal errors encountered -- cannot continue.
>>
>> The make file is attached. What is causing this error?
>>
>> Thanks
>>
>> Bob
>>
>>
>>
>> -----Original Message-----
>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent: Thursday, June 24, 2010 11:22 AM
>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>
>> Bob,
>>
>> I sent that email to John through the met_help ticketing system.
It was not
>> intended for you, sorry.
>>
>> We tested the PB2NCConfig file that you sent, and we found that it
did not
>> process the message types that are needed by point_stat. The
reason was that
>> the quality_mark_thresh setting was set to 1 instead of 2, and the
message
>> types
>> ADPSFC and ADPUPA were not processed as a result. Please try the
following
>> setting in your PB2NCConfig file:
>>
>> quality_mark_thresh = 2;
>>
>> Also, it is important to make sure that you are running the latest
MET
>> patches.
>> If you have not done so, please download and deploy the patches
from the MET
>> website at:
>>
>>
http://www.dtcenter.org/met/users/support/known_issues/METv2.0/index.php
>>
>> 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=38629 >
>>>
>>> John, you sent this to me - did you intend someone else?
>>>
>>> Bob
>>>
>>> -----Original Message-----
>>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>>> Sent: Wednesday, June 23, 2010 10:58 AM
>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>
>>> John,
>>>
>>> I'm trying to look into Bob's problems, but I'm struggling with
ncdump. Am
>>> I
>>> doing something wrong below where I attempt to view the variable
values?
>>>
>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -h
obs_2010062200.nc
>>> netcdf obs_2010062200 {
>>> dimensions:
>>> mxstr = 16 ;
>>> hdr_arr_len = 3 ;
>>> obs_arr_len = 5 ;
>>> nobs = UNLIMITED ; // (10004 currently)
>>> nhdr = 162 ;
>>> variables:
>>> float obs_arr(nobs, obs_arr_len) ;
>>> obs_arr:long_name = "array of observation values"
;
>>> obs_arr:_fill_value = -9999.f ;
>>> obs_arr:columns = "hdr_id gc lvl hgt ob" ;
>>> obs_arr:hdr_id_long_name = "index of matching
header data" ;
>>> obs_arr:gc_long_name = "grib code corresponding to
the
>>> observation type" ;
>>> obs_arr:lvl_long_name = "pressure level (hPa) or
>>> accumulation
>>> interval (h)" ;
>>> obs_arr:hgt_long_name = "height in meters above
sea level
>>> (msl)" ;
>>> obs_arr:ob_long_name = "observation value" ;
>>> char hdr_typ(nhdr, mxstr) ;
>>> hdr_typ:long_name = "message type" ;
>>> char hdr_sid(nhdr, mxstr) ;
>>> hdr_sid:long_name = "station identification" ;
>>> char hdr_vld(nhdr, mxstr) ;
>>> hdr_vld:long_name = "valid time" ;
>>> hdr_vld:units = "YYYYMMDD_HHMMSS" ;
>>> float hdr_arr(nhdr, hdr_arr_len) ;
>>> hdr_arr:long_name = "array of observation station
header
>>> values" ;
>>> hdr_arr:_fill_value = -9999.f ;
>>> hdr_arr:columns = "lat lon elv" ;
>>> hdr_arr:lat_long_name = "latitude" ;
>>> hdr_arr:lat_units = "degrees_north" ;
>>> hdr_arr:lon_long_name = "longitude" ;
>>> hdr_arr:lon_units = "degrees_east" ;
>>> hdr_arr:elv_long_name = "elevation" ;
>>> hdr_arr:elv_units = "meters above sea level (msl)"
;
>>>
>>> // global attributes:
>>> :FileOrigins = "File
>>> /h/data/global/WXQC/data/obs_prep/2010062200 generated
20100622_184241 UTC
>>> on
>>> host a31n09.afwa.af.mil by the PB2NC tool" ;
>>> }
>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_arr:long_name
>>> obs_2010062200.nc
>>> ncdump: NetCDF: Variable not found
>>> netcdf obs_2010062200 {
>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>>> obs_2010062200.nc
>>> ncdump: NetCDF: Variable not found
>>> netcdf obs_2010062200 {
>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>>> obs_2010062200.nc
>>> ncdump: NetCDF: Variable not found
>>> netcdf obs_2010062200 {
>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_typ:long_name
>>> obs_2010062200.nc
>>> ncdump: NetCDF: Variable not found
>>> netcdf obs_2010062200 {
>>>
>>>
>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>
>>>> I am getting the same 0 pair error with the attached ob file. I
built this
>>>> file using the attached config file. Is there surface
observations in my
>>>> file? Is the valid time of the file correct? I am using the
-dump option
>>>> but
>>>> I am not getting any file showing up in the directory I set for
this.
>>>>
>>>> Bob
>>>>
>>>> -----Original Message-----
>>>> From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
>>>> Sent: Tuesday, June 22, 2010 3:30 PM
>>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>
>>>> Bob,
>>>>
>>>> We looked more closely at the observation data that you included
in your
>>>> original email, and we noticed that they don't appear to be
surface
>>>> observations. When verifying 2m temp or 10m winds in point_stat,
only
>>>> observations of type ADPSFC or SFCSHP can be used (or ONLYSF for
either).
>>>>
>>>> Please change your point_stat config fcst_field setting to the
following
>>>> value:
>>>>
>>>> fcst_field[] = [ "33/Z10", "34/Z10"];
>>>>
>>>> also, please change the message type:
>>>>
>>>> message_type[] = [ "ADPSFC" ];
>>>>
>>>> Please let us know if you have any more questions.
>>>>
>>>> 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=38629 >
>>>>>
>>>>> Thanks for looking into this. I inadvertently left off my
config file
>>>>> name
>>>>> from the command string, hence it defaulted to profiler data. My
config
>>>>> file
>>>>> has "ANYSFC" which the documentation says includes ADPUPA (RAOB
data).
>>>>>
>>>>> For the observations, I am pulling the prepbufr files from
>>>>> ftp.ncep.noaa.gov
>>>>> in the /pub/data/nccf/com/gfs/prod/gdas.20100615 for the 15 Jun
files for
>>>>> the case below. An example of the filenames I am pulling are
>>>>> gdas2.t12z.prepbufr.nr. I rename these files to 20100615_12 in
this
>>>>> case.
>>>>> Since this file came out of the 20100615 directory I would
expect it to be
>>>>> valid for 20100615 at 12Z. Is there a way I can dump the header
of the
>>>>> prepbufr file to see the date that's in the file?
>>>>>
>>>>> Bob
>>>>>
>>>>>
>>>>> -----Original Message-----
>>>>> From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
>>>>> Sent: Friday, June 18, 2010 11:56 AM
>>>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>>
>>>>> Robert,
>>>>>
>>>>> Your observation file has a message type of PROFLR, but your
config file
>>>>> is
>>>>> specifying a message type of ADPSFC. This is one potential
problem.
>>>>> Another
>>>>> problem is that the model output appears to be valid at
20100615_12 and
>>>>> the
>>>>> observations are valid 24 hours later. You could change your
config file
>>>>> to
>>>>> open the valid time window to a larger value using beg_ds and
end_ds.
>>>>> Please
>>>>> let us know if you have any other questions.
>>>>>
>>>>> Thanks,
>>>>>
>>>>> Paul
>>>>>
>>>>>
>>>>>
>>>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN}
wrote:
>>>>>> Thu Jun 17 15:45:21 2010: Request 38629 was acted upon.
>>>>>> Transaction: Ticket created by craigr at offutt.af.mil
>>>>>> Queue: met_help
>>>>>> Subject: Point Stat Question
>>>>>> Owner: Nobody
>>>>>> Requestors: craigr at offutt.af.mil
>>>>>> Status: new
>>>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>>>
>>>>>>
>>>>>> I am running point stat on WRF winds but I am not getting any
stats. My
>>>>>> code to run this is:
>>>>>>
>>>>>>
>>>>>>
>>>>>> /sphome/qcteam/METv2.0/bin/point_stat
>>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>>> /h/data/global/WXQC/data/obs_pr
>>>>>>
>>>>>> ep/2010061512
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>> -outdir /h/data/global/WXQC/data/met/data -v 2
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> Attached is the config file along with the ob file and model
file.
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> The following is the output I get. Why the 0 pairs? Thanks
Bob Craig
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> GSL_RNG_TYPE=mt19937
>>>>>>
>>>>>> GSL_RNG_SEED=3818134092
>>>>>>
>>>>>> Forecast File:
>>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>>>
>>>>>> Climatology File: none
>>>>>>
>>>>>> Configuration File:
>>>>>> /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>>
>>>>>> Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>>>>>>
>>>>>>
>>>>>>
>>>>>> ----------------------------------------
>>>>>>
>>>>>>
>>>>>>
>>>>>> Reading records for UGRD/L10.
>>>>>>
>>>>>> For UGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>>
>>>>>>
>>>>>>
>>>>>> ----------------------------------------
>>>>>>
>>>>>>
>>>>>>
>>>>>> Reading records for VGRD/L10.
>>>>>>
>>>>>> For VGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>>
>>>>>>
>>>>>>
>>>>>> ----------------------------------------
>>>>>>
>>>>>>
>>>>>>
>>>>>> Searching 10473 observations from 170 PrepBufr messages.
>>>>>>
>>>>>>
>>>>>>
>>>>>> ----------------------------------------
>>>>>>
>>>>>>
>>>>>>
>>>>>> Processing UGRD/L10 versus UGRD/L10, for observation type
ADPSFC, over
>>>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>>>
>>>>>>
>>>>>>
>>>>>> ----------------------------------------
>>>>>>
>>>>>>
>>>>>>
>>>>>> Processing VGRD/L10 versus VGRD/L10, for observation type
ADPSFC, over
>>>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>>>
>>>>>>
>>>>>>
>>>>>> ----------------------------------------
>>>>>>
>>>>>>
>>>>>>
>>>>>> Output file:
>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
>>>>>> at
>>>>>>
>>>>>> Output file:
>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>>> c.txt
>>>>>>
>>>>>> Output file:
>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>>> s.txt
>>>>>>
>>>>>> Output file:
>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
>>>>>> t.txt
>>>>>>
>>>>>> Output file:
>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
>>>>>> 1l2.txt
>>>>>>
>>>>>> Output file:
>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
>>>>>> 1l2.txt
>>>>>>
>>
>>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
From: Paul Oldenburg
Time: Thu Jun 24 13:23:39 2010
Bob,
John tells me that recompiling BUFRLIB using the following commands
should fix
this problem:
xlc -c -DUNDERSCORE *.c
xlf -c -qextname *.f *.F
ar crv libbufr.a *.o
These are taken from this page:
http://www.dtcenter.org/met/users/support/online_tutorial/METv2.0/compilation/req_libs.php
And then try rebuilding MET.
Thanks,
Paul
Paul Oldenburg wrote:
> Bob,
>
> Can you please redirect the make output to a log file and send us
the entire
> make output? That will help us diagnose.
>
> 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=38629 >
>>
>> Gmake solved that problem - compiled with the following errors
(these were at the bottom of the listing):
>>
>> ld: 0711-224 WARNING: Duplicate symbol: .strlenFPCc
>> ld: 0711-345 Use the -bloadmap or -bnoquiet option to obtain more
information.
>> ld: 0711-317 ERROR: Undefined symbol: .openbf_
>> ld: 0711-317 ERROR: Undefined symbol: .datelen_
>> ld: 0711-317 ERROR: Undefined symbol: .ireadns_
>> ld: 0711-317 ERROR: Undefined symbol: .ufbint_
>> ld: 0711-317 ERROR: Undefined symbol: .ufbevn_
>> ld: 0711-317 ERROR: Undefined symbol: .ufbtab_
>> gmake[3]: *** [pb2nc] Error 8
>> gmake[3]: Leaving directory `/h/sphome/qcteam/METv2.0/src/pb2nc'
>> gmake[2]: *** [all] Error 2
>> gmake[2]: Leaving directory `/h/sphome/qcteam/METv2.0/src/pb2nc'
>> gmake[1]: *** [targets] Error 2
>> gmake[1]: Leaving directory `/h/sphome/qcteam/METv2.0'
>> gmake: *** [all] Error 2
>>
>> Should I worry about these - what is the fix?
>>
>> Thanks
>>
>> -----Original Message-----
>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent: Thursday, June 24, 2010 2:00 PM
>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>
>> Bob,
>>
>> These errors look like errors that are output by the IBM version of
make. MET
>> needs to be built with GNU make, which is called gmake on our IBM
machine.
>> Please look on your IBM machine to see if you have gmake installed
and use it
>> instead.
>>
>> 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=38629 >
>>>
>>> I will reset the quality threshold.
>>>
>>> I download the updates to MET and placed in the top level MET dir.
When I run
>>> a Make clean I get the following errors:
>>>
>>> [496 qcteam a31n09 /sphome/qcteam/METv2.0] make clean
>>> "Makefile", line 191: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 193: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 194: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 196: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 197: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 199: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 200: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 202: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 203: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 205: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 206: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 208: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 209: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 211: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 212: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 214: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 215: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 217: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 218: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 220: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> make: 1254-058 Fatal errors encountered -- cannot continue.
>>>
>>> The make file is attached. What is causing this error?
>>>
>>> Thanks
>>>
>>> Bob
>>>
>>>
>>>
>>> -----Original Message-----
>>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>>> Sent: Thursday, June 24, 2010 11:22 AM
>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>
>>> Bob,
>>>
>>> I sent that email to John through the met_help ticketing system.
It was not
>>> intended for you, sorry.
>>>
>>> We tested the PB2NCConfig file that you sent, and we found that it
did not
>>> process the message types that are needed by point_stat. The
reason was that
>>> the quality_mark_thresh setting was set to 1 instead of 2, and the
message
>>> types
>>> ADPSFC and ADPUPA were not processed as a result. Please try the
following
>>> setting in your PB2NCConfig file:
>>>
>>> quality_mark_thresh = 2;
>>>
>>> Also, it is important to make sure that you are running the latest
MET
>>> patches.
>>> If you have not done so, please download and deploy the patches
from the MET
>>> website at:
>>>
>>>
http://www.dtcenter.org/met/users/support/known_issues/METv2.0/index.php
>>>
>>> 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=38629 >
>>>>
>>>> John, you sent this to me - did you intend someone else?
>>>>
>>>> Bob
>>>>
>>>> -----Original Message-----
>>>> From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
>>>> Sent: Wednesday, June 23, 2010 10:58 AM
>>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>
>>>> John,
>>>>
>>>> I'm trying to look into Bob's problems, but I'm struggling with
ncdump. Am
>>>> I
>>>> doing something wrong below where I attempt to view the variable
values?
>>>>
>>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -h
obs_2010062200.nc
>>>> netcdf obs_2010062200 {
>>>> dimensions:
>>>> mxstr = 16 ;
>>>> hdr_arr_len = 3 ;
>>>> obs_arr_len = 5 ;
>>>> nobs = UNLIMITED ; // (10004 currently)
>>>> nhdr = 162 ;
>>>> variables:
>>>> float obs_arr(nobs, obs_arr_len) ;
>>>> obs_arr:long_name = "array of observation values"
;
>>>> obs_arr:_fill_value = -9999.f ;
>>>> obs_arr:columns = "hdr_id gc lvl hgt ob" ;
>>>> obs_arr:hdr_id_long_name = "index of matching
header data" ;
>>>> obs_arr:gc_long_name = "grib code corresponding
to the
>>>> observation type" ;
>>>> obs_arr:lvl_long_name = "pressure level (hPa) or
>>>> accumulation
>>>> interval (h)" ;
>>>> obs_arr:hgt_long_name = "height in meters above
sea level
>>>> (msl)" ;
>>>> obs_arr:ob_long_name = "observation value" ;
>>>> char hdr_typ(nhdr, mxstr) ;
>>>> hdr_typ:long_name = "message type" ;
>>>> char hdr_sid(nhdr, mxstr) ;
>>>> hdr_sid:long_name = "station identification" ;
>>>> char hdr_vld(nhdr, mxstr) ;
>>>> hdr_vld:long_name = "valid time" ;
>>>> hdr_vld:units = "YYYYMMDD_HHMMSS" ;
>>>> float hdr_arr(nhdr, hdr_arr_len) ;
>>>> hdr_arr:long_name = "array of observation station
header
>>>> values" ;
>>>> hdr_arr:_fill_value = -9999.f ;
>>>> hdr_arr:columns = "lat lon elv" ;
>>>> hdr_arr:lat_long_name = "latitude" ;
>>>> hdr_arr:lat_units = "degrees_north" ;
>>>> hdr_arr:lon_long_name = "longitude" ;
>>>> hdr_arr:lon_units = "degrees_east" ;
>>>> hdr_arr:elv_long_name = "elevation" ;
>>>> hdr_arr:elv_units = "meters above sea level
(msl)" ;
>>>>
>>>> // global attributes:
>>>> :FileOrigins = "File
>>>> /h/data/global/WXQC/data/obs_prep/2010062200 generated
20100622_184241 UTC
>>>> on
>>>> host a31n09.afwa.af.mil by the PB2NC tool" ;
>>>> }
>>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_arr:long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF: Variable not found
>>>> netcdf obs_2010062200 {
>>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF: Variable not found
>>>> netcdf obs_2010062200 {
>>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF: Variable not found
>>>> netcdf obs_2010062200 {
>>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_typ:long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF: Variable not found
>>>> netcdf obs_2010062200 {
>>>>
>>>>
>>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>>
>>>>> I am getting the same 0 pair error with the attached ob file. I
built this
>>>>> file using the attached config file. Is there surface
observations in my
>>>>> file? Is the valid time of the file correct? I am using the
-dump option
>>>>> but
>>>>> I am not getting any file showing up in the directory I set for
this.
>>>>>
>>>>> Bob
>>>>>
>>>>> -----Original Message-----
>>>>> From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
>>>>> Sent: Tuesday, June 22, 2010 3:30 PM
>>>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>>
>>>>> Bob,
>>>>>
>>>>> We looked more closely at the observation data that you included
in your
>>>>> original email, and we noticed that they don't appear to be
surface
>>>>> observations. When verifying 2m temp or 10m winds in
point_stat, only
>>>>> observations of type ADPSFC or SFCSHP can be used (or ONLYSF for
either).
>>>>>
>>>>> Please change your point_stat config fcst_field setting to the
following
>>>>> value:
>>>>>
>>>>> fcst_field[] = [ "33/Z10", "34/Z10"];
>>>>>
>>>>> also, please change the message type:
>>>>>
>>>>> message_type[] = [ "ADPSFC" ];
>>>>>
>>>>> Please let us know if you have any more questions.
>>>>>
>>>>> 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=38629 >
>>>>>>
>>>>>> Thanks for looking into this. I inadvertently left off my
config file
>>>>>> name
>>>>>> from the command string, hence it defaulted to profiler data.
My config
>>>>>> file
>>>>>> has "ANYSFC" which the documentation says includes ADPUPA (RAOB
data).
>>>>>>
>>>>>> For the observations, I am pulling the prepbufr files from
>>>>>> ftp.ncep.noaa.gov
>>>>>> in the /pub/data/nccf/com/gfs/prod/gdas.20100615 for the 15 Jun
files for
>>>>>> the case below. An example of the filenames I am pulling are
>>>>>> gdas2.t12z.prepbufr.nr. I rename these files to 20100615_12
in this
>>>>>> case.
>>>>>> Since this file came out of the 20100615 directory I would
expect it to be
>>>>>> valid for 20100615 at 12Z. Is there a way I can dump the
header of the
>>>>>> prepbufr file to see the date that's in the file?
>>>>>>
>>>>>> Bob
>>>>>>
>>>>>>
>>>>>> -----Original Message-----
>>>>>> From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
>>>>>> Sent: Friday, June 18, 2010 11:56 AM
>>>>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>>>
>>>>>> Robert,
>>>>>>
>>>>>> Your observation file has a message type of PROFLR, but your
config file
>>>>>> is
>>>>>> specifying a message type of ADPSFC. This is one potential
problem.
>>>>>> Another
>>>>>> problem is that the model output appears to be valid at
20100615_12 and
>>>>>> the
>>>>>> observations are valid 24 hours later. You could change your
config file
>>>>>> to
>>>>>> open the valid time window to a larger value using beg_ds and
end_ds.
>>>>>> Please
>>>>>> let us know if you have any other questions.
>>>>>>
>>>>>> Thanks,
>>>>>>
>>>>>> Paul
>>>>>>
>>>>>>
>>>>>>
>>>>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN}
wrote:
>>>>>>> Thu Jun 17 15:45:21 2010: Request 38629 was acted upon.
>>>>>>> Transaction: Ticket created by craigr at offutt.af.mil
>>>>>>> Queue: met_help
>>>>>>> Subject: Point Stat Question
>>>>>>> Owner: Nobody
>>>>>>> Requestors: craigr at offutt.af.mil
>>>>>>> Status: new
>>>>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>>>>
>>>>>>>
>>>>>>> I am running point stat on WRF winds but I am not getting any
stats. My
>>>>>>> code to run this is:
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> /sphome/qcteam/METv2.0/bin/point_stat
>>>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>>>> /h/data/global/WXQC/data/obs_pr
>>>>>>>
>>>>>>> ep/2010061512
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>>> -outdir /h/data/global/WXQC/data/met/data -v 2
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Attached is the config file along with the ob file and model
file.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> The following is the output I get. Why the 0 pairs? Thanks
Bob Craig
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> GSL_RNG_TYPE=mt19937
>>>>>>>
>>>>>>> GSL_RNG_SEED=3818134092
>>>>>>>
>>>>>>> Forecast File:
>>>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>>>>
>>>>>>> Climatology File: none
>>>>>>>
>>>>>>> Configuration File:
>>>>>>> /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>>>
>>>>>>> Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> ----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Reading records for UGRD/L10.
>>>>>>>
>>>>>>> For UGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> ----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Reading records for VGRD/L10.
>>>>>>>
>>>>>>> For VGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> ----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Searching 10473 observations from 170 PrepBufr messages.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> ----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Processing UGRD/L10 versus UGRD/L10, for observation type
ADPSFC, over
>>>>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> ----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Processing VGRD/L10 versus VGRD/L10, for observation type
ADPSFC, over
>>>>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> ----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
>>>>>>> at
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>>>> c.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>>>> s.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
>>>>>>> t.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
>>>>>>> 1l2.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
>>>>>>> 1l2.txt
>>>>>>>
>>>
>>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #38629] Point Stat Question
From: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Time: Thu Jun 24 13:31:06 2010
It looks like pointstat didn't compile, at least I didn't get a bin
file for it so I need to get past the previous error (previous email)
-----Original Message-----
From: RAL HelpDesk {for Craig, Robert J
Civ USAF AFWA 16 WS/WXN} [mailto:met_help at ucar.edu]
Sent: Thursday,
June 24, 2010 2:20 PM
To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Subject: RE: [rt.rap.ucar.edu #38629] Point Stat Question
Gmake
solved that problem - compiled with the following errors (these were
at the bottom of the listing):
ld: 0711-224 WARNING: Duplicate
symbol: .strlenFPCc
ld: 0711-345 Use the -bloadmap or -bnoquiet
option to obtain more information.
ld: 0711-317 ERROR: Undefined
symbol: .openbf_
ld: 0711-317 ERROR: Undefined symbol: .datelen_
ld:
0711-317 ERROR: Undefined symbol: .ireadns_
ld: 0711-317 ERROR:
Undefined symbol: .ufbint_
ld: 0711-317 ERROR: Undefined symbol:
.ufbevn_
ld: 0711-317 ERROR: Undefined symbol: .ufbtab_
gmake[3]:
*** [pb2nc] Error 8
gmake[3]: Leaving directory
`/h/sphome/qcteam/METv2.0/src/pb2nc'
gmake[2]: *** [all] Error 2
gmake[2]: Leaving directory `/h/sphome/qcteam/METv2.0/src/pb2nc'
gmake[1]: *** [targets] Error 2
gmake[1]: Leaving directory
`/h/sphome/qcteam/METv2.0'
gmake: *** [all] Error 2
Should I worry
about these - what is the fix?
Thanks
-----Original Message-----
From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
Sent: Thursday, June 24, 2010 2:00 PM
To: Craig, Robert J Civ USAF
AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat
Question
Bob,
These errors look like errors that are output by
the IBM version of make. MET
needs to be built with GNU make, which
is called gmake on our IBM machine.
Please look on your IBM machine
to see if you have gmake installed and use it
instead.
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=38629 >
>
> I will
reset the quality threshold.
>
> I download the updates to MET and
placed in the top level MET dir. When I run
> a Make clean I get the
following errors:
>
> [496 qcteam a31n09 /sphome/qcteam/METv2.0]
make clean
> "Makefile", line 191: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 193:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 194: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 196:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 197: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 199:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 200: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 202:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 203: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 205:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 206: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 208:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 209: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 211:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 212: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 214:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 215: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 217:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> "Makefile", line 218: make: 1254-055 Dependency line
needs colon or double
> colon operator.
> "Makefile", line 220:
make: 1254-055 Dependency line needs colon or double
> colon
operator.
> make: 1254-058 Fatal errors encountered -- cannot
continue.
>
> The make file is attached. What is causing this
error?
>
> Thanks
>
> Bob
>
>
>
> -----Original Message-----
> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent: Thursday, June 24, 2010 11:22 AM
> To: Craig, Robert J Civ
USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38629] Point
Stat Question
>
> Bob,
>
> I sent that email to John through the
met_help ticketing system. It was not
> intended for you, sorry.
>
> We tested the PB2NCConfig file that you sent, and we found that it
did not
> process the message types that are needed by point_stat.
The reason was that
> the quality_mark_thresh setting was set to 1
instead of 2, and the message
> types
> ADPSFC and ADPUPA were not
processed as a result. Please try the following
> setting in your
PB2NCConfig file:
>
> quality_mark_thresh = 2;
>
> Also, it is
important to make sure that you are running the latest MET
> patches.
> If you have not done so, please download and deploy the patches
from the MET
> website at:
>
>
http://www.dtcenter.org/met/users/support/known_issues/METv2.0/index.php
>
> 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=38629 >
>>
>>
John, you sent this to me - did you intend someone else?
>>
>> Bob
>>
>> -----Original Message-----
>> From: RAL HelpDesk {for Paul
Oldenburg} [mailto:met_help at ucar.edu]
>> Sent: Wednesday, June 23,
2010 10:58 AM
>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>
>>
John,
>>
>> I'm trying to look into Bob's problems, but I'm
struggling with ncdump. Am
>> I
>> doing something wrong below
where I attempt to view the variable values?
>>
>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -h
obs_2010062200.nc
>> netcdf obs_2010062200 {
>> dimensions:
>>
mxstr = 16 ;
>> hdr_arr_len = 3 ;
>> obs_arr_len = 5
;
>> nobs = UNLIMITED ; // (10004 currently)
>> nhdr
= 162 ;
>> variables:
>> float obs_arr(nobs, obs_arr_len) ;
>> obs_arr:long_name = "array of observation values" ;
>> obs_arr:_fill_value = -9999.f ;
>>
obs_arr:columns = "hdr_id gc lvl hgt ob" ;
>>
obs_arr:hdr_id_long_name = "index of matching header data" ;
>>
obs_arr:gc_long_name = "grib code corresponding to the
>> observation
type" ;
>> obs_arr:lvl_long_name = "pressure level
(hPa) or
>> accumulation
>> interval (h)" ;
>>
obs_arr:hgt_long_name = "height in meters above sea level
>> (msl)" ;
>> obs_arr:ob_long_name = "observation value" ;
>>
char hdr_typ(nhdr, mxstr) ;
>> hdr_typ:long_name =
"message type" ;
>> char hdr_sid(nhdr, mxstr) ;
>>
hdr_sid:long_name = "station identification" ;
>> char
hdr_vld(nhdr, mxstr) ;
>> hdr_vld:long_name = "valid
time" ;
>> hdr_vld:units = "YYYYMMDD_HHMMSS" ;
>>
float hdr_arr(nhdr, hdr_arr_len) ;
>>
hdr_arr:long_name = "array of observation station header
>> values" ;
>> hdr_arr:_fill_value = -9999.f ;
>>
hdr_arr:columns = "lat lon elv" ;
>>
hdr_arr:lat_long_name = "latitude" ;
>>
hdr_arr:lat_units = "degrees_north" ;
>>
hdr_arr:lon_long_name = "longitude" ;
>>
hdr_arr:lon_units = "degrees_east" ;
>>
hdr_arr:elv_long_name = "elevation" ;
>>
hdr_arr:elv_units = "meters above sea level (msl)" ;
>>
>> // global
attributes:
>> :FileOrigins = "File
>>
/h/data/global/WXQC/data/obs_prep/2010062200 generated 20100622_184241
UTC
>> on
>> host a31n09.afwa.af.mil by the PB2NC tool" ;
>> }
>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_arr:long_name
>> obs_2010062200.nc
>> ncdump: NetCDF: Variable
not found
>> netcdf obs_2010062200 {
>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>> obs_2010062200.nc
>> ncdump: NetCDF:
Variable not found
>> netcdf obs_2010062200 {
>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>> obs_2010062200.nc
>> ncdump: NetCDF:
Variable not found
>> netcdf obs_2010062200 {
>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_typ:long_name
>> obs_2010062200.nc
>> ncdump: NetCDF: Variable
not found
>> netcdf obs_2010062200 {
>>
>>
>> RAL HelpDesk {for
Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>
>>> I
am getting the same 0 pair error with the attached ob file. I built
this
>>> file using the attached config file. Is there surface
observations in my
>>> file? Is the valid time of the file correct?
I am using the -dump option
>>> but
>>> I am not getting any file
showing up in the directory I set for this.
>>>
>>> Bob
>>>
>>>
-----Original Message-----
>>> From: RAL HelpDesk {for Paul
Oldenburg} [mailto:met_help at ucar.edu]
>>> Sent: Tuesday, June 22,
2010 3:30 PM
>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>
>>>
Bob,
>>>
>>> We looked more closely at the observation data that you
included in your
>>> original email, and we noticed that they don't
appear to be surface
>>> observations. When verifying 2m temp or 10m
winds in point_stat, only
>>> observations of type ADPSFC or SFCSHP
can be used (or ONLYSF for either).
>>>
>>> Please change your
point_stat config fcst_field setting to the following
>>> value:
>>>
>>> fcst_field[] = [ "33/Z10", "34/Z10"];
>>>
>>> also, please
change the message type:
>>>
>>> message_type[] = [ "ADPSFC" ];
>>>
>>> Please let us know if you have any more questions.
>>>
>>>
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=38629 >
>>>>
>>>>
Thanks for looking into this. I inadvertently left off my config file
>>>> name
>>>> from the command string, hence it defaulted to
profiler data. My config
>>>> file
>>>> has "ANYSFC" which the
documentation says includes ADPUPA (RAOB data).
>>>>
>>>> For the
observations, I am pulling the prepbufr files from
>>>>
ftp.ncep.noaa.gov
>>>> in the
/pub/data/nccf/com/gfs/prod/gdas.20100615 for the 15 Jun files for
>>>> the case below. An example of the filenames I am pulling are
>>>> gdas2.t12z.prepbufr.nr. I rename these files to 20100615_12 in
this
>>>> case.
>>>> Since this file came out of the 20100615
directory I would expect it to be
>>>> valid for 20100615 at 12Z. Is
there a way I can dump the header of the
>>>> prepbufr file to see
the date that's in the file?
>>>>
>>>> Bob
>>>>
>>>>
>>>>
-----Original Message-----
>>>> From: RAL HelpDesk {for Paul
Oldenburg} [mailto:met_help at ucar.edu]
>>>> Sent: Friday, June 18,
2010 11:56 AM
>>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>>
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>
>>>>
Robert,
>>>>
>>>> Your observation file has a message type of
PROFLR, but your config file
>>>> is
>>>> specifying a message type
of ADPSFC. This is one potential problem.
>>>> Another
>>>> problem
is that the model output appears to be valid at 20100615_12 and
>>>>
the
>>>> observations are valid 24 hours later. You could change
your config file
>>>> to
>>>> open the valid time window to a larger
value using beg_ds and end_ds.
>>>> Please
>>>> let us know if you
have any other questions.
>>>>
>>>> Thanks,
>>>>
>>>> Paul
>>>>
>>>>
>>>>
>>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16
WS/WXN} wrote:
>>>>> Thu Jun 17 15:45:21 2010: Request 38629 was
acted upon.
>>>>> Transaction: Ticket created by craigr at offutt.af.mil
>>>>> Queue: met_help
>>>>> Subject: Point Stat Question
>>>>> Owner: Nobody
>>>>> Requestors: craigr at offutt.af.mil
>>>>> Status: new
>>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>>
>>>>>
>>>>> I am running point stat on WRF winds but I am not getting
any stats. My
>>>>> code to run this is:
>>>>>
>>>>>
>>>>>
>>>>>
/sphome/qcteam/METv2.0/bin/point_stat
>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>> /h/data/global/WXQC/data/obs_pr
>>>>>
>>>>> ep/2010061512
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>
-outdir /h/data/global/WXQC/data/met/data -v 2
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Attached is the config file along with the ob file
and model file.
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> The
following is the output I get. Why the 0 pairs? Thanks Bob Craig
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> GSL_RNG_TYPE=mt19937
>>>>>
>>>>> GSL_RNG_SEED=3818134092
>>>>>
>>>>> Forecast File:
>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>>
>>>>> Climatology File: none
>>>>>
>>>>> Configuration File:
>>>>> /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>
>>>>> Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>>>>>
>>>>>
>>>>>
>>>>> ----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>> Reading records for UGRD/L10.
>>>>>
>>>>>
For UGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>
>>>>>
>>>>>
>>>>> ----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>> Reading records for VGRD/L10.
>>>>>
>>>>> For
VGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>
>>>>>
>>>>>
>>>>> ----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>> Searching 10473 observations from 170 PrepBufr
messages.
>>>>>
>>>>>
>>>>>
>>>>>
----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>>
Processing UGRD/L10 versus UGRD/L10, for observation type ADPSFC, over
>>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>>
>>>>>
>>>>>
>>>>> ----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>> Processing VGRD/L10 versus VGRD/L10, for
observation type ADPSFC, over
>>>>> region FULL, for interpolation
method MAX(4), using 0 pairs.
>>>>>
>>>>>
>>>>>
>>>>>
----------------------------------------
>>>>>
>>>>>
>>>>>
>>>>>
Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
>>>>> at
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>> c.txt
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>> s.txt
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
>>>>> t.txt
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
>>>>> 1l2.txt
>>>>>
>>>>> Output file:
>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
>>>>> 1l2.txt
>>>>>
>>>
>>
>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
From: Paul Oldenburg
Time: Thu Jun 24 13:34:17 2010
Bob,
Please try to rebuild BUFRLIB using the following steps and then try a
clean
build of MET using gmake.
xlc -c -DUNDERSCORE *.c
xlf -c -qextname *.f *.F
ar crv libbufr.a *.o
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=38629 >
>
> It looks like pointstat didn't compile, at least I didn't get a bin
file for it so I need to get past the previous error (previous email)
>
>
>
> -----Original Message-----
> From: RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN}
[mailto:met_help at ucar.edu]
> Sent: Thursday, June 24, 2010 2:20 PM
> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
> Subject: RE: [rt.rap.ucar.edu #38629] Point Stat Question
>
> Gmake solved that problem - compiled with the following errors
(these were at the bottom of the listing):
>
> ld: 0711-224 WARNING: Duplicate symbol: .strlenFPCc
> ld: 0711-345 Use the -bloadmap or -bnoquiet option to obtain more
information.
> ld: 0711-317 ERROR: Undefined symbol: .openbf_
> ld: 0711-317 ERROR: Undefined symbol: .datelen_
> ld: 0711-317 ERROR: Undefined symbol: .ireadns_
> ld: 0711-317 ERROR: Undefined symbol: .ufbint_
> ld: 0711-317 ERROR: Undefined symbol: .ufbevn_
> ld: 0711-317 ERROR: Undefined symbol: .ufbtab_
> gmake[3]: *** [pb2nc] Error 8
> gmake[3]: Leaving directory `/h/sphome/qcteam/METv2.0/src/pb2nc'
> gmake[2]: *** [all] Error 2
> gmake[2]: Leaving directory `/h/sphome/qcteam/METv2.0/src/pb2nc'
> gmake[1]: *** [targets] Error 2
> gmake[1]: Leaving directory `/h/sphome/qcteam/METv2.0'
> gmake: *** [all] Error 2
>
> Should I worry about these - what is the fix?
>
> Thanks
>
> -----Original Message-----
> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
> Sent: Thursday, June 24, 2010 2:00 PM
> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>
> Bob,
>
> These errors look like errors that are output by the IBM version of
make. MET
> needs to be built with GNU make, which is called gmake on our IBM
machine.
> Please look on your IBM machine to see if you have gmake installed
and use it
> instead.
>
> 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=38629 >
>>
>> I will reset the quality threshold.
>>
>> I download the updates to MET and placed in the top level MET dir.
When I run
>> a Make clean I get the following errors:
>>
>> [496 qcteam a31n09 /sphome/qcteam/METv2.0] make clean
>> "Makefile", line 191: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 193: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 194: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 196: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 197: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 199: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 200: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 202: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 203: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 205: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 206: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 208: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 209: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 211: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 212: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 214: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 215: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 217: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 218: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> "Makefile", line 220: make: 1254-055 Dependency line needs colon or
double
>> colon operator.
>> make: 1254-058 Fatal errors encountered -- cannot continue.
>>
>> The make file is attached. What is causing this error?
>>
>> Thanks
>>
>> Bob
>>
>>
>>
>> -----Original Message-----
>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent: Thursday, June 24, 2010 11:22 AM
>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>
>> Bob,
>>
>> I sent that email to John through the met_help ticketing system.
It was not
>> intended for you, sorry.
>>
>> We tested the PB2NCConfig file that you sent, and we found that it
did not
>> process the message types that are needed by point_stat. The
reason was that
>> the quality_mark_thresh setting was set to 1 instead of 2, and the
message
>> types
>> ADPSFC and ADPUPA were not processed as a result. Please try the
following
>> setting in your PB2NCConfig file:
>>
>> quality_mark_thresh = 2;
>>
>> Also, it is important to make sure that you are running the latest
MET
>> patches.
>> If you have not done so, please download and deploy the patches
from the MET
>> website at:
>>
>>
http://www.dtcenter.org/met/users/support/known_issues/METv2.0/index.php
>>
>> 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=38629 >
>>>
>>> John, you sent this to me - did you intend someone else?
>>>
>>> Bob
>>>
>>> -----Original Message-----
>>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>>> Sent: Wednesday, June 23, 2010 10:58 AM
>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>
>>> John,
>>>
>>> I'm trying to look into Bob's problems, but I'm struggling with
ncdump. Am
>>> I
>>> doing something wrong below where I attempt to view the variable
values?
>>>
>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -h
obs_2010062200.nc
>>> netcdf obs_2010062200 {
>>> dimensions:
>>> mxstr = 16 ;
>>> hdr_arr_len = 3 ;
>>> obs_arr_len = 5 ;
>>> nobs = UNLIMITED ; // (10004 currently)
>>> nhdr = 162 ;
>>> variables:
>>> float obs_arr(nobs, obs_arr_len) ;
>>> obs_arr:long_name = "array of observation values"
;
>>> obs_arr:_fill_value = -9999.f ;
>>> obs_arr:columns = "hdr_id gc lvl hgt ob" ;
>>> obs_arr:hdr_id_long_name = "index of matching
header data" ;
>>> obs_arr:gc_long_name = "grib code corresponding to
the
>>> observation type" ;
>>> obs_arr:lvl_long_name = "pressure level (hPa) or
>>> accumulation
>>> interval (h)" ;
>>> obs_arr:hgt_long_name = "height in meters above
sea level
>>> (msl)" ;
>>> obs_arr:ob_long_name = "observation value" ;
>>> char hdr_typ(nhdr, mxstr) ;
>>> hdr_typ:long_name = "message type" ;
>>> char hdr_sid(nhdr, mxstr) ;
>>> hdr_sid:long_name = "station identification" ;
>>> char hdr_vld(nhdr, mxstr) ;
>>> hdr_vld:long_name = "valid time" ;
>>> hdr_vld:units = "YYYYMMDD_HHMMSS" ;
>>> float hdr_arr(nhdr, hdr_arr_len) ;
>>> hdr_arr:long_name = "array of observation station
header
>>> values" ;
>>> hdr_arr:_fill_value = -9999.f ;
>>> hdr_arr:columns = "lat lon elv" ;
>>> hdr_arr:lat_long_name = "latitude" ;
>>> hdr_arr:lat_units = "degrees_north" ;
>>> hdr_arr:lon_long_name = "longitude" ;
>>> hdr_arr:lon_units = "degrees_east" ;
>>> hdr_arr:elv_long_name = "elevation" ;
>>> hdr_arr:elv_units = "meters above sea level (msl)"
;
>>>
>>> // global attributes:
>>> :FileOrigins = "File
>>> /h/data/global/WXQC/data/obs_prep/2010062200 generated
20100622_184241 UTC
>>> on
>>> host a31n09.afwa.af.mil by the PB2NC tool" ;
>>> }
>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_arr:long_name
>>> obs_2010062200.nc
>>> ncdump: NetCDF: Variable not found
>>> netcdf obs_2010062200 {
>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>>> obs_2010062200.nc
>>> ncdump: NetCDF: Variable not found
>>> netcdf obs_2010062200 {
>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>>> obs_2010062200.nc
>>> ncdump: NetCDF: Variable not found
>>> netcdf obs_2010062200 {
>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_typ:long_name
>>> obs_2010062200.nc
>>> ncdump: NetCDF: Variable not found
>>> netcdf obs_2010062200 {
>>>
>>>
>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>
>>>> I am getting the same 0 pair error with the attached ob file. I
built this
>>>> file using the attached config file. Is there surface
observations in my
>>>> file? Is the valid time of the file correct? I am using the
-dump option
>>>> but
>>>> I am not getting any file showing up in the directory I set for
this.
>>>>
>>>> Bob
>>>>
>>>> -----Original Message-----
>>>> From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
>>>> Sent: Tuesday, June 22, 2010 3:30 PM
>>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>
>>>> Bob,
>>>>
>>>> We looked more closely at the observation data that you included
in your
>>>> original email, and we noticed that they don't appear to be
surface
>>>> observations. When verifying 2m temp or 10m winds in point_stat,
only
>>>> observations of type ADPSFC or SFCSHP can be used (or ONLYSF for
either).
>>>>
>>>> Please change your point_stat config fcst_field setting to the
following
>>>> value:
>>>>
>>>> fcst_field[] = [ "33/Z10", "34/Z10"];
>>>>
>>>> also, please change the message type:
>>>>
>>>> message_type[] = [ "ADPSFC" ];
>>>>
>>>> Please let us know if you have any more questions.
>>>>
>>>> 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=38629 >
>>>>>
>>>>> Thanks for looking into this. I inadvertently left off my
config file
>>>>> name
>>>>> from the command string, hence it defaulted to profiler data. My
config
>>>>> file
>>>>> has "ANYSFC" which the documentation says includes ADPUPA (RAOB
data).
>>>>>
>>>>> For the observations, I am pulling the prepbufr files from
>>>>> ftp.ncep.noaa.gov
>>>>> in the /pub/data/nccf/com/gfs/prod/gdas.20100615 for the 15 Jun
files for
>>>>> the case below. An example of the filenames I am pulling are
>>>>> gdas2.t12z.prepbufr.nr. I rename these files to 20100615_12 in
this
>>>>> case.
>>>>> Since this file came out of the 20100615 directory I would
expect it to be
>>>>> valid for 20100615 at 12Z. Is there a way I can dump the header
of the
>>>>> prepbufr file to see the date that's in the file?
>>>>>
>>>>> Bob
>>>>>
>>>>>
>>>>> -----Original Message-----
>>>>> From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
>>>>> Sent: Friday, June 18, 2010 11:56 AM
>>>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>>
>>>>> Robert,
>>>>>
>>>>> Your observation file has a message type of PROFLR, but your
config file
>>>>> is
>>>>> specifying a message type of ADPSFC. This is one potential
problem.
>>>>> Another
>>>>> problem is that the model output appears to be valid at
20100615_12 and
>>>>> the
>>>>> observations are valid 24 hours later. You could change your
config file
>>>>> to
>>>>> open the valid time window to a larger value using beg_ds and
end_ds.
>>>>> Please
>>>>> let us know if you have any other questions.
>>>>>
>>>>> Thanks,
>>>>>
>>>>> Paul
>>>>>
>>>>>
>>>>>
>>>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN}
wrote:
>>>>>> Thu Jun 17 15:45:21 2010: Request 38629 was acted upon.
>>>>>> Transaction: Ticket created by craigr at offutt.af.mil
>>>>>> Queue: met_help
>>>>>> Subject: Point Stat Question
>>>>>> Owner: Nobody
>>>>>> Requestors: craigr at offutt.af.mil
>>>>>> Status: new
>>>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>>>
>>>>>>
>>>>>> I am running point stat on WRF winds but I am not getting any
stats. My
>>>>>> code to run this is:
>>>>>>
>>>>>>
>>>>>>
>>>>>> /sphome/qcteam/METv2.0/bin/point_stat
>>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>>> /h/data/global/WXQC/data/obs_pr
>>>>>>
>>>>>> ep/2010061512
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>> -outdir /h/data/global/WXQC/data/met/data -v 2
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> Attached is the config file along with the ob file and model
file.
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> The following is the output I get. Why the 0 pairs? Thanks
Bob Craig
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> GSL_RNG_TYPE=mt19937
>>>>>>
>>>>>> GSL_RNG_SEED=3818134092
>>>>>>
>>>>>> Forecast File:
>>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>>>
>>>>>> Climatology File: none
>>>>>>
>>>>>> Configuration File:
>>>>>> /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>>
>>>>>> Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>>>>>>
>>>>>>
>>>>>>
>>>>>> ----------------------------------------
>>>>>>
>>>>>>
>>>>>>
>>>>>> Reading records for UGRD/L10.
>>>>>>
>>>>>> For UGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>>
>>>>>>
>>>>>>
>>>>>> ----------------------------------------
>>>>>>
>>>>>>
>>>>>>
>>>>>> Reading records for VGRD/L10.
>>>>>>
>>>>>> For VGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>>
>>>>>>
>>>>>>
>>>>>> ----------------------------------------
>>>>>>
>>>>>>
>>>>>>
>>>>>> Searching 10473 observations from 170 PrepBufr messages.
>>>>>>
>>>>>>
>>>>>>
>>>>>> ----------------------------------------
>>>>>>
>>>>>>
>>>>>>
>>>>>> Processing UGRD/L10 versus UGRD/L10, for observation type
ADPSFC, over
>>>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>>>
>>>>>>
>>>>>>
>>>>>> ----------------------------------------
>>>>>>
>>>>>>
>>>>>>
>>>>>> Processing VGRD/L10 versus VGRD/L10, for observation type
ADPSFC, over
>>>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>>>
>>>>>>
>>>>>>
>>>>>> ----------------------------------------
>>>>>>
>>>>>>
>>>>>>
>>>>>> Output file:
>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
>>>>>> at
>>>>>>
>>>>>> Output file:
>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>>> c.txt
>>>>>>
>>>>>> Output file:
>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>>> s.txt
>>>>>>
>>>>>> Output file:
>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
>>>>>> t.txt
>>>>>>
>>>>>> Output file:
>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
>>>>>> 1l2.txt
>>>>>>
>>>>>> Output file:
>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
>>>>>> 1l2.txt
>>>>>>
>>
>>
>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #38629] Point Stat Question
From: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Time: Thu Jun 24 13:46:52 2010
I am recompiling those files. On the xlf compile, compiling in the
bufrlib directory I am getting the following errors:
"addate.f", line
460.61: 1515-018 (S) Incorrect character found in source at line 460
in column 61. Hexadecimal value of character
is 00.
"addate.f",
line 460.62: 1515-018 (S) Incorrect character found in source at line
460 in column 62. Hexadecimal value of character
is 00.
"addate.f", line 460.63: 1515-018 (S) Incorrect character found in
source at line 460 in column 63. Hexadecimal value of character
is
00.
"addate.f", line 460.64: 1515-018 (S) Incorrect character found
in source at line 460 in column 64. Hexadecimal value of character
is 00.
"addate.f", line 460.65: 1515-018 (S) Incorrect character
found in source at line 460 in column 65. Hexadecimal value of
character
is 00.
"addate.f", line 460.66: 1515-018 (S) Incorrect
character found in source at line 460 in column 66. Hexadecimal value
of character
is 00.
"addate.f", line 460.67: 1515-018 (S) Incorrect
character found in source at line 460 in column 67. Hexadecimal value
of character
is 00.
"addate.f", line 460.68: 1515-018 (S) Incorrect
character found in source at line 460 in column 68. Hexadecimal value
of character
is 00.
"addate.f", line 460.69: 1515-018 (S) Incorrect
character found in source at line 460 in column 69. Hexadecimal value
of character
is 00.
"addate.f", line 460.70: 1515-018 (S) Incorrect
character found in source at line 460 in column 70. Hexadecimal value
of character
is 00.
"addate.f", line 460.71: 1515-018 (S) Incorrect
character found in source at line 460 in column 71. Hexadecimal value
of character
is 00.
"addate.f", line 460.72: 1515-018 (S) Incorrect
character found in source at line 460 in column 72. Hexadecimal value
of character
"addate.f", line 500.6: 1514-050 (S) Specification
statement is out of order. Statement is ignored.
The compile
appears to still be running. What would cause this?
-----Original Message-----
From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
Sent: Thursday, June 24, 2010 2:24 PM
To:
Craig, Robert J Civ USAF AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu
#38629] Point Stat Question
Bob,
John tells me that recompiling
BUFRLIB using the following commands should fix
this problem:
xlc
-c -DUNDERSCORE *.c
xlf -c -qextname *.f *.F
ar crv libbufr.a
*.o
These are taken from this page:
http://www.dtcenter.org/met/users/support/online_tutorial/METv2.0/compilation/req_libs.php
And then try rebuilding MET.
Thanks,
Paul
Paul Oldenburg
wrote:
> Bob,
>
> Can you please redirect the make output to a log
file and send us the entire
> make output? That will help us
diagnose.
>
> 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=38629 >
>>
>>
Gmake solved that problem - compiled with the following errors (these
were at the bottom of the listing):
>>
>> ld: 0711-224 WARNING:
Duplicate symbol: .strlenFPCc
>> ld: 0711-345 Use the -bloadmap or
-bnoquiet option to obtain more information.
>> ld: 0711-317 ERROR:
Undefined symbol: .openbf_
>> ld: 0711-317 ERROR: Undefined symbol:
.datelen_
>> ld: 0711-317 ERROR: Undefined symbol: .ireadns_
>> ld:
0711-317 ERROR: Undefined symbol: .ufbint_
>> ld: 0711-317 ERROR:
Undefined symbol: .ufbevn_
>> ld: 0711-317 ERROR: Undefined symbol:
.ufbtab_
>> gmake[3]: *** [pb2nc] Error 8
>> gmake[3]: Leaving
directory `/h/sphome/qcteam/METv2.0/src/pb2nc'
>> gmake[2]: *** [all]
Error 2
>> gmake[2]: Leaving directory
`/h/sphome/qcteam/METv2.0/src/pb2nc'
>> gmake[1]: *** [targets] Error
2
>> gmake[1]: Leaving directory `/h/sphome/qcteam/METv2.0'
>>
gmake: *** [all] Error 2
>>
>> Should I worry about these - what is
the fix?
>>
>> Thanks
>>
>> -----Original Message-----
>> From:
RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent:
Thursday, June 24, 2010 2:00 PM
>> To: Craig, Robert J Civ USAF AFWA
16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat
Question
>>
>> Bob,
>>
>> These errors look like errors that are
output by the IBM version of make. MET
>> needs to be built with GNU
make, which is called gmake on our IBM machine.
>> Please look on
your IBM machine to see if you have gmake installed and use it
>>
instead.
>>
>> 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=38629 >
>>>
>>> I
will reset the quality threshold.
>>>
>>> I download the updates to
MET and placed in the top level MET dir. When I run
>>> a Make clean
I get the following errors:
>>>
>>> [496 qcteam a31n09
/sphome/qcteam/METv2.0] make clean
>>> "Makefile", line 191: make:
1254-055 Dependency line needs colon or double
>>> colon operator.
>>> "Makefile", line 193: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 194: make: 1254-
055 Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 196: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 197: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 199: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 200: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 202: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 203: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 205: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 206: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 208: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 209: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 211: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 212: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 214: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 215: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 217: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 218: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 220: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> make: 1254-058 Fatal errors
encountered -- cannot continue.
>>>
>>> The make file is attached.
What is causing this error?
>>>
>>> Thanks
>>>
>>> Bob
>>>
>>>
>>>
>>> -----Original Message-----
>>> From: RAL HelpDesk {for Paul
Oldenburg} [mailto:met_help at ucar.edu]
>>> Sent: Thursday, June 24,
2010 11:22 AM
>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>
>>>
Bob,
>>>
>>> I sent that email to John through the met_help
ticketing system. It was not
>>> intended for you, sorry.
>>>
>>>
We tested the PB2NCConfig file that you sent, and we found that it did
not
>>> process the message types that are needed by point_stat. The
reason was that
>>> the quality_mark_thresh setting was set to 1
instead of 2, and the message
>>> types
>>> ADPSFC and ADPUPA were
not processed as a result. Please try the following
>>> setting in
your PB2NCConfig file:
>>>
>>> quality_mark_thresh = 2;
>>>
>>>
Also, it is important to make sure that you are running the latest MET
>>> patches.
>>> If you have not done so, please download and deploy
the patches from the MET
>>> website at:
>>>
>>>
http://www.dtcenter.org/met/users/support/known_issues/METv2.0/index.php
>>>
>>> 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=38629 >
>>>>
>>>>
John, you sent this to me - did you intend someone else?
>>>>
>>>>
Bob
>>>>
>>>> -----Original Message-----
>>>> From: RAL HelpDesk
{for Paul Oldenburg} [mailto:met_help at ucar.edu]
>>>> Sent: Wednesday,
June 23, 2010 10:58 AM
>>>> To: Craig, Robert J Civ USAF AFWA 16
WS/WXN
>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>
>>>> John,
>>>>
>>>> I'm trying to look into Bob's problems,
but I'm struggling with ncdump. Am
>>>> I
>>>> doing something
wrong below where I attempt to view the variable values?
>>>>
>>>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -h
obs_2010062200.nc
>>>> netcdf obs_2010062200 {
>>>> dimensions:
>>>> mxstr = 16 ;
>>>> hdr_arr_len = 3 ;
>>>>
obs_arr_len = 5 ;
>>>> nobs = UNLIMITED ; // (10004
currently)
>>>> nhdr = 162 ;
>>>> variables:
>>>>
float obs_arr(nobs, obs_arr_len) ;
>>>>
obs_arr:long_name = "array of observation values" ;
>>>>
obs_arr:_fill_value = -9999.f ;
>>>> obs_arr:columns
= "hdr_id gc lvl hgt ob" ;
>>>>
obs_arr:hdr_id_long_name = "index of matching header data" ;
>>>>
obs_arr:gc_long_name = "grib code corresponding to the
>>>>
observation type" ;
>>>> obs_arr:lvl_long_name =
"pressure level (hPa) or
>>>> accumulation
>>>> interval (h)" ;
>>>> obs_arr:hgt_long_name = "height in meters above
sea level
>>>> (msl)" ;
>>>> obs_arr:ob_long_name =
"observation value" ;
>>>> char hdr_typ(nhdr, mxstr) ;
>>>>
hdr_typ:long_name = "message type" ;
>>>> char hdr_sid(nhdr,
mxstr) ;
>>>> hdr_sid:long_name = "station
identification" ;
>>>> char hdr_vld(nhdr, mxstr) ;
>>>>
hdr_vld:long_name = "valid time" ;
>>>> hdr_vld:units
= "YYYYMMDD_HHMMSS" ;
>>>> float hdr_arr(nhdr, hdr_arr_len) ;
>>>> hdr_arr:long_name = "array of observation station
header
>>>> values" ;
>>>> hdr_arr:_fill_value =
-9999.f ;
>>>> hdr_arr:columns = "lat lon elv" ;
>>>> hdr_arr:lat_long_name = "latitude" ;
>>>>
hdr_arr:lat_units = "degrees_north" ;
>>>>
hdr_arr:lon_long_name = "longitude" ;
>>>>
hdr_arr:lon_units = "degrees_east" ;
>>>>
hdr_arr:elv_long_name = "elevation" ;
>>>>
hdr_arr:elv_units = "meters above sea level (msl)" ;
>>>>
>>>> //
global attributes:
>>>> :FileOrigins = "File
>>>>
/h/data/global/WXQC/data/obs_prep/2010062200 generated 20100622_184241
UTC
>>>> on
>>>> host a31n09.afwa.af.mil by the PB2NC tool" ;
>>>>
}
>>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_arr:long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF:
Variable not found
>>>> netcdf obs_2010062200 {
>>>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF:
Variable not found
>>>> netcdf obs_2010062200 {
>>>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF:
Variable not found
>>>> netcdf obs_2010062200 {
>>>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_typ:long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF:
Variable not found
>>>> netcdf obs_2010062200 {
>>>>
>>>>
>>>> RAL
HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>>>>
<URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>>
>>>>> I am getting the same 0 pair error with the attached ob file. I
built this
>>>>> file using the attached config file. Is there
surface observations in my
>>>>> file? Is the valid time of the file
correct? I am using the -dump option
>>>>> but
>>>>> I am not
getting any file showing up in the directory I set for this.
>>>>>
>>>>> Bob
>>>>>
>>>>> -----Original Message-----
>>>>> From: RAL
HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>>>>> Sent:
Tuesday, June 22, 2010 3:30 PM
>>>>> To: Craig, Robert J Civ USAF
AFWA 16 WS/WXN
>>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat
Question
>>>>>
>>>>> Bob,
>>>>>
>>>>> We looked more closely at
the observation data that you included in your
>>>>> original email,
and we noticed that they don't appear to be surface
>>>>>
observations. When verifying 2m temp or 10m winds in point_stat, only
>>>>> observations of type ADPSFC or SFCSHP can be used (or ONLYSF for
either).
>>>>>
>>>>> Please change your point_stat config fcst_field
setting to the following
>>>>> value:
>>>>>
>>>>> fcst_field[] = [
"33/Z10", "34/Z10"];
>>>>>
>>>>> also, please change the message
type:
>>>>>
>>>>> message_type[] = [ "ADPSFC" ];
>>>>>
>>>>>
Please let us know if you have any more questions.
>>>>>
>>>>>
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=38629 >
>>>>>>
>>>>>> Thanks for looking into this. I inadvertently left off my
config file
>>>>>> name
>>>>>> from the command string, hence it
defaulted to profiler data. My config
>>>>>> file
>>>>>> has
"ANYSFC" which the documentation says includes ADPUPA (RAOB data).
>>>>>>
>>>>>> For the observations, I am pulling the prepbufr files
from
>>>>>> ftp.ncep.noaa.gov
>>>>>> in the
/pub/data/nccf/com/gfs/prod/gdas.20100615 for the 15 Jun files for
>>>>>> the case below. An example of the filenames I am pulling are
>>>>>> gdas2.t12z.prepbufr.nr. I rename these files to 20100615_12
in this
>>>>>> case.
>>>>>> Since this file came out of the 20100615
directory I would expect it to be
>>>>>> valid for 20100615 at 12Z.
Is there a way I can dump the header of the
>>>>>> prepbufr file to
see the date that's in the file?
>>>>>>
>>>>>> Bob
>>>>>>
>>>>>>
>>>>>> -----Original Message-----
>>>>>> From: RAL HelpDesk {for Paul
Oldenburg} [mailto:met_help at ucar.edu]
>>>>>> Sent: Friday, June 18,
2010 11:56 AM
>>>>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>>>
>>>>>> Robert,
>>>>>>
>>>>>> Your observation file has a
message type of PROFLR, but your config file
>>>>>> is
>>>>>>
specifying a message type of ADPSFC. This is one potential problem.
>>>>>> Another
>>>>>> problem is that the model output appears to be
valid at 20100615_12 and
>>>>>> the
>>>>>> observations are valid 24
hours later. You could change your config file
>>>>>> to
>>>>>>
open the valid time window to a larger value using beg_ds and end_ds.
>>>>>> Please
>>>>>> let us know if you have any other questions.
>>>>>>
>>>>>> Thanks,
>>>>>>
>>>>>> Paul
>>>>>>
>>>>>>
>>>>>>
>>>>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN}
wrote:
>>>>>>> Thu Jun 17 15:45:21 2010: Request 38629 was acted
upon.
>>>>>>> Transaction: Ticket created by craigr at offutt.af.mil
>>>>>>> Queue: met_help
>>>>>>> Subject: Point Stat
Question
>>>>>>> Owner: Nobody
>>>>>>> Requestors:
craigr at offutt.af.mil
>>>>>>> Status: new
>>>>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>>>>
>>>>>>>
>>>>>>> I am running point stat on WRF winds but I am not
getting any stats. My
>>>>>>> code to run this is:
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> /sphome/qcteam/METv2.0/bin/point_stat
>>>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>>>> /h/data/global/WXQC/data/obs_pr
>>>>>>>
>>>>>>>
ep/2010061512 /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>>> -outdir /h/data/global/WXQC/data/met/data -v 2
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Attached is the config
file along with the ob file and model file.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> The following is the output I get. Why the
0 pairs? Thanks Bob Craig
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> GSL_RNG_TYPE=mt19937
>>>>>>>
>>>>>>>
GSL_RNG_SEED=3818134092
>>>>>>>
>>>>>>> Forecast File:
>>>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>>>>
>>>>>>> Climatology File: none
>>>>>>>
>>>>>>>
Configuration File:
>>>>>>>
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>>>
>>>>>>> Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Reading records for UGRD/L10.
>>>>>>>
>>>>>>> For UGRD/L10
found 1 forecast levels and 0 climatology levels.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> ----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Reading records for VGRD/L10.
>>>>>>>
>>>>>>> For VGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Searching 10473 observations from 170 PrepBufr messages.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Processing UGRD/L10 versus UGRD/L10, for observation type
ADPSFC, over
>>>>>>> region FULL, for interpolation method MAX(4),
using 0 pairs.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Processing VGRD/L10 versus VGRD/L10, for observation type
ADPSFC, over
>>>>>>> region FULL, for interpolation method MAX(4),
using 0 pairs.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
>>>>>>> at
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>>>> c.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>>>> s.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
>>>>>>> t.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
>>>>>>> 1l2.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
>>>>>>> 1l2.txt
>>>>>>>
>>>
>>
>
>
------------------------------------------------
Subject: Point Stat Question
From: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Time: Thu Jun 24 13:50:17 2010
Here is the complete listing of the xlf -c -qextname *.f *.F compile:
Bob
-----Original Message-----
From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
Sent: Thursday, June 24, 2010 2:24 PM
To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
Bob,
John tells me that recompiling BUFRLIB using the following commands
should fix
this problem:
xlc -c -DUNDERSCORE *.c
xlf -c -qextname *.f *.F
ar crv libbufr.a *.o
These are taken from this page:
http://www.dtcenter.org/met/users/support/online_tutorial/METv2.0/compilation/req_libs.php
And then try rebuilding MET.
Thanks,
Paul
Paul Oldenburg wrote:
> Bob,
>
> Can you please redirect the make output to a log file and send us
the entire
> make output? That will help us diagnose.
>
> 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=38629 >
>>
>> Gmake solved that problem - compiled with the following errors
(these were
>> at the bottom of the listing):
>>
>> ld: 0711-224 WARNING: Duplicate symbol: .strlenFPCc
>> ld: 0711-345 Use the -bloadmap or -bnoquiet option to obtain more
>> information.
>> ld: 0711-317 ERROR: Undefined symbol: .openbf_
>> ld: 0711-317 ERROR: Undefined symbol: .datelen_
>> ld: 0711-317 ERROR: Undefined symbol: .ireadns_
>> ld: 0711-317 ERROR: Undefined symbol: .ufbint_
>> ld: 0711-317 ERROR: Undefined symbol: .ufbevn_
>> ld: 0711-317 ERROR: Undefined symbol: .ufbtab_
>> gmake[3]: *** [pb2nc] Error 8
>> gmake[3]: Leaving directory `/h/sphome/qcteam/METv2.0/src/pb2nc'
>> gmake[2]: *** [all] Error 2
>> gmake[2]: Leaving directory `/h/sphome/qcteam/METv2.0/src/pb2nc'
>> gmake[1]: *** [targets] Error 2
>> gmake[1]: Leaving directory `/h/sphome/qcteam/METv2.0'
>> gmake: *** [all] Error 2
>>
>> Should I worry about these - what is the fix?
>>
>> Thanks
>>
>> -----Original Message-----
>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent: Thursday, June 24, 2010 2:00 PM
>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>
>> Bob,
>>
>> These errors look like errors that are output by the IBM version of
make.
>> MET
>> needs to be built with GNU make, which is called gmake on our IBM
machine.
>> Please look on your IBM machine to see if you have gmake installed
and use
>> it
>> instead.
>>
>> 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=38629 >
>>>
>>> I will reset the quality threshold.
>>>
>>> I download the updates to MET and placed in the top level MET dir.
When I
>>> run
>>> a Make clean I get the following errors:
>>>
>>> [496 qcteam a31n09 /sphome/qcteam/METv2.0] make clean
>>> "Makefile", line 191: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 193: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 194: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 196: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 197: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 199: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 200: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 202: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 203: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 205: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 206: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 208: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 209: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 211: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 212: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 214: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 215: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 217: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 218: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 220: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> make: 1254-058 Fatal errors encountered -- cannot continue.
>>>
>>> The make file is attached. What is causing this error?
>>>
>>> Thanks
>>>
>>> Bob
>>>
>>>
>>>
>>> -----Original Message-----
>>> From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>>> Sent: Thursday, June 24, 2010 11:22 AM
>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>
>>> Bob,
>>>
>>> I sent that email to John through the met_help ticketing system.
It was
>>> not
>>> intended for you, sorry.
>>>
>>> We tested the PB2NCConfig file that you sent, and we found that it
did not
>>> process the message types that are needed by point_stat. The
reason was
>>> that
>>> the quality_mark_thresh setting was set to 1 instead of 2, and the
message
>>> types
>>> ADPSFC and ADPUPA were not processed as a result. Please try the
>>> following
>>> setting in your PB2NCConfig file:
>>>
>>> quality_mark_thresh = 2;
>>>
>>> Also, it is important to make sure that you are running the latest
MET
>>> patches.
>>> If you have not done so, please download and deploy the patches
from the
>>> MET
>>> website at:
>>>
>>>
http://www.dtcenter.org/met/users/support/known_issues/METv2.0/index.php
>>>
>>> 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=38629 >
>>>>
>>>> John, you sent this to me - did you intend someone else?
>>>>
>>>> Bob
>>>>
>>>> -----Original Message-----
>>>> From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
>>>> Sent: Wednesday, June 23, 2010 10:58 AM
>>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>
>>>> John,
>>>>
>>>> I'm trying to look into Bob's problems, but I'm struggling with
ncdump.
>>>> Am
>>>> I
>>>> doing something wrong below where I attempt to view the variable
values?
>>>>
>>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -h
obs_2010062200.nc
>>>> netcdf obs_2010062200 {
>>>> dimensions:
>>>> mxstr = 16 ;
>>>> hdr_arr_len = 3 ;
>>>> obs_arr_len = 5 ;
>>>> nobs = UNLIMITED ; // (10004 currently)
>>>> nhdr = 162 ;
>>>> variables:
>>>> float obs_arr(nobs, obs_arr_len) ;
>>>> obs_arr:long_name = "array of observation values"
;
>>>> obs_arr:_fill_value = -9999.f ;
>>>> obs_arr:columns = "hdr_id gc lvl hgt ob" ;
>>>> obs_arr:hdr_id_long_name = "index of matching
header
>>>> data" ;
>>>> obs_arr:gc_long_name = "grib code corresponding
to the
>>>> observation type" ;
>>>> obs_arr:lvl_long_name = "pressure level (hPa) or
>>>> accumulation
>>>> interval (h)" ;
>>>> obs_arr:hgt_long_name = "height in meters above
sea level
>>>> (msl)" ;
>>>> obs_arr:ob_long_name = "observation value" ;
>>>> char hdr_typ(nhdr, mxstr) ;
>>>> hdr_typ:long_name = "message type" ;
>>>> char hdr_sid(nhdr, mxstr) ;
>>>> hdr_sid:long_name = "station identification" ;
>>>> char hdr_vld(nhdr, mxstr) ;
>>>> hdr_vld:long_name = "valid time" ;
>>>> hdr_vld:units = "YYYYMMDD_HHMMSS" ;
>>>> float hdr_arr(nhdr, hdr_arr_len) ;
>>>> hdr_arr:long_name = "array of observation station
header
>>>> values" ;
>>>> hdr_arr:_fill_value = -9999.f ;
>>>> hdr_arr:columns = "lat lon elv" ;
>>>> hdr_arr:lat_long_name = "latitude" ;
>>>> hdr_arr:lat_units = "degrees_north" ;
>>>> hdr_arr:lon_long_name = "longitude" ;
>>>> hdr_arr:lon_units = "degrees_east" ;
>>>> hdr_arr:elv_long_name = "elevation" ;
>>>> hdr_arr:elv_units = "meters above sea level
(msl)" ;
>>>>
>>>> // global attributes:
>>>> :FileOrigins = "File
>>>> /h/data/global/WXQC/data/obs_prep/2010062200 generated
20100622_184241
>>>> UTC
>>>> on
>>>> host a31n09.afwa.af.mil by the PB2NC tool" ;
>>>> }
>>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_arr:long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF: Variable not found
>>>> netcdf obs_2010062200 {
>>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
>>>> obs_arr:lvl_long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF: Variable not found
>>>> netcdf obs_2010062200 {
>>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
>>>> obs_arr:lvl_long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF: Variable not found
>>>> netcdf obs_2010062200 {
>>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_typ:long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF: Variable not found
>>>> netcdf obs_2010062200 {
>>>>
>>>>
>>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>>
>>>>> I am getting the same 0 pair error with the attached ob file. I
built
>>>>> this
>>>>> file using the attached config file. Is there surface
observations in
>>>>> my
>>>>> file? Is the valid time of the file correct? I am using the
-dump
>>>>> option
>>>>> but
>>>>> I am not getting any file showing up in the directory I set for
this.
>>>>>
>>>>> Bob
>>>>>
>>>>> -----Original Message-----
>>>>> From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
>>>>> Sent: Tuesday, June 22, 2010 3:30 PM
>>>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>>
>>>>> Bob,
>>>>>
>>>>> We looked more closely at the observation data that you included
in your
>>>>> original email, and we noticed that they don't appear to be
surface
>>>>> observations. When verifying 2m temp or 10m winds in
point_stat, only
>>>>> observations of type ADPSFC or SFCSHP can be used (or ONLYSF for
>>>>> either).
>>>>>
>>>>> Please change your point_stat config fcst_field setting to the
following
>>>>> value:
>>>>>
>>>>> fcst_field[] = [ "33/Z10", "34/Z10"];
>>>>>
>>>>> also, please change the message type:
>>>>>
>>>>> message_type[] = [ "ADPSFC" ];
>>>>>
>>>>> Please let us know if you have any more questions.
>>>>>
>>>>> 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=38629 >
>>>>>>
>>>>>> Thanks for looking into this. I inadvertently left off my
config file
>>>>>> name
>>>>>> from the command string, hence it defaulted to profiler data.
My config
>>>>>> file
>>>>>> has "ANYSFC" which the documentation says includes ADPUPA (RAOB
data).
>>>>>>
>>>>>> For the observations, I am pulling the prepbufr files from
>>>>>> ftp.ncep.noaa.gov
>>>>>> in the /pub/data/nccf/com/gfs/prod/gdas.20100615 for the 15 Jun
files
>>>>>> for
>>>>>> the case below. An example of the filenames I am pulling are
>>>>>> gdas2.t12z.prepbufr.nr. I rename these files to 20100615_12
in this
>>>>>> case.
>>>>>> Since this file came out of the 20100615 directory I would
expect it to
>>>>>> be
>>>>>> valid for 20100615 at 12Z. Is there a way I can dump the
header of the
>>>>>> prepbufr file to see the date that's in the file?
>>>>>>
>>>>>> Bob
>>>>>>
>>>>>>
>>>>>> -----Original Message-----
>>>>>> From: RAL HelpDesk {for Paul Oldenburg}
[mailto:met_help at ucar.edu]
>>>>>> Sent: Friday, June 18, 2010 11:56 AM
>>>>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>>>
>>>>>> Robert,
>>>>>>
>>>>>> Your observation file has a message type of PROFLR, but your
config
>>>>>> file
>>>>>> is
>>>>>> specifying a message type of ADPSFC. This is one potential
problem.
>>>>>> Another
>>>>>> problem is that the model output appears to be valid at
20100615_12 and
>>>>>> the
>>>>>> observations are valid 24 hours later. You could change your
config
>>>>>> file
>>>>>> to
>>>>>> open the valid time window to a larger value using beg_ds and
end_ds.
>>>>>> Please
>>>>>> let us know if you have any other questions.
>>>>>>
>>>>>> Thanks,
>>>>>>
>>>>>> Paul
>>>>>>
>>>>>>
>>>>>>
>>>>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN}
wrote:
>>>>>>> Thu Jun 17 15:45:21 2010: Request 38629 was acted upon.
>>>>>>> Transaction: Ticket created by craigr at offutt.af.mil
>>>>>>> Queue: met_help
>>>>>>> Subject: Point Stat Question
>>>>>>> Owner: Nobody
>>>>>>> Requestors: craigr at offutt.af.mil
>>>>>>> Status: new
>>>>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629
>>>>>>> >
>>>>>>>
>>>>>>>
>>>>>>> I am running point stat on WRF winds but I am not getting any
stats.
>>>>>>> My
>>>>>>> code to run this is:
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> /sphome/qcteam/METv2.0/bin/point_stat
>>>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>>>> /h/data/global/WXQC/data/obs_pr
>>>>>>>
>>>>>>> ep/2010061512
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>>> -outdir /h/data/global/WXQC/data/met/data -v 2
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Attached is the config file along with the ob file and model
file.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> The following is the output I get. Why the 0 pairs? Thanks
Bob
>>>>>>> Craig
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> GSL_RNG_TYPE=mt19937
>>>>>>>
>>>>>>> GSL_RNG_SEED=3818134092
>>>>>>>
>>>>>>> Forecast File:
>>>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>>>>
>>>>>>> Climatology File: none
>>>>>>>
>>>>>>> Configuration File:
>>>>>>> /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>>>
>>>>>>> Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> ----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Reading records for UGRD/L10.
>>>>>>>
>>>>>>> For UGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> ----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Reading records for VGRD/L10.
>>>>>>>
>>>>>>> For VGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> ----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Searching 10473 observations from 170 PrepBufr messages.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> ----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Processing UGRD/L10 versus UGRD/L10, for observation type
ADPSFC, over
>>>>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> ----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Processing VGRD/L10 versus VGRD/L10, for observation type
ADPSFC, over
>>>>>>> region FULL, for interpolation method MAX(4), using 0 pairs.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> ----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
>>>>>>> at
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>>>> c.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>>>> s.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
>>>>>>> t.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
>>>>>>> 1l2.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
>>>>>>> 1l2.txt
>>>>>>>
>>>
>>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #38629] Point Stat Question
From: Craig, Robert J Civ USAF AFWA 16 WS/WXN
Time: Thu Jun 24 15:45:59 2010
I found that the addate.f file somehow got corrupted. I downloaded a
new version and bufrlib compiled as it should and MET compiled as it
should.
Thanks for your help.
Bob
-----Original Message-----
From: RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
Sent: Thursday, June 24, 2010 2:24 PM
To: Craig, Robert J Civ USAF
AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat
Question
Bob,
John tells me that recompiling BUFRLIB using the
following commands should fix
this problem:
xlc -c -DUNDERSCORE
*.c
xlf -c -qextname *.f *.F
ar crv libbufr.a *.o
These are
taken from this page:
http://www.dtcenter.org/met/users/support/online_tutorial/METv2.0/compilation/req_libs.php
And then try rebuilding MET.
Thanks,
Paul
Paul Oldenburg
wrote:
> Bob,
>
> Can you please redirect the make output to a log
file and send us the entire
> make output? That will help us
diagnose.
>
> 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=38629 >
>>
>>
Gmake solved that problem - compiled with the following errors (these
were at the bottom of the listing):
>>
>> ld: 0711-224 WARNING:
Duplicate symbol: .strlenFPCc
>> ld: 0711-345 Use the -bloadmap or
-bnoquiet option to obtain more information.
>> ld: 0711-317 ERROR:
Undefined symbol: .openbf_
>> ld: 0711-317 ERROR: Undefined symbol:
.datelen_
>> ld: 0711-317 ERROR: Undefined symbol: .ireadns_
>> ld:
0711-317 ERROR: Undefined symbol: .ufbint_
>> ld: 0711-317 ERROR:
Undefined symbol: .ufbevn_
>> ld: 0711-317 ERROR: Undefined symbol:
.ufbtab_
>> gmake[3]: *** [pb2nc] Error 8
>> gmake[3]: Leaving
directory `/h/sphome/qcteam/METv2.0/src/pb2nc'
>> gmake[2]: *** [all]
Error 2
>> gmake[2]: Leaving directory
`/h/sphome/qcteam/METv2.0/src/pb2nc'
>> gmake[1]: *** [targets] Error
2
>> gmake[1]: Leaving directory `/h/sphome/qcteam/METv2.0'
>>
gmake: *** [all] Error 2
>>
>> Should I worry about these - what is
the fix?
>>
>> Thanks
>>
>> -----Original Message-----
>> From:
RAL HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>> Sent:
Thursday, June 24, 2010 2:00 PM
>> To: Craig, Robert J Civ USAF AFWA
16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat
Question
>>
>> Bob,
>>
>> These errors look like errors that are
output by the IBM version of make. MET
>> needs to be built with GNU
make, which is called gmake on our IBM machine.
>> Please look on
your IBM machine to see if you have gmake installed and use it
>>
instead.
>>
>> 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=38629 >
>>>
>>> I
will reset the quality threshold.
>>>
>>> I download the updates to
MET and placed in the top level MET dir. When I run
>>> a Make clean
I get the following errors:
>>>
>>> [496 qcteam a31n09
/sphome/qcteam/METv2.0] make clean
>>> "Makefile", line 191: make:
1254-055 Dependency line needs colon or double
>>> colon operator.
>>> "Makefile", line 193: make: 1254-055 Dependency line needs colon
or double
>>> colon operator.
>>> "Makefile", line 194: make: 1254-
055 Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 196: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 197: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 199: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 200: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 202: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 203: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 205: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 206: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 208: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 209: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 211: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 212: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 214: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 215: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 217: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> "Makefile", line 218: make: 1254-055
Dependency line needs colon or double
>>> colon operator.
>>>
"Makefile", line 220: make: 1254-055 Dependency line needs colon or
double
>>> colon operator.
>>> make: 1254-058 Fatal errors
encountered -- cannot continue.
>>>
>>> The make file is attached.
What is causing this error?
>>>
>>> Thanks
>>>
>>> Bob
>>>
>>>
>>>
>>> -----Original Message-----
>>> From: RAL HelpDesk {for Paul
Oldenburg} [mailto:met_help at ucar.edu]
>>> Sent: Thursday, June 24,
2010 11:22 AM
>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>
Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>
>>>
Bob,
>>>
>>> I sent that email to John through the met_help
ticketing system. It was not
>>> intended for you, sorry.
>>>
>>>
We tested the PB2NCConfig file that you sent, and we found that it did
not
>>> process the message types that are needed by point_stat. The
reason was that
>>> the quality_mark_thresh setting was set to 1
instead of 2, and the message
>>> types
>>> ADPSFC and ADPUPA were
not processed as a result. Please try the following
>>> setting in
your PB2NCConfig file:
>>>
>>> quality_mark_thresh = 2;
>>>
>>>
Also, it is important to make sure that you are running the latest MET
>>> patches.
>>> If you have not done so, please download and deploy
the patches from the MET
>>> website at:
>>>
>>>
http://www.dtcenter.org/met/users/support/known_issues/METv2.0/index.php
>>>
>>> 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=38629 >
>>>>
>>>>
John, you sent this to me - did you intend someone else?
>>>>
>>>>
Bob
>>>>
>>>> -----Original Message-----
>>>> From: RAL HelpDesk
{for Paul Oldenburg} [mailto:met_help at ucar.edu]
>>>> Sent: Wednesday,
June 23, 2010 10:58 AM
>>>> To: Craig, Robert J Civ USAF AFWA 16
WS/WXN
>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>
>>>> John,
>>>>
>>>> I'm trying to look into Bob's problems,
but I'm struggling with ncdump. Am
>>>> I
>>>> doing something
wrong below where I attempt to view the variable values?
>>>>
>>>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -h
obs_2010062200.nc
>>>> netcdf obs_2010062200 {
>>>> dimensions:
>>>> mxstr = 16 ;
>>>> hdr_arr_len = 3 ;
>>>>
obs_arr_len = 5 ;
>>>> nobs = UNLIMITED ; // (10004
currently)
>>>> nhdr = 162 ;
>>>> variables:
>>>>
float obs_arr(nobs, obs_arr_len) ;
>>>>
obs_arr:long_name = "array of observation values" ;
>>>>
obs_arr:_fill_value = -9999.f ;
>>>> obs_arr:columns
= "hdr_id gc lvl hgt ob" ;
>>>>
obs_arr:hdr_id_long_name = "index of matching header data" ;
>>>>
obs_arr:gc_long_name = "grib code corresponding to the
>>>>
observation type" ;
>>>> obs_arr:lvl_long_name =
"pressure level (hPa) or
>>>> accumulation
>>>> interval (h)" ;
>>>> obs_arr:hgt_long_name = "height in meters above
sea level
>>>> (msl)" ;
>>>> obs_arr:ob_long_name =
"observation value" ;
>>>> char hdr_typ(nhdr, mxstr) ;
>>>>
hdr_typ:long_name = "message type" ;
>>>> char hdr_sid(nhdr,
mxstr) ;
>>>> hdr_sid:long_name = "station
identification" ;
>>>> char hdr_vld(nhdr, mxstr) ;
>>>>
hdr_vld:long_name = "valid time" ;
>>>> hdr_vld:units
= "YYYYMMDD_HHMMSS" ;
>>>> float hdr_arr(nhdr, hdr_arr_len) ;
>>>> hdr_arr:long_name = "array of observation station
header
>>>> values" ;
>>>> hdr_arr:_fill_value =
-9999.f ;
>>>> hdr_arr:columns = "lat lon elv" ;
>>>> hdr_arr:lat_long_name = "latitude" ;
>>>>
hdr_arr:lat_units = "degrees_north" ;
>>>>
hdr_arr:lon_long_name = "longitude" ;
>>>>
hdr_arr:lon_units = "degrees_east" ;
>>>>
hdr_arr:elv_long_name = "elevation" ;
>>>>
hdr_arr:elv_units = "meters above sea level (msl)" ;
>>>>
>>>> //
global attributes:
>>>> :FileOrigins = "File
>>>>
/h/data/global/WXQC/data/obs_prep/2010062200 generated 20100622_184241
UTC
>>>> on
>>>> host a31n09.afwa.af.mil by the PB2NC tool" ;
>>>>
}
>>>> pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_arr:long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF:
Variable not found
>>>> netcdf obs_2010062200 {
>>>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF:
Variable not found
>>>> netcdf obs_2010062200 {
>>>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
obs_arr:lvl_long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF:
Variable not found
>>>> netcdf obs_2010062200 {
>>>>
pgoldenb at kemosabe:~/met_help/20100618.craig$ ncdump -v
hdr_typ:long_name
>>>> obs_2010062200.nc
>>>> ncdump: NetCDF:
Variable not found
>>>> netcdf obs_2010062200 {
>>>>
>>>>
>>>> RAL
HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN} wrote:
>>>>>
<URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>>
>>>>> I am getting the same 0 pair error with the attached ob file. I
built this
>>>>> file using the attached config file. Is there
surface observations in my
>>>>> file? Is the valid time of the file
correct? I am using the -dump option
>>>>> but
>>>>> I am not
getting any file showing up in the directory I set for this.
>>>>>
>>>>> Bob
>>>>>
>>>>> -----Original Message-----
>>>>> From: RAL
HelpDesk {for Paul Oldenburg} [mailto:met_help at ucar.edu]
>>>>> Sent:
Tuesday, June 22, 2010 3:30 PM
>>>>> To: Craig, Robert J Civ USAF
AFWA 16 WS/WXN
>>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat
Question
>>>>>
>>>>> Bob,
>>>>>
>>>>> We looked more closely at
the observation data that you included in your
>>>>> original email,
and we noticed that they don't appear to be surface
>>>>>
observations. When verifying 2m temp or 10m winds in point_stat, only
>>>>> observations of type ADPSFC or SFCSHP can be used (or ONLYSF for
either).
>>>>>
>>>>> Please change your point_stat config fcst_field
setting to the following
>>>>> value:
>>>>>
>>>>> fcst_field[] = [
"33/Z10", "34/Z10"];
>>>>>
>>>>> also, please change the message
type:
>>>>>
>>>>> message_type[] = [ "ADPSFC" ];
>>>>>
>>>>>
Please let us know if you have any more questions.
>>>>>
>>>>>
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=38629 >
>>>>>>
>>>>>> Thanks for looking into this. I inadvertently left off my
config file
>>>>>> name
>>>>>> from the command string, hence it
defaulted to profiler data. My config
>>>>>> file
>>>>>> has
"ANYSFC" which the documentation says includes ADPUPA (RAOB data).
>>>>>>
>>>>>> For the observations, I am pulling the prepbufr files
from
>>>>>> ftp.ncep.noaa.gov
>>>>>> in the
/pub/data/nccf/com/gfs/prod/gdas.20100615 for the 15 Jun files for
>>>>>> the case below. An example of the filenames I am pulling are
>>>>>> gdas2.t12z.prepbufr.nr. I rename these files to 20100615_12
in this
>>>>>> case.
>>>>>> Since this file came out of the 20100615
directory I would expect it to be
>>>>>> valid for 20100615 at 12Z.
Is there a way I can dump the header of the
>>>>>> prepbufr file to
see the date that's in the file?
>>>>>>
>>>>>> Bob
>>>>>>
>>>>>>
>>>>>> -----Original Message-----
>>>>>> From: RAL HelpDesk {for Paul
Oldenburg} [mailto:met_help at ucar.edu]
>>>>>> Sent: Friday, June 18,
2010 11:56 AM
>>>>>> To: Craig, Robert J Civ USAF AFWA 16 WS/WXN
>>>>>> Subject: Re: [rt.rap.ucar.edu #38629] Point Stat Question
>>>>>>
>>>>>> Robert,
>>>>>>
>>>>>> Your observation file has a
message type of PROFLR, but your config file
>>>>>> is
>>>>>>
specifying a message type of ADPSFC. This is one potential problem.
>>>>>> Another
>>>>>> problem is that the model output appears to be
valid at 20100615_12 and
>>>>>> the
>>>>>> observations are valid 24
hours later. You could change your config file
>>>>>> to
>>>>>>
open the valid time window to a larger value using beg_ds and end_ds.
>>>>>> Please
>>>>>> let us know if you have any other questions.
>>>>>>
>>>>>> Thanks,
>>>>>>
>>>>>> Paul
>>>>>>
>>>>>>
>>>>>>
>>>>>> RAL HelpDesk {for Craig, Robert J Civ USAF AFWA 16 WS/WXN}
wrote:
>>>>>>> Thu Jun 17 15:45:21 2010: Request 38629 was acted
upon.
>>>>>>> Transaction: Ticket created by craigr at offutt.af.mil
>>>>>>> Queue: met_help
>>>>>>> Subject: Point Stat
Question
>>>>>>> Owner: Nobody
>>>>>>> Requestors:
craigr at offutt.af.mil
>>>>>>> Status: new
>>>>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=38629 >
>>>>>>>
>>>>>>>
>>>>>>> I am running point stat on WRF winds but I am not
getting any stats. My
>>>>>>> code to run this is:
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> /sphome/qcteam/METv2.0/bin/point_stat
>>>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>>>> /h/data/global/WXQC/data/obs_pr
>>>>>>>
>>>>>>>
ep/2010061512 /sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>>> -outdir /h/data/global/WXQC/data/met/data -v 2
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Attached is the config
file along with the ob file and model file.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> The following is the output I get. Why the
0 pairs? Thanks Bob Craig
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> GSL_RNG_TYPE=mt19937
>>>>>>>
>>>>>>>
GSL_RNG_SEED=3818134092
>>>>>>>
>>>>>>> Forecast File:
>>>>>>>
/gds/data/global/WXQC/data/wrf/t2/qc_100615_us057g1011t02b060000600
>>>>>>>
>>>>>>> Climatology File: none
>>>>>>>
>>>>>>>
Configuration File:
>>>>>>>
/sphome/qcteam/METv2.0/data/config/point_stat_config_WRF
>>>>>>>
>>>>>>> Observation File: /h/data/global/WXQC/data/obs_prep/2010061512
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Reading records for UGRD/L10.
>>>>>>>
>>>>>>> For UGRD/L10
found 1 forecast levels and 0 climatology levels.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> ----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Reading records for VGRD/L10.
>>>>>>>
>>>>>>> For VGRD/L10 found 1 forecast levels and 0 climatology levels.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Searching 10473 observations from 170 PrepBufr messages.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Processing UGRD/L10 versus UGRD/L10, for observation type
ADPSFC, over
>>>>>>> region FULL, for interpolation method MAX(4),
using 0 pairs.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Processing VGRD/L10 versus VGRD/L10, for observation type
ADPSFC, over
>>>>>>> region FULL, for interpolation method MAX(4),
using 0 pairs.
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
----------------------------------------
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V.st
>>>>>>> at
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>>>> c.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_ct
>>>>>>> s.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_cn
>>>>>>> t.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_sl
>>>>>>> 1l2.txt
>>>>>>>
>>>>>>> Output file:
>>>>>>>
/h/data/global/WXQC/data/met/data/point_stat_WRF_060000L_20100615_120000V_vl
>>>>>>> 1l2.txt
>>>>>>>
>>>
>>
>
>
------------------------------------------------
More information about the Met_help
mailing list