[Met_help] [rt.rap.ucar.edu #58307] History for netcdf config file
John Halley Gotway via RT
met_help at ucar.edu
Mon Oct 8 08:47:38 MDT 2012
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
hi, I try to run point-stat tool whit fcst_file in netcdf format but no run
and i think that the problem is in my config file can you show me an
example for that file?
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #58307] netcdf config file
From: John Halley Gotway
Time: Mon Sep 17 07:12:55 2012
Jose,
The Point-Stat configuration file is described in section 4.3.2 of the
MET User's Guide:
www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v4.0.1.pdf
Sample Point-Stat configuration files can also be found in the
released code (I'm assuming METv4.0, here):
METv4.0/data/config/PointStatConfig_default
METv4.0/scripts/config/PointStatConfig
And this README file describes the config file entries - search for
"NetCDF" in that file:
METv4.0/data/config/README
You may also find it helpful to look at the online tutorial:
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.0/point_stat/index.php
However, I see that none of those examples show using an input NetCDF
file. But here's an example of that for Grid-Stat:
METv4.0/scripts/config/GridStatConfig_APCP_12
It's pretty simple, for gridded NetCDF files, the "name" entry is set
to the NetCDF variable name. And the "level" specifies how to extract
the 2D gridded field from that variable. See that README
file. One word of warning, if you're verifying a NetCDF forecast file
in Point-Stat, you'll need to define the "obs" dictionary to specify
which observations should be used for the verification. In
Point-Stat, the observations are always specified using GRIB
conventions.
Give that a shot, and feel free to send us your config file and input
data for Point-Stat if you continue to experience problems.
Hope that helps.
John Halley Gotway
met_help at ucar.edu
On 09/14/2012 04:02 PM, Josue Don via RT wrote:
>
> Fri Sep 14 16:02:49 2012: Request 58307 was acted upon.
> Transaction: Ticket created by eljoz89 at gmail.com
> Queue: met_help
> Subject: netcdf config file
> Owner: Nobody
> Requestors: eljoz89 at gmail.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
>
>
> hi, I try to run point-stat tool whit fcst_file in netcdf format but
no run
> and i think that the problem is in my config file can you show me an
> example for that file?
>
------------------------------------------------
Subject: netcdf config file
From: Josue Don
Time: Tue Sep 18 09:09:36 2012
thanks for answer I try all your advices but I can't do that poin-stat
tool
works I'm using for inputs:
fcst_file: wrf out file in netcdf format in levels zigma and
pressure.(i
try to use pcp_combine but don't works)
obs_file: http://rda.ucar.edu/datasets/ds337.0/#description netcdf
suffice
and http://rda.ucar.edu/datasets/ds336.0/#description in netcdf metar
(i
try to use madis2nc tool but don't works)
config_file: attach it.
the error message says:
DEBUG 1: Default Config File:
/home/xozue/MET/METv4.0/data/config/PointStatConfig_default
DEBUG 1: User Config File: proyecto/config/PointStatConfig_tutorial
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
Met2dDataFile object of type "FileType_NcPinterp".
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of
type "FileType_NcPinterp".
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of
type "FileType_NcPinterp".
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of
type "FileType_NcPinterp".
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of
type "FileType_NcPinterp".
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of
type "FileType_NcPinterp".
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=2307545168
DEBUG 1: Forecast File: ../datos/wrfout_d01_2012-07-03_PLEVELS
DEBUG 1: Climatology File: none
DEBUG 1: Observation File: ../datos/Surface_METAR_20120703_0000.nc
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for SPFHP500.
ERROR :
ERROR : process_fcst_climo_files() -> no fields matching SPFHP500
found in
file: ../datos/wrfout_d01_2012-07-03_PLEVELS
ERROR :
i never call this SPFHP500 var.
i think that the problem is in the config file and i don't understand
the
format for obs dictionary and fcst for netcdf can you expand the
explication whit a example like tutorial page for point-stat tool.
thanks for your time.
On 17 September 2012 08:12, John Halley Gotway via RT
<met_help at ucar.edu>wrote:
> Jose,
>
> The Point-Stat configuration file is described in section 4.3.2 of
the MET
> User's Guide:
>
www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v4.0.1.pdf
>
> Sample Point-Stat configuration files can also be found in the
released
> code (I'm assuming METv4.0, here):
> METv4.0/data/config/PointStatConfig_default
> METv4.0/scripts/config/PointStatConfig
>
> And this README file describes the config file entries - search for
> "NetCDF" in that file:
> METv4.0/data/config/README
>
> You may also find it helpful to look at the online tutorial:
>
>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.0/point_stat/index.php
>
> However, I see that none of those examples show using an input
NetCDF
> file. But here's an example of that for Grid-Stat:
> METv4.0/scripts/config/GridStatConfig_APCP_12
>
> It's pretty simple, for gridded NetCDF files, the "name" entry is
set to
> the NetCDF variable name. And the "level" specifies how to extract
the 2D
> gridded field from that variable. See that README
> file. One word of warning, if you're verifying a NetCDF forecast
file in
> Point-Stat, you'll need to define the "obs" dictionary to specify
which
> observations should be used for the verification. In
> Point-Stat, the observations are always specified using GRIB
conventions.
>
> Give that a shot, and feel free to send us your config file and
input data
> for Point-Stat if you continue to experience problems.
>
> Hope that helps.
>
> John Halley Gotway
> met_help at ucar.edu
>
> On 09/14/2012 04:02 PM, Josue Don via RT wrote:
> >
> > Fri Sep 14 16:02:49 2012: Request 58307 was acted upon.
> > Transaction: Ticket created by eljoz89 at gmail.com
> > Queue: met_help
> > Subject: netcdf config file
> > Owner: Nobody
> > Requestors: eljoz89 at gmail.com
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
> >
> >
> > hi, I try to run point-stat tool whit fcst_file in netcdf format
but no
> run
> > and i think that the problem is in my config file can you show me
an
> > example for that file?
> >
>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #58307] netcdf config file
From: John Halley Gotway
Time: Tue Sep 18 09:24:40 2012
Jose,
It would probably be easiest for you to just send me your data files,
and I'll run it here. I'll need:
- What version of MET are you running?
- Forecast file: datos/wrfout_d01_2012-07-03_PLEVELS
- Observation file: datos/Surface_METAR_20120703_0000.nc
- Config file: proyecto/config/PointStatConfig_tutorial
Please post them to our anonymous ftp site, following the instructions
listed here:
http://www.dtcenter.org/met/users/support/met_help.php#ftp
Thanks,
John
On 09/18/2012 09:09 AM, Josue Don via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
>
> thanks for answer I try all your advices but I can't do that poin-
stat tool
> works I'm using for inputs:
> fcst_file: wrf out file in netcdf format in levels zigma and
pressure.(i
> try to use pcp_combine but don't works)
> obs_file: http://rda.ucar.edu/datasets/ds337.0/#description netcdf
suffice
> and http://rda.ucar.edu/datasets/ds336.0/#description in netcdf
metar (i
> try to use madis2nc tool but don't works)
> config_file: attach it.
>
> the error message says:
> DEBUG 1: Default Config File:
> /home/xozue/MET/METv4.0/data/config/PointStatConfig_default
> DEBUG 1: User Config File: proyecto/config/PointStatConfig_tutorial
> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
> Met2dDataFile object of type "FileType_NcPinterp".
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> type "FileType_NcPinterp".
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> type "FileType_NcPinterp".
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> type "FileType_NcPinterp".
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> type "FileType_NcPinterp".
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> type "FileType_NcPinterp".
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=2307545168
> DEBUG 1: Forecast File: ../datos/wrfout_d01_2012-07-03_PLEVELS
> DEBUG 1: Climatology File: none
> DEBUG 1: Observation File: ../datos/Surface_METAR_20120703_0000.nc
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for SPFHP500.
> ERROR :
> ERROR : process_fcst_climo_files() -> no fields matching SPFHP500
found in
> file: ../datos/wrfout_d01_2012-07-03_PLEVELS
> ERROR :
>
> i never call this SPFHP500 var.
>
> i think that the problem is in the config file and i don't
understand the
> format for obs dictionary and fcst for netcdf can you expand the
> explication whit a example like tutorial page for point-stat tool.
>
> thanks for your time.
>
> On 17 September 2012 08:12, John Halley Gotway via RT
<met_help at ucar.edu>wrote:
>
>> Jose,
>>
>> The Point-Stat configuration file is described in section 4.3.2 of
the MET
>> User's Guide:
>>
www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v4.0.1.pdf
>>
>> Sample Point-Stat configuration files can also be found in the
released
>> code (I'm assuming METv4.0, here):
>> METv4.0/data/config/PointStatConfig_default
>> METv4.0/scripts/config/PointStatConfig
>>
>> And this README file describes the config file entries - search for
>> "NetCDF" in that file:
>> METv4.0/data/config/README
>>
>> You may also find it helpful to look at the online tutorial:
>>
>>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.0/point_stat/index.php
>>
>> However, I see that none of those examples show using an input
NetCDF
>> file. But here's an example of that for Grid-Stat:
>> METv4.0/scripts/config/GridStatConfig_APCP_12
>>
>> It's pretty simple, for gridded NetCDF files, the "name" entry is
set to
>> the NetCDF variable name. And the "level" specifies how to extract
the 2D
>> gridded field from that variable. See that README
>> file. One word of warning, if you're verifying a NetCDF forecast
file in
>> Point-Stat, you'll need to define the "obs" dictionary to specify
which
>> observations should be used for the verification. In
>> Point-Stat, the observations are always specified using GRIB
conventions.
>>
>> Give that a shot, and feel free to send us your config file and
input data
>> for Point-Stat if you continue to experience problems.
>>
>> Hope that helps.
>>
>> John Halley Gotway
>> met_help at ucar.edu
>>
>> On 09/14/2012 04:02 PM, Josue Don via RT wrote:
>>>
>>> Fri Sep 14 16:02:49 2012: Request 58307 was acted upon.
>>> Transaction: Ticket created by eljoz89 at gmail.com
>>> Queue: met_help
>>> Subject: netcdf config file
>>> Owner: Nobody
>>> Requestors: eljoz89 at gmail.com
>>> Status: new
>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
>>>
>>>
>>> hi, I try to run point-stat tool whit fcst_file in netcdf format
but no
>> run
>>> and i think that the problem is in my config file can you show me
an
>>> example for that file?
>>>
>>
>>
>>
------------------------------------------------
Subject: netcdf config file
From: Josue Don
Time: Tue Sep 18 23:28:48 2012
hi, i upload the files in don_data folder:
fsct: PLEVELS -> post proces the SLEVEL file whit p_interp tool (WRF
post
proces tools) (72 hours)
SLEVELS -> direct to WRF v.3.1 output (72 hours)
obs: ds336.0 -> Surface_METAR_20120703_0000.nc (24 hours)
Surface_Synoptic_20120703_0000.nc.bz2 (24 hours)
http://rda.ucar.edu/datasets/ds336.0/#description
ds337.0 -> netcdf.gdas.2012070300.nc (6 hours)
http://rda.ucar.edu/datasets/ds337.0/#description
I try whit all this files.
On 18 September 2012 10:24, John Halley Gotway via RT
<met_help at ucar.edu>wrote:
> Jose,
>
> It would probably be easiest for you to just send me your data
files, and
> I'll run it here. I'll need:
> - What version of MET are you running?
> - Forecast file: datos/wrfout_d01_2012-07-03_PLEVELS
> - Observation file: datos/Surface_METAR_20120703_0000.nc
> - Config file: proyecto/config/PointStatConfig_tutorial
>
> Please post them to our anonymous ftp site, following the
instructions
> listed here:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> Thanks,
> John
>
> On 09/18/2012 09:09 AM, Josue Don via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
> >
> > thanks for answer I try all your advices but I can't do that poin-
stat
> tool
> > works I'm using for inputs:
> > fcst_file: wrf out file in netcdf format in levels zigma and
pressure.(i
> > try to use pcp_combine but don't works)
> > obs_file: http://rda.ucar.edu/datasets/ds337.0/#description netcdf
> suffice
> > and http://rda.ucar.edu/datasets/ds336.0/#description in netcdf
metar (i
> > try to use madis2nc tool but don't works)
> > config_file: attach it.
> >
> > the error message says:
> > DEBUG 1: Default Config File:
> > /home/xozue/MET/METv4.0/data/config/PointStatConfig_default
> > DEBUG 1: User Config File:
proyecto/config/PointStatConfig_tutorial
> > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
> > Met2dDataFile object of type "FileType_NcPinterp".
> > DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> > type "FileType_NcPinterp".
> > DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> > type "FileType_NcPinterp".
> > DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> > type "FileType_NcPinterp".
> > DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> > type "FileType_NcPinterp".
> > DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> > type "FileType_NcPinterp".
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=2307545168
> > DEBUG 1: Forecast File: ../datos/wrfout_d01_2012-07-03_PLEVELS
> > DEBUG 1: Climatology File: none
> > DEBUG 1: Observation File: ../datos/Surface_METAR_20120703_0000.nc
> > DEBUG 2:
> > DEBUG 2:
> >
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Reading data for SPFHP500.
> > ERROR :
> > ERROR : process_fcst_climo_files() -> no fields matching SPFHP500
found
> in
> > file: ../datos/wrfout_d01_2012-07-03_PLEVELS
> > ERROR :
> >
> > i never call this SPFHP500 var.
> >
> > i think that the problem is in the config file and i don't
understand the
> > format for obs dictionary and fcst for netcdf can you expand the
> > explication whit a example like tutorial page for point-stat tool.
> >
> > thanks for your time.
> >
> > On 17 September 2012 08:12, John Halley Gotway via RT
<met_help at ucar.edu
> >wrote:
> >
> >> Jose,
> >>
> >> The Point-Stat configuration file is described in section 4.3.2
of the
> MET
> >> User's Guide:
> >>
>
www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v4.0.1.pdf
> >>
> >> Sample Point-Stat configuration files can also be found in the
released
> >> code (I'm assuming METv4.0, here):
> >> METv4.0/data/config/PointStatConfig_default
> >> METv4.0/scripts/config/PointStatConfig
> >>
> >> And this README file describes the config file entries - search
for
> >> "NetCDF" in that file:
> >> METv4.0/data/config/README
> >>
> >> You may also find it helpful to look at the online tutorial:
> >>
> >>
>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.0/point_stat/index.php
> >>
> >> However, I see that none of those examples show using an input
NetCDF
> >> file. But here's an example of that for Grid-Stat:
> >> METv4.0/scripts/config/GridStatConfig_APCP_12
> >>
> >> It's pretty simple, for gridded NetCDF files, the "name" entry is
set to
> >> the NetCDF variable name. And the "level" specifies how to
extract the
> 2D
> >> gridded field from that variable. See that README
> >> file. One word of warning, if you're verifying a NetCDF forecast
file
> in
> >> Point-Stat, you'll need to define the "obs" dictionary to specify
which
> >> observations should be used for the verification. In
> >> Point-Stat, the observations are always specified using GRIB
> conventions.
> >>
> >> Give that a shot, and feel free to send us your config file and
input
> data
> >> for Point-Stat if you continue to experience problems.
> >>
> >> Hope that helps.
> >>
> >> John Halley Gotway
> >> met_help at ucar.edu
> >>
> >> On 09/14/2012 04:02 PM, Josue Don via RT wrote:
> >>>
> >>> Fri Sep 14 16:02:49 2012: Request 58307 was acted upon.
> >>> Transaction: Ticket created by eljoz89 at gmail.com
> >>> Queue: met_help
> >>> Subject: netcdf config file
> >>> Owner: Nobody
> >>> Requestors: eljoz89 at gmail.com
> >>> Status: new
> >>> Ticket <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
> >>>
> >>>
> >>> hi, I try to run point-stat tool whit fcst_file in netcdf format
but no
> >> run
> >>> and i think that the problem is in my config file can you show
me an
> >>> example for that file?
> >>>
> >>
> >>
> >>
>
>
------------------------------------------------
Subject: netcdf config file
From: John Halley Gotway
Time: Wed Sep 19 08:43:04 2012
Jose,
I retrieved your data files and ran them with MET. I started just by
using the plot_data_plane tool to plot some data and make sure it
looks OK:
# Plot RAINC for time 11
METv4.0/bin/plot_data_plane wrfout_d01_2012-07-03_PLEVELS
rainc_time11.ps 'name="RAINC" ; level="(10,*,*)";' -v 10
# Plot SKINTEMP for time 11
METv4.0/bin/plot_data_plane wrfout_d01_2012-07-03_PLEVELS
skintemp_time11.ps 'name="SKINTEMP" ; level="(10,*,*)";' -v 10
The output images look fine and I see the your domain is primarily
over Mexico.
Next, by looking in your Point-Stat config file, I see that you want
to verify the precipitation in the "RAINC(0,*,*)" variable. There are
a couple problems with that:
- First, there are no observations of precipitation in the GDAS
point observation file (netcdf.gdas.2012070300.nc).
- Second, while there are observations of precipitation in the
METAR file (Surface_METAR_20120703_0000.nc), that file is not in a
format that MET can handle. However, I have made a note of this
data source and we will consider adding support for it in future
releases of MET. For now, to use those observations of precip, you'd
need to reformat them in a way that MET can handle.
To get you going, I switched over from verifying precip to verifying
surface temperature. Here I'm verifying against observations of
surface temperature from the GDAS file with message type ADPSFC
(means surface). You can plot the location of those observations on
your domain using the following command:
# Plot location of TMP for ADPSFC over the model domain
METv4.0/bin/plot_point_obs \
netcdf.gdas.2012070300.nc \
point_obs_tmp_adpsfc.ps \
-gc 11 -msg_typ ADPSFC \
-data_file wrfout_d01_2012-07-03_PLEVELS
# Next, run Point-Stat to compare the model output to the
observations:
METv4.0/bin/point_stat \
wrfout_d01_2012-07-03_PLEVELS \
netcdf.gdas.2012070300.nc \
PointStatConfig_TMP_ADPSFC \
-outdir out \
-v 2
I've attached the config file I used. Hope this helps getting you
going in the right direction.
John
On 09/18/2012 11:28 PM, Josue Don via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
>
> hi, i upload the files in don_data folder:
>
> fsct: PLEVELS -> post proces the SLEVEL file whit p_interp tool (WRF
post
> proces tools) (72 hours)
> SLEVELS -> direct to WRF v.3.1 output (72 hours)
>
> obs: ds336.0 -> Surface_METAR_20120703_0000.nc (24 hours)
> Surface_Synoptic_20120703_0000.nc.bz2 (24 hours)
> http://rda.ucar.edu/datasets/ds336.0/#description
>
> ds337.0 -> netcdf.gdas.2012070300.nc (6 hours)
> http://rda.ucar.edu/datasets/ds337.0/#description
>
> I try whit all this files.
>
>
> On 18 September 2012 10:24, John Halley Gotway via RT
<met_help at ucar.edu>wrote:
>
>> Jose,
>>
>> It would probably be easiest for you to just send me your data
files, and
>> I'll run it here. I'll need:
>> - What version of MET are you running?
>> - Forecast file: datos/wrfout_d01_2012-07-03_PLEVELS
>> - Observation file: datos/Surface_METAR_20120703_0000.nc
>> - Config file: proyecto/config/PointStatConfig_tutorial
>>
>> Please post them to our anonymous ftp site, following the
instructions
>> listed here:
>> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>>
>> Thanks,
>> John
>>
>> On 09/18/2012 09:09 AM, Josue Don via RT wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
>>>
>>> thanks for answer I try all your advices but I can't do that poin-
stat
>> tool
>>> works I'm using for inputs:
>>> fcst_file: wrf out file in netcdf format in levels zigma and
pressure.(i
>>> try to use pcp_combine but don't works)
>>> obs_file: http://rda.ucar.edu/datasets/ds337.0/#description netcdf
>> suffice
>>> and http://rda.ucar.edu/datasets/ds336.0/#description in netcdf
metar (i
>>> try to use madis2nc tool but don't works)
>>> config_file: attach it.
>>>
>>> the error message says:
>>> DEBUG 1: Default Config File:
>>> /home/xozue/MET/METv4.0/data/config/PointStatConfig_default
>>> DEBUG 1: User Config File:
proyecto/config/PointStatConfig_tutorial
>>> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
>>> Met2dDataFile object of type "FileType_NcPinterp".
>>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
>>> type "FileType_NcPinterp".
>>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
>>> type "FileType_NcPinterp".
>>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
>>> type "FileType_NcPinterp".
>>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
>>> type "FileType_NcPinterp".
>>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
>>> type "FileType_NcPinterp".
>>> GSL_RNG_TYPE=mt19937
>>> GSL_RNG_SEED=2307545168
>>> DEBUG 1: Forecast File: ../datos/wrfout_d01_2012-07-03_PLEVELS
>>> DEBUG 1: Climatology File: none
>>> DEBUG 1: Observation File: ../datos/Surface_METAR_20120703_0000.nc
>>> DEBUG 2:
>>> DEBUG 2:
>>>
>>
--------------------------------------------------------------------------------
>>> DEBUG 2:
>>> DEBUG 2: Reading data for SPFHP500.
>>> ERROR :
>>> ERROR : process_fcst_climo_files() -> no fields matching SPFHP500
found
>> in
>>> file: ../datos/wrfout_d01_2012-07-03_PLEVELS
>>> ERROR :
>>>
>>> i never call this SPFHP500 var.
>>>
>>> i think that the problem is in the config file and i don't
understand the
>>> format for obs dictionary and fcst for netcdf can you expand the
>>> explication whit a example like tutorial page for point-stat tool.
>>>
>>> thanks for your time.
>>>
>>> On 17 September 2012 08:12, John Halley Gotway via RT
<met_help at ucar.edu
>>> wrote:
>>>
>>>> Jose,
>>>>
>>>> The Point-Stat configuration file is described in section 4.3.2
of the
>> MET
>>>> User's Guide:
>>>>
>>
www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v4.0.1.pdf
>>>>
>>>> Sample Point-Stat configuration files can also be found in the
released
>>>> code (I'm assuming METv4.0, here):
>>>> METv4.0/data/config/PointStatConfig_default
>>>> METv4.0/scripts/config/PointStatConfig
>>>>
>>>> And this README file describes the config file entries - search
for
>>>> "NetCDF" in that file:
>>>> METv4.0/data/config/README
>>>>
>>>> You may also find it helpful to look at the online tutorial:
>>>>
>>>>
>>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.0/point_stat/index.php
>>>>
>>>> However, I see that none of those examples show using an input
NetCDF
>>>> file. But here's an example of that for Grid-Stat:
>>>> METv4.0/scripts/config/GridStatConfig_APCP_12
>>>>
>>>> It's pretty simple, for gridded NetCDF files, the "name" entry is
set to
>>>> the NetCDF variable name. And the "level" specifies how to
extract the
>> 2D
>>>> gridded field from that variable. See that README
>>>> file. One word of warning, if you're verifying a NetCDF forecast
file
>> in
>>>> Point-Stat, you'll need to define the "obs" dictionary to specify
which
>>>> observations should be used for the verification. In
>>>> Point-Stat, the observations are always specified using GRIB
>> conventions.
>>>>
>>>> Give that a shot, and feel free to send us your config file and
input
>> data
>>>> for Point-Stat if you continue to experience problems.
>>>>
>>>> Hope that helps.
>>>>
>>>> John Halley Gotway
>>>> met_help at ucar.edu
>>>>
>>>> On 09/14/2012 04:02 PM, Josue Don via RT wrote:
>>>>>
>>>>> Fri Sep 14 16:02:49 2012: Request 58307 was acted upon.
>>>>> Transaction: Ticket created by eljoz89 at gmail.com
>>>>> Queue: met_help
>>>>> Subject: netcdf config file
>>>>> Owner: Nobody
>>>>> Requestors: eljoz89 at gmail.com
>>>>> Status: new
>>>>> Ticket <URL:
>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
>>>>>
>>>>>
>>>>> hi, I try to run point-stat tool whit fcst_file in netcdf format
but no
>>>> run
>>>>> and i think that the problem is in my config file can you show
me an
>>>>> example for that file?
>>>>>
>>>>
>>>>
>>>>
>>
>>
------------------------------------------------
Subject: netcdf config file
From: John Halley Gotway
Time: Wed Sep 19 08:43:04 2012
////////////////////////////////////////////////////////////////////////////////
//
// Point-Stat configuration file.
//
// For additional information, see the MET_BASE/data/config/README
file.
//
////////////////////////////////////////////////////////////////////////////////
//
// Output model name to be written
//
model = "WRF";
////////////////////////////////////////////////////////////////////////////////
//
// Forecast and observation fields to be verified
//
fcst = {
wind_thresh = [ NA ];
field = [
{
name = "SKINTEMP";
level = [ "(0,*,*)" ];
cat_thresh = [ >273 ];
}
];
};
obs = {
message_type = "ADPSFC";
field = [
{
name = "TMP";
level = [ "Z0" ];
cat_thresh = [ >273 ];
}
];
};
////////////////////////////////////////////////////////////////////////////////
//
// Point observation time window
//
obs_window = {
beg = -5400;
end = 5400;
}
////////////////////////////////////////////////////////////////////////////////
//
// Verification masking regions
//
mask = {
grid = [ "FULL" ];
poly = [];
sid = "";
};
////////////////////////////////////////////////////////////////////////////////
//
// Confidence interval settings
//
ci_alpha = [ 0.05 ];
boot = {
interval = PCTILE;
rep_prop = 1.0;
n_rep = 1000;
rng = "mt19937";
seed = "";
};
////////////////////////////////////////////////////////////////////////////////
//
// Interpolation methods
//
interp = {
field = FCST;
vld_thresh = 1.0;
type = [
{
method = UW_MEAN;
width = 1;
}
];
};
////////////////////////////////////////////////////////////////////////////////
//
// Statistical output types
//
output_flag = {
fho = BOTH;
ctc = BOTH;
cts = BOTH;
mctc = BOTH;
mcts = BOTH;
cnt = BOTH;
sl1l2 = BOTH;
sal1l2 = BOTH;
vl1l2 = BOTH;
val1l2 = BOTH;
pct = BOTH;
pstd = BOTH;
pjc = BOTH;
prc = BOTH;
mpr = BOTH;
};
//
// NetCDF matched pairs output file
//
nc_pairs_flag = TRUE;
////////////////////////////////////////////////////////////////////////////////
//duplicate_flag = NONE;
rank_corr_flag = TRUE;
tmp_dir = "/tmp";
output_prefix = "";
version = "V4.0";
////////////////////////////////////////////////////////////////////////////////
------------------------------------------------
Subject: netcdf config file
From: Josue Don
Time: Wed Sep 19 14:05:05 2012
thanks, It was really helpful.
I have some questions:
- Is it the same if using wrfout_d01_2012-07-03_SLEVELS file Zigma
levels
of pressure? or
Do I need to process it with "p_interp"?
- How I can know the variables in a "netcdf" file and equivalent in
format
"Grid1" to generate the values of "obs" and "FCST"?
- Could you tell me where to get the table of equivalence between
variables
"netcdf" and "GRIB"?
- Is there a tool that allows me to change the file format "metar"
(in Surface_METAR_20120703_0000.nc) to a compatible format and what is
the
right format?
On 19 September 2012 09:43, John Halley Gotway via RT
<met_help at ucar.edu>wrote:
> Jose,
>
> I retrieved your data files and ran them with MET. I started just
by
> using the plot_data_plane tool to plot some data and make sure it
looks OK:
>
> # Plot RAINC for time 11
> METv4.0/bin/plot_data_plane wrfout_d01_2012-07-03_PLEVELS
rainc_time11.ps'name="RAINC" ; level="(10,*,*)";' -v 10
>
> # Plot SKINTEMP for time 11
> METv4.0/bin/plot_data_plane wrfout_d01_2012-07-03_PLEVELS
> skintemp_time11.ps 'name="SKINTEMP" ; level="(10,*,*)";' -v 10
>
> The output images look fine and I see the your domain is primarily
over
> Mexico.
>
> Next, by looking in your Point-Stat config file, I see that you want
to
> verify the precipitation in the "RAINC(0,*,*)" variable. There are
a
> couple problems with that:
> - First, there are no observations of precipitation in the GDAS
point
> observation file (netcdf.gdas.2012070300.nc).
> - Second, while there are observations of precipitation in the
METAR
> file (Surface_METAR_20120703_0000.nc), that file is not in a format
that
> MET can handle. However, I have made a note of this
> data source and we will consider adding support for it in future
releases
> of MET. For now, to use those observations of precip, you'd need to
> reformat them in a way that MET can handle.
>
> To get you going, I switched over from verifying precip to verifying
> surface temperature. Here I'm verifying against observations of
surface
> temperature from the GDAS file with message type ADPSFC
> (means surface). You can plot the location of those observations on
your
> domain using the following command:
>
> # Plot location of TMP for ADPSFC over the model domain
> METv4.0/bin/plot_point_obs \
> netcdf.gdas.2012070300.nc \
> point_obs_tmp_adpsfc.ps \
> -gc 11 -msg_typ ADPSFC \
> -data_file wrfout_d01_2012-07-03_PLEVELS
>
> # Next, run Point-Stat to compare the model output to the
observations:
> METv4.0/bin/point_stat \
> wrfout_d01_2012-07-03_PLEVELS \
> netcdf.gdas.2012070300.nc \
> PointStatConfig_TMP_ADPSFC \
> -outdir out \
> -v 2
>
> I've attached the config file I used. Hope this helps getting you
going
> in the right direction.
>
> John
>
> On 09/18/2012 11:28 PM, Josue Don via RT wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
> >
> > hi, i upload the files in don_data folder:
> >
> > fsct: PLEVELS -> post proces the SLEVEL file whit p_interp tool
(WRF post
> > proces tools) (72 hours)
> > SLEVELS -> direct to WRF v.3.1 output (72 hours)
> >
> > obs: ds336.0 -> Surface_METAR_20120703_0000.nc (24 hours)
> > Surface_Synoptic_20120703_0000.nc.bz2 (24 hours)
> > http://rda.ucar.edu/datasets/ds336.0/#description
> >
> > ds337.0 -> netcdf.gdas.2012070300.nc (6 hours)
> > http://rda.ucar.edu/datasets/ds337.0/#description
> >
> > I try whit all this files.
> >
> >
> > On 18 September 2012 10:24, John Halley Gotway via RT
<met_help at ucar.edu
> >wrote:
> >
> >> Jose,
> >>
> >> It would probably be easiest for you to just send me your data
files,
> and
> >> I'll run it here. I'll need:
> >> - What version of MET are you running?
> >> - Forecast file: datos/wrfout_d01_2012-07-03_PLEVELS
> >> - Observation file: datos/Surface_METAR_20120703_0000.nc
> >> - Config file: proyecto/config/PointStatConfig_tutorial
> >>
> >> Please post them to our anonymous ftp site, following the
instructions
> >> listed here:
> >> http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >>
> >> Thanks,
> >> John
> >>
> >> On 09/18/2012 09:09 AM, Josue Don via RT wrote:
> >>>
> >>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
> >>>
> >>> thanks for answer I try all your advices but I can't do that
poin-stat
> >> tool
> >>> works I'm using for inputs:
> >>> fcst_file: wrf out file in netcdf format in levels zigma and
> pressure.(i
> >>> try to use pcp_combine but don't works)
> >>> obs_file: http://rda.ucar.edu/datasets/ds337.0/#description
netcdf
> >> suffice
> >>> and http://rda.ucar.edu/datasets/ds336.0/#description in netcdf
metar
> (i
> >>> try to use madis2nc tool but don't works)
> >>> config_file: attach it.
> >>>
> >>> the error message says:
> >>> DEBUG 1: Default Config File:
> >>> /home/xozue/MET/METv4.0/data/config/PointStatConfig_default
> >>> DEBUG 1: User Config File:
proyecto/config/PointStatConfig_tutorial
> >>> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
> >>> Met2dDataFile object of type "FileType_NcPinterp".
> >>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> of
> >>> type "FileType_NcPinterp".
> >>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> of
> >>> type "FileType_NcPinterp".
> >>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> of
> >>> type "FileType_NcPinterp".
> >>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> of
> >>> type "FileType_NcPinterp".
> >>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> of
> >>> type "FileType_NcPinterp".
> >>> GSL_RNG_TYPE=mt19937
> >>> GSL_RNG_SEED=2307545168
> >>> DEBUG 1: Forecast File: ../datos/wrfout_d01_2012-07-03_PLEVELS
> >>> DEBUG 1: Climatology File: none
> >>> DEBUG 1: Observation File:
../datos/Surface_METAR_20120703_0000.nc
> >>> DEBUG 2:
> >>> DEBUG 2:
> >>>
> >>
>
--------------------------------------------------------------------------------
> >>> DEBUG 2:
> >>> DEBUG 2: Reading data for SPFHP500.
> >>> ERROR :
> >>> ERROR : process_fcst_climo_files() -> no fields matching
SPFHP500
> found
> >> in
> >>> file: ../datos/wrfout_d01_2012-07-03_PLEVELS
> >>> ERROR :
> >>>
> >>> i never call this SPFHP500 var.
> >>>
> >>> i think that the problem is in the config file and i don't
understand
> the
> >>> format for obs dictionary and fcst for netcdf can you expand the
> >>> explication whit a example like tutorial page for point-stat
tool.
> >>>
> >>> thanks for your time.
> >>>
> >>> On 17 September 2012 08:12, John Halley Gotway via RT <
> met_help at ucar.edu
> >>> wrote:
> >>>
> >>>> Jose,
> >>>>
> >>>> The Point-Stat configuration file is described in section 4.3.2
of the
> >> MET
> >>>> User's Guide:
> >>>>
> >>
www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v4.0.1.pdf
> >>>>
> >>>> Sample Point-Stat configuration files can also be found in the
> released
> >>>> code (I'm assuming METv4.0, here):
> >>>> METv4.0/data/config/PointStatConfig_default
> >>>> METv4.0/scripts/config/PointStatConfig
> >>>>
> >>>> And this README file describes the config file entries - search
for
> >>>> "NetCDF" in that file:
> >>>> METv4.0/data/config/README
> >>>>
> >>>> You may also find it helpful to look at the online tutorial:
> >>>>
> >>>>
> >>
>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.0/point_stat/index.php
> >>>>
> >>>> However, I see that none of those examples show using an input
NetCDF
> >>>> file. But here's an example of that for Grid-Stat:
> >>>> METv4.0/scripts/config/GridStatConfig_APCP_12
> >>>>
> >>>> It's pretty simple, for gridded NetCDF files, the "name" entry
is set
> to
> >>>> the NetCDF variable name. And the "level" specifies how to
extract
> the
> >> 2D
> >>>> gridded field from that variable. See that README
> >>>> file. One word of warning, if you're verifying a NetCDF
forecast file
> >> in
> >>>> Point-Stat, you'll need to define the "obs" dictionary to
specify
> which
> >>>> observations should be used for the verification. In
> >>>> Point-Stat, the observations are always specified using GRIB
> >> conventions.
> >>>>
> >>>> Give that a shot, and feel free to send us your config file and
input
> >> data
> >>>> for Point-Stat if you continue to experience problems.
> >>>>
> >>>> Hope that helps.
> >>>>
> >>>> John Halley Gotway
> >>>> met_help at ucar.edu
> >>>>
> >>>> On 09/14/2012 04:02 PM, Josue Don via RT wrote:
> >>>>>
> >>>>> Fri Sep 14 16:02:49 2012: Request 58307 was acted upon.
> >>>>> Transaction: Ticket created by eljoz89 at gmail.com
> >>>>> Queue: met_help
> >>>>> Subject: netcdf config file
> >>>>> Owner: Nobody
> >>>>> Requestors: eljoz89 at gmail.com
> >>>>> Status: new
> >>>>> Ticket <URL:
> >> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
> >>>>>
> >>>>>
> >>>>> hi, I try to run point-stat tool whit fcst_file in netcdf
format but
> no
> >>>> run
> >>>>> and i think that the problem is in my config file can you show
me an
> >>>>> example for that file?
> >>>>>
> >>>>
> >>>>
> >>>>
> >>
> >>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #58307] netcdf config file
From: John Halley Gotway
Time: Mon Sep 24 07:15:34 2012
Jose,
I apologize for the delay in getting back to you on this. Let me
start by saying that we typically recommend that users make use of the
Unified Post-Processor instead of the p_interp tool:
http://www.dtcenter.org/wrf-nmm/users/overview/upp_overview.php
The output of UPP is in GRIB format, for which MET was originally
designed. In addition, the copygb utility allows you to easily regrid
GRIB file, which is necessary when trying to do grid-to-grid
verification. No such regridder is available for the NetCDF p_interp
output. Also, MET can only handle p_interp variables that are defined
as mass points - meaning that it will not read the wind
fields which are defined using the staggered dimensions. Generally
speaking, MET is easier to use with GRIB data than with the NetCDF
output of p_interp.
Here are answers to your specific questions:
- Is it the same if using wrfout_d01_2012-07-03_SLEVELS file Zigma
levels of pressure? or Do I need to process it with "p_interp"?
You need to post-process the direct WRF output files prior to reading
them in MET. Above, I've recommended that you use UPP to do this.
Both UPP and p_interp interpolate the output from model sigma
levels to pressure levels. The verifying point observations are
defined on pressure and/or height levels. That's why the model output
must be placed onto pressure and/or height levels - the forecast
and observations must be on comparable vertical levels before they can
be compared. UPP also de-staggers the horizontal grid, so MET can
handle the wind fields when post-processing with UPP.
- How I can know the variables in a "netcdf" file and equivalent in
format "Grid1" to generate the values of "obs" and "FCST"?
- Could you tell me where to get the table of equivalence between
variables "netcdf" and "GRIB"?
There is no chart listing a correspondence from one to the other. But
it should be pretty easy to look at the GRIB table and the NetCDF
header and figure out which goes with which:
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
Please pay close attention to the units. MET is not doing any unit
conversions, and instead we assume that all the variables match the
GRIB1 convention. It's your responsibility to make sure that
the forecast and observations are in the same units.
- Is there a tool that allows me to change the file format "metar" (in
Surface_METAR_20120703_0000.nc) to a compatible format and what is the
right format?
No such tool currently exists. We've made note of this data source
and hope to add support for it in the next release of MET. But for
now, the easiest way to use those obs would be to write a script
to dump the contents of that NetCDF file into the 10-column ASCII
format accepted by the MET ascii2nc tool. Then, run it through
ascii2nc to get it into a form that MET expects, and use it for the
verification in Point-Stat.
Hope that helps.
Thanks,
John
On 09/19/2012 02:05 PM, Josue Don via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
>
> thanks, It was really helpful.
> I have some questions:
>
> - Is it the same if using wrfout_d01_2012-07-03_SLEVELS file Zigma
levels
> of pressure? or
> Do I need to process it with "p_interp"?
> - How I can know the variables in a "netcdf" file and equivalent in
format
> "Grid1" to generate the values of "obs" and "FCST"?
> - Could you tell me where to get the table of equivalence between
variables
> "netcdf" and "GRIB"?
> - Is there a tool that allows me to change the file format "metar"
> (in Surface_METAR_20120703_0000.nc) to a compatible format and what
is the
> right format?
>
>
> On 19 September 2012 09:43, John Halley Gotway via RT
<met_help at ucar.edu>wrote:
>
>> Jose,
>>
>> I retrieved your data files and ran them with MET. I started just
by
>> using the plot_data_plane tool to plot some data and make sure it
looks OK:
>>
>> # Plot RAINC for time 11
>> METv4.0/bin/plot_data_plane wrfout_d01_2012-07-03_PLEVELS
rainc_time11.ps'name="RAINC" ; level="(10,*,*)";' -v 10
>>
>> # Plot SKINTEMP for time 11
>> METv4.0/bin/plot_data_plane wrfout_d01_2012-07-03_PLEVELS
>> skintemp_time11.ps 'name="SKINTEMP" ; level="(10,*,*)";' -v 10
>>
>> The output images look fine and I see the your domain is primarily
over
>> Mexico.
>>
>> Next, by looking in your Point-Stat config file, I see that you
want to
>> verify the precipitation in the "RAINC(0,*,*)" variable. There are
a
>> couple problems with that:
>> - First, there are no observations of precipitation in the
GDAS point
>> observation file (netcdf.gdas.2012070300.nc).
>> - Second, while there are observations of precipitation in the
METAR
>> file (Surface_METAR_20120703_0000.nc), that file is not in a format
that
>> MET can handle. However, I have made a note of this
>> data source and we will consider adding support for it in future
releases
>> of MET. For now, to use those observations of precip, you'd need
to
>> reformat them in a way that MET can handle.
>>
>> To get you going, I switched over from verifying precip to
verifying
>> surface temperature. Here I'm verifying against observations of
surface
>> temperature from the GDAS file with message type ADPSFC
>> (means surface). You can plot the location of those observations
on your
>> domain using the following command:
>>
>> # Plot location of TMP for ADPSFC over the model domain
>> METv4.0/bin/plot_point_obs \
>> netcdf.gdas.2012070300.nc \
>> point_obs_tmp_adpsfc.ps \
>> -gc 11 -msg_typ ADPSFC \
>> -data_file wrfout_d01_2012-07-03_PLEVELS
>>
>> # Next, run Point-Stat to compare the model output to the
observations:
>> METv4.0/bin/point_stat \
>> wrfout_d01_2012-07-03_PLEVELS \
>> netcdf.gdas.2012070300.nc \
>> PointStatConfig_TMP_ADPSFC \
>> -outdir out \
>> -v 2
>>
>> I've attached the config file I used. Hope this helps getting you
going
>> in the right direction.
>>
>> John
>>
>> On 09/18/2012 11:28 PM, Josue Don via RT wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
>>>
>>> hi, i upload the files in don_data folder:
>>>
>>> fsct: PLEVELS -> post proces the SLEVEL file whit p_interp tool
(WRF post
>>> proces tools) (72 hours)
>>> SLEVELS -> direct to WRF v.3.1 output (72 hours)
>>>
>>> obs: ds336.0 -> Surface_METAR_20120703_0000.nc (24 hours)
>>> Surface_Synoptic_20120703_0000.nc.bz2 (24 hours)
>>> http://rda.ucar.edu/datasets/ds336.0/#description
>>>
>>> ds337.0 -> netcdf.gdas.2012070300.nc (6 hours)
>>> http://rda.ucar.edu/datasets/ds337.0/#description
>>>
>>> I try whit all this files.
>>>
>>>
>>> On 18 September 2012 10:24, John Halley Gotway via RT
<met_help at ucar.edu
>>> wrote:
>>>
>>>> Jose,
>>>>
>>>> It would probably be easiest for you to just send me your data
files,
>> and
>>>> I'll run it here. I'll need:
>>>> - What version of MET are you running?
>>>> - Forecast file: datos/wrfout_d01_2012-07-03_PLEVELS
>>>> - Observation file: datos/Surface_METAR_20120703_0000.nc
>>>> - Config file: proyecto/config/PointStatConfig_tutorial
>>>>
>>>> Please post them to our anonymous ftp site, following the
instructions
>>>> listed here:
>>>> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>>>>
>>>> Thanks,
>>>> John
>>>>
>>>> On 09/18/2012 09:09 AM, Josue Don via RT wrote:
>>>>>
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
>>>>>
>>>>> thanks for answer I try all your advices but I can't do that
poin-stat
>>>> tool
>>>>> works I'm using for inputs:
>>>>> fcst_file: wrf out file in netcdf format in levels zigma and
>> pressure.(i
>>>>> try to use pcp_combine but don't works)
>>>>> obs_file: http://rda.ucar.edu/datasets/ds337.0/#description
netcdf
>>>> suffice
>>>>> and http://rda.ucar.edu/datasets/ds336.0/#description in netcdf
metar
>> (i
>>>>> try to use madis2nc tool but don't works)
>>>>> config_file: attach it.
>>>>>
>>>>> the error message says:
>>>>> DEBUG 1: Default Config File:
>>>>> /home/xozue/MET/METv4.0/data/config/PointStatConfig_default
>>>>> DEBUG 1: User Config File:
proyecto/config/PointStatConfig_tutorial
>>>>> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
>>>>> Met2dDataFile object of type "FileType_NcPinterp".
>>>>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
>> of
>>>>> type "FileType_NcPinterp".
>>>>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
>> of
>>>>> type "FileType_NcPinterp".
>>>>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
>> of
>>>>> type "FileType_NcPinterp".
>>>>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
>> of
>>>>> type "FileType_NcPinterp".
>>>>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
>> of
>>>>> type "FileType_NcPinterp".
>>>>> GSL_RNG_TYPE=mt19937
>>>>> GSL_RNG_SEED=2307545168
>>>>> DEBUG 1: Forecast File: ../datos/wrfout_d01_2012-07-03_PLEVELS
>>>>> DEBUG 1: Climatology File: none
>>>>> DEBUG 1: Observation File:
../datos/Surface_METAR_20120703_0000.nc
>>>>> DEBUG 2:
>>>>> DEBUG 2:
>>>>>
>>>>
>>
--------------------------------------------------------------------------------
>>>>> DEBUG 2:
>>>>> DEBUG 2: Reading data for SPFHP500.
>>>>> ERROR :
>>>>> ERROR : process_fcst_climo_files() -> no fields matching
SPFHP500
>> found
>>>> in
>>>>> file: ../datos/wrfout_d01_2012-07-03_PLEVELS
>>>>> ERROR :
>>>>>
>>>>> i never call this SPFHP500 var.
>>>>>
>>>>> i think that the problem is in the config file and i don't
understand
>> the
>>>>> format for obs dictionary and fcst for netcdf can you expand the
>>>>> explication whit a example like tutorial page for point-stat
tool.
>>>>>
>>>>> thanks for your time.
>>>>>
>>>>> On 17 September 2012 08:12, John Halley Gotway via RT <
>> met_help at ucar.edu
>>>>> wrote:
>>>>>
>>>>>> Jose,
>>>>>>
>>>>>> The Point-Stat configuration file is described in section 4.3.2
of the
>>>> MET
>>>>>> User's Guide:
>>>>>>
>>>>
www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v4.0.1.pdf
>>>>>>
>>>>>> Sample Point-Stat configuration files can also be found in the
>> released
>>>>>> code (I'm assuming METv4.0, here):
>>>>>> METv4.0/data/config/PointStatConfig_default
>>>>>> METv4.0/scripts/config/PointStatConfig
>>>>>>
>>>>>> And this README file describes the config file entries - search
for
>>>>>> "NetCDF" in that file:
>>>>>> METv4.0/data/config/README
>>>>>>
>>>>>> You may also find it helpful to look at the online tutorial:
>>>>>>
>>>>>>
>>>>
>>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.0/point_stat/index.php
>>>>>>
>>>>>> However, I see that none of those examples show using an input
NetCDF
>>>>>> file. But here's an example of that for Grid-Stat:
>>>>>> METv4.0/scripts/config/GridStatConfig_APCP_12
>>>>>>
>>>>>> It's pretty simple, for gridded NetCDF files, the "name" entry
is set
>> to
>>>>>> the NetCDF variable name. And the "level" specifies how to
extract
>> the
>>>> 2D
>>>>>> gridded field from that variable. See that README
>>>>>> file. One word of warning, if you're verifying a NetCDF
forecast file
>>>> in
>>>>>> Point-Stat, you'll need to define the "obs" dictionary to
specify
>> which
>>>>>> observations should be used for the verification. In
>>>>>> Point-Stat, the observations are always specified using GRIB
>>>> conventions.
>>>>>>
>>>>>> Give that a shot, and feel free to send us your config file and
input
>>>> data
>>>>>> for Point-Stat if you continue to experience problems.
>>>>>>
>>>>>> Hope that helps.
>>>>>>
>>>>>> John Halley Gotway
>>>>>> met_help at ucar.edu
>>>>>>
>>>>>> On 09/14/2012 04:02 PM, Josue Don via RT wrote:
>>>>>>>
>>>>>>> Fri Sep 14 16:02:49 2012: Request 58307 was acted upon.
>>>>>>> Transaction: Ticket created by eljoz89 at gmail.com
>>>>>>> Queue: met_help
>>>>>>> Subject: netcdf config file
>>>>>>> Owner: Nobody
>>>>>>> Requestors: eljoz89 at gmail.com
>>>>>>> Status: new
>>>>>>> Ticket <URL:
>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=58307 >
>>>>>>>
>>>>>>>
>>>>>>> hi, I try to run point-stat tool whit fcst_file in netcdf
format but
>> no
>>>>>> run
>>>>>>> and i think that the problem is in my config file can you show
me an
>>>>>>> example for that file?
>>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>
>>>>
>>
>>
------------------------------------------------
More information about the Met_help
mailing list