[Met_help] [rt.rap.ucar.edu #62156] History for Issue using Mode tool
John Halley Gotway via RT
met_help at ucar.edu
Tue Aug 27 10:38:25 MDT 2013
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello,
I am attempting to use the Mode tool to analyze model performance regarding
precipitation forecast. After using the pcp_combine tool to get the
appropriate accumulation variable (6hr), I tried to run the Mode tool as
below:
wrf at srbwks35369-02-> bin/mode \
/data/wrf/wrf_mode/pcp_comb_out/wrf_pcp_d02_2005_01_01_12_6hr.grb \
/data/wrf/precip_obs/ST4.2005010112.06h \
00_for_frank/mode_config \
-outdir /data/wrf/wrf_mode/ \
-v 2
The first input is the file that was created using pcp_combine, and the
second input is the Stage IV observed precipitation data from the EOL data
archive. This results in the following error:
ERROR: GCInfo::set_gcinfo_grib() -> each GRIB code specified must be
followed by an accumulation, level, or presssure level indicator "APCP_A6".
The config file is setup with the comparison variables written as such:
fcst_field = "APCP_A6";
obs_field = "APCP_A6";
I'm pretty sure I have those correct, but the mode tool won't accept that.
To determine if the relevant data was present, I imported each of the files
into IDV and was able to plot the accumulations for each, but I did notice
in looking at the details of each file in IDV that the pcp_combine-created
file's field is named "APCP" with the level "A6", while the Stage IV obs
file's field is also "APCP," but with no designated level in the metadata
that I can find. Metadata for both files is attached.
Thank you for your help,
--
Zachary Ebenstein
University of Michigan
College of Engineering
Dept. Atmospheric, Oceanic and Space Sciences
Climate Physics (B.S.E), Dec. 2011
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #62156] Issue using Mode tool
From: John Halley Gotway
Time: Tue Jul 09 15:19:09 2013
Zachary,
I see that you're having problems getting the configuration correct
for the MODE tool. And since you're setting "fcst_field" and
"obs_field", I assume you're using METv3.1. If that's not correct,
please let me know.
First, let me tell you that there were two major issues discovered in
MODE that have been corrected in METv4.1. They are described in the
"Existing Tools" section of the METv4.1 release notes:
http://www.dtcenter.org/met/users/support/release_notes/METv4.1_release_notes.php
So if you're able to switch to the newer version, I'd suggest doing
so. One warning though - the configuration files have been updated a
lot since METv3.1.
Regarding the problem at hand, you're using MODE to compare a NetCDF
file (wrf_pcp_d02_2005_01_01_12_6hr.grb) to a GRIB file
(ST4.2005010112.06h). First, naming the NetCDF output of pcp_combine
with
a ".grb" extension was a poor choice. The MET tools assume that files
ending in .grb are GRIB and try to process them as such. There's a
way to override that behavior in METv4.1 but not in METv3.1.
I'd suggest using a ".nc" extension.
Also, the way the field information is specified differs depending on
the file type. For GRIB, we use the GRIB code abbreviation followed
by the accumulation interval to be used. For NetCDF, we use
the NetCDF variable name followed by the dimensions to be used. So
please try using the following:
fcst_field = "APCP_06(*,*)"; # for NetCDF
obs_field = "APCP/A6"; # for GRIB
Please give that a shot and let me know how it goes.
Also, I was a bit confused about the GRIB metadata you sent. That
looks like NetCDF metadata, but I'm assuming you converted the GRIB
file to NetCDF to see what it contains. I'd suggest using the
"wgrib" utility to inventory the contents of a GRIB file directly.
Hope that helps.
Thanks,
John Halley Gotway
met_help at ucar.edu
On 07/09/2013 12:52 PM, Zachary Ebenstein via RT wrote:
>
> Tue Jul 09 12:52:36 2013: Request 62156 was acted upon.
> Transaction: Ticket created by ebenszac at umich.edu
> Queue: met_help
> Subject: Issue using Mode tool
> Owner: Nobody
> Requestors: ebenszac at umich.edu
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62156 >
>
>
> Hello,
>
> I am attempting to use the Mode tool to analyze model performance
regarding
> precipitation forecast. After using the pcp_combine tool to get the
> appropriate accumulation variable (6hr), I tried to run the Mode
tool as
> below:
>
> wrf at srbwks35369-02-> bin/mode \
> /data/wrf/wrf_mode/pcp_comb_out/wrf_pcp_d02_2005_01_01_12_6hr.grb \
> /data/wrf/precip_obs/ST4.2005010112.06h \
> 00_for_frank/mode_config \
> -outdir /data/wrf/wrf_mode/ \
> -v 2
>
> The first input is the file that was created using pcp_combine, and
the
> second input is the Stage IV observed precipitation data from the
EOL data
> archive. This results in the following error:
>
> ERROR: GCInfo::set_gcinfo_grib() -> each GRIB code specified must be
> followed by an accumulation, level, or presssure level indicator
"APCP_A6".
>
> The config file is setup with the comparison variables written as
such:
>
> fcst_field = "APCP_A6";
> obs_field = "APCP_A6";
>
> I'm pretty sure I have those correct, but the mode tool won't accept
that.
> To determine if the relevant data was present, I imported each of
the files
> into IDV and was able to plot the accumulations for each, but I did
notice
> in looking at the details of each file in IDV that the pcp_combine-
created
> file's field is named "APCP" with the level "A6", while the Stage IV
obs
> file's field is also "APCP," but with no designated level in the
metadata
> that I can find. Metadata for both files is attached.
>
> Thank you for your help,
>
------------------------------------------------
Subject: Issue using Mode tool
From: Zachary Ebenstein
Time: Fri Jul 12 11:23:19 2013
John,
I am using METv3, and unfortunately I am not able to switch to METv4
at
this time. As for the naming of the pcp_combine output, I had
mistakenly
thought the output was still in GRIB format. I have produced new
output
with the .nc file extension, as well as with no extension, with the
config
changed to include these fields:
fcst_field = "APCP_A6(*,*)";
obs_field = "APCP_A6";
Unfortunately, I have still had the same error as a result. The
metadata
for the observation file looks like this:
wrf at srbwks35369-02-> wgrib /data/wrf/precip_obs/ST4.2005010112.06h
1:0:d=04123112:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=18:P2=24:TimeU=1:sfc:18-
24hr
acc:NAve=0
Zac
On Tue, Jul 9, 2013 at 5:19 PM, John Halley Gotway via RT
<met_help at ucar.edu
> wrote:
> Zachary,
>
> I see that you're having problems getting the configuration correct
for
> the MODE tool. And since you're setting "fcst_field" and
"obs_field", I
> assume you're using METv3.1. If that's not correct,
> please let me know.
>
> First, let me tell you that there were two major issues discovered
in MODE
> that have been corrected in METv4.1. They are described in the
"Existing
> Tools" section of the METv4.1 release notes:
>
>
http://www.dtcenter.org/met/users/support/release_notes/METv4.1_release_notes.php
>
> So if you're able to switch to the newer version, I'd suggest doing
so.
> One warning though - the configuration files have been updated a
lot since
> METv3.1.
>
> Regarding the problem at hand, you're using MODE to compare a NetCDF
file
> (wrf_pcp_d02_2005_01_01_12_6hr.grb) to a GRIB file
(ST4.2005010112.06h).
> First, naming the NetCDF output of pcp_combine with
> a ".grb" extension was a poor choice. The MET tools assume that
files
> ending in .grb are GRIB and try to process them as such. There's a
way to
> override that behavior in METv4.1 but not in METv3.1.
> I'd suggest using a ".nc" extension.
>
> Also, the way the field information is specified differs depending
on the
> file type. For GRIB, we use the GRIB code abbreviation followed by
the
> accumulation interval to be used. For NetCDF, we use
> the NetCDF variable name followed by the dimensions to be used. So
please
> try using the following:
>
> fcst_field = "APCP_06(*,*)"; # for NetCDF
> obs_field = "APCP/A6"; # for GRIB
>
> Please give that a shot and let me know how it goes.
>
> Also, I was a bit confused about the GRIB metadata you sent. That
looks
> like NetCDF metadata, but I'm assuming you converted the GRIB file
to
> NetCDF to see what it contains. I'd suggest using the
> "wgrib" utility to inventory the contents of a GRIB file directly.
>
> Hope that helps.
>
> Thanks,
> John Halley Gotway
> met_help at ucar.edu
>
>
>
> On 07/09/2013 12:52 PM, Zachary Ebenstein via RT wrote:
> >
> > Tue Jul 09 12:52:36 2013: Request 62156 was acted upon.
> > Transaction: Ticket created by ebenszac at umich.edu
> > Queue: met_help
> > Subject: Issue using Mode tool
> > Owner: Nobody
> > Requestors: ebenszac at umich.edu
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62156 >
> >
> >
> > Hello,
> >
> > I am attempting to use the Mode tool to analyze model performance
> regarding
> > precipitation forecast. After using the pcp_combine tool to get
the
> > appropriate accumulation variable (6hr), I tried to run the Mode
tool as
> > below:
> >
> > wrf at srbwks35369-02-> bin/mode \
> > /data/wrf/wrf_mode/pcp_comb_out/wrf_pcp_d02_2005_01_01_12_6hr.grb
\
> > /data/wrf/precip_obs/ST4.2005010112.06h \
> > 00_for_frank/mode_config \
> > -outdir /data/wrf/wrf_mode/ \
> > -v 2
> >
> > The first input is the file that was created using pcp_combine,
and the
> > second input is the Stage IV observed precipitation data from the
EOL
> data
> > archive. This results in the following error:
> >
> > ERROR: GCInfo::set_gcinfo_grib() -> each GRIB code specified must
be
> > followed by an accumulation, level, or presssure level indicator
> "APCP_A6".
> >
> > The config file is setup with the comparison variables written as
such:
> >
> > fcst_field = "APCP_A6";
> > obs_field = "APCP_A6";
> >
> > I'm pretty sure I have those correct, but the mode tool won't
accept
> that.
> > To determine if the relevant data was present, I imported each of
the
> files
> > into IDV and was able to plot the accumulations for each, but I
did
> notice
> > in looking at the details of each file in IDV that the
> pcp_combine-created
> > file's field is named "APCP" with the level "A6", while the Stage
IV obs
> > file's field is also "APCP," but with no designated level in the
metadata
> > that I can find. Metadata for both files is attached.
> >
> > Thank you for your help,
> >
>
>
--
Zachary Ebenstein
University of Michigan
College of Engineering
Dept. Atmospheric, Oceanic and Space Sciences
Climate Physics (B.S.E), Dec. 2011
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #62156] Issue using Mode tool
From: John Halley Gotway
Time: Fri Jul 12 11:38:18 2013
Zac,
OK, probably the easiest thing to do now is for you to send me a
sample forecast file, observation file, and MODE configuration file.
Also send me the exact command line you're using to run MODE and
then cut-and-paste the output you get printed to the screen when you
run it.
Rather than sending those large file via email, just post them to our
anonymous ftp site following these instructions:
http://www.dtcenter.org/met/users/support/met_help.php#ftp
I'll run them through METv3.0 here and try to replicate the behavior
you're seeing. Then I can debug and figure out what's going on.
One question - have you updated to the latest set of METv3.0 patches
or are you using METv3.0 as originally released? Here's the known
issues page, including bug fixes, for that release:
http://www.dtcenter.org/met/users/support/known_issues/METv3.0/index.php
Please let me know when you've posted that data so I can go grab it
from ftp.
Thanks,
John
On 07/12/2013 11:23 AM, Zachary Ebenstein via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62156 >
>
> John,
>
> I am using METv3, and unfortunately I am not able to switch to METv4
at
> this time. As for the naming of the pcp_combine output, I had
mistakenly
> thought the output was still in GRIB format. I have produced new
output
> with the .nc file extension, as well as with no extension, with the
config
> changed to include these fields:
>
> fcst_field = "APCP_A6(*,*)";
> obs_field = "APCP_A6";
>
>
> Unfortunately, I have still had the same error as a result. The
metadata
> for the observation file looks like this:
>
> wrf at srbwks35369-02-> wgrib /data/wrf/precip_obs/ST4.2005010112.06h
>
1:0:d=04123112:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=18:P2=24:TimeU=1:sfc:18-
24hr
> acc:NAve=0
>
> Zac
>
>
>
> On Tue, Jul 9, 2013 at 5:19 PM, John Halley Gotway via RT
<met_help at ucar.edu
>> wrote:
>
>> Zachary,
>>
>> I see that you're having problems getting the configuration correct
for
>> the MODE tool. And since you're setting "fcst_field" and
"obs_field", I
>> assume you're using METv3.1. If that's not correct,
>> please let me know.
>>
>> First, let me tell you that there were two major issues discovered
in MODE
>> that have been corrected in METv4.1. They are described in the
"Existing
>> Tools" section of the METv4.1 release notes:
>>
>>
http://www.dtcenter.org/met/users/support/release_notes/METv4.1_release_notes.php
>>
>> So if you're able to switch to the newer version, I'd suggest doing
so.
>> One warning though - the configuration files have been updated a
lot since
>> METv3.1.
>>
>> Regarding the problem at hand, you're using MODE to compare a
NetCDF file
>> (wrf_pcp_d02_2005_01_01_12_6hr.grb) to a GRIB file
(ST4.2005010112.06h).
>> First, naming the NetCDF output of pcp_combine with
>> a ".grb" extension was a poor choice. The MET tools assume that
files
>> ending in .grb are GRIB and try to process them as such. There's a
way to
>> override that behavior in METv4.1 but not in METv3.1.
>> I'd suggest using a ".nc" extension.
>>
>> Also, the way the field information is specified differs depending
on the
>> file type. For GRIB, we use the GRIB code abbreviation followed by
the
>> accumulation interval to be used. For NetCDF, we use
>> the NetCDF variable name followed by the dimensions to be used. So
please
>> try using the following:
>>
>> fcst_field = "APCP_06(*,*)"; # for NetCDF
>> obs_field = "APCP/A6"; # for GRIB
>>
>> Please give that a shot and let me know how it goes.
>>
>> Also, I was a bit confused about the GRIB metadata you sent. That
looks
>> like NetCDF metadata, but I'm assuming you converted the GRIB file
to
>> NetCDF to see what it contains. I'd suggest using the
>> "wgrib" utility to inventory the contents of a GRIB file directly.
>>
>> Hope that helps.
>>
>> Thanks,
>> John Halley Gotway
>> met_help at ucar.edu
>>
>>
>>
>> On 07/09/2013 12:52 PM, Zachary Ebenstein via RT wrote:
>>>
>>> Tue Jul 09 12:52:36 2013: Request 62156 was acted upon.
>>> Transaction: Ticket created by ebenszac at umich.edu
>>> Queue: met_help
>>> Subject: Issue using Mode tool
>>> Owner: Nobody
>>> Requestors: ebenszac at umich.edu
>>> Status: new
>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62156 >
>>>
>>>
>>> Hello,
>>>
>>> I am attempting to use the Mode tool to analyze model performance
>> regarding
>>> precipitation forecast. After using the pcp_combine tool to get
the
>>> appropriate accumulation variable (6hr), I tried to run the Mode
tool as
>>> below:
>>>
>>> wrf at srbwks35369-02-> bin/mode \
>>> /data/wrf/wrf_mode/pcp_comb_out/wrf_pcp_d02_2005_01_01_12_6hr.grb
\
>>> /data/wrf/precip_obs/ST4.2005010112.06h \
>>> 00_for_frank/mode_config \
>>> -outdir /data/wrf/wrf_mode/ \
>>> -v 2
>>>
>>> The first input is the file that was created using pcp_combine,
and the
>>> second input is the Stage IV observed precipitation data from the
EOL
>> data
>>> archive. This results in the following error:
>>>
>>> ERROR: GCInfo::set_gcinfo_grib() -> each GRIB code specified must
be
>>> followed by an accumulation, level, or presssure level indicator
>> "APCP_A6".
>>>
>>> The config file is setup with the comparison variables written as
such:
>>>
>>> fcst_field = "APCP_A6";
>>> obs_field = "APCP_A6";
>>>
>>> I'm pretty sure I have those correct, but the mode tool won't
accept
>> that.
>>> To determine if the relevant data was present, I imported each of
the
>> files
>>> into IDV and was able to plot the accumulations for each, but I
did
>> notice
>>> in looking at the details of each file in IDV that the
>> pcp_combine-created
>>> file's field is named "APCP" with the level "A6", while the Stage
IV obs
>>> file's field is also "APCP," but with no designated level in the
metadata
>>> that I can find. Metadata for both files is attached.
>>>
>>> Thank you for your help,
>>>
>>
>>
>
>
------------------------------------------------
Subject: Issue using Mode tool
From: Zachary Ebenstein
Time: Fri Jul 12 12:45:55 2013
John,
I have posted the data to the files to the ftp site as follows:
forecast file: wrf_pcp_d02_2005_01_01_12_6hr.nc
obs file: ST4.2005010112.06h
config file: mode_config
The following is the copied and pasted command line and output:
wrf at srbwks35369-02-> bin/mode \
? /data/wrf/wrf_mode/pcp_comb_out/wrf_pcp_d02_2005_01_01_12_6hr.nc \
? /data/wrf/precip_obs/ST4.2005010112.06h \
? 00_for_frank/mode_config \
? outdir /data/wrf/wrf_mode/ \
? -v 2
ERROR: GCInfo::set_gcinfo_grib() -> each GRIB code specified must be
followed by an accumulation, level, or presssure level indicator
"APCP_A6".
Also, as far as I can tell, I am running the original METv3.0. I was
not
the person that installed the program, but I believe it was installed
shortly after v3.0 was released and hasn't been updated since. I know
it
has been used to successfully produce the desired output in the past,
using
input files that were produced in the same way (fcst file) or
retrieved
from the same source (obs file), though again I was not the one that
did it.
Zac
On Fri, Jul 12, 2013 at 1:38 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Zac,
>
> OK, probably the easiest thing to do now is for you to send me a
sample
> forecast file, observation file, and MODE configuration file. Also
send me
> the exact command line you're using to run MODE and
> then cut-and-paste the output you get printed to the screen when you
run
> it.
>
> Rather than sending those large file via email, just post them to
our
> anonymous ftp site following these instructions:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> I'll run them through METv3.0 here and try to replicate the behavior
> you're seeing. Then I can debug and figure out what's going on.
>
> One question - have you updated to the latest set of METv3.0 patches
or
> are you using METv3.0 as originally released? Here's the known
issues
> page, including bug fixes, for that release:
>
>
http://www.dtcenter.org/met/users/support/known_issues/METv3.0/index.php
>
> Please let me know when you've posted that data so I can go grab it
from
> ftp.
>
> Thanks,
> John
>
> On 07/12/2013 11:23 AM, Zachary Ebenstein via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62156 >
> >
> > John,
> >
> > I am using METv3, and unfortunately I am not able to switch to
METv4 at
> > this time. As for the naming of the pcp_combine output, I had
mistakenly
> > thought the output was still in GRIB format. I have produced new
output
> > with the .nc file extension, as well as with no extension, with
the
> config
> > changed to include these fields:
> >
> > fcst_field = "APCP_A6(*,*)";
> > obs_field = "APCP_A6";
> >
> >
> > Unfortunately, I have still had the same error as a result. The
metadata
> > for the observation file looks like this:
> >
> > wrf at srbwks35369-02-> wgrib /data/wrf/precip_obs/ST4.2005010112.06h
> >
>
1:0:d=04123112:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=18:P2=24:TimeU=1:sfc:18-
24hr
> > acc:NAve=0
> >
> > Zac
> >
> >
> >
> > On Tue, Jul 9, 2013 at 5:19 PM, John Halley Gotway via RT <
> met_help at ucar.edu
> >> wrote:
> >
> >> Zachary,
> >>
> >> I see that you're having problems getting the configuration
correct for
> >> the MODE tool. And since you're setting "fcst_field" and
"obs_field", I
> >> assume you're using METv3.1. If that's not correct,
> >> please let me know.
> >>
> >> First, let me tell you that there were two major issues
discovered in
> MODE
> >> that have been corrected in METv4.1. They are described in the
> "Existing
> >> Tools" section of the METv4.1 release notes:
> >>
> >>
>
http://www.dtcenter.org/met/users/support/release_notes/METv4.1_release_notes.php
> >>
> >> So if you're able to switch to the newer version, I'd suggest
doing so.
> >> One warning though - the configuration files have been updated
a lot
> since
> >> METv3.1.
> >>
> >> Regarding the problem at hand, you're using MODE to compare a
NetCDF
> file
> >> (wrf_pcp_d02_2005_01_01_12_6hr.grb) to a GRIB file
(ST4.2005010112.06h).
> >> First, naming the NetCDF output of pcp_combine with
> >> a ".grb" extension was a poor choice. The MET tools assume that
files
> >> ending in .grb are GRIB and try to process them as such. There's
a way
> to
> >> override that behavior in METv4.1 but not in METv3.1.
> >> I'd suggest using a ".nc" extension.
> >>
> >> Also, the way the field information is specified differs
depending on
> the
> >> file type. For GRIB, we use the GRIB code abbreviation followed
by the
> >> accumulation interval to be used. For NetCDF, we use
> >> the NetCDF variable name followed by the dimensions to be used.
So
> please
> >> try using the following:
> >>
> >> fcst_field = "APCP_06(*,*)"; # for NetCDF
> >> obs_field = "APCP/A6"; # for GRIB
> >>
> >> Please give that a shot and let me know how it goes.
> >>
> >> Also, I was a bit confused about the GRIB metadata you sent.
That looks
> >> like NetCDF metadata, but I'm assuming you converted the GRIB
file to
> >> NetCDF to see what it contains. I'd suggest using the
> >> "wgrib" utility to inventory the contents of a GRIB file
directly.
> >>
> >> Hope that helps.
> >>
> >> Thanks,
> >> John Halley Gotway
> >> met_help at ucar.edu
> >>
> >>
> >>
> >> On 07/09/2013 12:52 PM, Zachary Ebenstein via RT wrote:
> >>>
> >>> Tue Jul 09 12:52:36 2013: Request 62156 was acted upon.
> >>> Transaction: Ticket created by ebenszac at umich.edu
> >>> Queue: met_help
> >>> Subject: Issue using Mode tool
> >>> Owner: Nobody
> >>> Requestors: ebenszac at umich.edu
> >>> Status: new
> >>> Ticket <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62156 >
> >>>
> >>>
> >>> Hello,
> >>>
> >>> I am attempting to use the Mode tool to analyze model
performance
> >> regarding
> >>> precipitation forecast. After using the pcp_combine tool to get
the
> >>> appropriate accumulation variable (6hr), I tried to run the Mode
tool
> as
> >>> below:
> >>>
> >>> wrf at srbwks35369-02-> bin/mode \
> >>>
/data/wrf/wrf_mode/pcp_comb_out/wrf_pcp_d02_2005_01_01_12_6hr.grb \
> >>> /data/wrf/precip_obs/ST4.2005010112.06h \
> >>> 00_for_frank/mode_config \
> >>> -outdir /data/wrf/wrf_mode/ \
> >>> -v 2
> >>>
> >>> The first input is the file that was created using pcp_combine,
and the
> >>> second input is the Stage IV observed precipitation data from
the EOL
> >> data
> >>> archive. This results in the following error:
> >>>
> >>> ERROR: GCInfo::set_gcinfo_grib() -> each GRIB code specified
must be
> >>> followed by an accumulation, level, or presssure level indicator
> >> "APCP_A6".
> >>>
> >>> The config file is setup with the comparison variables written
as such:
> >>>
> >>> fcst_field = "APCP_A6";
> >>> obs_field = "APCP_A6";
> >>>
> >>> I'm pretty sure I have those correct, but the mode tool won't
accept
> >> that.
> >>> To determine if the relevant data was present, I imported each
of the
> >> files
> >>> into IDV and was able to plot the accumulations for each, but I
did
> >> notice
> >>> in looking at the details of each file in IDV that the
> >> pcp_combine-created
> >>> file's field is named "APCP" with the level "A6", while the
Stage IV
> obs
> >>> file's field is also "APCP," but with no designated level in the
> metadata
> >>> that I can find. Metadata for both files is attached.
> >>>
> >>> Thank you for your help,
> >>>
> >>
> >>
> >
> >
>
>
--
Zachary Ebenstein
University of Michigan
College of Engineering
Dept. Atmospheric, Oceanic and Space Sciences
Climate Physics (B.S.E), Dec. 2011
------------------------------------------------
Subject: Issue using Mode tool
From: John Halley Gotway
Time: Fri Jul 12 15:54:36 2013
Zach,
Using the following settings in the config file will get you past that
error:
fcst_field = "APCP_06(*,*)";
obs_field = "APCP/A6";
The fcst_field selects data from a NetCDF file while the obs_field
selects a file from a GRIB file.
However, once you're past that error, you'll get this one:
ERROR: process_fcst_obs_files() -> The forecast and observation
grids do not match.
You'll need to put the forecast and observation files onto a common
grid. You can easily regrid GRIB files using the copygb utility. So
we want to regrid the StageIV observation data to your model's
domain. Here's the current grid information:
- Forecast file:
ncdump -h wrf_pcp_d02_2005_01_01_12_6hr.nc
:Projection = "Lambert Conformal" ;
:scale_lat_1 = "50.000000" ;
:scale_lat_2 = "40.000000" ;
:lat_pin = "34.147000" ;
:lon_pin = "-99.751000" ;
:x_pin = "0.000000" ;
:y_pin = "0.000000" ;
:lon_orient = "-100.008000" ;
:d_km = "12.000000" ;
:r_km = "6367.470000" ;
:nx = "234" ;
:ny = "189 grid_points" ;
- Observation file:
wgrib -d 1 -V ST4.2005010112.06h
polar stereo: Lat1 23.117000 Long1 -119.023000 Orient -105.000000
north pole (1121 x 881) Dx 4763 Dy 4763 scan 64 mode 8
Here's how to regrid the StageIV file to the model grid. This is
described in the MET online tutorial:
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.0/copygb/run3.php
copygb -i3 -xg"255 3 234 189 34147 -99751 8 -100008 12000 12000 0
64 50000 40000" ST4.2005010112.06h ST4.2005010112.06h_regrid
Since this is precip, I used the -i3 option to select budget
interpolation. Then you can use the output of copygb as the
observation file for MODE. I've attached the first page of the
resulting
postscript output from MODE.
Hope that helps.
Thanks,
John
On 07/12/2013 12:45 PM, Zachary Ebenstein via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62156 >
>
> John,
>
> I have posted the data to the files to the ftp site as follows:
>
> forecast file: wrf_pcp_d02_2005_01_01_12_6hr.nc
> obs file: ST4.2005010112.06h
> config file: mode_config
>
> The following is the copied and pasted command line and output:
>
> wrf at srbwks35369-02-> bin/mode \
> ? /data/wrf/wrf_mode/pcp_comb_out/wrf_pcp_d02_2005_01_01_12_6hr.nc \
> ? /data/wrf/precip_obs/ST4.2005010112.06h \
> ? 00_for_frank/mode_config \
> ? outdir /data/wrf/wrf_mode/ \
> ? -v 2
>
>
> ERROR: GCInfo::set_gcinfo_grib() -> each GRIB code specified must be
> followed by an accumulation, level, or presssure level indicator
"APCP_A6".
>
> Also, as far as I can tell, I am running the original METv3.0. I was
not
> the person that installed the program, but I believe it was
installed
> shortly after v3.0 was released and hasn't been updated since. I
know it
> has been used to successfully produce the desired output in the
past, using
> input files that were produced in the same way (fcst file) or
retrieved
> from the same source (obs file), though again I was not the one that
did it.
>
> Zac
>
>
>
> On Fri, Jul 12, 2013 at 1:38 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
>> Zac,
>>
>> OK, probably the easiest thing to do now is for you to send me a
sample
>> forecast file, observation file, and MODE configuration file. Also
send me
>> the exact command line you're using to run MODE and
>> then cut-and-paste the output you get printed to the screen when
you run
>> it.
>>
>> Rather than sending those large file via email, just post them to
our
>> anonymous ftp site following these instructions:
>> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>>
>> I'll run them through METv3.0 here and try to replicate the
behavior
>> you're seeing. Then I can debug and figure out what's going on.
>>
>> One question - have you updated to the latest set of METv3.0
patches or
>> are you using METv3.0 as originally released? Here's the known
issues
>> page, including bug fixes, for that release:
>>
>>
http://www.dtcenter.org/met/users/support/known_issues/METv3.0/index.php
>>
>> Please let me know when you've posted that data so I can go grab it
from
>> ftp.
>>
>> Thanks,
>> John
>>
>> On 07/12/2013 11:23 AM, Zachary Ebenstein via RT wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62156 >
>>>
>>> John,
>>>
>>> I am using METv3, and unfortunately I am not able to switch to
METv4 at
>>> this time. As for the naming of the pcp_combine output, I had
mistakenly
>>> thought the output was still in GRIB format. I have produced new
output
>>> with the .nc file extension, as well as with no extension, with
the
>> config
>>> changed to include these fields:
>>>
>>> fcst_field = "APCP_A6(*,*)";
>>> obs_field = "APCP_A6";
>>>
>>>
>>> Unfortunately, I have still had the same error as a result. The
metadata
>>> for the observation file looks like this:
>>>
>>> wrf at srbwks35369-02-> wgrib /data/wrf/precip_obs/ST4.2005010112.06h
>>>
>>
1:0:d=04123112:APCP:kpds5=61:kpds6=1:kpds7=0:TR=4:P1=18:P2=24:TimeU=1:sfc:18-
24hr
>>> acc:NAve=0
>>>
>>> Zac
>>>
>>>
>>>
>>> On Tue, Jul 9, 2013 at 5:19 PM, John Halley Gotway via RT <
>> met_help at ucar.edu
>>>> wrote:
>>>
>>>> Zachary,
>>>>
>>>> I see that you're having problems getting the configuration
correct for
>>>> the MODE tool. And since you're setting "fcst_field" and
"obs_field", I
>>>> assume you're using METv3.1. If that's not correct,
>>>> please let me know.
>>>>
>>>> First, let me tell you that there were two major issues
discovered in
>> MODE
>>>> that have been corrected in METv4.1. They are described in the
>> "Existing
>>>> Tools" section of the METv4.1 release notes:
>>>>
>>>>
>>
http://www.dtcenter.org/met/users/support/release_notes/METv4.1_release_notes.php
>>>>
>>>> So if you're able to switch to the newer version, I'd suggest
doing so.
>>>> One warning though - the configuration files have been updated
a lot
>> since
>>>> METv3.1.
>>>>
>>>> Regarding the problem at hand, you're using MODE to compare a
NetCDF
>> file
>>>> (wrf_pcp_d02_2005_01_01_12_6hr.grb) to a GRIB file
(ST4.2005010112.06h).
>>>> First, naming the NetCDF output of pcp_combine with
>>>> a ".grb" extension was a poor choice. The MET tools assume that
files
>>>> ending in .grb are GRIB and try to process them as such. There's
a way
>> to
>>>> override that behavior in METv4.1 but not in METv3.1.
>>>> I'd suggest using a ".nc" extension.
>>>>
>>>> Also, the way the field information is specified differs
depending on
>> the
>>>> file type. For GRIB, we use the GRIB code abbreviation followed
by the
>>>> accumulation interval to be used. For NetCDF, we use
>>>> the NetCDF variable name followed by the dimensions to be used.
So
>> please
>>>> try using the following:
>>>>
>>>> fcst_field = "APCP_06(*,*)"; # for NetCDF
>>>> obs_field = "APCP/A6"; # for GRIB
>>>>
>>>> Please give that a shot and let me know how it goes.
>>>>
>>>> Also, I was a bit confused about the GRIB metadata you sent.
That looks
>>>> like NetCDF metadata, but I'm assuming you converted the GRIB
file to
>>>> NetCDF to see what it contains. I'd suggest using the
>>>> "wgrib" utility to inventory the contents of a GRIB file
directly.
>>>>
>>>> Hope that helps.
>>>>
>>>> Thanks,
>>>> John Halley Gotway
>>>> met_help at ucar.edu
>>>>
>>>>
>>>>
>>>> On 07/09/2013 12:52 PM, Zachary Ebenstein via RT wrote:
>>>>>
>>>>> Tue Jul 09 12:52:36 2013: Request 62156 was acted upon.
>>>>> Transaction: Ticket created by ebenszac at umich.edu
>>>>> Queue: met_help
>>>>> Subject: Issue using Mode tool
>>>>> Owner: Nobody
>>>>> Requestors: ebenszac at umich.edu
>>>>> Status: new
>>>>> Ticket <URL:
>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62156 >
>>>>>
>>>>>
>>>>> Hello,
>>>>>
>>>>> I am attempting to use the Mode tool to analyze model
performance
>>>> regarding
>>>>> precipitation forecast. After using the pcp_combine tool to get
the
>>>>> appropriate accumulation variable (6hr), I tried to run the Mode
tool
>> as
>>>>> below:
>>>>>
>>>>> wrf at srbwks35369-02-> bin/mode \
>>>>>
/data/wrf/wrf_mode/pcp_comb_out/wrf_pcp_d02_2005_01_01_12_6hr.grb \
>>>>> /data/wrf/precip_obs/ST4.2005010112.06h \
>>>>> 00_for_frank/mode_config \
>>>>> -outdir /data/wrf/wrf_mode/ \
>>>>> -v 2
>>>>>
>>>>> The first input is the file that was created using pcp_combine,
and the
>>>>> second input is the Stage IV observed precipitation data from
the EOL
>>>> data
>>>>> archive. This results in the following error:
>>>>>
>>>>> ERROR: GCInfo::set_gcinfo_grib() -> each GRIB code specified
must be
>>>>> followed by an accumulation, level, or presssure level indicator
>>>> "APCP_A6".
>>>>>
>>>>> The config file is setup with the comparison variables written
as such:
>>>>>
>>>>> fcst_field = "APCP_A6";
>>>>> obs_field = "APCP_A6";
>>>>>
>>>>> I'm pretty sure I have those correct, but the mode tool won't
accept
>>>> that.
>>>>> To determine if the relevant data was present, I imported each
of the
>>>> files
>>>>> into IDV and was able to plot the accumulations for each, but I
did
>>>> notice
>>>>> in looking at the details of each file in IDV that the
>>>> pcp_combine-created
>>>>> file's field is named "APCP" with the level "A6", while the
Stage IV
>> obs
>>>>> file's field is also "APCP," but with no designated level in the
>> metadata
>>>>> that I can find. Metadata for both files is attached.
>>>>>
>>>>> Thank you for your help,
>>>>>
>>>>
>>>>
>>>
>>>
>>
>>
>
>
------------------------------------------------
More information about the Met_help
mailing list