[Met_help] [rt.rap.ucar.edu #40804] History for netcdf file for grid_stat (grid code error)
RAL HelpDesk {for Paul Oldenburg}
met_help at ucar.edu
Mon Oct 11 15:28:32 MDT 2010
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear Sir,
I want to use grid_stat tool to verify 2m temperature. Both (observation and
forecast) files are in netcdf format. when i am defining fcst_field, it is
asking for grid code. I am not using any grib file for this test. What
should I do ?
Thanks and Regards,
--
Hashmi Fatima
Research Scholar
Centre for Atmospheric Sciences
Indian Institute of Technology Delhi
Hauz Khas
New Delhi - 110016
INDIA
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #40804] netcdf file for grid_stat (grid code error)
From: Paul Oldenburg
Time: Thu Sep 16 08:57:42 2010
Hashmi,
Can you please send us an example of your model and obs data along
with your grid_stat config file? If the data is too
large to attach to an email, please upload it to our FTP site:
ftp ftp.rap.ucar.edu
username = anonymous
password = "your email address"
cd incoming/irap/met_help/fatima_data
put "your files"
Thanks,
Paul
RAL HelpDesk {for Hashmi Fatima} wrote:
> Thu Sep 16 01:44:32 2010: Request 40804 was acted upon.
> Transaction: Ticket created by hashmifatima at gmail.com
> Queue: met_help
> Subject: netcdf file for grid_stat (grid code error)
> Owner: Nobody
> Requestors: hashmifatima at gmail.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>
>
> Dear Sir,
>
> I want to use grid_stat tool to verify 2m temperature. Both
(observation and
> forecast) files are in netcdf format. when i am defining fcst_field,
it is
> asking for grid code. I am not using any grib file for this test.
What
> should I do ?
>
> Thanks and Regards,
------------------------------------------------
Subject: netcdf file for grid_stat (grid code error)
From: Hashmi Fatima
Time: Mon Sep 20 00:15:16 2010
Dear Sir,
I have put my files at the FTP site told by you. Please check.
One more thing, i want to know that do i need to give equal number of
threshold values as fcst_field/obs_field? Suppose I there is only one
fcst_field/ob_field, but I want to check it for various threshold
values.
In some other case, I have given 5 fcst_field/obs_field and
corresponding 5
threshold (fcst_thresh/obs_thresh) values. But I did not get anything
in
the output. I have got netcdf file in the output but nothing in the
STAT and
.txt files.
Thanks and Regards
Hashmi Fatima
On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for Paul Oldenburg} <
met_help at ucar.edu> wrote:
> Hashmi,
>
> Can you please send us an example of your model and obs data along
with
> your grid_stat config file? If the data is too
> large to attach to an email, please upload it to our FTP site:
>
> ftp ftp.rap.ucar.edu
> username = anonymous
> password = "your email address"
> cd incoming/irap/met_help/fatima_data
> put "your files"
>
> Thanks,
>
> Paul
>
> RAL HelpDesk {for Hashmi Fatima} wrote:
> > Thu Sep 16 01:44:32 2010: Request 40804 was acted upon.
> > Transaction: Ticket created by hashmifatima at gmail.com
> > Queue: met_help
> > Subject: netcdf file for grid_stat (grid code error)
> > Owner: Nobody
> > Requestors: hashmifatima at gmail.com
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >
> >
> > Dear Sir,
> >
> > I want to use grid_stat tool to verify 2m temperature. Both
(observation
> and
> > forecast) files are in netcdf format. when i am defining
fcst_field, it
> is
> > asking for grid code. I am not using any grib file for this test.
What
> > should I do ?
> >
> > Thanks and Regards,
>
>
>
--
Hashmi Fatima
Research Scholar
Centre for Atmospheric Sciences
Indian Institute of Technology Delhi
Hauz Khas
New Delhi - 110016
INDIA
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #40804] netcdf file for grid_stat (grid code error)
From: John Halley Gotway
Time: Mon Sep 20 15:27:03 2010
Hashmi,
Paul will be out of the office for the next couple of weeks, but I'll
be happy to help you with this issue.
There are a couple of points I'd like to make:
First, the NetCDF files you sent us are the NetCDF output of WRF. MET
does not directly read the NetCDF output of WRF. Instead, for
METv2.0, your WRF output must first be run through the
WRF-PostProcessor (WPP). Details of running WPP can be found in
chapter 7 of the WRF user's guide:
http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
WPP writes out GRIB data, which MET can read.
Second, once you've run your data through the WRF PostProcessor, you
can select multiple thresholds for Grid-Stat to use. You do so by
setting up the Grid-Stat configuration file something like this:
fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
obs_field[] = [];
fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240 ge260 ge280
ge300" ];
obs_thresh[] = [];
The settings above tell Grid-Stat to verify temperature at 2-meters as
well as temperature at 500mb. For 2-meter temp, it'll use the 2
thresholds listed (ge273 ge283) and for temp at 500mb, it'll use
the 6 thresholds listed (ge200.0 ge220 ge240 ge260 ge280 ge300).
So the real task is running your WRF output files through WPP.
Lastly, I wanted to tell you that in METv3.0 (to be released shortly),
we've added support for reading the NetCDF output of the pinterp tool
(which only processes ARW output). So once METv3.0 is
available, you can use either the GRIB output of WPP or the NetCDF
output of the pinterp tool.
Hope that helps.
Thanks,
John Halley Gotway
met_help at ucar.edu
RAL HelpDesk {for Hashmi Fatima} wrote:
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>
> Dear Sir,
>
> I have put my files at the FTP site told by you. Please check.
>
> One more thing, i want to know that do i need to give equal number
of
> threshold values as fcst_field/obs_field? Suppose I there is only
one
> fcst_field/ob_field, but I want to check it for various threshold
values.
>
> In some other case, I have given 5 fcst_field/obs_field and
corresponding 5
> threshold (fcst_thresh/obs_thresh) values. But I did not get
anything in
> the output. I have got netcdf file in the output but nothing in the
STAT and
> .txt files.
>
> Thanks and Regards
> Hashmi Fatima
>
> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for Paul Oldenburg} <
> met_help at ucar.edu> wrote:
>
>> Hashmi,
>>
>> Can you please send us an example of your model and obs data along
with
>> your grid_stat config file? If the data is too
>> large to attach to an email, please upload it to our FTP site:
>>
>> ftp ftp.rap.ucar.edu
>> username = anonymous
>> password = "your email address"
>> cd incoming/irap/met_help/fatima_data
>> put "your files"
>>
>> Thanks,
>>
>> Paul
>>
>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted upon.
>>> Transaction: Ticket created by hashmifatima at gmail.com
>>> Queue: met_help
>>> Subject: netcdf file for grid_stat (grid code error)
>>> Owner: Nobody
>>> Requestors: hashmifatima at gmail.com
>>> Status: new
>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>
>>>
>>> Dear Sir,
>>>
>>> I want to use grid_stat tool to verify 2m temperature. Both
(observation
>> and
>>> forecast) files are in netcdf format. when i am defining
fcst_field, it
>> is
>>> asking for grid code. I am not using any grib file for this test.
What
>>> should I do ?
>>>
>>> Thanks and Regards,
>>
>>
>
>
------------------------------------------------
Subject: netcdf file for grid_stat (grid code error)
From: Hashmi Fatima
Time: Mon Sep 27 04:35:55 2010
Dear Sir,
Thanks for your help.Its working now.
I have run wavelet tool for computation of RH. Now i want to use the
output
of wavelet for StatAnalysis. but its not working. I am sending you
configuration file of wavelet and StatAnalysis. Could you please
check.
One more thing I want to know that in the wavelet output, why does it
take
that much of area. I have not assigned any area for the verification
except
full. I have masked the whole grid. I am attaching that output too as
Wavelet_output.png.
Thanks and Regards
Hashmi fatima
On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for John Halley Gotway}
<
met_help at ucar.edu> wrote:
> Hashmi,
>
> Paul will be out of the office for the next couple of weeks, but
I'll be
> happy to help you with this issue.
>
> There are a couple of points I'd like to make:
>
> First, the NetCDF files you sent us are the NetCDF output of WRF.
MET does
> not directly read the NetCDF output of WRF. Instead, for METv2.0,
your WRF
> output must first be run through the
> WRF-PostProcessor (WPP). Details of running WPP can be found in
chapter 7
> of the WRF user's guide:
>
> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
> WPP writes out GRIB data, which MET can read.
>
> Second, once you've run your data through the WRF PostProcessor, you
can
> select multiple thresholds for Grid-Stat to use. You do so by
setting up
> the Grid-Stat configuration file something like this:
>
> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
> obs_field[] = [];
> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240 ge260 ge280
ge300" ];
> obs_thresh[] = [];
>
> The settings above tell Grid-Stat to verify temperature at 2-meters
as well
> as temperature at 500mb. For 2-meter temp, it'll use the 2
thresholds
> listed (ge273 ge283) and for temp at 500mb, it'll use
> the 6 thresholds listed (ge200.0 ge220 ge240 ge260 ge280 ge300).
>
> So the real task is running your WRF output files through WPP.
>
> Lastly, I wanted to tell you that in METv3.0 (to be released
shortly),
> we've added support for reading the NetCDF output of the pinterp
tool (which
> only processes ARW output). So once METv3.0 is
> available, you can use either the GRIB output of WPP or the NetCDF
output
> of the pinterp tool.
>
> Hope that helps.
>
> Thanks,
> John Halley Gotway
> met_help at ucar.edu
>
>
> RAL HelpDesk {for Hashmi Fatima} wrote:
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >
> > Dear Sir,
> >
> > I have put my files at the FTP site told by you. Please check.
> >
> > One more thing, i want to know that do i need to give equal number
of
> > threshold values as fcst_field/obs_field? Suppose I there is only
one
> > fcst_field/ob_field, but I want to check it for various threshold
values.
> >
> > In some other case, I have given 5 fcst_field/obs_field and
corresponding
> 5
> > threshold (fcst_thresh/obs_thresh) values. But I did not get
anything in
> > the output. I have got netcdf file in the output but nothing in
the STAT
> and
> > .txt files.
> >
> > Thanks and Regards
> > Hashmi Fatima
> >
> > On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for Paul Oldenburg}
<
> > met_help at ucar.edu> wrote:
> >
> >> Hashmi,
> >>
> >> Can you please send us an example of your model and obs data
along with
> >> your grid_stat config file? If the data is too
> >> large to attach to an email, please upload it to our FTP site:
> >>
> >> ftp ftp.rap.ucar.edu
> >> username = anonymous
> >> password = "your email address"
> >> cd incoming/irap/met_help/fatima_data
> >> put "your files"
> >>
> >> Thanks,
> >>
> >> Paul
> >>
> >> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>> Thu Sep 16 01:44:32 2010: Request 40804 was acted upon.
> >>> Transaction: Ticket created by hashmifatima at gmail.com
> >>> Queue: met_help
> >>> Subject: netcdf file for grid_stat (grid code error)
> >>> Owner: Nobody
> >>> Requestors: hashmifatima at gmail.com
> >>> Status: new
> >>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
> >>>
> >>>
> >>> Dear Sir,
> >>>
> >>> I want to use grid_stat tool to verify 2m temperature. Both
> (observation
> >> and
> >>> forecast) files are in netcdf format. when i am defining
fcst_field, it
> >> is
> >>> asking for grid code. I am not using any grib file for this
test. What
> >>> should I do ?
> >>>
> >>> Thanks and Regards,
> >>
> >>
> >
> >
>
>
--
Hashmi Fatima
Research Scholar
Centre for Atmospheric Sciences
Indian Institute of Technology Delhi
Hauz Khas
New Delhi - 110016
INDIA
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #40804] netcdf file for grid_stat (grid code error)
From: Paul Oldenburg
Time: Mon Sep 27 10:55:29 2010
Hashmi,
I (Paul) have returned to my met_help responsibilities and I have
followed your interactions with John. In order to
help you run stat_analysis, I would like to have you send me some of
your data via FTP. First, can you give us the WPP
output that you fed into wavelet_stat? Also, can you give us your
wavelet_stat output? Finally, please copy the error
message that you get from stat_analysis in your response. That
information will be very helpful.
Concerning the masking region selected by wavelet_stat, the entire
grid cannot be processed because wavelet_stat can
process only dyadic regions - that is square regions whose sides are a
power of 2 grid points (like 256 or 2048). By
default, wavelet_stat should try to guess the best regions to use in
your domain. If you send us your WPP output, we
will be able to help you configure the best masking region for your
case.
We look forward to receiving your data. If you have any questions,
please let me know.
Thanks,
Paul
RAL HelpDesk {for Hashmi Fatima} wrote:
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>
> Dear Sir,
>
> Thanks for your help.Its working now.
>
> I have run wavelet tool for computation of RH. Now i want to use the
output
> of wavelet for StatAnalysis. but its not working. I am sending you
> configuration file of wavelet and StatAnalysis. Could you please
check.
>
> One more thing I want to know that in the wavelet output, why does
it take
> that much of area. I have not assigned any area for the verification
except
> full. I have masked the whole grid. I am attaching that output too
as
> Wavelet_output.png.
>
> Thanks and Regards
> Hashmi fatima
>
>
> On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for John Halley
Gotway} <
> met_help at ucar.edu> wrote:
>
>> Hashmi,
>>
>> Paul will be out of the office for the next couple of weeks, but
I'll be
>> happy to help you with this issue.
>>
>> There are a couple of points I'd like to make:
>>
>> First, the NetCDF files you sent us are the NetCDF output of WRF.
MET does
>> not directly read the NetCDF output of WRF. Instead, for METv2.0,
your WRF
>> output must first be run through the
>> WRF-PostProcessor (WPP). Details of running WPP can be found in
chapter 7
>> of the WRF user's guide:
>>
>> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
>> WPP writes out GRIB data, which MET can read.
>>
>> Second, once you've run your data through the WRF PostProcessor,
you can
>> select multiple thresholds for Grid-Stat to use. You do so by
setting up
>> the Grid-Stat configuration file something like this:
>>
>> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
>> obs_field[] = [];
>> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240 ge260 ge280
ge300" ];
>> obs_thresh[] = [];
>>
>> The settings above tell Grid-Stat to verify temperature at 2-meters
as well
>> as temperature at 500mb. For 2-meter temp, it'll use the 2
thresholds
>> listed (ge273 ge283) and for temp at 500mb, it'll use
>> the 6 thresholds listed (ge200.0 ge220 ge240 ge260 ge280 ge300).
>>
>> So the real task is running your WRF output files through WPP.
>>
>> Lastly, I wanted to tell you that in METv3.0 (to be released
shortly),
>> we've added support for reading the NetCDF output of the pinterp
tool (which
>> only processes ARW output). So once METv3.0 is
>> available, you can use either the GRIB output of WPP or the NetCDF
output
>> of the pinterp tool.
>>
>> Hope that helps.
>>
>> Thanks,
>> John Halley Gotway
>> met_help at ucar.edu
>>
>>
>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>
>>> Dear Sir,
>>>
>>> I have put my files at the FTP site told by you. Please check.
>>>
>>> One more thing, i want to know that do i need to give equal number
of
>>> threshold values as fcst_field/obs_field? Suppose I there is only
one
>>> fcst_field/ob_field, but I want to check it for various threshold
values.
>>>
>>> In some other case, I have given 5 fcst_field/obs_field and
corresponding
>> 5
>>> threshold (fcst_thresh/obs_thresh) values. But I did not get
anything in
>>> the output. I have got netcdf file in the output but nothing in
the STAT
>> and
>>> .txt files.
>>>
>>> Thanks and Regards
>>> Hashmi Fatima
>>>
>>> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for Paul Oldenburg}
<
>>> met_help at ucar.edu> wrote:
>>>
>>>> Hashmi,
>>>>
>>>> Can you please send us an example of your model and obs data
along with
>>>> your grid_stat config file? If the data is too
>>>> large to attach to an email, please upload it to our FTP site:
>>>>
>>>> ftp ftp.rap.ucar.edu
>>>> username = anonymous
>>>> password = "your email address"
>>>> cd incoming/irap/met_help/fatima_data
>>>> put "your files"
>>>>
>>>> Thanks,
>>>>
>>>> Paul
>>>>
>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted upon.
>>>>> Transaction: Ticket created by hashmifatima at gmail.com
>>>>> Queue: met_help
>>>>> Subject: netcdf file for grid_stat (grid code error)
>>>>> Owner: Nobody
>>>>> Requestors: hashmifatima at gmail.com
>>>>> Status: new
>>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
>>>>>
>>>>>
>>>>> Dear Sir,
>>>>>
>>>>> I want to use grid_stat tool to verify 2m temperature. Both
>> (observation
>>>> and
>>>>> forecast) files are in netcdf format. when i am defining
fcst_field, it
>>>> is
>>>>> asking for grid code. I am not using any grib file for this
test. What
>>>>> should I do ?
>>>>>
>>>>> Thanks and Regards,
>>>>
>>>
>>
>
>
>
>
------------------------------------------------------------------------
>
------------------------------------------------
Subject: netcdf file for grid_stat (grid code error)
From: Hashmi Fatima
Time: Tue Sep 28 06:09:09 2010
Dear Paul
I tried to upload my files at the FTP server you have told me in
previous
mails. But it says " Permission denied on server".
Regards
Hashmi Fatima
On Mon, Sep 27, 2010 at 10:25 PM, RAL HelpDesk {for Paul Oldenburg} <
met_help at ucar.edu> wrote:
> Hashmi,
>
> I (Paul) have returned to my met_help responsibilities and I have
followed
> your interactions with John. In order to
> help you run stat_analysis, I would like to have you send me some of
your
> data via FTP. First, can you give us the WPP
> output that you fed into wavelet_stat? Also, can you give us your
> wavelet_stat output? Finally, please copy the error
> message that you get from stat_analysis in your response. That
information
> will be very helpful.
>
> Concerning the masking region selected by wavelet_stat, the entire
grid
> cannot be processed because wavelet_stat can
> process only dyadic regions - that is square regions whose sides are
a
> power of 2 grid points (like 256 or 2048). By
> default, wavelet_stat should try to guess the best regions to use in
your
> domain. If you send us your WPP output, we
> will be able to help you configure the best masking region for your
case.
>
> We look forward to receiving your data. If you have any questions,
please
> let me know.
>
> Thanks,
>
> Paul
>
>
> RAL HelpDesk {for Hashmi Fatima} wrote:
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >
> > Dear Sir,
> >
> > Thanks for your help.Its working now.
> >
> > I have run wavelet tool for computation of RH. Now i want to use
the
> output
> > of wavelet for StatAnalysis. but its not working. I am sending you
> > configuration file of wavelet and StatAnalysis. Could you please
check.
> >
> > One more thing I want to know that in the wavelet output, why does
it
> take
> > that much of area. I have not assigned any area for the
verification
> except
> > full. I have masked the whole grid. I am attaching that output too
as
> > Wavelet_output.png.
> >
> > Thanks and Regards
> > Hashmi fatima
> >
> >
> > On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for John Halley
Gotway} <
> > met_help at ucar.edu> wrote:
> >
> >> Hashmi,
> >>
> >> Paul will be out of the office for the next couple of weeks, but
I'll be
> >> happy to help you with this issue.
> >>
> >> There are a couple of points I'd like to make:
> >>
> >> First, the NetCDF files you sent us are the NetCDF output of WRF.
MET
> does
> >> not directly read the NetCDF output of WRF. Instead, for
METv2.0, your
> WRF
> >> output must first be run through the
> >> WRF-PostProcessor (WPP). Details of running WPP can be found in
chapter
> 7
> >> of the WRF user's guide:
> >>
> >>
> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
> >> WPP writes out GRIB data, which MET can read.
> >>
> >> Second, once you've run your data through the WRF PostProcessor,
you can
> >> select multiple thresholds for Grid-Stat to use. You do so by
setting
> up
> >> the Grid-Stat configuration file something like this:
> >>
> >> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
> >> obs_field[] = [];
> >> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240 ge260 ge280
ge300"
> ];
> >> obs_thresh[] = [];
> >>
> >> The settings above tell Grid-Stat to verify temperature at 2-
meters as
> well
> >> as temperature at 500mb. For 2-meter temp, it'll use the 2
thresholds
> >> listed (ge273 ge283) and for temp at 500mb, it'll use
> >> the 6 thresholds listed (ge200.0 ge220 ge240 ge260 ge280 ge300).
> >>
> >> So the real task is running your WRF output files through WPP.
> >>
> >> Lastly, I wanted to tell you that in METv3.0 (to be released
shortly),
> >> we've added support for reading the NetCDF output of the pinterp
tool
> (which
> >> only processes ARW output). So once METv3.0 is
> >> available, you can use either the GRIB output of WPP or the
NetCDF
> output
> >> of the pinterp tool.
> >>
> >> Hope that helps.
> >>
> >> Thanks,
> >> John Halley Gotway
> >> met_help at ucar.edu
> >>
> >>
> >> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >>>
> >>> Dear Sir,
> >>>
> >>> I have put my files at the FTP site told by you. Please check.
> >>>
> >>> One more thing, i want to know that do i need to give equal
number of
> >>> threshold values as fcst_field/obs_field? Suppose I there is
only one
> >>> fcst_field/ob_field, but I want to check it for various
threshold
> values.
> >>>
> >>> In some other case, I have given 5 fcst_field/obs_field and
> corresponding
> >> 5
> >>> threshold (fcst_thresh/obs_thresh) values. But I did not get
anything
> in
> >>> the output. I have got netcdf file in the output but nothing in
the
> STAT
> >> and
> >>> .txt files.
> >>>
> >>> Thanks and Regards
> >>> Hashmi Fatima
> >>>
> >>> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for Paul
Oldenburg} <
> >>> met_help at ucar.edu> wrote:
> >>>
> >>>> Hashmi,
> >>>>
> >>>> Can you please send us an example of your model and obs data
along
> with
> >>>> your grid_stat config file? If the data is too
> >>>> large to attach to an email, please upload it to our FTP site:
> >>>>
> >>>> ftp ftp.rap.ucar.edu
> >>>> username = anonymous
> >>>> password = "your email address"
> >>>> cd incoming/irap/met_help/fatima_data
> >>>> put "your files"
> >>>>
> >>>> Thanks,
> >>>>
> >>>> Paul
> >>>>
> >>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted upon.
> >>>>> Transaction: Ticket created by hashmifatima at gmail.com
> >>>>> Queue: met_help
> >>>>> Subject: netcdf file for grid_stat (grid code error)
> >>>>> Owner: Nobody
> >>>>> Requestors: hashmifatima at gmail.com
> >>>>> Status: new
> >>>>> Ticket <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
> >>>>>
> >>>>>
> >>>>> Dear Sir,
> >>>>>
> >>>>> I want to use grid_stat tool to verify 2m temperature. Both
> >> (observation
> >>>> and
> >>>>> forecast) files are in netcdf format. when i am defining
fcst_field,
> it
> >>>> is
> >>>>> asking for grid code. I am not using any grib file for this
test.
> What
> >>>>> should I do ?
> >>>>>
> >>>>> Thanks and Regards,
> >>>>
> >>>
> >>
> >
> >
> >
> >
------------------------------------------------------------------------
> >
>
>
>
--
Hashmi Fatima
Research Scholar
Centre for Atmospheric Sciences
Indian Institute of Technology Delhi
Hauz Khas
New Delhi - 110016
INDIA
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #40804] netcdf file for grid_stat (grid code error)
From: Paul Oldenburg
Time: Tue Sep 28 10:55:37 2010
Hashmi,
We are investigating a potential problem with our FTP site. I will
let you know when it is fixed, and then we can
proceed with addressing your wavelet_stat question.
Thanks,
Paul
RAL HelpDesk {for Hashmi Fatima} wrote:
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>
> Dear Paul
>
> I tried to upload my files at the FTP server you have told me in
previous
> mails. But it says " Permission denied on server".
>
> Regards
> Hashmi Fatima
>
> On Mon, Sep 27, 2010 at 10:25 PM, RAL HelpDesk {for Paul Oldenburg}
<
> met_help at ucar.edu> wrote:
>
>> Hashmi,
>>
>> I (Paul) have returned to my met_help responsibilities and I have
followed
>> your interactions with John. In order to
>> help you run stat_analysis, I would like to have you send me some
of your
>> data via FTP. First, can you give us the WPP
>> output that you fed into wavelet_stat? Also, can you give us your
>> wavelet_stat output? Finally, please copy the error
>> message that you get from stat_analysis in your response. That
information
>> will be very helpful.
>>
>> Concerning the masking region selected by wavelet_stat, the entire
grid
>> cannot be processed because wavelet_stat can
>> process only dyadic regions - that is square regions whose sides
are a
>> power of 2 grid points (like 256 or 2048). By
>> default, wavelet_stat should try to guess the best regions to use
in your
>> domain. If you send us your WPP output, we
>> will be able to help you configure the best masking region for your
case.
>>
>> We look forward to receiving your data. If you have any questions,
please
>> let me know.
>>
>> Thanks,
>>
>> Paul
>>
>>
>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>
>>> Dear Sir,
>>>
>>> Thanks for your help.Its working now.
>>>
>>> I have run wavelet tool for computation of RH. Now i want to use
the
>> output
>>> of wavelet for StatAnalysis. but its not working. I am sending you
>>> configuration file of wavelet and StatAnalysis. Could you please
check.
>>>
>>> One more thing I want to know that in the wavelet output, why does
it
>> take
>>> that much of area. I have not assigned any area for the
verification
>> except
>>> full. I have masked the whole grid. I am attaching that output too
as
>>> Wavelet_output.png.
>>>
>>> Thanks and Regards
>>> Hashmi fatima
>>>
>>>
>>> On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for John Halley
Gotway} <
>>> met_help at ucar.edu> wrote:
>>>
>>>> Hashmi,
>>>>
>>>> Paul will be out of the office for the next couple of weeks, but
I'll be
>>>> happy to help you with this issue.
>>>>
>>>> There are a couple of points I'd like to make:
>>>>
>>>> First, the NetCDF files you sent us are the NetCDF output of WRF.
MET
>> does
>>>> not directly read the NetCDF output of WRF. Instead, for
METv2.0, your
>> WRF
>>>> output must first be run through the
>>>> WRF-PostProcessor (WPP). Details of running WPP can be found in
chapter
>> 7
>>>> of the WRF user's guide:
>>>>
>>>>
>> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
>>>> WPP writes out GRIB data, which MET can read.
>>>>
>>>> Second, once you've run your data through the WRF PostProcessor,
you can
>>>> select multiple thresholds for Grid-Stat to use. You do so by
setting
>> up
>>>> the Grid-Stat configuration file something like this:
>>>>
>>>> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
>>>> obs_field[] = [];
>>>> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240 ge260 ge280
ge300"
>> ];
>>>> obs_thresh[] = [];
>>>>
>>>> The settings above tell Grid-Stat to verify temperature at 2-
meters as
>> well
>>>> as temperature at 500mb. For 2-meter temp, it'll use the 2
thresholds
>>>> listed (ge273 ge283) and for temp at 500mb, it'll use
>>>> the 6 thresholds listed (ge200.0 ge220 ge240 ge260 ge280 ge300).
>>>>
>>>> So the real task is running your WRF output files through WPP.
>>>>
>>>> Lastly, I wanted to tell you that in METv3.0 (to be released
shortly),
>>>> we've added support for reading the NetCDF output of the pinterp
tool
>> (which
>>>> only processes ARW output). So once METv3.0 is
>>>> available, you can use either the GRIB output of WPP or the
NetCDF
>> output
>>>> of the pinterp tool.
>>>>
>>>> Hope that helps.
>>>>
>>>> Thanks,
>>>> John Halley Gotway
>>>> met_help at ucar.edu
>>>>
>>>>
>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>
>>>>> Dear Sir,
>>>>>
>>>>> I have put my files at the FTP site told by you. Please check.
>>>>>
>>>>> One more thing, i want to know that do i need to give equal
number of
>>>>> threshold values as fcst_field/obs_field? Suppose I there is
only one
>>>>> fcst_field/ob_field, but I want to check it for various
threshold
>> values.
>>>>> In some other case, I have given 5 fcst_field/obs_field and
>> corresponding
>>>> 5
>>>>> threshold (fcst_thresh/obs_thresh) values. But I did not get
anything
>> in
>>>>> the output. I have got netcdf file in the output but nothing in
the
>> STAT
>>>> and
>>>>> .txt files.
>>>>>
>>>>> Thanks and Regards
>>>>> Hashmi Fatima
>>>>>
>>>>> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for Paul
Oldenburg} <
>>>>> met_help at ucar.edu> wrote:
>>>>>
>>>>>> Hashmi,
>>>>>>
>>>>>> Can you please send us an example of your model and obs data
along
>> with
>>>>>> your grid_stat config file? If the data is too
>>>>>> large to attach to an email, please upload it to our FTP site:
>>>>>>
>>>>>> ftp ftp.rap.ucar.edu
>>>>>> username = anonymous
>>>>>> password = "your email address"
>>>>>> cd incoming/irap/met_help/fatima_data
>>>>>> put "your files"
>>>>>>
>>>>>> Thanks,
>>>>>>
>>>>>> Paul
>>>>>>
>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted upon.
>>>>>>> Transaction: Ticket created by hashmifatima at gmail.com
>>>>>>> Queue: met_help
>>>>>>> Subject: netcdf file for grid_stat (grid code error)
>>>>>>> Owner: Nobody
>>>>>>> Requestors: hashmifatima at gmail.com
>>>>>>> Status: new
>>>>>>> Ticket <URL:
>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
>>>>>>>
>>>>>>> Dear Sir,
>>>>>>>
>>>>>>> I want to use grid_stat tool to verify 2m temperature. Both
>>>> (observation
>>>>>> and
>>>>>>> forecast) files are in netcdf format. when i am defining
fcst_field,
>> it
>>>>>> is
>>>>>>> asking for grid code. I am not using any grib file for this
test.
>> What
>>>>>>> should I do ?
>>>>>>>
>>>>>>> Thanks and Regards,
>>>
>>>
>>>
------------------------------------------------------------------------
>>>
>>
>>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #40804] netcdf file for grid_stat (grid code error)
From: Paul Oldenburg
Time: Tue Sep 28 12:52:16 2010
Hashmi,
Please try to invoke your ftp client using the -p command line
argument, like so:
ftp -p ftp.rap.ucar.edu
username = anonymous
password = "your email address"
cd incoming/irap/met_help/fatima_data
put "your files"
Let us know how this approach goes.
Thanks,
Paul
Paul Oldenburg wrote:
> Hashmi,
>
> We are investigating a potential problem with our FTP site. I will
let you know when it is fixed, and then we can
> proceed with addressing your wavelet_stat question.
>
> Thanks,
>
> Paul
>
> RAL HelpDesk {for Hashmi Fatima} wrote:
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>
>> Dear Paul
>>
>> I tried to upload my files at the FTP server you have told me in
previous
>> mails. But it says " Permission denied on server".
>>
>> Regards
>> Hashmi Fatima
>>
>> On Mon, Sep 27, 2010 at 10:25 PM, RAL HelpDesk {for Paul Oldenburg}
<
>> met_help at ucar.edu> wrote:
>>
>>> Hashmi,
>>>
>>> I (Paul) have returned to my met_help responsibilities and I have
followed
>>> your interactions with John. In order to
>>> help you run stat_analysis, I would like to have you send me some
of your
>>> data via FTP. First, can you give us the WPP
>>> output that you fed into wavelet_stat? Also, can you give us your
>>> wavelet_stat output? Finally, please copy the error
>>> message that you get from stat_analysis in your response. That
information
>>> will be very helpful.
>>>
>>> Concerning the masking region selected by wavelet_stat, the entire
grid
>>> cannot be processed because wavelet_stat can
>>> process only dyadic regions - that is square regions whose sides
are a
>>> power of 2 grid points (like 256 or 2048). By
>>> default, wavelet_stat should try to guess the best regions to use
in your
>>> domain. If you send us your WPP output, we
>>> will be able to help you configure the best masking region for
your case.
>>>
>>> We look forward to receiving your data. If you have any
questions, please
>>> let me know.
>>>
>>> Thanks,
>>>
>>> Paul
>>>
>>>
>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>
>>>> Dear Sir,
>>>>
>>>> Thanks for your help.Its working now.
>>>>
>>>> I have run wavelet tool for computation of RH. Now i want to use
the
>>> output
>>>> of wavelet for StatAnalysis. but its not working. I am sending
you
>>>> configuration file of wavelet and StatAnalysis. Could you please
check.
>>>>
>>>> One more thing I want to know that in the wavelet output, why
does it
>>> take
>>>> that much of area. I have not assigned any area for the
verification
>>> except
>>>> full. I have masked the whole grid. I am attaching that output
too as
>>>> Wavelet_output.png.
>>>>
>>>> Thanks and Regards
>>>> Hashmi fatima
>>>>
>>>>
>>>> On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for John Halley
Gotway} <
>>>> met_help at ucar.edu> wrote:
>>>>
>>>>> Hashmi,
>>>>>
>>>>> Paul will be out of the office for the next couple of weeks, but
I'll be
>>>>> happy to help you with this issue.
>>>>>
>>>>> There are a couple of points I'd like to make:
>>>>>
>>>>> First, the NetCDF files you sent us are the NetCDF output of
WRF. MET
>>> does
>>>>> not directly read the NetCDF output of WRF. Instead, for
METv2.0, your
>>> WRF
>>>>> output must first be run through the
>>>>> WRF-PostProcessor (WPP). Details of running WPP can be found in
chapter
>>> 7
>>>>> of the WRF user's guide:
>>>>>
>>>>>
>>> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
>>>>> WPP writes out GRIB data, which MET can read.
>>>>>
>>>>> Second, once you've run your data through the WRF PostProcessor,
you can
>>>>> select multiple thresholds for Grid-Stat to use. You do so by
setting
>>> up
>>>>> the Grid-Stat configuration file something like this:
>>>>>
>>>>> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
>>>>> obs_field[] = [];
>>>>> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240 ge260
ge280 ge300"
>>> ];
>>>>> obs_thresh[] = [];
>>>>>
>>>>> The settings above tell Grid-Stat to verify temperature at 2-
meters as
>>> well
>>>>> as temperature at 500mb. For 2-meter temp, it'll use the 2
thresholds
>>>>> listed (ge273 ge283) and for temp at 500mb, it'll use
>>>>> the 6 thresholds listed (ge200.0 ge220 ge240 ge260 ge280 ge300).
>>>>>
>>>>> So the real task is running your WRF output files through WPP.
>>>>>
>>>>> Lastly, I wanted to tell you that in METv3.0 (to be released
shortly),
>>>>> we've added support for reading the NetCDF output of the pinterp
tool
>>> (which
>>>>> only processes ARW output). So once METv3.0 is
>>>>> available, you can use either the GRIB output of WPP or the
NetCDF
>>> output
>>>>> of the pinterp tool.
>>>>>
>>>>> Hope that helps.
>>>>>
>>>>> Thanks,
>>>>> John Halley Gotway
>>>>> met_help at ucar.edu
>>>>>
>>>>>
>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>>
>>>>>> Dear Sir,
>>>>>>
>>>>>> I have put my files at the FTP site told by you. Please check.
>>>>>>
>>>>>> One more thing, i want to know that do i need to give equal
number of
>>>>>> threshold values as fcst_field/obs_field? Suppose I there is
only one
>>>>>> fcst_field/ob_field, but I want to check it for various
threshold
>>> values.
>>>>>> In some other case, I have given 5 fcst_field/obs_field and
>>> corresponding
>>>>> 5
>>>>>> threshold (fcst_thresh/obs_thresh) values. But I did not get
anything
>>> in
>>>>>> the output. I have got netcdf file in the output but nothing in
the
>>> STAT
>>>>> and
>>>>>> .txt files.
>>>>>>
>>>>>> Thanks and Regards
>>>>>> Hashmi Fatima
>>>>>>
>>>>>> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for Paul
Oldenburg} <
>>>>>> met_help at ucar.edu> wrote:
>>>>>>
>>>>>>> Hashmi,
>>>>>>>
>>>>>>> Can you please send us an example of your model and obs data
along
>>> with
>>>>>>> your grid_stat config file? If the data is too
>>>>>>> large to attach to an email, please upload it to our FTP site:
>>>>>>>
>>>>>>> ftp ftp.rap.ucar.edu
>>>>>>> username = anonymous
>>>>>>> password = "your email address"
>>>>>>> cd incoming/irap/met_help/fatima_data
>>>>>>> put "your files"
>>>>>>>
>>>>>>> Thanks,
>>>>>>>
>>>>>>> Paul
>>>>>>>
>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted upon.
>>>>>>>> Transaction: Ticket created by hashmifatima at gmail.com
>>>>>>>> Queue: met_help
>>>>>>>> Subject: netcdf file for grid_stat (grid code error)
>>>>>>>> Owner: Nobody
>>>>>>>> Requestors: hashmifatima at gmail.com
>>>>>>>> Status: new
>>>>>>>> Ticket <URL:
>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
>>>>>>>> Dear Sir,
>>>>>>>>
>>>>>>>> I want to use grid_stat tool to verify 2m temperature. Both
>>>>> (observation
>>>>>>> and
>>>>>>>> forecast) files are in netcdf format. when i am defining
fcst_field,
>>> it
>>>>>>> is
>>>>>>>> asking for grid code. I am not using any grib file for this
test.
>>> What
>>>>>>>> should I do ?
>>>>>>>>
>>>>>>>> Thanks and Regards,
>>>>
>>>>
------------------------------------------------------------------------
>>>>
>>>
>>
>
>
------------------------------------------------
Subject: netcdf file for grid_stat (grid code error)
From: Hashmi Fatima
Time: Wed Sep 29 00:09:55 2010
Dear Paul
ftp -p ftp.rap.ucar.edu is not working. It is giving me error
message.
ftp -p ftp.rap.ucar.edu
ftp: p: unknown option
Usage: ftp [-v] [-d] [-i] [-n] [-g] [-k realm] [-f] [-x] [-u] [-t]
[host]
Regards
Hashmi Fatima
On Wed, Sep 29, 2010 at 12:22 AM, RAL HelpDesk {for Paul Oldenburg} <
met_help at ucar.edu> wrote:
> Hashmi,
>
> Please try to invoke your ftp client using the -p command line
argument,
> like so:
>
> ftp -p ftp.rap.ucar.edu
> username = anonymous
> password = "your email address"
> cd incoming/irap/met_help/fatima_data
> put "your files"
>
> Let us know how this approach goes.
>
> Thanks,
>
> Paul
>
>
> Paul Oldenburg wrote:
> > Hashmi,
> >
> > We are investigating a potential problem with our FTP site. I
will let
> you know when it is fixed, and then we can
> > proceed with addressing your wavelet_stat question.
> >
> > Thanks,
> >
> > Paul
> >
> > RAL HelpDesk {for Hashmi Fatima} wrote:
> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >>
> >> Dear Paul
> >>
> >> I tried to upload my files at the FTP server you have told me in
> previous
> >> mails. But it says " Permission denied on server".
> >>
> >> Regards
> >> Hashmi Fatima
> >>
> >> On Mon, Sep 27, 2010 at 10:25 PM, RAL HelpDesk {for Paul
Oldenburg} <
> >> met_help at ucar.edu> wrote:
> >>
> >>> Hashmi,
> >>>
> >>> I (Paul) have returned to my met_help responsibilities and I
have
> followed
> >>> your interactions with John. In order to
> >>> help you run stat_analysis, I would like to have you send me
some of
> your
> >>> data via FTP. First, can you give us the WPP
> >>> output that you fed into wavelet_stat? Also, can you give us
your
> >>> wavelet_stat output? Finally, please copy the error
> >>> message that you get from stat_analysis in your response. That
> information
> >>> will be very helpful.
> >>>
> >>> Concerning the masking region selected by wavelet_stat, the
entire grid
> >>> cannot be processed because wavelet_stat can
> >>> process only dyadic regions - that is square regions whose sides
are a
> >>> power of 2 grid points (like 256 or 2048). By
> >>> default, wavelet_stat should try to guess the best regions to
use in
> your
> >>> domain. If you send us your WPP output, we
> >>> will be able to help you configure the best masking region for
your
> case.
> >>>
> >>> We look forward to receiving your data. If you have any
questions,
> please
> >>> let me know.
> >>>
> >>> Thanks,
> >>>
> >>> Paul
> >>>
> >>>
> >>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >>>>
> >>>> Dear Sir,
> >>>>
> >>>> Thanks for your help.Its working now.
> >>>>
> >>>> I have run wavelet tool for computation of RH. Now i want to
use the
> >>> output
> >>>> of wavelet for StatAnalysis. but its not working. I am sending
you
> >>>> configuration file of wavelet and StatAnalysis. Could you
please
> check.
> >>>>
> >>>> One more thing I want to know that in the wavelet output, why
does it
> >>> take
> >>>> that much of area. I have not assigned any area for the
verification
> >>> except
> >>>> full. I have masked the whole grid. I am attaching that output
too as
> >>>> Wavelet_output.png.
> >>>>
> >>>> Thanks and Regards
> >>>> Hashmi fatima
> >>>>
> >>>>
> >>>> On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for John Halley
Gotway}
> <
> >>>> met_help at ucar.edu> wrote:
> >>>>
> >>>>> Hashmi,
> >>>>>
> >>>>> Paul will be out of the office for the next couple of weeks,
but I'll
> be
> >>>>> happy to help you with this issue.
> >>>>>
> >>>>> There are a couple of points I'd like to make:
> >>>>>
> >>>>> First, the NetCDF files you sent us are the NetCDF output of
WRF.
> MET
> >>> does
> >>>>> not directly read the NetCDF output of WRF. Instead, for
METv2.0,
> your
> >>> WRF
> >>>>> output must first be run through the
> >>>>> WRF-PostProcessor (WPP). Details of running WPP can be found
in
> chapter
> >>> 7
> >>>>> of the WRF user's guide:
> >>>>>
> >>>>>
> >>>
> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
> >>>>> WPP writes out GRIB data, which MET can read.
> >>>>>
> >>>>> Second, once you've run your data through the WRF
PostProcessor, you
> can
> >>>>> select multiple thresholds for Grid-Stat to use. You do so by
> setting
> >>> up
> >>>>> the Grid-Stat configuration file something like this:
> >>>>>
> >>>>> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
> >>>>> obs_field[] = [];
> >>>>> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240 ge260
ge280
> ge300"
> >>> ];
> >>>>> obs_thresh[] = [];
> >>>>>
> >>>>> The settings above tell Grid-Stat to verify temperature at 2-
meters
> as
> >>> well
> >>>>> as temperature at 500mb. For 2-meter temp, it'll use the 2
> thresholds
> >>>>> listed (ge273 ge283) and for temp at 500mb, it'll use
> >>>>> the 6 thresholds listed (ge200.0 ge220 ge240 ge260 ge280
ge300).
> >>>>>
> >>>>> So the real task is running your WRF output files through WPP.
> >>>>>
> >>>>> Lastly, I wanted to tell you that in METv3.0 (to be released
> shortly),
> >>>>> we've added support for reading the NetCDF output of the
pinterp tool
> >>> (which
> >>>>> only processes ARW output). So once METv3.0 is
> >>>>> available, you can use either the GRIB output of WPP or the
NetCDF
> >>> output
> >>>>> of the pinterp tool.
> >>>>>
> >>>>> Hope that helps.
> >>>>>
> >>>>> Thanks,
> >>>>> John Halley Gotway
> >>>>> met_help at ucar.edu
> >>>>>
> >>>>>
> >>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804
>
> >>>>>>
> >>>>>> Dear Sir,
> >>>>>>
> >>>>>> I have put my files at the FTP site told by you. Please
check.
> >>>>>>
> >>>>>> One more thing, i want to know that do i need to give equal
number
> of
> >>>>>> threshold values as fcst_field/obs_field? Suppose I there is
only
> one
> >>>>>> fcst_field/ob_field, but I want to check it for various
threshold
> >>> values.
> >>>>>> In some other case, I have given 5 fcst_field/obs_field and
> >>> corresponding
> >>>>> 5
> >>>>>> threshold (fcst_thresh/obs_thresh) values. But I did not get
> anything
> >>> in
> >>>>>> the output. I have got netcdf file in the output but nothing
in the
> >>> STAT
> >>>>> and
> >>>>>> .txt files.
> >>>>>>
> >>>>>> Thanks and Regards
> >>>>>> Hashmi Fatima
> >>>>>>
> >>>>>> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for Paul
Oldenburg} <
> >>>>>> met_help at ucar.edu> wrote:
> >>>>>>
> >>>>>>> Hashmi,
> >>>>>>>
> >>>>>>> Can you please send us an example of your model and obs data
along
> >>> with
> >>>>>>> your grid_stat config file? If the data is too
> >>>>>>> large to attach to an email, please upload it to our FTP
site:
> >>>>>>>
> >>>>>>> ftp ftp.rap.ucar.edu
> >>>>>>> username = anonymous
> >>>>>>> password = "your email address"
> >>>>>>> cd incoming/irap/met_help/fatima_data
> >>>>>>> put "your files"
> >>>>>>>
> >>>>>>> Thanks,
> >>>>>>>
> >>>>>>> Paul
> >>>>>>>
> >>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted upon.
> >>>>>>>> Transaction: Ticket created by hashmifatima at gmail.com
> >>>>>>>> Queue: met_help
> >>>>>>>> Subject: netcdf file for grid_stat (grid code error)
> >>>>>>>> Owner: Nobody
> >>>>>>>> Requestors: hashmifatima at gmail.com
> >>>>>>>> Status: new
> >>>>>>>> Ticket <URL:
> >>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
> >>>>>>>> Dear Sir,
> >>>>>>>>
> >>>>>>>> I want to use grid_stat tool to verify 2m temperature. Both
> >>>>> (observation
> >>>>>>> and
> >>>>>>>> forecast) files are in netcdf format. when i am defining
> fcst_field,
> >>> it
> >>>>>>> is
> >>>>>>>> asking for grid code. I am not using any grib file for this
test.
> >>> What
> >>>>>>>> should I do ?
> >>>>>>>>
> >>>>>>>> Thanks and Regards,
> >>>>
> >>>>
>
------------------------------------------------------------------------
> >>>>
> >>>
> >>
> >
> >
>
>
>
--
Hashmi Fatima
Research Scholar
Centre for Atmospheric Sciences
Indian Institute of Technology Delhi
Hauz Khas
New Delhi - 110016
INDIA
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #40804] netcdf file for grid_stat (grid code error)
From: Paul Oldenburg
Time: Wed Sep 29 09:11:34 2010
Hashmi,
We apparently have a new FTP server configuration that we are still
getting used to. Can you please refer to the
instructions on the following website for uploading data and let us
know if you can get it to work:
RAL HelpDesk {for Hashmi Fatima} wrote:
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>
> Dear Paul
>
> ftp -p ftp.rap.ucar.edu is not working. It is giving me error
message.
>
> ftp -p ftp.rap.ucar.edu
> ftp: p: unknown option
> Usage: ftp [-v] [-d] [-i] [-n] [-g] [-k realm] [-f] [-x] [-u] [-t]
[host]
>
> Regards
> Hashmi Fatima
>
>
>
> On Wed, Sep 29, 2010 at 12:22 AM, RAL HelpDesk {for Paul Oldenburg}
<
> met_help at ucar.edu> wrote:
>
>> Hashmi,
>>
>> Please try to invoke your ftp client using the -p command line
argument,
>> like so:
>>
>> ftp -p ftp.rap.ucar.edu
>> username = anonymous
>> password = "your email address"
>> cd incoming/irap/met_help/fatima_data
>> put "your files"
>>
>> Let us know how this approach goes.
>>
>> Thanks,
>>
>> Paul
>>
>>
>> Paul Oldenburg wrote:
>>> Hashmi,
>>>
>>> We are investigating a potential problem with our FTP site. I
will let
>> you know when it is fixed, and then we can
>>> proceed with addressing your wavelet_stat question.
>>>
>>> Thanks,
>>>
>>> Paul
>>>
>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>
>>>> Dear Paul
>>>>
>>>> I tried to upload my files at the FTP server you have told me in
>> previous
>>>> mails. But it says " Permission denied on server".
>>>>
>>>> Regards
>>>> Hashmi Fatima
>>>>
>>>> On Mon, Sep 27, 2010 at 10:25 PM, RAL HelpDesk {for Paul
Oldenburg} <
>>>> met_help at ucar.edu> wrote:
>>>>
>>>>> Hashmi,
>>>>>
>>>>> I (Paul) have returned to my met_help responsibilities and I
have
>> followed
>>>>> your interactions with John. In order to
>>>>> help you run stat_analysis, I would like to have you send me
some of
>> your
>>>>> data via FTP. First, can you give us the WPP
>>>>> output that you fed into wavelet_stat? Also, can you give us
your
>>>>> wavelet_stat output? Finally, please copy the error
>>>>> message that you get from stat_analysis in your response. That
>> information
>>>>> will be very helpful.
>>>>>
>>>>> Concerning the masking region selected by wavelet_stat, the
entire grid
>>>>> cannot be processed because wavelet_stat can
>>>>> process only dyadic regions - that is square regions whose sides
are a
>>>>> power of 2 grid points (like 256 or 2048). By
>>>>> default, wavelet_stat should try to guess the best regions to
use in
>> your
>>>>> domain. If you send us your WPP output, we
>>>>> will be able to help you configure the best masking region for
your
>> case.
>>>>> We look forward to receiving your data. If you have any
questions,
>> please
>>>>> let me know.
>>>>>
>>>>> Thanks,
>>>>>
>>>>> Paul
>>>>>
>>>>>
>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>>
>>>>>> Dear Sir,
>>>>>>
>>>>>> Thanks for your help.Its working now.
>>>>>>
>>>>>> I have run wavelet tool for computation of RH. Now i want to
use the
>>>>> output
>>>>>> of wavelet for StatAnalysis. but its not working. I am sending
you
>>>>>> configuration file of wavelet and StatAnalysis. Could you
please
>> check.
>>>>>> One more thing I want to know that in the wavelet output, why
does it
>>>>> take
>>>>>> that much of area. I have not assigned any area for the
verification
>>>>> except
>>>>>> full. I have masked the whole grid. I am attaching that output
too as
>>>>>> Wavelet_output.png.
>>>>>>
>>>>>> Thanks and Regards
>>>>>> Hashmi fatima
>>>>>>
>>>>>>
>>>>>> On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for John Halley
Gotway}
>> <
>>>>>> met_help at ucar.edu> wrote:
>>>>>>
>>>>>>> Hashmi,
>>>>>>>
>>>>>>> Paul will be out of the office for the next couple of weeks,
but I'll
>> be
>>>>>>> happy to help you with this issue.
>>>>>>>
>>>>>>> There are a couple of points I'd like to make:
>>>>>>>
>>>>>>> First, the NetCDF files you sent us are the NetCDF output of
WRF.
>> MET
>>>>> does
>>>>>>> not directly read the NetCDF output of WRF. Instead, for
METv2.0,
>> your
>>>>> WRF
>>>>>>> output must first be run through the
>>>>>>> WRF-PostProcessor (WPP). Details of running WPP can be found
in
>> chapter
>>>>> 7
>>>>>>> of the WRF user's guide:
>>>>>>>
>>>>>>>
>> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
>>>>>>> WPP writes out GRIB data, which MET can read.
>>>>>>>
>>>>>>> Second, once you've run your data through the WRF
PostProcessor, you
>> can
>>>>>>> select multiple thresholds for Grid-Stat to use. You do so by
>> setting
>>>>> up
>>>>>>> the Grid-Stat configuration file something like this:
>>>>>>>
>>>>>>> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
>>>>>>> obs_field[] = [];
>>>>>>> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240 ge260
ge280
>> ge300"
>>>>> ];
>>>>>>> obs_thresh[] = [];
>>>>>>>
>>>>>>> The settings above tell Grid-Stat to verify temperature at 2-
meters
>> as
>>>>> well
>>>>>>> as temperature at 500mb. For 2-meter temp, it'll use the 2
>> thresholds
>>>>>>> listed (ge273 ge283) and for temp at 500mb, it'll use
>>>>>>> the 6 thresholds listed (ge200.0 ge220 ge240 ge260 ge280
ge300).
>>>>>>>
>>>>>>> So the real task is running your WRF output files through WPP.
>>>>>>>
>>>>>>> Lastly, I wanted to tell you that in METv3.0 (to be released
>> shortly),
>>>>>>> we've added support for reading the NetCDF output of the
pinterp tool
>>>>> (which
>>>>>>> only processes ARW output). So once METv3.0 is
>>>>>>> available, you can use either the GRIB output of WPP or the
NetCDF
>>>>> output
>>>>>>> of the pinterp tool.
>>>>>>>
>>>>>>> Hope that helps.
>>>>>>>
>>>>>>> Thanks,
>>>>>>> John Halley Gotway
>>>>>>> met_help at ucar.edu
>>>>>>>
>>>>>>>
>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804
>
>>>>>>>>
>>>>>>>> Dear Sir,
>>>>>>>>
>>>>>>>> I have put my files at the FTP site told by you. Please
check.
>>>>>>>>
>>>>>>>> One more thing, i want to know that do i need to give equal
number
>> of
>>>>>>>> threshold values as fcst_field/obs_field? Suppose I there is
only
>> one
>>>>>>>> fcst_field/ob_field, but I want to check it for various
threshold
>>>>> values.
>>>>>>>> In some other case, I have given 5 fcst_field/obs_field and
>>>>> corresponding
>>>>>>> 5
>>>>>>>> threshold (fcst_thresh/obs_thresh) values. But I did not get
>> anything
>>>>> in
>>>>>>>> the output. I have got netcdf file in the output but nothing
in the
>>>>> STAT
>>>>>>> and
>>>>>>>> .txt files.
>>>>>>>>
>>>>>>>> Thanks and Regards
>>>>>>>> Hashmi Fatima
>>>>>>>>
>>>>>>>> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for Paul
Oldenburg} <
>>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>>
>>>>>>>>> Hashmi,
>>>>>>>>>
>>>>>>>>> Can you please send us an example of your model and obs data
along
>>>>> with
>>>>>>>>> your grid_stat config file? If the data is too
>>>>>>>>> large to attach to an email, please upload it to our FTP
site:
>>>>>>>>>
>>>>>>>>> ftp ftp.rap.ucar.edu
>>>>>>>>> username = anonymous
>>>>>>>>> password = "your email address"
>>>>>>>>> cd incoming/irap/met_help/fatima_data
>>>>>>>>> put "your files"
>>>>>>>>>
>>>>>>>>> Thanks,
>>>>>>>>>
>>>>>>>>> Paul
>>>>>>>>>
>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted upon.
>>>>>>>>>> Transaction: Ticket created by hashmifatima at gmail.com
>>>>>>>>>> Queue: met_help
>>>>>>>>>> Subject: netcdf file for grid_stat (grid code error)
>>>>>>>>>> Owner: Nobody
>>>>>>>>>> Requestors: hashmifatima at gmail.com
>>>>>>>>>> Status: new
>>>>>>>>>> Ticket <URL:
>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
>>>>>>>>>> Dear Sir,
>>>>>>>>>>
>>>>>>>>>> I want to use grid_stat tool to verify 2m temperature. Both
>>>>>>> (observation
>>>>>>>>> and
>>>>>>>>>> forecast) files are in netcdf format. when i am defining
>> fcst_field,
>>>>> it
>>>>>>>>> is
>>>>>>>>>> asking for grid code. I am not using any grib file for this
test.
>>>>> What
>>>>>>>>>> should I do ?
>>>>>>>>>>
>>>>>>>>>> Thanks and Regards,
>>>>>>
>>
------------------------------------------------------------------------
>>>
>>
>>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #40804] netcdf file for grid_stat (grid code error)
From: Paul Oldenburg
Time: Wed Sep 29 09:19:13 2010
Hashmi,
Here is the website, sorry:
http://www.dtcenter.org/met/users/support/met_help.php#ftp
If you cannot get this to work, please let us know and we will have to
make other arrangements for you to send us your
data. Sorry for all the trouble!
Thanks,
Paul
RAL HelpDesk {for Paul Oldenburg} wrote:
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>
> Hashmi,
>
> We apparently have a new FTP server configuration that we are still
getting used to. Can you please refer to the
> instructions on the following website for uploading data and let us
know if you can get it to work:
>
>
>
> RAL HelpDesk {for Hashmi Fatima} wrote:
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>
>> Dear Paul
>>
>> ftp -p ftp.rap.ucar.edu is not working. It is giving me error
message.
>>
>> ftp -p ftp.rap.ucar.edu
>> ftp: p: unknown option
>> Usage: ftp [-v] [-d] [-i] [-n] [-g] [-k realm] [-f] [-x] [-u] [-t]
[host]
>>
>> Regards
>> Hashmi Fatima
>>
>>
>>
>> On Wed, Sep 29, 2010 at 12:22 AM, RAL HelpDesk {for Paul Oldenburg}
<
>> met_help at ucar.edu> wrote:
>>
>>> Hashmi,
>>>
>>> Please try to invoke your ftp client using the -p command line
argument,
>>> like so:
>>>
>>> ftp -p ftp.rap.ucar.edu
>>> username = anonymous
>>> password = "your email address"
>>> cd incoming/irap/met_help/fatima_data
>>> put "your files"
>>>
>>> Let us know how this approach goes.
>>>
>>> Thanks,
>>>
>>> Paul
>>>
>>>
>>> Paul Oldenburg wrote:
>>>> Hashmi,
>>>>
>>>> We are investigating a potential problem with our FTP site. I
will let
>>> you know when it is fixed, and then we can
>>>> proceed with addressing your wavelet_stat question.
>>>>
>>>> Thanks,
>>>>
>>>> Paul
>>>>
>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>
>>>>> Dear Paul
>>>>>
>>>>> I tried to upload my files at the FTP server you have told me in
>>> previous
>>>>> mails. But it says " Permission denied on server".
>>>>>
>>>>> Regards
>>>>> Hashmi Fatima
>>>>>
>>>>> On Mon, Sep 27, 2010 at 10:25 PM, RAL HelpDesk {for Paul
Oldenburg} <
>>>>> met_help at ucar.edu> wrote:
>>>>>
>>>>>> Hashmi,
>>>>>>
>>>>>> I (Paul) have returned to my met_help responsibilities and I
have
>>> followed
>>>>>> your interactions with John. In order to
>>>>>> help you run stat_analysis, I would like to have you send me
some of
>>> your
>>>>>> data via FTP. First, can you give us the WPP
>>>>>> output that you fed into wavelet_stat? Also, can you give us
your
>>>>>> wavelet_stat output? Finally, please copy the error
>>>>>> message that you get from stat_analysis in your response. That
>>> information
>>>>>> will be very helpful.
>>>>>>
>>>>>> Concerning the masking region selected by wavelet_stat, the
entire grid
>>>>>> cannot be processed because wavelet_stat can
>>>>>> process only dyadic regions - that is square regions whose
sides are a
>>>>>> power of 2 grid points (like 256 or 2048). By
>>>>>> default, wavelet_stat should try to guess the best regions to
use in
>>> your
>>>>>> domain. If you send us your WPP output, we
>>>>>> will be able to help you configure the best masking region for
your
>>> case.
>>>>>> We look forward to receiving your data. If you have any
questions,
>>> please
>>>>>> let me know.
>>>>>>
>>>>>> Thanks,
>>>>>>
>>>>>> Paul
>>>>>>
>>>>>>
>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804
>
>>>>>>>
>>>>>>> Dear Sir,
>>>>>>>
>>>>>>> Thanks for your help.Its working now.
>>>>>>>
>>>>>>> I have run wavelet tool for computation of RH. Now i want to
use the
>>>>>> output
>>>>>>> of wavelet for StatAnalysis. but its not working. I am sending
you
>>>>>>> configuration file of wavelet and StatAnalysis. Could you
please
>>> check.
>>>>>>> One more thing I want to know that in the wavelet output, why
does it
>>>>>> take
>>>>>>> that much of area. I have not assigned any area for the
verification
>>>>>> except
>>>>>>> full. I have masked the whole grid. I am attaching that output
too as
>>>>>>> Wavelet_output.png.
>>>>>>>
>>>>>>> Thanks and Regards
>>>>>>> Hashmi fatima
>>>>>>>
>>>>>>>
>>>>>>> On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for John Halley
Gotway}
>>> <
>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>
>>>>>>>> Hashmi,
>>>>>>>>
>>>>>>>> Paul will be out of the office for the next couple of weeks,
but I'll
>>> be
>>>>>>>> happy to help you with this issue.
>>>>>>>>
>>>>>>>> There are a couple of points I'd like to make:
>>>>>>>>
>>>>>>>> First, the NetCDF files you sent us are the NetCDF output of
WRF.
>>> MET
>>>>>> does
>>>>>>>> not directly read the NetCDF output of WRF. Instead, for
METv2.0,
>>> your
>>>>>> WRF
>>>>>>>> output must first be run through the
>>>>>>>> WRF-PostProcessor (WPP). Details of running WPP can be found
in
>>> chapter
>>>>>> 7
>>>>>>>> of the WRF user's guide:
>>>>>>>>
>>>>>>>>
>>> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
>>>>>>>> WPP writes out GRIB data, which MET can read.
>>>>>>>>
>>>>>>>> Second, once you've run your data through the WRF
PostProcessor, you
>>> can
>>>>>>>> select multiple thresholds for Grid-Stat to use. You do so
by
>>> setting
>>>>>> up
>>>>>>>> the Grid-Stat configuration file something like this:
>>>>>>>>
>>>>>>>> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
>>>>>>>> obs_field[] = [];
>>>>>>>> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240 ge260
ge280
>>> ge300"
>>>>>> ];
>>>>>>>> obs_thresh[] = [];
>>>>>>>>
>>>>>>>> The settings above tell Grid-Stat to verify temperature at 2-
meters
>>> as
>>>>>> well
>>>>>>>> as temperature at 500mb. For 2-meter temp, it'll use the 2
>>> thresholds
>>>>>>>> listed (ge273 ge283) and for temp at 500mb, it'll use
>>>>>>>> the 6 thresholds listed (ge200.0 ge220 ge240 ge260 ge280
ge300).
>>>>>>>>
>>>>>>>> So the real task is running your WRF output files through
WPP.
>>>>>>>>
>>>>>>>> Lastly, I wanted to tell you that in METv3.0 (to be released
>>> shortly),
>>>>>>>> we've added support for reading the NetCDF output of the
pinterp tool
>>>>>> (which
>>>>>>>> only processes ARW output). So once METv3.0 is
>>>>>>>> available, you can use either the GRIB output of WPP or the
NetCDF
>>>>>> output
>>>>>>>> of the pinterp tool.
>>>>>>>>
>>>>>>>> Hope that helps.
>>>>>>>>
>>>>>>>> Thanks,
>>>>>>>> John Halley Gotway
>>>>>>>> met_help at ucar.edu
>>>>>>>>
>>>>>>>>
>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>>>>>
>>>>>>>>> Dear Sir,
>>>>>>>>>
>>>>>>>>> I have put my files at the FTP site told by you. Please
check.
>>>>>>>>>
>>>>>>>>> One more thing, i want to know that do i need to give equal
number
>>> of
>>>>>>>>> threshold values as fcst_field/obs_field? Suppose I there is
only
>>> one
>>>>>>>>> fcst_field/ob_field, but I want to check it for various
threshold
>>>>>> values.
>>>>>>>>> In some other case, I have given 5 fcst_field/obs_field and
>>>>>> corresponding
>>>>>>>> 5
>>>>>>>>> threshold (fcst_thresh/obs_thresh) values. But I did not
get
>>> anything
>>>>>> in
>>>>>>>>> the output. I have got netcdf file in the output but nothing
in the
>>>>>> STAT
>>>>>>>> and
>>>>>>>>> .txt files.
>>>>>>>>>
>>>>>>>>> Thanks and Regards
>>>>>>>>> Hashmi Fatima
>>>>>>>>>
>>>>>>>>> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for Paul
Oldenburg} <
>>>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>>>
>>>>>>>>>> Hashmi,
>>>>>>>>>>
>>>>>>>>>> Can you please send us an example of your model and obs
data along
>>>>>> with
>>>>>>>>>> your grid_stat config file? If the data is too
>>>>>>>>>> large to attach to an email, please upload it to our FTP
site:
>>>>>>>>>>
>>>>>>>>>> ftp ftp.rap.ucar.edu
>>>>>>>>>> username = anonymous
>>>>>>>>>> password = "your email address"
>>>>>>>>>> cd incoming/irap/met_help/fatima_data
>>>>>>>>>> put "your files"
>>>>>>>>>>
>>>>>>>>>> Thanks,
>>>>>>>>>>
>>>>>>>>>> Paul
>>>>>>>>>>
>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted upon.
>>>>>>>>>>> Transaction: Ticket created by hashmifatima at gmail.com
>>>>>>>>>>> Queue: met_help
>>>>>>>>>>> Subject: netcdf file for grid_stat (grid code error)
>>>>>>>>>>> Owner: Nobody
>>>>>>>>>>> Requestors: hashmifatima at gmail.com
>>>>>>>>>>> Status: new
>>>>>>>>>>> Ticket <URL:
>>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>
>>>>>>>>>>> I want to use grid_stat tool to verify 2m temperature.
Both
>>>>>>>> (observation
>>>>>>>>>> and
>>>>>>>>>>> forecast) files are in netcdf format. when i am defining
>>> fcst_field,
>>>>>> it
>>>>>>>>>> is
>>>>>>>>>>> asking for grid code. I am not using any grib file for
this test.
>>>>>> What
>>>>>>>>>>> should I do ?
>>>>>>>>>>>
>>>>>>>>>>> Thanks and Regards,
>>>
------------------------------------------------------------------------
>>>
>>
>
------------------------------------------------
Subject: netcdf file for grid_stat (grid code error)
From: Hashmi Fatima
Time: Thu Sep 30 23:54:02 2010
Dear Paul,
Sorry for the delay. I have put my files at ftp site in my folder. I
have
done the verufication through wavelet for 1 month data. Here I am
sending
few files.
*input files (WPP output):*
gfs.t00z.pgrbf00 (analysis file)
gfs.t00z.pgrbf24 (forecast file)
*wavelet output file (stat analysis input) :*
wavelet_stat_000000L_20100920_000000V.nc
wavelet_stat_000000L_20100920_000000V.ps
wavelet_stat_000000L_20100920_000000V.stat
wavelet_stat_000000L_20100920_000000V_isc.txt
The error message which I am getting after running stat_analysis tool
after
using wavelet output is in file *error_msg.txt*.
Thanks and Regards
Hashmi Fatima
On Wed, Sep 29, 2010 at 8:49 PM, RAL HelpDesk {for Paul Oldenburg} <
met_help at ucar.edu> wrote:
> Hashmi,
>
> Here is the website, sorry:
>
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> If you cannot get this to work, please let us know and we will have
to make
> other arrangements for you to send us your
> data. Sorry for all the trouble!
>
> Thanks,
>
> Paul
>
>
> RAL HelpDesk {for Paul Oldenburg} wrote:
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >
> > Hashmi,
> >
> > We apparently have a new FTP server configuration that we are
still
> getting used to. Can you please refer to the
> > instructions on the following website for uploading data and let
us know
> if you can get it to work:
> >
> >
> >
> > RAL HelpDesk {for Hashmi Fatima} wrote:
> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >>
> >> Dear Paul
> >>
> >> ftp -p ftp.rap.ucar.edu is not working. It is giving me error
message.
> >>
> >> ftp -p ftp.rap.ucar.edu
> >> ftp: p: unknown option
> >> Usage: ftp [-v] [-d] [-i] [-n] [-g] [-k realm] [-f] [-x] [-u] [-
t]
> [host]
> >>
> >> Regards
> >> Hashmi Fatima
> >>
> >>
> >>
> >> On Wed, Sep 29, 2010 at 12:22 AM, RAL HelpDesk {for Paul
Oldenburg} <
> >> met_help at ucar.edu> wrote:
> >>
> >>> Hashmi,
> >>>
> >>> Please try to invoke your ftp client using the -p command line
> argument,
> >>> like so:
> >>>
> >>> ftp -p ftp.rap.ucar.edu
> >>> username = anonymous
> >>> password = "your email address"
> >>> cd incoming/irap/met_help/fatima_data
> >>> put "your files"
> >>>
> >>> Let us know how this approach goes.
> >>>
> >>> Thanks,
> >>>
> >>> Paul
> >>>
> >>>
> >>> Paul Oldenburg wrote:
> >>>> Hashmi,
> >>>>
> >>>> We are investigating a potential problem with our FTP site. I
will
> let
> >>> you know when it is fixed, and then we can
> >>>> proceed with addressing your wavelet_stat question.
> >>>>
> >>>> Thanks,
> >>>>
> >>>> Paul
> >>>>
> >>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804
>
> >>>>>
> >>>>> Dear Paul
> >>>>>
> >>>>> I tried to upload my files at the FTP server you have told me
in
> >>> previous
> >>>>> mails. But it says " Permission denied on server".
> >>>>>
> >>>>> Regards
> >>>>> Hashmi Fatima
> >>>>>
> >>>>> On Mon, Sep 27, 2010 at 10:25 PM, RAL HelpDesk {for Paul
Oldenburg} <
> >>>>> met_help at ucar.edu> wrote:
> >>>>>
> >>>>>> Hashmi,
> >>>>>>
> >>>>>> I (Paul) have returned to my met_help responsibilities and I
have
> >>> followed
> >>>>>> your interactions with John. In order to
> >>>>>> help you run stat_analysis, I would like to have you send me
some of
> >>> your
> >>>>>> data via FTP. First, can you give us the WPP
> >>>>>> output that you fed into wavelet_stat? Also, can you give us
your
> >>>>>> wavelet_stat output? Finally, please copy the error
> >>>>>> message that you get from stat_analysis in your response.
That
> >>> information
> >>>>>> will be very helpful.
> >>>>>>
> >>>>>> Concerning the masking region selected by wavelet_stat, the
entire
> grid
> >>>>>> cannot be processed because wavelet_stat can
> >>>>>> process only dyadic regions - that is square regions whose
sides are
> a
> >>>>>> power of 2 grid points (like 256 or 2048). By
> >>>>>> default, wavelet_stat should try to guess the best regions to
use in
> >>> your
> >>>>>> domain. If you send us your WPP output, we
> >>>>>> will be able to help you configure the best masking region
for your
> >>> case.
> >>>>>> We look forward to receiving your data. If you have any
questions,
> >>> please
> >>>>>> let me know.
> >>>>>>
> >>>>>> Thanks,
> >>>>>>
> >>>>>> Paul
> >>>>>>
> >>>>>>
> >>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >>>>>>>
> >>>>>>> Dear Sir,
> >>>>>>>
> >>>>>>> Thanks for your help.Its working now.
> >>>>>>>
> >>>>>>> I have run wavelet tool for computation of RH. Now i want to
use
> the
> >>>>>> output
> >>>>>>> of wavelet for StatAnalysis. but its not working. I am
sending you
> >>>>>>> configuration file of wavelet and StatAnalysis. Could you
please
> >>> check.
> >>>>>>> One more thing I want to know that in the wavelet output,
why does
> it
> >>>>>> take
> >>>>>>> that much of area. I have not assigned any area for the
> verification
> >>>>>> except
> >>>>>>> full. I have masked the whole grid. I am attaching that
output too
> as
> >>>>>>> Wavelet_output.png.
> >>>>>>>
> >>>>>>> Thanks and Regards
> >>>>>>> Hashmi fatima
> >>>>>>>
> >>>>>>>
> >>>>>>> On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for John
Halley
> Gotway}
> >>> <
> >>>>>>> met_help at ucar.edu> wrote:
> >>>>>>>
> >>>>>>>> Hashmi,
> >>>>>>>>
> >>>>>>>> Paul will be out of the office for the next couple of
weeks, but
> I'll
> >>> be
> >>>>>>>> happy to help you with this issue.
> >>>>>>>>
> >>>>>>>> There are a couple of points I'd like to make:
> >>>>>>>>
> >>>>>>>> First, the NetCDF files you sent us are the NetCDF output
of WRF.
> >>> MET
> >>>>>> does
> >>>>>>>> not directly read the NetCDF output of WRF. Instead, for
METv2.0,
> >>> your
> >>>>>> WRF
> >>>>>>>> output must first be run through the
> >>>>>>>> WRF-PostProcessor (WPP). Details of running WPP can be
found in
> >>> chapter
> >>>>>> 7
> >>>>>>>> of the WRF user's guide:
> >>>>>>>>
> >>>>>>>>
> >>>
> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
> >>>>>>>> WPP writes out GRIB data, which MET can read.
> >>>>>>>>
> >>>>>>>> Second, once you've run your data through the WRF
PostProcessor,
> you
> >>> can
> >>>>>>>> select multiple thresholds for Grid-Stat to use. You do so
by
> >>> setting
> >>>>>> up
> >>>>>>>> the Grid-Stat configuration file something like this:
> >>>>>>>>
> >>>>>>>> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
> >>>>>>>> obs_field[] = [];
> >>>>>>>> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240 ge260
ge280
> >>> ge300"
> >>>>>> ];
> >>>>>>>> obs_thresh[] = [];
> >>>>>>>>
> >>>>>>>> The settings above tell Grid-Stat to verify temperature at
> 2-meters
> >>> as
> >>>>>> well
> >>>>>>>> as temperature at 500mb. For 2-meter temp, it'll use the 2
> >>> thresholds
> >>>>>>>> listed (ge273 ge283) and for temp at 500mb, it'll use
> >>>>>>>> the 6 thresholds listed (ge200.0 ge220 ge240 ge260 ge280
ge300).
> >>>>>>>>
> >>>>>>>> So the real task is running your WRF output files through
WPP.
> >>>>>>>>
> >>>>>>>> Lastly, I wanted to tell you that in METv3.0 (to be
released
> >>> shortly),
> >>>>>>>> we've added support for reading the NetCDF output of the
pinterp
> tool
> >>>>>> (which
> >>>>>>>> only processes ARW output). So once METv3.0 is
> >>>>>>>> available, you can use either the GRIB output of WPP or the
NetCDF
> >>>>>> output
> >>>>>>>> of the pinterp tool.
> >>>>>>>>
> >>>>>>>> Hope that helps.
> >>>>>>>>
> >>>>>>>> Thanks,
> >>>>>>>> John Halley Gotway
> >>>>>>>> met_help at ucar.edu
> >>>>>>>>
> >>>>>>>>
> >>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >>>>>>>>>
> >>>>>>>>> Dear Sir,
> >>>>>>>>>
> >>>>>>>>> I have put my files at the FTP site told by you. Please
check.
> >>>>>>>>>
> >>>>>>>>> One more thing, i want to know that do i need to give
equal
> number
> >>> of
> >>>>>>>>> threshold values as fcst_field/obs_field? Suppose I there
is only
> >>> one
> >>>>>>>>> fcst_field/ob_field, but I want to check it for various
threshold
> >>>>>> values.
> >>>>>>>>> In some other case, I have given 5 fcst_field/obs_field
and
> >>>>>> corresponding
> >>>>>>>> 5
> >>>>>>>>> threshold (fcst_thresh/obs_thresh) values. But I did not
get
> >>> anything
> >>>>>> in
> >>>>>>>>> the output. I have got netcdf file in the output but
nothing in
> the
> >>>>>> STAT
> >>>>>>>> and
> >>>>>>>>> .txt files.
> >>>>>>>>>
> >>>>>>>>> Thanks and Regards
> >>>>>>>>> Hashmi Fatima
> >>>>>>>>>
> >>>>>>>>> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for Paul
> Oldenburg} <
> >>>>>>>>> met_help at ucar.edu> wrote:
> >>>>>>>>>
> >>>>>>>>>> Hashmi,
> >>>>>>>>>>
> >>>>>>>>>> Can you please send us an example of your model and obs
data
> along
> >>>>>> with
> >>>>>>>>>> your grid_stat config file? If the data is too
> >>>>>>>>>> large to attach to an email, please upload it to our FTP
site:
> >>>>>>>>>>
> >>>>>>>>>> ftp ftp.rap.ucar.edu
> >>>>>>>>>> username = anonymous
> >>>>>>>>>> password = "your email address"
> >>>>>>>>>> cd incoming/irap/met_help/fatima_data
> >>>>>>>>>> put "your files"
> >>>>>>>>>>
> >>>>>>>>>> Thanks,
> >>>>>>>>>>
> >>>>>>>>>> Paul
> >>>>>>>>>>
> >>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>>>>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted upon.
> >>>>>>>>>>> Transaction: Ticket created by hashmifatima at gmail.com
> >>>>>>>>>>> Queue: met_help
> >>>>>>>>>>> Subject: netcdf file for grid_stat (grid code
error)
> >>>>>>>>>>> Owner: Nobody
> >>>>>>>>>>> Requestors: hashmifatima at gmail.com
> >>>>>>>>>>> Status: new
> >>>>>>>>>>> Ticket <URL:
> >>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
> >>>>>>>>>>> Dear Sir,
> >>>>>>>>>>>
> >>>>>>>>>>> I want to use grid_stat tool to verify 2m temperature.
Both
> >>>>>>>> (observation
> >>>>>>>>>> and
> >>>>>>>>>>> forecast) files are in netcdf format. when i am defining
> >>> fcst_field,
> >>>>>> it
> >>>>>>>>>> is
> >>>>>>>>>>> asking for grid code. I am not using any grib file for
this
> test.
> >>>>>> What
> >>>>>>>>>>> should I do ?
> >>>>>>>>>>>
> >>>>>>>>>>> Thanks and Regards,
> >>>
>
------------------------------------------------------------------------
> >>>
> >>
> >
>
>
>
--
Hashmi Fatima
Research Scholar
Centre for Atmospheric Sciences
Indian Institute of Technology Delhi
Hauz Khas
New Delhi - 110016
INDIA
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #40804] netcdf file for grid_stat (grid code error)
From: John Halley Gotway
Time: Fri Oct 01 15:19:08 2010
Hashmi,
Paul's out of the office, so hopefully I can provide some advice on
this issue.
It looks like you're using the STATAnalysis config file that's
included with the MET test scripts. This is really only meant as an
example to work with the test data. You'll need to define the
analysis jobs explicitly to match the output data you have.
Since you're running Wavelet-Stat, I assume you may want to aggregate
the output of Wavelet-Stat across multiple cases. Please keep in mind
that Stat-Analysis can be run on the command line OR using
a config file. When you're trying to set up your analysis job at the
beginning, I'd suggest using the command line option. For example,
here's the output I get for the following command:
+++++++++++++++++++++++++++++++++++++++
bin/stat_analysis -lookin wavelet_stat_000000L_20100920_000000V.stat
-job aggregate -line_type ISC -vx_mask TILE_TOT -fcst_thresh
">=50.000" -dump_row aggr_ISC_50.stat -v 2
STAT Lines read = 330
STAT Lines retained = 330
Processing Job 1: -job aggregate -vx_mask TILE_TOT -fcst_thresh
>=50.000 -line_type ISC -dump_row aggr_ISC_50.stat
Creating STAT output file "aggr_ISC_50.stat"
JOB_LIST: -job aggregate -vx_mask TILE_TOT -fcst_thresh >=50.000
-line_type ISC -dump_row aggr_ISC_50.stat
COL_NAME: TOTAL TILE_DIM TILE_XLL TILE_YLL NSCALE ISCALE MSE ISC
FENERGY2 OENERGY2 BASER FBIAS
ISC: 131072 NA NA NA 8 0 0.09870
0.63219 0.16414 0.15513 0.15513 1.05808
131072 NA NA NA 8 1 0.01761
0.40938 0.00935 0.00921 0.15513 1.05808
131072 NA NA NA 8 2 0.01881
0.36913 0.01064 0.01070 0.15513 1.05808
131072 NA NA NA 8 3 0.02372
0.20445 0.01519 0.01605 0.15513 1.05808
131072 NA NA NA 8 4 0.01836
0.38422 0.01856 0.01899 0.15513 1.05808
131072 NA NA NA 8 5 0.01390
0.53381 0.01933 0.01968 0.15513 1.05808
131072 NA NA NA 8 6 0.00524
0.82426 0.02334 0.02096 0.15513 1.05808
131072 NA NA NA 8 7 0.00082
0.97250 0.02947 0.02577 0.15513 1.05808
131072 NA NA NA 8 8 0.00011
0.99631 0.00854 0.00764 0.15513 1.05808
131072 NA NA NA 8 9 0.00013
0.99564 0.02972 0.02614 0.15513 1.05808
Job 1 used 10 out of 330 STAT lines.
+++++++++++++++++++++++++++++++++++++++
Since you just sent us one Wavelet-Stat output file, that's all we
could test with. But in this example, I'm aggregating together output
for the ">=50.000" threshold. I'm guessing you'll want to
rerun this type of job for each of the thresholds you used in Wavelet-
Stat. When you're first setting up jobs, it's a good idea to use the
"-dump_row" option. That'll write all of the STAT lines the
job actually used to the file name you specify. That way, you can
make sure that STAT-Analysis is operating on the exact subset of data
that you expect.
Hope that helps clarify.
John Halley Gotway
met_help at ucar.edu
RAL HelpDesk {for Hashmi Fatima} wrote:
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>
> Dear Paul,
>
> Sorry for the delay. I have put my files at ftp site in my folder. I
have
> done the verufication through wavelet for 1 month data. Here I am
sending
> few files.
>
> *input files (WPP output):*
> gfs.t00z.pgrbf00 (analysis file)
> gfs.t00z.pgrbf24 (forecast file)
>
> *wavelet output file (stat analysis input) :*
> wavelet_stat_000000L_20100920_000000V.nc
> wavelet_stat_000000L_20100920_000000V.ps
> wavelet_stat_000000L_20100920_000000V.stat
> wavelet_stat_000000L_20100920_000000V_isc.txt
>
> The error message which I am getting after running stat_analysis
tool after
> using wavelet output is in file *error_msg.txt*.
>
> Thanks and Regards
> Hashmi Fatima
>
> On Wed, Sep 29, 2010 at 8:49 PM, RAL HelpDesk {for Paul Oldenburg} <
> met_help at ucar.edu> wrote:
>
>> Hashmi,
>>
>> Here is the website, sorry:
>>
>> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>>
>> If you cannot get this to work, please let us know and we will have
to make
>> other arrangements for you to send us your
>> data. Sorry for all the trouble!
>>
>> Thanks,
>>
>> Paul
>>
>>
>> RAL HelpDesk {for Paul Oldenburg} wrote:
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>
>>> Hashmi,
>>>
>>> We apparently have a new FTP server configuration that we are
still
>> getting used to. Can you please refer to the
>>> instructions on the following website for uploading data and let
us know
>> if you can get it to work:
>>>
>>>
>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>
>>>> Dear Paul
>>>>
>>>> ftp -p ftp.rap.ucar.edu is not working. It is giving me error
message.
>>>>
>>>> ftp -p ftp.rap.ucar.edu
>>>> ftp: p: unknown option
>>>> Usage: ftp [-v] [-d] [-i] [-n] [-g] [-k realm] [-f] [-x] [-u] [-
t]
>> [host]
>>>> Regards
>>>> Hashmi Fatima
>>>>
>>>>
>>>>
>>>> On Wed, Sep 29, 2010 at 12:22 AM, RAL HelpDesk {for Paul
Oldenburg} <
>>>> met_help at ucar.edu> wrote:
>>>>
>>>>> Hashmi,
>>>>>
>>>>> Please try to invoke your ftp client using the -p command line
>> argument,
>>>>> like so:
>>>>>
>>>>> ftp -p ftp.rap.ucar.edu
>>>>> username = anonymous
>>>>> password = "your email address"
>>>>> cd incoming/irap/met_help/fatima_data
>>>>> put "your files"
>>>>>
>>>>> Let us know how this approach goes.
>>>>>
>>>>> Thanks,
>>>>>
>>>>> Paul
>>>>>
>>>>>
>>>>> Paul Oldenburg wrote:
>>>>>> Hashmi,
>>>>>>
>>>>>> We are investigating a potential problem with our FTP site. I
will
>> let
>>>>> you know when it is fixed, and then we can
>>>>>> proceed with addressing your wavelet_stat question.
>>>>>>
>>>>>> Thanks,
>>>>>>
>>>>>> Paul
>>>>>>
>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804
>
>>>>>>>
>>>>>>> Dear Paul
>>>>>>>
>>>>>>> I tried to upload my files at the FTP server you have told me
in
>>>>> previous
>>>>>>> mails. But it says " Permission denied on server".
>>>>>>>
>>>>>>> Regards
>>>>>>> Hashmi Fatima
>>>>>>>
>>>>>>> On Mon, Sep 27, 2010 at 10:25 PM, RAL HelpDesk {for Paul
Oldenburg} <
>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>
>>>>>>>> Hashmi,
>>>>>>>>
>>>>>>>> I (Paul) have returned to my met_help responsibilities and I
have
>>>>> followed
>>>>>>>> your interactions with John. In order to
>>>>>>>> help you run stat_analysis, I would like to have you send me
some of
>>>>> your
>>>>>>>> data via FTP. First, can you give us the WPP
>>>>>>>> output that you fed into wavelet_stat? Also, can you give us
your
>>>>>>>> wavelet_stat output? Finally, please copy the error
>>>>>>>> message that you get from stat_analysis in your response.
That
>>>>> information
>>>>>>>> will be very helpful.
>>>>>>>>
>>>>>>>> Concerning the masking region selected by wavelet_stat, the
entire
>> grid
>>>>>>>> cannot be processed because wavelet_stat can
>>>>>>>> process only dyadic regions - that is square regions whose
sides are
>> a
>>>>>>>> power of 2 grid points (like 256 or 2048). By
>>>>>>>> default, wavelet_stat should try to guess the best regions to
use in
>>>>> your
>>>>>>>> domain. If you send us your WPP output, we
>>>>>>>> will be able to help you configure the best masking region
for your
>>>>> case.
>>>>>>>> We look forward to receiving your data. If you have any
questions,
>>>>> please
>>>>>>>> let me know.
>>>>>>>>
>>>>>>>> Thanks,
>>>>>>>>
>>>>>>>> Paul
>>>>>>>>
>>>>>>>>
>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>>>>>
>>>>>>>>> Dear Sir,
>>>>>>>>>
>>>>>>>>> Thanks for your help.Its working now.
>>>>>>>>>
>>>>>>>>> I have run wavelet tool for computation of RH. Now i want to
use
>> the
>>>>>>>> output
>>>>>>>>> of wavelet for StatAnalysis. but its not working. I am
sending you
>>>>>>>>> configuration file of wavelet and StatAnalysis. Could you
please
>>>>> check.
>>>>>>>>> One more thing I want to know that in the wavelet output,
why does
>> it
>>>>>>>> take
>>>>>>>>> that much of area. I have not assigned any area for the
>> verification
>>>>>>>> except
>>>>>>>>> full. I have masked the whole grid. I am attaching that
output too
>> as
>>>>>>>>> Wavelet_output.png.
>>>>>>>>>
>>>>>>>>> Thanks and Regards
>>>>>>>>> Hashmi fatima
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for John
Halley
>> Gotway}
>>>>> <
>>>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>>>
>>>>>>>>>> Hashmi,
>>>>>>>>>>
>>>>>>>>>> Paul will be out of the office for the next couple of
weeks, but
>> I'll
>>>>> be
>>>>>>>>>> happy to help you with this issue.
>>>>>>>>>>
>>>>>>>>>> There are a couple of points I'd like to make:
>>>>>>>>>>
>>>>>>>>>> First, the NetCDF files you sent us are the NetCDF output
of WRF.
>>>>> MET
>>>>>>>> does
>>>>>>>>>> not directly read the NetCDF output of WRF. Instead, for
METv2.0,
>>>>> your
>>>>>>>> WRF
>>>>>>>>>> output must first be run through the
>>>>>>>>>> WRF-PostProcessor (WPP). Details of running WPP can be
found in
>>>>> chapter
>>>>>>>> 7
>>>>>>>>>> of the WRF user's guide:
>>>>>>>>>>
>>>>>>>>>>
>> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
>>>>>>>>>> WPP writes out GRIB data, which MET can read.
>>>>>>>>>>
>>>>>>>>>> Second, once you've run your data through the WRF
PostProcessor,
>> you
>>>>> can
>>>>>>>>>> select multiple thresholds for Grid-Stat to use. You do so
by
>>>>> setting
>>>>>>>> up
>>>>>>>>>> the Grid-Stat configuration file something like this:
>>>>>>>>>>
>>>>>>>>>> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
>>>>>>>>>> obs_field[] = [];
>>>>>>>>>> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240 ge260
ge280
>>>>> ge300"
>>>>>>>> ];
>>>>>>>>>> obs_thresh[] = [];
>>>>>>>>>>
>>>>>>>>>> The settings above tell Grid-Stat to verify temperature at
>> 2-meters
>>>>> as
>>>>>>>> well
>>>>>>>>>> as temperature at 500mb. For 2-meter temp, it'll use the 2
>>>>> thresholds
>>>>>>>>>> listed (ge273 ge283) and for temp at 500mb, it'll use
>>>>>>>>>> the 6 thresholds listed (ge200.0 ge220 ge240 ge260 ge280
ge300).
>>>>>>>>>>
>>>>>>>>>> So the real task is running your WRF output files through
WPP.
>>>>>>>>>>
>>>>>>>>>> Lastly, I wanted to tell you that in METv3.0 (to be
released
>>>>> shortly),
>>>>>>>>>> we've added support for reading the NetCDF output of the
pinterp
>> tool
>>>>>>>> (which
>>>>>>>>>> only processes ARW output). So once METv3.0 is
>>>>>>>>>> available, you can use either the GRIB output of WPP or the
NetCDF
>>>>>>>> output
>>>>>>>>>> of the pinterp tool.
>>>>>>>>>>
>>>>>>>>>> Hope that helps.
>>>>>>>>>>
>>>>>>>>>> Thanks,
>>>>>>>>>> John Halley Gotway
>>>>>>>>>> met_help at ucar.edu
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>>>>>>>
>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>
>>>>>>>>>>> I have put my files at the FTP site told by you. Please
check.
>>>>>>>>>>>
>>>>>>>>>>> One more thing, i want to know that do i need to give
equal
>> number
>>>>> of
>>>>>>>>>>> threshold values as fcst_field/obs_field? Suppose I there
is only
>>>>> one
>>>>>>>>>>> fcst_field/ob_field, but I want to check it for various
threshold
>>>>>>>> values.
>>>>>>>>>>> In some other case, I have given 5 fcst_field/obs_field
and
>>>>>>>> corresponding
>>>>>>>>>> 5
>>>>>>>>>>> threshold (fcst_thresh/obs_thresh) values. But I did not
get
>>>>> anything
>>>>>>>> in
>>>>>>>>>>> the output. I have got netcdf file in the output but
nothing in
>> the
>>>>>>>> STAT
>>>>>>>>>> and
>>>>>>>>>>> .txt files.
>>>>>>>>>>>
>>>>>>>>>>> Thanks and Regards
>>>>>>>>>>> Hashmi Fatima
>>>>>>>>>>>
>>>>>>>>>>> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for Paul
>> Oldenburg} <
>>>>>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>>>>>
>>>>>>>>>>>> Hashmi,
>>>>>>>>>>>>
>>>>>>>>>>>> Can you please send us an example of your model and obs
data
>> along
>>>>>>>> with
>>>>>>>>>>>> your grid_stat config file? If the data is too
>>>>>>>>>>>> large to attach to an email, please upload it to our FTP
site:
>>>>>>>>>>>>
>>>>>>>>>>>> ftp ftp.rap.ucar.edu
>>>>>>>>>>>> username = anonymous
>>>>>>>>>>>> password = "your email address"
>>>>>>>>>>>> cd incoming/irap/met_help/fatima_data
>>>>>>>>>>>> put "your files"
>>>>>>>>>>>>
>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>
>>>>>>>>>>>> Paul
>>>>>>>>>>>>
>>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>>>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted upon.
>>>>>>>>>>>>> Transaction: Ticket created by hashmifatima at gmail.com
>>>>>>>>>>>>> Queue: met_help
>>>>>>>>>>>>> Subject: netcdf file for grid_stat (grid code
error)
>>>>>>>>>>>>> Owner: Nobody
>>>>>>>>>>>>> Requestors: hashmifatima at gmail.com
>>>>>>>>>>>>> Status: new
>>>>>>>>>>>>> Ticket <URL:
>>>>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
>>>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>>>
>>>>>>>>>>>>> I want to use grid_stat tool to verify 2m temperature.
Both
>>>>>>>>>> (observation
>>>>>>>>>>>> and
>>>>>>>>>>>>> forecast) files are in netcdf format. when i am defining
>>>>> fcst_field,
>>>>>>>> it
>>>>>>>>>>>> is
>>>>>>>>>>>>> asking for grid code. I am not using any grib file for
this
>> test.
>>>>>>>> What
>>>>>>>>>>>>> should I do ?
>>>>>>>>>>>>>
>>>>>>>>>>>>> Thanks and Regards,
>>
------------------------------------------------------------------------
>>
>>
>
>
------------------------------------------------
Subject: netcdf file for grid_stat (grid code error)
From: Hashmi Fatima
Time: Tue Oct 05 05:04:48 2010
Dear Sir
I have configured STATconfig file for the previous case.
But now I am trying STATanalysis with grid_stat output for HGT at
P500. But
its again giving some error message. I have uploaded the grid_stat
output
files to your FTP server, which are input for STATanalysis case (all
grid_stat_HGT_P500-240000L_20101001_00000V_* files). STATconfig file
is
STATAnalysisConfig_HGT. I have configurered the STATconfig file
according to
the input. but still unable to get the output.
The error message which I am getting is srored in the file
grid_output_stat_analysis_error.txt.
Please check and let me know about the error.
Thanks and Regards
Hashmi Fatima
On Sat, Oct 2, 2010 at 2:49 AM, RAL HelpDesk {for John Halley Gotway}
<
met_help at ucar.edu> wrote:
> Hashmi,
>
> Paul's out of the office, so hopefully I can provide some advice on
this
> issue.
>
> It looks like you're using the STATAnalysis config file that's
included
> with the MET test scripts. This is really only meant as an example
to work
> with the test data. You'll need to define the
> analysis jobs explicitly to match the output data you have.
>
> Since you're running Wavelet-Stat, I assume you may want to
aggregate the
> output of Wavelet-Stat across multiple cases. Please keep in mind
that
> Stat-Analysis can be run on the command line OR using
> a config file. When you're trying to set up your analysis job at
the
> beginning, I'd suggest using the command line option. For example,
here's
> the output I get for the following command:
>
> +++++++++++++++++++++++++++++++++++++++
> bin/stat_analysis -lookin wavelet_stat_000000L_20100920_000000V.stat
-job
> aggregate -line_type ISC -vx_mask TILE_TOT -fcst_thresh ">=50.000"
-dump_row
> aggr_ISC_50.stat -v 2
> STAT Lines read = 330
> STAT Lines retained = 330
>
> Processing Job 1: -job aggregate -vx_mask TILE_TOT -fcst_thresh
>=50.000
> -line_type ISC -dump_row aggr_ISC_50.stat
> Creating STAT output file "aggr_ISC_50.stat"
> JOB_LIST: -job aggregate -vx_mask TILE_TOT -fcst_thresh
>=50.000
> -line_type ISC -dump_row aggr_ISC_50.stat
> COL_NAME: TOTAL TILE_DIM TILE_XLL TILE_YLL NSCALE ISCALE MSE
ISC
> FENERGY2 OENERGY2 BASER FBIAS
> ISC: 131072 NA NA NA 8 0 0.09870
0.63219
> 0.16414 0.15513 0.15513 1.05808
> 131072 NA NA NA 8 1 0.01761
0.40938
> 0.00935 0.00921 0.15513 1.05808
> 131072 NA NA NA 8 2 0.01881
0.36913
> 0.01064 0.01070 0.15513 1.05808
> 131072 NA NA NA 8 3 0.02372
0.20445
> 0.01519 0.01605 0.15513 1.05808
> 131072 NA NA NA 8 4 0.01836
0.38422
> 0.01856 0.01899 0.15513 1.05808
> 131072 NA NA NA 8 5 0.01390
0.53381
> 0.01933 0.01968 0.15513 1.05808
> 131072 NA NA NA 8 6 0.00524
0.82426
> 0.02334 0.02096 0.15513 1.05808
> 131072 NA NA NA 8 7 0.00082
0.97250
> 0.02947 0.02577 0.15513 1.05808
> 131072 NA NA NA 8 8 0.00011
0.99631
> 0.00854 0.00764 0.15513 1.05808
> 131072 NA NA NA 8 9 0.00013
0.99564
> 0.02972 0.02614 0.15513 1.05808
>
> Job 1 used 10 out of 330 STAT lines.
> +++++++++++++++++++++++++++++++++++++++
>
> Since you just sent us one Wavelet-Stat output file, that's all we
could
> test with. But in this example, I'm aggregating together output for
the
> ">=50.000" threshold. I'm guessing you'll want to
> rerun this type of job for each of the thresholds you used in
Wavelet-Stat.
> When you're first setting up jobs, it's a good idea to use the "-
dump_row"
> option. That'll write all of the STAT lines the
> job actually used to the file name you specify. That way, you can
make
> sure that STAT-Analysis is operating on the exact subset of data
that you
> expect.
>
> Hope that helps clarify.
>
> John Halley Gotway
> met_help at ucar.edu
>
> RAL HelpDesk {for Hashmi Fatima} wrote:
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >
> > Dear Paul,
> >
> > Sorry for the delay. I have put my files at ftp site in my folder.
I have
> > done the verufication through wavelet for 1 month data. Here I am
sending
> > few files.
> >
> > *input files (WPP output):*
> > gfs.t00z.pgrbf00 (analysis file)
> > gfs.t00z.pgrbf24 (forecast file)
> >
> > *wavelet output file (stat analysis input) :*
> > wavelet_stat_000000L_20100920_000000V.nc
> > wavelet_stat_000000L_20100920_000000V.ps
> > wavelet_stat_000000L_20100920_000000V.stat
> > wavelet_stat_000000L_20100920_000000V_isc.txt
> >
> > The error message which I am getting after running stat_analysis
tool
> after
> > using wavelet output is in file *error_msg.txt*.
> >
> > Thanks and Regards
> > Hashmi Fatima
> >
> > On Wed, Sep 29, 2010 at 8:49 PM, RAL HelpDesk {for Paul Oldenburg}
<
> > met_help at ucar.edu> wrote:
> >
> >> Hashmi,
> >>
> >> Here is the website, sorry:
> >>
> >> http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >>
> >> If you cannot get this to work, please let us know and we will
have to
> make
> >> other arrangements for you to send us your
> >> data. Sorry for all the trouble!
> >>
> >> Thanks,
> >>
> >> Paul
> >>
> >>
> >> RAL HelpDesk {for Paul Oldenburg} wrote:
> >>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >>>
> >>> Hashmi,
> >>>
> >>> We apparently have a new FTP server configuration that we are
still
> >> getting used to. Can you please refer to the
> >>> instructions on the following website for uploading data and let
us
> know
> >> if you can get it to work:
> >>>
> >>>
> >>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >>>>
> >>>> Dear Paul
> >>>>
> >>>> ftp -p ftp.rap.ucar.edu is not working. It is giving me error
> message.
> >>>>
> >>>> ftp -p ftp.rap.ucar.edu
> >>>> ftp: p: unknown option
> >>>> Usage: ftp [-v] [-d] [-i] [-n] [-g] [-k realm] [-f] [-x] [-u]
[-t]
> >> [host]
> >>>> Regards
> >>>> Hashmi Fatima
> >>>>
> >>>>
> >>>>
> >>>> On Wed, Sep 29, 2010 at 12:22 AM, RAL HelpDesk {for Paul
Oldenburg} <
> >>>> met_help at ucar.edu> wrote:
> >>>>
> >>>>> Hashmi,
> >>>>>
> >>>>> Please try to invoke your ftp client using the -p command line
> >> argument,
> >>>>> like so:
> >>>>>
> >>>>> ftp -p ftp.rap.ucar.edu
> >>>>> username = anonymous
> >>>>> password = "your email address"
> >>>>> cd incoming/irap/met_help/fatima_data
> >>>>> put "your files"
> >>>>>
> >>>>> Let us know how this approach goes.
> >>>>>
> >>>>> Thanks,
> >>>>>
> >>>>> Paul
> >>>>>
> >>>>>
> >>>>> Paul Oldenburg wrote:
> >>>>>> Hashmi,
> >>>>>>
> >>>>>> We are investigating a potential problem with our FTP site.
I will
> >> let
> >>>>> you know when it is fixed, and then we can
> >>>>>> proceed with addressing your wavelet_stat question.
> >>>>>>
> >>>>>> Thanks,
> >>>>>>
> >>>>>> Paul
> >>>>>>
> >>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >>>>>>>
> >>>>>>> Dear Paul
> >>>>>>>
> >>>>>>> I tried to upload my files at the FTP server you have told
me in
> >>>>> previous
> >>>>>>> mails. But it says " Permission denied on server".
> >>>>>>>
> >>>>>>> Regards
> >>>>>>> Hashmi Fatima
> >>>>>>>
> >>>>>>> On Mon, Sep 27, 2010 at 10:25 PM, RAL HelpDesk {for Paul
Oldenburg}
> <
> >>>>>>> met_help at ucar.edu> wrote:
> >>>>>>>
> >>>>>>>> Hashmi,
> >>>>>>>>
> >>>>>>>> I (Paul) have returned to my met_help responsibilities and
I have
> >>>>> followed
> >>>>>>>> your interactions with John. In order to
> >>>>>>>> help you run stat_analysis, I would like to have you send
me some
> of
> >>>>> your
> >>>>>>>> data via FTP. First, can you give us the WPP
> >>>>>>>> output that you fed into wavelet_stat? Also, can you give
us your
> >>>>>>>> wavelet_stat output? Finally, please copy the error
> >>>>>>>> message that you get from stat_analysis in your response.
That
> >>>>> information
> >>>>>>>> will be very helpful.
> >>>>>>>>
> >>>>>>>> Concerning the masking region selected by wavelet_stat, the
entire
> >> grid
> >>>>>>>> cannot be processed because wavelet_stat can
> >>>>>>>> process only dyadic regions - that is square regions whose
sides
> are
> >> a
> >>>>>>>> power of 2 grid points (like 256 or 2048). By
> >>>>>>>> default, wavelet_stat should try to guess the best regions
to use
> in
> >>>>> your
> >>>>>>>> domain. If you send us your WPP output, we
> >>>>>>>> will be able to help you configure the best masking region
for
> your
> >>>>> case.
> >>>>>>>> We look forward to receiving your data. If you have any
> questions,
> >>>>> please
> >>>>>>>> let me know.
> >>>>>>>>
> >>>>>>>> Thanks,
> >>>>>>>>
> >>>>>>>> Paul
> >>>>>>>>
> >>>>>>>>
> >>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >>>>>>>>>
> >>>>>>>>> Dear Sir,
> >>>>>>>>>
> >>>>>>>>> Thanks for your help.Its working now.
> >>>>>>>>>
> >>>>>>>>> I have run wavelet tool for computation of RH. Now i want
to use
> >> the
> >>>>>>>> output
> >>>>>>>>> of wavelet for StatAnalysis. but its not working. I am
sending
> you
> >>>>>>>>> configuration file of wavelet and StatAnalysis. Could you
please
> >>>>> check.
> >>>>>>>>> One more thing I want to know that in the wavelet output,
why
> does
> >> it
> >>>>>>>> take
> >>>>>>>>> that much of area. I have not assigned any area for the
> >> verification
> >>>>>>>> except
> >>>>>>>>> full. I have masked the whole grid. I am attaching that
output
> too
> >> as
> >>>>>>>>> Wavelet_output.png.
> >>>>>>>>>
> >>>>>>>>> Thanks and Regards
> >>>>>>>>> Hashmi fatima
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>> On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for John
Halley
> >> Gotway}
> >>>>> <
> >>>>>>>>> met_help at ucar.edu> wrote:
> >>>>>>>>>
> >>>>>>>>>> Hashmi,
> >>>>>>>>>>
> >>>>>>>>>> Paul will be out of the office for the next couple of
weeks, but
> >> I'll
> >>>>> be
> >>>>>>>>>> happy to help you with this issue.
> >>>>>>>>>>
> >>>>>>>>>> There are a couple of points I'd like to make:
> >>>>>>>>>>
> >>>>>>>>>> First, the NetCDF files you sent us are the NetCDF output
of
> WRF.
> >>>>> MET
> >>>>>>>> does
> >>>>>>>>>> not directly read the NetCDF output of WRF. Instead, for
> METv2.0,
> >>>>> your
> >>>>>>>> WRF
> >>>>>>>>>> output must first be run through the
> >>>>>>>>>> WRF-PostProcessor (WPP). Details of running WPP can be
found in
> >>>>> chapter
> >>>>>>>> 7
> >>>>>>>>>> of the WRF user's guide:
> >>>>>>>>>>
> >>>>>>>>>>
> >>
> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
> >>>>>>>>>> WPP writes out GRIB data, which MET can read.
> >>>>>>>>>>
> >>>>>>>>>> Second, once you've run your data through the WRF
PostProcessor,
> >> you
> >>>>> can
> >>>>>>>>>> select multiple thresholds for Grid-Stat to use. You do
so by
> >>>>> setting
> >>>>>>>> up
> >>>>>>>>>> the Grid-Stat configuration file something like this:
> >>>>>>>>>>
> >>>>>>>>>> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
> >>>>>>>>>> obs_field[] = [];
> >>>>>>>>>> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240
ge260
> ge280
> >>>>> ge300"
> >>>>>>>> ];
> >>>>>>>>>> obs_thresh[] = [];
> >>>>>>>>>>
> >>>>>>>>>> The settings above tell Grid-Stat to verify temperature
at
> >> 2-meters
> >>>>> as
> >>>>>>>> well
> >>>>>>>>>> as temperature at 500mb. For 2-meter temp, it'll use the
2
> >>>>> thresholds
> >>>>>>>>>> listed (ge273 ge283) and for temp at 500mb, it'll use
> >>>>>>>>>> the 6 thresholds listed (ge200.0 ge220 ge240 ge260 ge280
ge300).
> >>>>>>>>>>
> >>>>>>>>>> So the real task is running your WRF output files through
WPP.
> >>>>>>>>>>
> >>>>>>>>>> Lastly, I wanted to tell you that in METv3.0 (to be
released
> >>>>> shortly),
> >>>>>>>>>> we've added support for reading the NetCDF output of the
pinterp
> >> tool
> >>>>>>>> (which
> >>>>>>>>>> only processes ARW output). So once METv3.0 is
> >>>>>>>>>> available, you can use either the GRIB output of WPP or
the
> NetCDF
> >>>>>>>> output
> >>>>>>>>>> of the pinterp tool.
> >>>>>>>>>>
> >>>>>>>>>> Hope that helps.
> >>>>>>>>>>
> >>>>>>>>>> Thanks,
> >>>>>>>>>> John Halley Gotway
> >>>>>>>>>> met_help at ucar.edu
> >>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
> >>>>>>>>>>>
> >>>>>>>>>>> Dear Sir,
> >>>>>>>>>>>
> >>>>>>>>>>> I have put my files at the FTP site told by you. Please
check.
> >>>>>>>>>>>
> >>>>>>>>>>> One more thing, i want to know that do i need to give
equal
> >> number
> >>>>> of
> >>>>>>>>>>> threshold values as fcst_field/obs_field? Suppose I
there is
> only
> >>>>> one
> >>>>>>>>>>> fcst_field/ob_field, but I want to check it for various
> threshold
> >>>>>>>> values.
> >>>>>>>>>>> In some other case, I have given 5 fcst_field/obs_field
and
> >>>>>>>> corresponding
> >>>>>>>>>> 5
> >>>>>>>>>>> threshold (fcst_thresh/obs_thresh) values. But I did
not get
> >>>>> anything
> >>>>>>>> in
> >>>>>>>>>>> the output. I have got netcdf file in the output but
nothing in
> >> the
> >>>>>>>> STAT
> >>>>>>>>>> and
> >>>>>>>>>>> .txt files.
> >>>>>>>>>>>
> >>>>>>>>>>> Thanks and Regards
> >>>>>>>>>>> Hashmi Fatima
> >>>>>>>>>>>
> >>>>>>>>>>> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for Paul
> >> Oldenburg} <
> >>>>>>>>>>> met_help at ucar.edu> wrote:
> >>>>>>>>>>>
> >>>>>>>>>>>> Hashmi,
> >>>>>>>>>>>>
> >>>>>>>>>>>> Can you please send us an example of your model and obs
data
> >> along
> >>>>>>>> with
> >>>>>>>>>>>> your grid_stat config file? If the data is too
> >>>>>>>>>>>> large to attach to an email, please upload it to our
FTP site:
> >>>>>>>>>>>>
> >>>>>>>>>>>> ftp ftp.rap.ucar.edu
> >>>>>>>>>>>> username = anonymous
> >>>>>>>>>>>> password = "your email address"
> >>>>>>>>>>>> cd incoming/irap/met_help/fatima_data
> >>>>>>>>>>>> put "your files"
> >>>>>>>>>>>>
> >>>>>>>>>>>> Thanks,
> >>>>>>>>>>>>
> >>>>>>>>>>>> Paul
> >>>>>>>>>>>>
> >>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>>>>>>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted
upon.
> >>>>>>>>>>>>> Transaction: Ticket created by hashmifatima at gmail.com
> >>>>>>>>>>>>> Queue: met_help
> >>>>>>>>>>>>> Subject: netcdf file for grid_stat (grid code
error)
> >>>>>>>>>>>>> Owner: Nobody
> >>>>>>>>>>>>> Requestors: hashmifatima at gmail.com
> >>>>>>>>>>>>> Status: new
> >>>>>>>>>>>>> Ticket <URL:
> >>>>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
> >>>>>>>>>>>>> Dear Sir,
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> I want to use grid_stat tool to verify 2m temperature.
Both
> >>>>>>>>>> (observation
> >>>>>>>>>>>> and
> >>>>>>>>>>>>> forecast) files are in netcdf format. when i am
defining
> >>>>> fcst_field,
> >>>>>>>> it
> >>>>>>>>>>>> is
> >>>>>>>>>>>>> asking for grid code. I am not using any grib file for
this
> >> test.
> >>>>>>>> What
> >>>>>>>>>>>>> should I do ?
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> Thanks and Regards,
> >>
------------------------------------------------------------------------
> >>
> >>
> >
> >
>
>
--
Hashmi Fatima
Research Scholar
Centre for Atmospheric Sciences
Indian Institute of Technology Delhi
Hauz Khas
New Delhi - 110016
INDIA
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #40804] netcdf file for grid_stat (grid code error)
From: Paul Oldenburg
Time: Tue Oct 05 14:44:40 2010
Hashmi,
We inspected your stat_analysis jobs and saw that your criteria
resulted in no MET output data being found. The reason
for this is that you are trying to select data in line types SL1L2 and
CNT, yet you are also specifying thresholds in
the same job. Line types SL1L2 and CNT contain statistics that are
not divided among thresholds, because they are
continuous (as opposed to categorical or thresholded). Please look
carefully at your MET output and you will see that
these lines have a value of NA in the fcst_thresh line. If you have
any questions, please let us know.
Thanks,
Paul
RAL HelpDesk {for Hashmi Fatima} wrote:
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>
> Dear Sir
>
> I have configured STATconfig file for the previous case.
>
> But now I am trying STATanalysis with grid_stat output for HGT at
P500. But
> its again giving some error message. I have uploaded the grid_stat
output
> files to your FTP server, which are input for STATanalysis case (all
> grid_stat_HGT_P500-240000L_20101001_00000V_* files). STATconfig file
is
> STATAnalysisConfig_HGT. I have configurered the STATconfig file
according to
> the input. but still unable to get the output.
> The error message which I am getting is srored in the file
> grid_output_stat_analysis_error.txt.
>
> Please check and let me know about the error.
>
> Thanks and Regards
> Hashmi Fatima
>
>
>
> On Sat, Oct 2, 2010 at 2:49 AM, RAL HelpDesk {for John Halley
Gotway} <
> met_help at ucar.edu> wrote:
>
>> Hashmi,
>>
>> Paul's out of the office, so hopefully I can provide some advice on
this
>> issue.
>>
>> It looks like you're using the STATAnalysis config file that's
included
>> with the MET test scripts. This is really only meant as an example
to work
>> with the test data. You'll need to define the
>> analysis jobs explicitly to match the output data you have.
>>
>> Since you're running Wavelet-Stat, I assume you may want to
aggregate the
>> output of Wavelet-Stat across multiple cases. Please keep in mind
that
>> Stat-Analysis can be run on the command line OR using
>> a config file. When you're trying to set up your analysis job at
the
>> beginning, I'd suggest using the command line option. For example,
here's
>> the output I get for the following command:
>>
>> +++++++++++++++++++++++++++++++++++++++
>> bin/stat_analysis -lookin
wavelet_stat_000000L_20100920_000000V.stat -job
>> aggregate -line_type ISC -vx_mask TILE_TOT -fcst_thresh ">=50.000"
-dump_row
>> aggr_ISC_50.stat -v 2
>> STAT Lines read = 330
>> STAT Lines retained = 330
>>
>> Processing Job 1: -job aggregate -vx_mask TILE_TOT -fcst_thresh
>=50.000
>> -line_type ISC -dump_row aggr_ISC_50.stat
>> Creating STAT output file "aggr_ISC_50.stat"
>> JOB_LIST: -job aggregate -vx_mask TILE_TOT -fcst_thresh
>=50.000
>> -line_type ISC -dump_row aggr_ISC_50.stat
>> COL_NAME: TOTAL TILE_DIM TILE_XLL TILE_YLL NSCALE ISCALE MSE
ISC
>> FENERGY2 OENERGY2 BASER FBIAS
>> ISC: 131072 NA NA NA 8 0 0.09870
0.63219
>> 0.16414 0.15513 0.15513 1.05808
>> 131072 NA NA NA 8 1 0.01761
0.40938
>> 0.00935 0.00921 0.15513 1.05808
>> 131072 NA NA NA 8 2 0.01881
0.36913
>> 0.01064 0.01070 0.15513 1.05808
>> 131072 NA NA NA 8 3 0.02372
0.20445
>> 0.01519 0.01605 0.15513 1.05808
>> 131072 NA NA NA 8 4 0.01836
0.38422
>> 0.01856 0.01899 0.15513 1.05808
>> 131072 NA NA NA 8 5 0.01390
0.53381
>> 0.01933 0.01968 0.15513 1.05808
>> 131072 NA NA NA 8 6 0.00524
0.82426
>> 0.02334 0.02096 0.15513 1.05808
>> 131072 NA NA NA 8 7 0.00082
0.97250
>> 0.02947 0.02577 0.15513 1.05808
>> 131072 NA NA NA 8 8 0.00011
0.99631
>> 0.00854 0.00764 0.15513 1.05808
>> 131072 NA NA NA 8 9 0.00013
0.99564
>> 0.02972 0.02614 0.15513 1.05808
>>
>> Job 1 used 10 out of 330 STAT lines.
>> +++++++++++++++++++++++++++++++++++++++
>>
>> Since you just sent us one Wavelet-Stat output file, that's all we
could
>> test with. But in this example, I'm aggregating together output
for the
>> ">=50.000" threshold. I'm guessing you'll want to
>> rerun this type of job for each of the thresholds you used in
Wavelet-Stat.
>> When you're first setting up jobs, it's a good idea to use the "-
dump_row"
>> option. That'll write all of the STAT lines the
>> job actually used to the file name you specify. That way, you can
make
>> sure that STAT-Analysis is operating on the exact subset of data
that you
>> expect.
>>
>> Hope that helps clarify.
>>
>> John Halley Gotway
>> met_help at ucar.edu
>>
>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>
>>> Dear Paul,
>>>
>>> Sorry for the delay. I have put my files at ftp site in my folder.
I have
>>> done the verufication through wavelet for 1 month data. Here I am
sending
>>> few files.
>>>
>>> *input files (WPP output):*
>>> gfs.t00z.pgrbf00 (analysis file)
>>> gfs.t00z.pgrbf24 (forecast file)
>>>
>>> *wavelet output file (stat analysis input) :*
>>> wavelet_stat_000000L_20100920_000000V.nc
>>> wavelet_stat_000000L_20100920_000000V.ps
>>> wavelet_stat_000000L_20100920_000000V.stat
>>> wavelet_stat_000000L_20100920_000000V_isc.txt
>>>
>>> The error message which I am getting after running stat_analysis
tool
>> after
>>> using wavelet output is in file *error_msg.txt*.
>>>
>>> Thanks and Regards
>>> Hashmi Fatima
>>>
>>> On Wed, Sep 29, 2010 at 8:49 PM, RAL HelpDesk {for Paul Oldenburg}
<
>>> met_help at ucar.edu> wrote:
>>>
>>>> Hashmi,
>>>>
>>>> Here is the website, sorry:
>>>>
>>>> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>>>>
>>>> If you cannot get this to work, please let us know and we will
have to
>> make
>>>> other arrangements for you to send us your
>>>> data. Sorry for all the trouble!
>>>>
>>>> Thanks,
>>>>
>>>> Paul
>>>>
>>>>
>>>> RAL HelpDesk {for Paul Oldenburg} wrote:
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>
>>>>> Hashmi,
>>>>>
>>>>> We apparently have a new FTP server configuration that we are
still
>>>> getting used to. Can you please refer to the
>>>>> instructions on the following website for uploading data and let
us
>> know
>>>> if you can get it to work:
>>>>>
>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>>
>>>>>> Dear Paul
>>>>>>
>>>>>> ftp -p ftp.rap.ucar.edu is not working. It is giving me error
>> message.
>>>>>> ftp -p ftp.rap.ucar.edu
>>>>>> ftp: p: unknown option
>>>>>> Usage: ftp [-v] [-d] [-i] [-n] [-g] [-k realm] [-f] [-x] [-u]
[-t]
>>>> [host]
>>>>>> Regards
>>>>>> Hashmi Fatima
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Wed, Sep 29, 2010 at 12:22 AM, RAL HelpDesk {for Paul
Oldenburg} <
>>>>>> met_help at ucar.edu> wrote:
>>>>>>
>>>>>>> Hashmi,
>>>>>>>
>>>>>>> Please try to invoke your ftp client using the -p command line
>>>> argument,
>>>>>>> like so:
>>>>>>>
>>>>>>> ftp -p ftp.rap.ucar.edu
>>>>>>> username = anonymous
>>>>>>> password = "your email address"
>>>>>>> cd incoming/irap/met_help/fatima_data
>>>>>>> put "your files"
>>>>>>>
>>>>>>> Let us know how this approach goes.
>>>>>>>
>>>>>>> Thanks,
>>>>>>>
>>>>>>> Paul
>>>>>>>
>>>>>>>
>>>>>>> Paul Oldenburg wrote:
>>>>>>>> Hashmi,
>>>>>>>>
>>>>>>>> We are investigating a potential problem with our FTP site.
I will
>>>> let
>>>>>>> you know when it is fixed, and then we can
>>>>>>>> proceed with addressing your wavelet_stat question.
>>>>>>>>
>>>>>>>> Thanks,
>>>>>>>>
>>>>>>>> Paul
>>>>>>>>
>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>>>>>
>>>>>>>>> Dear Paul
>>>>>>>>>
>>>>>>>>> I tried to upload my files at the FTP server you have told
me in
>>>>>>> previous
>>>>>>>>> mails. But it says " Permission denied on server".
>>>>>>>>>
>>>>>>>>> Regards
>>>>>>>>> Hashmi Fatima
>>>>>>>>>
>>>>>>>>> On Mon, Sep 27, 2010 at 10:25 PM, RAL HelpDesk {for Paul
Oldenburg}
>> <
>>>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>>>
>>>>>>>>>> Hashmi,
>>>>>>>>>>
>>>>>>>>>> I (Paul) have returned to my met_help responsibilities and
I have
>>>>>>> followed
>>>>>>>>>> your interactions with John. In order to
>>>>>>>>>> help you run stat_analysis, I would like to have you send
me some
>> of
>>>>>>> your
>>>>>>>>>> data via FTP. First, can you give us the WPP
>>>>>>>>>> output that you fed into wavelet_stat? Also, can you give
us your
>>>>>>>>>> wavelet_stat output? Finally, please copy the error
>>>>>>>>>> message that you get from stat_analysis in your response.
That
>>>>>>> information
>>>>>>>>>> will be very helpful.
>>>>>>>>>>
>>>>>>>>>> Concerning the masking region selected by wavelet_stat, the
entire
>>>> grid
>>>>>>>>>> cannot be processed because wavelet_stat can
>>>>>>>>>> process only dyadic regions - that is square regions whose
sides
>> are
>>>> a
>>>>>>>>>> power of 2 grid points (like 256 or 2048). By
>>>>>>>>>> default, wavelet_stat should try to guess the best regions
to use
>> in
>>>>>>> your
>>>>>>>>>> domain. If you send us your WPP output, we
>>>>>>>>>> will be able to help you configure the best masking region
for
>> your
>>>>>>> case.
>>>>>>>>>> We look forward to receiving your data. If you have any
>> questions,
>>>>>>> please
>>>>>>>>>> let me know.
>>>>>>>>>>
>>>>>>>>>> Thanks,
>>>>>>>>>>
>>>>>>>>>> Paul
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>>>>>>>
>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>
>>>>>>>>>>> Thanks for your help.Its working now.
>>>>>>>>>>>
>>>>>>>>>>> I have run wavelet tool for computation of RH. Now i want
to use
>>>> the
>>>>>>>>>> output
>>>>>>>>>>> of wavelet for StatAnalysis. but its not working. I am
sending
>> you
>>>>>>>>>>> configuration file of wavelet and StatAnalysis. Could you
please
>>>>>>> check.
>>>>>>>>>>> One more thing I want to know that in the wavelet output,
why
>> does
>>>> it
>>>>>>>>>> take
>>>>>>>>>>> that much of area. I have not assigned any area for the
>>>> verification
>>>>>>>>>> except
>>>>>>>>>>> full. I have masked the whole grid. I am attaching that
output
>> too
>>>> as
>>>>>>>>>>> Wavelet_output.png.
>>>>>>>>>>>
>>>>>>>>>>> Thanks and Regards
>>>>>>>>>>> Hashmi fatima
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for John
Halley
>>>> Gotway}
>>>>>>> <
>>>>>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>>>>>
>>>>>>>>>>>> Hashmi,
>>>>>>>>>>>>
>>>>>>>>>>>> Paul will be out of the office for the next couple of
weeks, but
>>>> I'll
>>>>>>> be
>>>>>>>>>>>> happy to help you with this issue.
>>>>>>>>>>>>
>>>>>>>>>>>> There are a couple of points I'd like to make:
>>>>>>>>>>>>
>>>>>>>>>>>> First, the NetCDF files you sent us are the NetCDF output
of
>> WRF.
>>>>>>> MET
>>>>>>>>>> does
>>>>>>>>>>>> not directly read the NetCDF output of WRF. Instead, for
>> METv2.0,
>>>>>>> your
>>>>>>>>>> WRF
>>>>>>>>>>>> output must first be run through the
>>>>>>>>>>>> WRF-PostProcessor (WPP). Details of running WPP can be
found in
>>>>>>> chapter
>>>>>>>>>> 7
>>>>>>>>>>>> of the WRF user's guide:
>>>>>>>>>>>>
>>>>>>>>>>>>
>> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
>>>>>>>>>>>> WPP writes out GRIB data, which MET can read.
>>>>>>>>>>>>
>>>>>>>>>>>> Second, once you've run your data through the WRF
PostProcessor,
>>>> you
>>>>>>> can
>>>>>>>>>>>> select multiple thresholds for Grid-Stat to use. You do
so by
>>>>>>> setting
>>>>>>>>>> up
>>>>>>>>>>>> the Grid-Stat configuration file something like this:
>>>>>>>>>>>>
>>>>>>>>>>>> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
>>>>>>>>>>>> obs_field[] = [];
>>>>>>>>>>>> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240
ge260
>> ge280
>>>>>>> ge300"
>>>>>>>>>> ];
>>>>>>>>>>>> obs_thresh[] = [];
>>>>>>>>>>>>
>>>>>>>>>>>> The settings above tell Grid-Stat to verify temperature
at
>>>> 2-meters
>>>>>>> as
>>>>>>>>>> well
>>>>>>>>>>>> as temperature at 500mb. For 2-meter temp, it'll use the
2
>>>>>>> thresholds
>>>>>>>>>>>> listed (ge273 ge283) and for temp at 500mb, it'll use
>>>>>>>>>>>> the 6 thresholds listed (ge200.0 ge220 ge240 ge260 ge280
ge300).
>>>>>>>>>>>>
>>>>>>>>>>>> So the real task is running your WRF output files through
WPP.
>>>>>>>>>>>>
>>>>>>>>>>>> Lastly, I wanted to tell you that in METv3.0 (to be
released
>>>>>>> shortly),
>>>>>>>>>>>> we've added support for reading the NetCDF output of the
pinterp
>>>> tool
>>>>>>>>>> (which
>>>>>>>>>>>> only processes ARW output). So once METv3.0 is
>>>>>>>>>>>> available, you can use either the GRIB output of WPP or
the
>> NetCDF
>>>>>>>>>> output
>>>>>>>>>>>> of the pinterp tool.
>>>>>>>>>>>>
>>>>>>>>>>>> Hope that helps.
>>>>>>>>>>>>
>>>>>>>>>>>> Thanks,
>>>>>>>>>>>> John Halley Gotway
>>>>>>>>>>>> met_help at ucar.edu
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
>>>>>>>>>>>>>
>>>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>>>
>>>>>>>>>>>>> I have put my files at the FTP site told by you. Please
check.
>>>>>>>>>>>>>
>>>>>>>>>>>>> One more thing, i want to know that do i need to give
equal
>>>> number
>>>>>>> of
>>>>>>>>>>>>> threshold values as fcst_field/obs_field? Suppose I
there is
>> only
>>>>>>> one
>>>>>>>>>>>>> fcst_field/ob_field, but I want to check it for various
>> threshold
>>>>>>>>>> values.
>>>>>>>>>>>>> In some other case, I have given 5 fcst_field/obs_field
and
>>>>>>>>>> corresponding
>>>>>>>>>>>> 5
>>>>>>>>>>>>> threshold (fcst_thresh/obs_thresh) values. But I did
not get
>>>>>>> anything
>>>>>>>>>> in
>>>>>>>>>>>>> the output. I have got netcdf file in the output but
nothing in
>>>> the
>>>>>>>>>> STAT
>>>>>>>>>>>> and
>>>>>>>>>>>>> .txt files.
>>>>>>>>>>>>>
>>>>>>>>>>>>> Thanks and Regards
>>>>>>>>>>>>> Hashmi Fatima
>>>>>>>>>>>>>
>>>>>>>>>>>>> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for Paul
>>>> Oldenburg} <
>>>>>>>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>>>>>>>
>>>>>>>>>>>>>> Hashmi,
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Can you please send us an example of your model and obs
data
>>>> along
>>>>>>>>>> with
>>>>>>>>>>>>>> your grid_stat config file? If the data is too
>>>>>>>>>>>>>> large to attach to an email, please upload it to our
FTP site:
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> ftp ftp.rap.ucar.edu
>>>>>>>>>>>>>> username = anonymous
>>>>>>>>>>>>>> password = "your email address"
>>>>>>>>>>>>>> cd incoming/irap/met_help/fatima_data
>>>>>>>>>>>>>> put "your files"
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Paul
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>>>>>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted
upon.
>>>>>>>>>>>>>>> Transaction: Ticket created by hashmifatima at gmail.com
>>>>>>>>>>>>>>> Queue: met_help
>>>>>>>>>>>>>>> Subject: netcdf file for grid_stat (grid code
error)
>>>>>>>>>>>>>>> Owner: Nobody
>>>>>>>>>>>>>>> Requestors: hashmifatima at gmail.com
>>>>>>>>>>>>>>> Status: new
>>>>>>>>>>>>>>> Ticket <URL:
>>>>>>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
>>>>>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> I want to use grid_stat tool to verify 2m temperature.
Both
>>>>>>>>>>>> (observation
>>>>>>>>>>>>>> and
>>>>>>>>>>>>>>> forecast) files are in netcdf format. when i am
defining
>>>>>>> fcst_field,
>>>>>>>>>> it
>>>>>>>>>>>>>> is
>>>>>>>>>>>>>>> asking for grid code. I am not using any grib file for
this
>>>> test.
>>>>>>>>>> What
>>>>>>>>>>>>>>> should I do ?
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Thanks and Regards,
>>>>
------------------------------------------------------------------------
>>>>
>>>>
>>>
>>
>
>
------------------------------------------------
Subject: netcdf file for grid_stat (grid code error)
From: Hashmi Fatima
Time: Wed Oct 06 03:34:37 2010
Dear Sir,
Thanks for this information.
Could you please tell me that what line_type should i use to get
output in
stat_analysis?
Regards
Hashmi Fatima
On Wed, Oct 6, 2010 at 2:14 AM, RAL HelpDesk {for Paul Oldenburg} <
met_help at ucar.edu> wrote:
> Hashmi,
>
> We inspected your stat_analysis jobs and saw that your criteria
resulted in
> no MET output data being found. The reason
> for this is that you are trying to select data in line types SL1L2
and CNT,
> yet you are also specifying thresholds in
> the same job. Line types SL1L2 and CNT contain statistics that are
not
> divided among thresholds, because they are
> continuous (as opposed to categorical or thresholded). Please look
> carefully at your MET output and you will see that
> these lines have a value of NA in the fcst_thresh line. If you have
any
> questions, please let us know.
>
> Thanks,
>
> Paul
>
>
> RAL HelpDesk {for Hashmi Fatima} wrote:
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >
> > Dear Sir
> >
> > I have configured STATconfig file for the previous case.
> >
> > But now I am trying STATanalysis with grid_stat output for HGT at
P500.
> But
> > its again giving some error message. I have uploaded the grid_stat
output
> > files to your FTP server, which are input for STATanalysis case
(all
> > grid_stat_HGT_P500-240000L_20101001_00000V_* files). STATconfig
file is
> > STATAnalysisConfig_HGT. I have configurered the STATconfig file
according
> to
> > the input. but still unable to get the output.
> > The error message which I am getting is srored in the file
> > grid_output_stat_analysis_error.txt.
> >
> > Please check and let me know about the error.
> >
> > Thanks and Regards
> > Hashmi Fatima
> >
> >
> >
> > On Sat, Oct 2, 2010 at 2:49 AM, RAL HelpDesk {for John Halley
Gotway} <
> > met_help at ucar.edu> wrote:
> >
> >> Hashmi,
> >>
> >> Paul's out of the office, so hopefully I can provide some advice
on this
> >> issue.
> >>
> >> It looks like you're using the STATAnalysis config file that's
included
> >> with the MET test scripts. This is really only meant as an
example to
> work
> >> with the test data. You'll need to define the
> >> analysis jobs explicitly to match the output data you have.
> >>
> >> Since you're running Wavelet-Stat, I assume you may want to
aggregate
> the
> >> output of Wavelet-Stat across multiple cases. Please keep in
mind that
> >> Stat-Analysis can be run on the command line OR using
> >> a config file. When you're trying to set up your analysis job at
the
> >> beginning, I'd suggest using the command line option. For
example,
> here's
> >> the output I get for the following command:
> >>
> >> +++++++++++++++++++++++++++++++++++++++
> >> bin/stat_analysis -lookin
wavelet_stat_000000L_20100920_000000V.stat
> -job
> >> aggregate -line_type ISC -vx_mask TILE_TOT -fcst_thresh
">=50.000"
> -dump_row
> >> aggr_ISC_50.stat -v 2
> >> STAT Lines read = 330
> >> STAT Lines retained = 330
> >>
> >> Processing Job 1: -job aggregate -vx_mask TILE_TOT -fcst_thresh
>=50.000
> >> -line_type ISC -dump_row aggr_ISC_50.stat
> >> Creating STAT output file "aggr_ISC_50.stat"
> >> JOB_LIST: -job aggregate -vx_mask TILE_TOT -fcst_thresh
>=50.000
> >> -line_type ISC -dump_row aggr_ISC_50.stat
> >> COL_NAME: TOTAL TILE_DIM TILE_XLL TILE_YLL NSCALE ISCALE MSE
ISC
> >> FENERGY2 OENERGY2 BASER FBIAS
> >> ISC: 131072 NA NA NA 8 0 0.09870
0.63219
> >> 0.16414 0.15513 0.15513 1.05808
> >> 131072 NA NA NA 8 1 0.01761
0.40938
> >> 0.00935 0.00921 0.15513 1.05808
> >> 131072 NA NA NA 8 2 0.01881
0.36913
> >> 0.01064 0.01070 0.15513 1.05808
> >> 131072 NA NA NA 8 3 0.02372
0.20445
> >> 0.01519 0.01605 0.15513 1.05808
> >> 131072 NA NA NA 8 4 0.01836
0.38422
> >> 0.01856 0.01899 0.15513 1.05808
> >> 131072 NA NA NA 8 5 0.01390
0.53381
> >> 0.01933 0.01968 0.15513 1.05808
> >> 131072 NA NA NA 8 6 0.00524
0.82426
> >> 0.02334 0.02096 0.15513 1.05808
> >> 131072 NA NA NA 8 7 0.00082
0.97250
> >> 0.02947 0.02577 0.15513 1.05808
> >> 131072 NA NA NA 8 8 0.00011
0.99631
> >> 0.00854 0.00764 0.15513 1.05808
> >> 131072 NA NA NA 8 9 0.00013
0.99564
> >> 0.02972 0.02614 0.15513 1.05808
> >>
> >> Job 1 used 10 out of 330 STAT lines.
> >> +++++++++++++++++++++++++++++++++++++++
> >>
> >> Since you just sent us one Wavelet-Stat output file, that's all
we could
> >> test with. But in this example, I'm aggregating together output
for the
> >> ">=50.000" threshold. I'm guessing you'll want to
> >> rerun this type of job for each of the thresholds you used in
> Wavelet-Stat.
> >> When you're first setting up jobs, it's a good idea to use the
> "-dump_row"
> >> option. That'll write all of the STAT lines the
> >> job actually used to the file name you specify. That way, you
can make
> >> sure that STAT-Analysis is operating on the exact subset of data
that
> you
> >> expect.
> >>
> >> Hope that helps clarify.
> >>
> >> John Halley Gotway
> >> met_help at ucar.edu
> >>
> >> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >>>
> >>> Dear Paul,
> >>>
> >>> Sorry for the delay. I have put my files at ftp site in my
folder. I
> have
> >>> done the verufication through wavelet for 1 month data. Here I
am
> sending
> >>> few files.
> >>>
> >>> *input files (WPP output):*
> >>> gfs.t00z.pgrbf00 (analysis file)
> >>> gfs.t00z.pgrbf24 (forecast file)
> >>>
> >>> *wavelet output file (stat analysis input) :*
> >>> wavelet_stat_000000L_20100920_000000V.nc
> >>> wavelet_stat_000000L_20100920_000000V.ps
> >>> wavelet_stat_000000L_20100920_000000V.stat
> >>> wavelet_stat_000000L_20100920_000000V_isc.txt
> >>>
> >>> The error message which I am getting after running stat_analysis
tool
> >> after
> >>> using wavelet output is in file *error_msg.txt*.
> >>>
> >>> Thanks and Regards
> >>> Hashmi Fatima
> >>>
> >>> On Wed, Sep 29, 2010 at 8:49 PM, RAL HelpDesk {for Paul
Oldenburg} <
> >>> met_help at ucar.edu> wrote:
> >>>
> >>>> Hashmi,
> >>>>
> >>>> Here is the website, sorry:
> >>>>
> >>>> http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >>>>
> >>>> If you cannot get this to work, please let us know and we will
have to
> >> make
> >>>> other arrangements for you to send us your
> >>>> data. Sorry for all the trouble!
> >>>>
> >>>> Thanks,
> >>>>
> >>>> Paul
> >>>>
> >>>>
> >>>> RAL HelpDesk {for Paul Oldenburg} wrote:
> >>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804
>
> >>>>>
> >>>>> Hashmi,
> >>>>>
> >>>>> We apparently have a new FTP server configuration that we are
still
> >>>> getting used to. Can you please refer to the
> >>>>> instructions on the following website for uploading data and
let us
> >> know
> >>>> if you can get it to work:
> >>>>>
> >>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804
>
> >>>>>>
> >>>>>> Dear Paul
> >>>>>>
> >>>>>> ftp -p ftp.rap.ucar.edu is not working. It is giving me
error
> >> message.
> >>>>>> ftp -p ftp.rap.ucar.edu
> >>>>>> ftp: p: unknown option
> >>>>>> Usage: ftp [-v] [-d] [-i] [-n] [-g] [-k realm] [-f] [-x] [-u]
[-t]
> >>>> [host]
> >>>>>> Regards
> >>>>>> Hashmi Fatima
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>> On Wed, Sep 29, 2010 at 12:22 AM, RAL HelpDesk {for Paul
Oldenburg}
> <
> >>>>>> met_help at ucar.edu> wrote:
> >>>>>>
> >>>>>>> Hashmi,
> >>>>>>>
> >>>>>>> Please try to invoke your ftp client using the -p command
line
> >>>> argument,
> >>>>>>> like so:
> >>>>>>>
> >>>>>>> ftp -p ftp.rap.ucar.edu
> >>>>>>> username = anonymous
> >>>>>>> password = "your email address"
> >>>>>>> cd incoming/irap/met_help/fatima_data
> >>>>>>> put "your files"
> >>>>>>>
> >>>>>>> Let us know how this approach goes.
> >>>>>>>
> >>>>>>> Thanks,
> >>>>>>>
> >>>>>>> Paul
> >>>>>>>
> >>>>>>>
> >>>>>>> Paul Oldenburg wrote:
> >>>>>>>> Hashmi,
> >>>>>>>>
> >>>>>>>> We are investigating a potential problem with our FTP site.
I
> will
> >>>> let
> >>>>>>> you know when it is fixed, and then we can
> >>>>>>>> proceed with addressing your wavelet_stat question.
> >>>>>>>>
> >>>>>>>> Thanks,
> >>>>>>>>
> >>>>>>>> Paul
> >>>>>>>>
> >>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >>>>>>>>>
> >>>>>>>>> Dear Paul
> >>>>>>>>>
> >>>>>>>>> I tried to upload my files at the FTP server you have told
me in
> >>>>>>> previous
> >>>>>>>>> mails. But it says " Permission denied on server".
> >>>>>>>>>
> >>>>>>>>> Regards
> >>>>>>>>> Hashmi Fatima
> >>>>>>>>>
> >>>>>>>>> On Mon, Sep 27, 2010 at 10:25 PM, RAL HelpDesk {for Paul
> Oldenburg}
> >> <
> >>>>>>>>> met_help at ucar.edu> wrote:
> >>>>>>>>>
> >>>>>>>>>> Hashmi,
> >>>>>>>>>>
> >>>>>>>>>> I (Paul) have returned to my met_help responsibilities
and I
> have
> >>>>>>> followed
> >>>>>>>>>> your interactions with John. In order to
> >>>>>>>>>> help you run stat_analysis, I would like to have you send
me
> some
> >> of
> >>>>>>> your
> >>>>>>>>>> data via FTP. First, can you give us the WPP
> >>>>>>>>>> output that you fed into wavelet_stat? Also, can you
give us
> your
> >>>>>>>>>> wavelet_stat output? Finally, please copy the error
> >>>>>>>>>> message that you get from stat_analysis in your response.
That
> >>>>>>> information
> >>>>>>>>>> will be very helpful.
> >>>>>>>>>>
> >>>>>>>>>> Concerning the masking region selected by wavelet_stat,
the
> entire
> >>>> grid
> >>>>>>>>>> cannot be processed because wavelet_stat can
> >>>>>>>>>> process only dyadic regions - that is square regions
whose sides
> >> are
> >>>> a
> >>>>>>>>>> power of 2 grid points (like 256 or 2048). By
> >>>>>>>>>> default, wavelet_stat should try to guess the best
regions to
> use
> >> in
> >>>>>>> your
> >>>>>>>>>> domain. If you send us your WPP output, we
> >>>>>>>>>> will be able to help you configure the best masking
region for
> >> your
> >>>>>>> case.
> >>>>>>>>>> We look forward to receiving your data. If you have any
> >> questions,
> >>>>>>> please
> >>>>>>>>>> let me know.
> >>>>>>>>>>
> >>>>>>>>>> Thanks,
> >>>>>>>>>>
> >>>>>>>>>> Paul
> >>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
> >>>>>>>>>>>
> >>>>>>>>>>> Dear Sir,
> >>>>>>>>>>>
> >>>>>>>>>>> Thanks for your help.Its working now.
> >>>>>>>>>>>
> >>>>>>>>>>> I have run wavelet tool for computation of RH. Now i
want to
> use
> >>>> the
> >>>>>>>>>> output
> >>>>>>>>>>> of wavelet for StatAnalysis. but its not working. I am
sending
> >> you
> >>>>>>>>>>> configuration file of wavelet and StatAnalysis. Could
you
> please
> >>>>>>> check.
> >>>>>>>>>>> One more thing I want to know that in the wavelet
output, why
> >> does
> >>>> it
> >>>>>>>>>> take
> >>>>>>>>>>> that much of area. I have not assigned any area for the
> >>>> verification
> >>>>>>>>>> except
> >>>>>>>>>>> full. I have masked the whole grid. I am attaching that
output
> >> too
> >>>> as
> >>>>>>>>>>> Wavelet_output.png.
> >>>>>>>>>>>
> >>>>>>>>>>> Thanks and Regards
> >>>>>>>>>>> Hashmi fatima
> >>>>>>>>>>>
> >>>>>>>>>>>
> >>>>>>>>>>> On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for John
Halley
> >>>> Gotway}
> >>>>>>> <
> >>>>>>>>>>> met_help at ucar.edu> wrote:
> >>>>>>>>>>>
> >>>>>>>>>>>> Hashmi,
> >>>>>>>>>>>>
> >>>>>>>>>>>> Paul will be out of the office for the next couple of
weeks,
> but
> >>>> I'll
> >>>>>>> be
> >>>>>>>>>>>> happy to help you with this issue.
> >>>>>>>>>>>>
> >>>>>>>>>>>> There are a couple of points I'd like to make:
> >>>>>>>>>>>>
> >>>>>>>>>>>> First, the NetCDF files you sent us are the NetCDF
output of
> >> WRF.
> >>>>>>> MET
> >>>>>>>>>> does
> >>>>>>>>>>>> not directly read the NetCDF output of WRF. Instead,
for
> >> METv2.0,
> >>>>>>> your
> >>>>>>>>>> WRF
> >>>>>>>>>>>> output must first be run through the
> >>>>>>>>>>>> WRF-PostProcessor (WPP). Details of running WPP can be
found
> in
> >>>>>>> chapter
> >>>>>>>>>> 7
> >>>>>>>>>>>> of the WRF user's guide:
> >>>>>>>>>>>>
> >>>>>>>>>>>>
> >>
> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
> >>>>>>>>>>>> WPP writes out GRIB data, which MET can read.
> >>>>>>>>>>>>
> >>>>>>>>>>>> Second, once you've run your data through the WRF
> PostProcessor,
> >>>> you
> >>>>>>> can
> >>>>>>>>>>>> select multiple thresholds for Grid-Stat to use. You
do so by
> >>>>>>> setting
> >>>>>>>>>> up
> >>>>>>>>>>>> the Grid-Stat configuration file something like this:
> >>>>>>>>>>>>
> >>>>>>>>>>>> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
> >>>>>>>>>>>> obs_field[] = [];
> >>>>>>>>>>>> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240
ge260
> >> ge280
> >>>>>>> ge300"
> >>>>>>>>>> ];
> >>>>>>>>>>>> obs_thresh[] = [];
> >>>>>>>>>>>>
> >>>>>>>>>>>> The settings above tell Grid-Stat to verify temperature
at
> >>>> 2-meters
> >>>>>>> as
> >>>>>>>>>> well
> >>>>>>>>>>>> as temperature at 500mb. For 2-meter temp, it'll use
the 2
> >>>>>>> thresholds
> >>>>>>>>>>>> listed (ge273 ge283) and for temp at 500mb, it'll use
> >>>>>>>>>>>> the 6 thresholds listed (ge200.0 ge220 ge240 ge260
ge280
> ge300).
> >>>>>>>>>>>>
> >>>>>>>>>>>> So the real task is running your WRF output files
through WPP.
> >>>>>>>>>>>>
> >>>>>>>>>>>> Lastly, I wanted to tell you that in METv3.0 (to be
released
> >>>>>>> shortly),
> >>>>>>>>>>>> we've added support for reading the NetCDF output of
the
> pinterp
> >>>> tool
> >>>>>>>>>> (which
> >>>>>>>>>>>> only processes ARW output). So once METv3.0 is
> >>>>>>>>>>>> available, you can use either the GRIB output of WPP or
the
> >> NetCDF
> >>>>>>>>>> output
> >>>>>>>>>>>> of the pinterp tool.
> >>>>>>>>>>>>
> >>>>>>>>>>>> Hope that helps.
> >>>>>>>>>>>>
> >>>>>>>>>>>> Thanks,
> >>>>>>>>>>>> John Halley Gotway
> >>>>>>>>>>>> met_help at ucar.edu
> >>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>>>>>>>> <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> Dear Sir,
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> I have put my files at the FTP site told by you.
Please
> check.
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> One more thing, i want to know that do i need to give
equal
> >>>> number
> >>>>>>> of
> >>>>>>>>>>>>> threshold values as fcst_field/obs_field? Suppose I
there is
> >> only
> >>>>>>> one
> >>>>>>>>>>>>> fcst_field/ob_field, but I want to check it for
various
> >> threshold
> >>>>>>>>>> values.
> >>>>>>>>>>>>> In some other case, I have given 5
fcst_field/obs_field and
> >>>>>>>>>> corresponding
> >>>>>>>>>>>> 5
> >>>>>>>>>>>>> threshold (fcst_thresh/obs_thresh) values. But I did
not get
> >>>>>>> anything
> >>>>>>>>>> in
> >>>>>>>>>>>>> the output. I have got netcdf file in the output but
nothing
> in
> >>>> the
> >>>>>>>>>> STAT
> >>>>>>>>>>>> and
> >>>>>>>>>>>>> .txt files.
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> Thanks and Regards
> >>>>>>>>>>>>> Hashmi Fatima
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for
Paul
> >>>> Oldenburg} <
> >>>>>>>>>>>>> met_help at ucar.edu> wrote:
> >>>>>>>>>>>>>
> >>>>>>>>>>>>>> Hashmi,
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> Can you please send us an example of your model and
obs data
> >>>> along
> >>>>>>>>>> with
> >>>>>>>>>>>>>> your grid_stat config file? If the data is too
> >>>>>>>>>>>>>> large to attach to an email, please upload it to our
FTP
> site:
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> ftp ftp.rap.ucar.edu
> >>>>>>>>>>>>>> username = anonymous
> >>>>>>>>>>>>>> password = "your email address"
> >>>>>>>>>>>>>> cd incoming/irap/met_help/fatima_data
> >>>>>>>>>>>>>> put "your files"
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> Thanks,
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> Paul
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>>>>>>>>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted
upon.
> >>>>>>>>>>>>>>> Transaction: Ticket created by
hashmifatima at gmail.com
> >>>>>>>>>>>>>>> Queue: met_help
> >>>>>>>>>>>>>>> Subject: netcdf file for grid_stat (grid code
error)
> >>>>>>>>>>>>>>> Owner: Nobody
> >>>>>>>>>>>>>>> Requestors: hashmifatima at gmail.com
> >>>>>>>>>>>>>>> Status: new
> >>>>>>>>>>>>>>> Ticket <URL:
> >>>>>>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
> >>>>>>>>>>>>>>> Dear Sir,
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>> I want to use grid_stat tool to verify 2m
temperature. Both
> >>>>>>>>>>>> (observation
> >>>>>>>>>>>>>> and
> >>>>>>>>>>>>>>> forecast) files are in netcdf format. when i am
defining
> >>>>>>> fcst_field,
> >>>>>>>>>> it
> >>>>>>>>>>>>>> is
> >>>>>>>>>>>>>>> asking for grid code. I am not using any grib file
for this
> >>>> test.
> >>>>>>>>>> What
> >>>>>>>>>>>>>>> should I do ?
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>> Thanks and Regards,
> >>>>
--------------------------------------------------------------------
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #40804] netcdf file for grid_stat (grid code error)
From: Paul Oldenburg
Time: Mon Oct 11 13:02:39 2010
Hashmi,
We are puzzled by your question. Which verification statistics are
you interested in? If you would like stat_analysis
to aggregate continuous statistics using certain criteria, try a
command line job like this:
$ stat_analysis -lookin . -job aggregate_stat -line_type SL1L2
-out_line_type CNT -vx_mask FULL -fcst_lead 240000
This job reads SL1L2 lines from all stat files in the current
directory and computes continuous statistics for the FULL
domain data with lead time 240000. SL1L2 lines contain partial sum
data information. The output of this command is
continuous statistic values like mean error (ME).
If you want to aggregate categorical or thresholded statistics, try a
command line job like this:
$ stat_analysis -lookin . -job aggregate_stat -line_type CTC
-out_line_type CTS -vx_mask FULL -fcst_lead 240000
-fcst_thresh ">=6000.000"
The output of this command is categorical statistic values like
probability of detection (PODY) for the threshold
>=6000.000.
Note that neither of the above two commands uses a configuration file.
It is often simplest and easiest to start
building stat_analysis jobs on the command line, so that you can test
the criteria. If you have any questions, please
let us know.
Thanks,
Paul
On 10/06/2010 03:34 AM, RAL HelpDesk {for Hashmi Fatima} wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>
> Dear Sir,
>
> Thanks for this information.
>
> Could you please tell me that what line_type should i use to get
output in
> stat_analysis?
>
> Regards
> Hashmi Fatima
>
> On Wed, Oct 6, 2010 at 2:14 AM, RAL HelpDesk {for Paul Oldenburg} <
> met_help at ucar.edu> wrote:
>
>> Hashmi,
>>
>> We inspected your stat_analysis jobs and saw that your criteria
resulted in
>> no MET output data being found. The reason
>> for this is that you are trying to select data in line types SL1L2
and CNT,
>> yet you are also specifying thresholds in
>> the same job. Line types SL1L2 and CNT contain statistics that are
not
>> divided among thresholds, because they are
>> continuous (as opposed to categorical or thresholded). Please look
>> carefully at your MET output and you will see that
>> these lines have a value of NA in the fcst_thresh line. If you
have any
>> questions, please let us know.
>>
>> Thanks,
>>
>> Paul
>>
>>
>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>
>>> Dear Sir
>>>
>>> I have configured STATconfig file for the previous case.
>>>
>>> But now I am trying STATanalysis with grid_stat output for HGT at
P500.
>> But
>>> its again giving some error message. I have uploaded the grid_stat
output
>>> files to your FTP server, which are input for STATanalysis case
(all
>>> grid_stat_HGT_P500-240000L_20101001_00000V_* files). STATconfig
file is
>>> STATAnalysisConfig_HGT. I have configurered the STATconfig file
according
>> to
>>> the input. but still unable to get the output.
>>> The error message which I am getting is srored in the file
>>> grid_output_stat_analysis_error.txt.
>>>
>>> Please check and let me know about the error.
>>>
>>> Thanks and Regards
>>> Hashmi Fatima
>>>
>>>
>>>
>>> On Sat, Oct 2, 2010 at 2:49 AM, RAL HelpDesk {for John Halley
Gotway} <
>>> met_help at ucar.edu> wrote:
>>>
>>>> Hashmi,
>>>>
>>>> Paul's out of the office, so hopefully I can provide some advice
on this
>>>> issue.
>>>>
>>>> It looks like you're using the STATAnalysis config file that's
included
>>>> with the MET test scripts. This is really only meant as an
example to
>> work
>>>> with the test data. You'll need to define the
>>>> analysis jobs explicitly to match the output data you have.
>>>>
>>>> Since you're running Wavelet-Stat, I assume you may want to
aggregate
>> the
>>>> output of Wavelet-Stat across multiple cases. Please keep in
mind that
>>>> Stat-Analysis can be run on the command line OR using
>>>> a config file. When you're trying to set up your analysis job at
the
>>>> beginning, I'd suggest using the command line option. For
example,
>> here's
>>>> the output I get for the following command:
>>>>
>>>> +++++++++++++++++++++++++++++++++++++++
>>>> bin/stat_analysis -lookin
wavelet_stat_000000L_20100920_000000V.stat
>> -job
>>>> aggregate -line_type ISC -vx_mask TILE_TOT -fcst_thresh
">=50.000"
>> -dump_row
>>>> aggr_ISC_50.stat -v 2
>>>> STAT Lines read = 330
>>>> STAT Lines retained = 330
>>>>
>>>> Processing Job 1: -job aggregate -vx_mask TILE_TOT -fcst_thresh
>=50.000
>>>> -line_type ISC -dump_row aggr_ISC_50.stat
>>>> Creating STAT output file "aggr_ISC_50.stat"
>>>> JOB_LIST: -job aggregate -vx_mask TILE_TOT -fcst_thresh
>=50.000
>>>> -line_type ISC -dump_row aggr_ISC_50.stat
>>>> COL_NAME: TOTAL TILE_DIM TILE_XLL TILE_YLL NSCALE ISCALE MSE
ISC
>>>> FENERGY2 OENERGY2 BASER FBIAS
>>>> ISC: 131072 NA NA NA 8 0 0.09870
0.63219
>>>> 0.16414 0.15513 0.15513 1.05808
>>>> 131072 NA NA NA 8 1 0.01761
0.40938
>>>> 0.00935 0.00921 0.15513 1.05808
>>>> 131072 NA NA NA 8 2 0.01881
0.36913
>>>> 0.01064 0.01070 0.15513 1.05808
>>>> 131072 NA NA NA 8 3 0.02372
0.20445
>>>> 0.01519 0.01605 0.15513 1.05808
>>>> 131072 NA NA NA 8 4 0.01836
0.38422
>>>> 0.01856 0.01899 0.15513 1.05808
>>>> 131072 NA NA NA 8 5 0.01390
0.53381
>>>> 0.01933 0.01968 0.15513 1.05808
>>>> 131072 NA NA NA 8 6 0.00524
0.82426
>>>> 0.02334 0.02096 0.15513 1.05808
>>>> 131072 NA NA NA 8 7 0.00082
0.97250
>>>> 0.02947 0.02577 0.15513 1.05808
>>>> 131072 NA NA NA 8 8 0.00011
0.99631
>>>> 0.00854 0.00764 0.15513 1.05808
>>>> 131072 NA NA NA 8 9 0.00013
0.99564
>>>> 0.02972 0.02614 0.15513 1.05808
>>>>
>>>> Job 1 used 10 out of 330 STAT lines.
>>>> +++++++++++++++++++++++++++++++++++++++
>>>>
>>>> Since you just sent us one Wavelet-Stat output file, that's all
we could
>>>> test with. But in this example, I'm aggregating together output
for the
>>>> ">=50.000" threshold. I'm guessing you'll want to
>>>> rerun this type of job for each of the thresholds you used in
>> Wavelet-Stat.
>>>> When you're first setting up jobs, it's a good idea to use the
>> "-dump_row"
>>>> option. That'll write all of the STAT lines the
>>>> job actually used to the file name you specify. That way, you
can make
>>>> sure that STAT-Analysis is operating on the exact subset of data
that
>> you
>>>> expect.
>>>>
>>>> Hope that helps clarify.
>>>>
>>>> John Halley Gotway
>>>> met_help at ucar.edu
>>>>
>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>
>>>>> Dear Paul,
>>>>>
>>>>> Sorry for the delay. I have put my files at ftp site in my
folder. I
>> have
>>>>> done the verufication through wavelet for 1 month data. Here I
am
>> sending
>>>>> few files.
>>>>>
>>>>> *input files (WPP output):*
>>>>> gfs.t00z.pgrbf00 (analysis file)
>>>>> gfs.t00z.pgrbf24 (forecast file)
>>>>>
>>>>> *wavelet output file (stat analysis input) :*
>>>>> wavelet_stat_000000L_20100920_000000V.nc
>>>>> wavelet_stat_000000L_20100920_000000V.ps
>>>>> wavelet_stat_000000L_20100920_000000V.stat
>>>>> wavelet_stat_000000L_20100920_000000V_isc.txt
>>>>>
>>>>> The error message which I am getting after running stat_analysis
tool
>>>> after
>>>>> using wavelet output is in file *error_msg.txt*.
>>>>>
>>>>> Thanks and Regards
>>>>> Hashmi Fatima
>>>>>
>>>>> On Wed, Sep 29, 2010 at 8:49 PM, RAL HelpDesk {for Paul
Oldenburg} <
>>>>> met_help at ucar.edu> wrote:
>>>>>
>>>>>> Hashmi,
>>>>>>
>>>>>> Here is the website, sorry:
>>>>>>
>>>>>> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>>>>>>
>>>>>> If you cannot get this to work, please let us know and we will
have to
>>>> make
>>>>>> other arrangements for you to send us your
>>>>>> data. Sorry for all the trouble!
>>>>>>
>>>>>> Thanks,
>>>>>>
>>>>>> Paul
>>>>>>
>>>>>>
>>>>>> RAL HelpDesk {for Paul Oldenburg} wrote:
>>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804
>
>>>>>>>
>>>>>>> Hashmi,
>>>>>>>
>>>>>>> We apparently have a new FTP server configuration that we are
still
>>>>>> getting used to. Can you please refer to the
>>>>>>> instructions on the following website for uploading data and
let us
>>>> know
>>>>>> if you can get it to work:
>>>>>>>
>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804
>
>>>>>>>>
>>>>>>>> Dear Paul
>>>>>>>>
>>>>>>>> ftp -p ftp.rap.ucar.edu is not working. It is giving me
error
>>>> message.
>>>>>>>> ftp -p ftp.rap.ucar.edu
>>>>>>>> ftp: p: unknown option
>>>>>>>> Usage: ftp [-v] [-d] [-i] [-n] [-g] [-k realm] [-f] [-x] [-u]
[-t]
>>>>>> [host]
>>>>>>>> Regards
>>>>>>>> Hashmi Fatima
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> On Wed, Sep 29, 2010 at 12:22 AM, RAL HelpDesk {for Paul
Oldenburg}
>> <
>>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>>
>>>>>>>>> Hashmi,
>>>>>>>>>
>>>>>>>>> Please try to invoke your ftp client using the -p command
line
>>>>>> argument,
>>>>>>>>> like so:
>>>>>>>>>
>>>>>>>>> ftp -p ftp.rap.ucar.edu
>>>>>>>>> username = anonymous
>>>>>>>>> password = "your email address"
>>>>>>>>> cd incoming/irap/met_help/fatima_data
>>>>>>>>> put "your files"
>>>>>>>>>
>>>>>>>>> Let us know how this approach goes.
>>>>>>>>>
>>>>>>>>> Thanks,
>>>>>>>>>
>>>>>>>>> Paul
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> Paul Oldenburg wrote:
>>>>>>>>>> Hashmi,
>>>>>>>>>>
>>>>>>>>>> We are investigating a potential problem with our FTP site.
I
>> will
>>>>>> let
>>>>>>>>> you know when it is fixed, and then we can
>>>>>>>>>> proceed with addressing your wavelet_stat question.
>>>>>>>>>>
>>>>>>>>>> Thanks,
>>>>>>>>>>
>>>>>>>>>> Paul
>>>>>>>>>>
>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>>>>>>>
>>>>>>>>>>> Dear Paul
>>>>>>>>>>>
>>>>>>>>>>> I tried to upload my files at the FTP server you have told
me in
>>>>>>>>> previous
>>>>>>>>>>> mails. But it says " Permission denied on server".
>>>>>>>>>>>
>>>>>>>>>>> Regards
>>>>>>>>>>> Hashmi Fatima
>>>>>>>>>>>
>>>>>>>>>>> On Mon, Sep 27, 2010 at 10:25 PM, RAL HelpDesk {for Paul
>> Oldenburg}
>>>> <
>>>>>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>>>>>
>>>>>>>>>>>> Hashmi,
>>>>>>>>>>>>
>>>>>>>>>>>> I (Paul) have returned to my met_help responsibilities
and I
>> have
>>>>>>>>> followed
>>>>>>>>>>>> your interactions with John. In order to
>>>>>>>>>>>> help you run stat_analysis, I would like to have you send
me
>> some
>>>> of
>>>>>>>>> your
>>>>>>>>>>>> data via FTP. First, can you give us the WPP
>>>>>>>>>>>> output that you fed into wavelet_stat? Also, can you
give us
>> your
>>>>>>>>>>>> wavelet_stat output? Finally, please copy the error
>>>>>>>>>>>> message that you get from stat_analysis in your response.
That
>>>>>>>>> information
>>>>>>>>>>>> will be very helpful.
>>>>>>>>>>>>
>>>>>>>>>>>> Concerning the masking region selected by wavelet_stat,
the
>> entire
>>>>>> grid
>>>>>>>>>>>> cannot be processed because wavelet_stat can
>>>>>>>>>>>> process only dyadic regions - that is square regions
whose sides
>>>> are
>>>>>> a
>>>>>>>>>>>> power of 2 grid points (like 256 or 2048). By
>>>>>>>>>>>> default, wavelet_stat should try to guess the best
regions to
>> use
>>>> in
>>>>>>>>> your
>>>>>>>>>>>> domain. If you send us your WPP output, we
>>>>>>>>>>>> will be able to help you configure the best masking
region for
>>>> your
>>>>>>>>> case.
>>>>>>>>>>>> We look forward to receiving your data. If you have any
>>>> questions,
>>>>>>>>> please
>>>>>>>>>>>> let me know.
>>>>>>>>>>>>
>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>
>>>>>>>>>>>> Paul
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
>>>>>>>>>>>>>
>>>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>>>
>>>>>>>>>>>>> Thanks for your help.Its working now.
>>>>>>>>>>>>>
>>>>>>>>>>>>> I have run wavelet tool for computation of RH. Now i
want to
>> use
>>>>>> the
>>>>>>>>>>>> output
>>>>>>>>>>>>> of wavelet for StatAnalysis. but its not working. I am
sending
>>>> you
>>>>>>>>>>>>> configuration file of wavelet and StatAnalysis. Could
you
>> please
>>>>>>>>> check.
>>>>>>>>>>>>> One more thing I want to know that in the wavelet
output, why
>>>> does
>>>>>> it
>>>>>>>>>>>> take
>>>>>>>>>>>>> that much of area. I have not assigned any area for the
>>>>>> verification
>>>>>>>>>>>> except
>>>>>>>>>>>>> full. I have masked the whole grid. I am attaching that
output
>>>> too
>>>>>> as
>>>>>>>>>>>>> Wavelet_output.png.
>>>>>>>>>>>>>
>>>>>>>>>>>>> Thanks and Regards
>>>>>>>>>>>>> Hashmi fatima
>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>> On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for John
Halley
>>>>>> Gotway}
>>>>>>>>> <
>>>>>>>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>>>>>>>
>>>>>>>>>>>>>> Hashmi,
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Paul will be out of the office for the next couple of
weeks,
>> but
>>>>>> I'll
>>>>>>>>> be
>>>>>>>>>>>>>> happy to help you with this issue.
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> There are a couple of points I'd like to make:
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> First, the NetCDF files you sent us are the NetCDF
output of
>>>> WRF.
>>>>>>>>> MET
>>>>>>>>>>>> does
>>>>>>>>>>>>>> not directly read the NetCDF output of WRF. Instead,
for
>>>> METv2.0,
>>>>>>>>> your
>>>>>>>>>>>> WRF
>>>>>>>>>>>>>> output must first be run through the
>>>>>>>>>>>>>> WRF-PostProcessor (WPP). Details of running WPP can be
found
>> in
>>>>>>>>> chapter
>>>>>>>>>>>> 7
>>>>>>>>>>>>>> of the WRF user's guide:
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>
>> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
>>>>>>>>>>>>>> WPP writes out GRIB data, which MET can read.
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Second, once you've run your data through the WRF
>> PostProcessor,
>>>>>> you
>>>>>>>>> can
>>>>>>>>>>>>>> select multiple thresholds for Grid-Stat to use. You
do so by
>>>>>>>>> setting
>>>>>>>>>>>> up
>>>>>>>>>>>>>> the Grid-Stat configuration file something like this:
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
>>>>>>>>>>>>>> obs_field[] = [];
>>>>>>>>>>>>>> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240
ge260
>>>> ge280
>>>>>>>>> ge300"
>>>>>>>>>>>> ];
>>>>>>>>>>>>>> obs_thresh[] = [];
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> The settings above tell Grid-Stat to verify temperature
at
>>>>>> 2-meters
>>>>>>>>> as
>>>>>>>>>>>> well
>>>>>>>>>>>>>> as temperature at 500mb. For 2-meter temp, it'll use
the 2
>>>>>>>>> thresholds
>>>>>>>>>>>>>> listed (ge273 ge283) and for temp at 500mb, it'll use
>>>>>>>>>>>>>> the 6 thresholds listed (ge200.0 ge220 ge240 ge260
ge280
>> ge300).
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> So the real task is running your WRF output files
through WPP.
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Lastly, I wanted to tell you that in METv3.0 (to be
released
>>>>>>>>> shortly),
>>>>>>>>>>>>>> we've added support for reading the NetCDF output of
the
>> pinterp
>>>>>> tool
>>>>>>>>>>>> (which
>>>>>>>>>>>>>> only processes ARW output). So once METv3.0 is
>>>>>>>>>>>>>> available, you can use either the GRIB output of WPP or
the
>>>> NetCDF
>>>>>>>>>>>> output
>>>>>>>>>>>>>> of the pinterp tool.
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Hope that helps.
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>>> John Halley Gotway
>>>>>>>>>>>>>> met_help at ucar.edu
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>>>>>>> <URL:
>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> I have put my files at the FTP site told by you.
Please
>> check.
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> One more thing, i want to know that do i need to give
equal
>>>>>> number
>>>>>>>>> of
>>>>>>>>>>>>>>> threshold values as fcst_field/obs_field? Suppose I
there is
>>>> only
>>>>>>>>> one
>>>>>>>>>>>>>>> fcst_field/ob_field, but I want to check it for
various
>>>> threshold
>>>>>>>>>>>> values.
>>>>>>>>>>>>>>> In some other case, I have given 5
fcst_field/obs_field and
>>>>>>>>>>>> corresponding
>>>>>>>>>>>>>> 5
>>>>>>>>>>>>>>> threshold (fcst_thresh/obs_thresh) values. But I did
not get
>>>>>>>>> anything
>>>>>>>>>>>> in
>>>>>>>>>>>>>>> the output. I have got netcdf file in the output but
nothing
>> in
>>>>>> the
>>>>>>>>>>>> STAT
>>>>>>>>>>>>>> and
>>>>>>>>>>>>>>> .txt files.
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Thanks and Regards
>>>>>>>>>>>>>>> Hashmi Fatima
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for
Paul
>>>>>> Oldenburg} <
>>>>>>>>>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Hashmi,
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Can you please send us an example of your model and
obs data
>>>>>> along
>>>>>>>>>>>> with
>>>>>>>>>>>>>>>> your grid_stat config file? If the data is too
>>>>>>>>>>>>>>>> large to attach to an email, please upload it to our
FTP
>> site:
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> ftp ftp.rap.ucar.edu
>>>>>>>>>>>>>>>> username = anonymous
>>>>>>>>>>>>>>>> password = "your email address"
>>>>>>>>>>>>>>>> cd incoming/irap/met_help/fatima_data
>>>>>>>>>>>>>>>> put "your files"
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Paul
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>>>>>>>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted
upon.
>>>>>>>>>>>>>>>>> Transaction: Ticket created by
hashmifatima at gmail.com
>>>>>>>>>>>>>>>>> Queue: met_help
>>>>>>>>>>>>>>>>> Subject: netcdf file for grid_stat (grid code
error)
>>>>>>>>>>>>>>>>> Owner: Nobody
>>>>>>>>>>>>>>>>> Requestors: hashmifatima at gmail.com
>>>>>>>>>>>>>>>>> Status: new
>>>>>>>>>>>>>>>>> Ticket <URL:
>>>>>>>>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
>>>>>>>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> I want to use grid_stat tool to verify 2m
temperature. Both
>>>>>>>>>>>>>> (observation
>>>>>>>>>>>>>>>> and
>>>>>>>>>>>>>>>>> forecast) files are in netcdf format. when i am
defining
>>>>>>>>> fcst_field,
>>>>>>>>>>>> it
>>>>>>>>>>>>>>>> is
>>>>>>>>>>>>>>>>> asking for grid code. I am not using any grib file
for this
>>>>>> test.
>>>>>>>>>>>> What
>>>>>>>>>>>>>>>>> should I do ?
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> Thanks and Regards,
>>>>>>
--------------------------------------------------------------------
>>
>
------------------------------------------------
Subject: netcdf file for grid_stat (grid code error)
From: Hashmi Fatima
Time: Mon Oct 11 13:12:19 2010
Thanks
I am able to get the output by reading few lines but not SL1L2. It
says "no
valid data".
Regards
On Tue, Oct 12, 2010 at 12:32 AM, RAL HelpDesk {for Paul Oldenburg} <
met_help at ucar.edu> wrote:
> Hashmi,
>
> We are puzzled by your question. Which verification statistics are
you
> interested in? If you would like stat_analysis
> to aggregate continuous statistics using certain criteria, try a
command
> line job like this:
>
> $ stat_analysis -lookin . -job aggregate_stat -line_type SL1L2
> -out_line_type CNT -vx_mask FULL -fcst_lead 240000
>
> This job reads SL1L2 lines from all stat files in the current
directory and
> computes continuous statistics for the FULL
> domain data with lead time 240000. SL1L2 lines contain partial sum
data
> information. The output of this command is
> continuous statistic values like mean error (ME).
>
> If you want to aggregate categorical or thresholded statistics, try
a
> command line job like this:
>
> $ stat_analysis -lookin . -job aggregate_stat -line_type CTC
-out_line_type
> CTS -vx_mask FULL -fcst_lead 240000
> -fcst_thresh ">=6000.000"
>
> The output of this command is categorical statistic values like
probability
> of detection (PODY) for the threshold
> >=6000.000.
>
> Note that neither of the above two commands uses a configuration
file. It
> is often simplest and easiest to start
> building stat_analysis jobs on the command line, so that you can
test the
> criteria. If you have any questions, please
> let us know.
>
> Thanks,
>
> Paul
>
>
>
>
> On 10/06/2010 03:34 AM, RAL HelpDesk {for Hashmi Fatima} wrote:
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >
> > Dear Sir,
> >
> > Thanks for this information.
> >
> > Could you please tell me that what line_type should i use to get
output
> in
> > stat_analysis?
> >
> > Regards
> > Hashmi Fatima
> >
> > On Wed, Oct 6, 2010 at 2:14 AM, RAL HelpDesk {for Paul Oldenburg}
<
> > met_help at ucar.edu> wrote:
> >
> >> Hashmi,
> >>
> >> We inspected your stat_analysis jobs and saw that your criteria
resulted
> in
> >> no MET output data being found. The reason
> >> for this is that you are trying to select data in line types
SL1L2 and
> CNT,
> >> yet you are also specifying thresholds in
> >> the same job. Line types SL1L2 and CNT contain statistics that
are not
> >> divided among thresholds, because they are
> >> continuous (as opposed to categorical or thresholded). Please
look
> >> carefully at your MET output and you will see that
> >> these lines have a value of NA in the fcst_thresh line. If you
have any
> >> questions, please let us know.
> >>
> >> Thanks,
> >>
> >> Paul
> >>
> >>
> >> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >>>
> >>> Dear Sir
> >>>
> >>> I have configured STATconfig file for the previous case.
> >>>
> >>> But now I am trying STATanalysis with grid_stat output for HGT
at P500.
> >> But
> >>> its again giving some error message. I have uploaded the
grid_stat
> output
> >>> files to your FTP server, which are input for STATanalysis case
(all
> >>> grid_stat_HGT_P500-240000L_20101001_00000V_* files). STATconfig
file is
> >>> STATAnalysisConfig_HGT. I have configurered the STATconfig file
> according
> >> to
> >>> the input. but still unable to get the output.
> >>> The error message which I am getting is srored in the file
> >>> grid_output_stat_analysis_error.txt.
> >>>
> >>> Please check and let me know about the error.
> >>>
> >>> Thanks and Regards
> >>> Hashmi Fatima
> >>>
> >>>
> >>>
> >>> On Sat, Oct 2, 2010 at 2:49 AM, RAL HelpDesk {for John Halley
Gotway} <
> >>> met_help at ucar.edu> wrote:
> >>>
> >>>> Hashmi,
> >>>>
> >>>> Paul's out of the office, so hopefully I can provide some
advice on
> this
> >>>> issue.
> >>>>
> >>>> It looks like you're using the STATAnalysis config file that's
> included
> >>>> with the MET test scripts. This is really only meant as an
example to
> >> work
> >>>> with the test data. You'll need to define the
> >>>> analysis jobs explicitly to match the output data you have.
> >>>>
> >>>> Since you're running Wavelet-Stat, I assume you may want to
aggregate
> >> the
> >>>> output of Wavelet-Stat across multiple cases. Please keep in
mind
> that
> >>>> Stat-Analysis can be run on the command line OR using
> >>>> a config file. When you're trying to set up your analysis job
at the
> >>>> beginning, I'd suggest using the command line option. For
example,
> >> here's
> >>>> the output I get for the following command:
> >>>>
> >>>> +++++++++++++++++++++++++++++++++++++++
> >>>> bin/stat_analysis -lookin
wavelet_stat_000000L_20100920_000000V.stat
> >> -job
> >>>> aggregate -line_type ISC -vx_mask TILE_TOT -fcst_thresh
">=50.000"
> >> -dump_row
> >>>> aggr_ISC_50.stat -v 2
> >>>> STAT Lines read = 330
> >>>> STAT Lines retained = 330
> >>>>
> >>>> Processing Job 1: -job aggregate -vx_mask TILE_TOT -fcst_thresh
> >=50.000
> >>>> -line_type ISC -dump_row aggr_ISC_50.stat
> >>>> Creating STAT output file "aggr_ISC_50.stat"
> >>>> JOB_LIST: -job aggregate -vx_mask TILE_TOT -fcst_thresh
>=50.000
> >>>> -line_type ISC -dump_row aggr_ISC_50.stat
> >>>> COL_NAME: TOTAL TILE_DIM TILE_XLL TILE_YLL NSCALE ISCALE MSE
ISC
> >>>> FENERGY2 OENERGY2 BASER FBIAS
> >>>> ISC: 131072 NA NA NA 8 0
0.09870
> 0.63219
> >>>> 0.16414 0.15513 0.15513 1.05808
> >>>> 131072 NA NA NA 8 1
0.01761
> 0.40938
> >>>> 0.00935 0.00921 0.15513 1.05808
> >>>> 131072 NA NA NA 8 2
0.01881
> 0.36913
> >>>> 0.01064 0.01070 0.15513 1.05808
> >>>> 131072 NA NA NA 8 3
0.02372
> 0.20445
> >>>> 0.01519 0.01605 0.15513 1.05808
> >>>> 131072 NA NA NA 8 4
0.01836
> 0.38422
> >>>> 0.01856 0.01899 0.15513 1.05808
> >>>> 131072 NA NA NA 8 5
0.01390
> 0.53381
> >>>> 0.01933 0.01968 0.15513 1.05808
> >>>> 131072 NA NA NA 8 6
0.00524
> 0.82426
> >>>> 0.02334 0.02096 0.15513 1.05808
> >>>> 131072 NA NA NA 8 7
0.00082
> 0.97250
> >>>> 0.02947 0.02577 0.15513 1.05808
> >>>> 131072 NA NA NA 8 8
0.00011
> 0.99631
> >>>> 0.00854 0.00764 0.15513 1.05808
> >>>> 131072 NA NA NA 8 9
0.00013
> 0.99564
> >>>> 0.02972 0.02614 0.15513 1.05808
> >>>>
> >>>> Job 1 used 10 out of 330 STAT lines.
> >>>> +++++++++++++++++++++++++++++++++++++++
> >>>>
> >>>> Since you just sent us one Wavelet-Stat output file, that's all
we
> could
> >>>> test with. But in this example, I'm aggregating together
output for
> the
> >>>> ">=50.000" threshold. I'm guessing you'll want to
> >>>> rerun this type of job for each of the thresholds you used in
> >> Wavelet-Stat.
> >>>> When you're first setting up jobs, it's a good idea to use the
> >> "-dump_row"
> >>>> option. That'll write all of the STAT lines the
> >>>> job actually used to the file name you specify. That way, you
can
> make
> >>>> sure that STAT-Analysis is operating on the exact subset of
data that
> >> you
> >>>> expect.
> >>>>
> >>>> Hope that helps clarify.
> >>>>
> >>>> John Halley Gotway
> >>>> met_help at ucar.edu
> >>>>
> >>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804
>
> >>>>>
> >>>>> Dear Paul,
> >>>>>
> >>>>> Sorry for the delay. I have put my files at ftp site in my
folder. I
> >> have
> >>>>> done the verufication through wavelet for 1 month data. Here I
am
> >> sending
> >>>>> few files.
> >>>>>
> >>>>> *input files (WPP output):*
> >>>>> gfs.t00z.pgrbf00 (analysis file)
> >>>>> gfs.t00z.pgrbf24 (forecast file)
> >>>>>
> >>>>> *wavelet output file (stat analysis input) :*
> >>>>> wavelet_stat_000000L_20100920_000000V.nc
> >>>>> wavelet_stat_000000L_20100920_000000V.ps
> >>>>> wavelet_stat_000000L_20100920_000000V.stat
> >>>>> wavelet_stat_000000L_20100920_000000V_isc.txt
> >>>>>
> >>>>> The error message which I am getting after running
stat_analysis tool
> >>>> after
> >>>>> using wavelet output is in file *error_msg.txt*.
> >>>>>
> >>>>> Thanks and Regards
> >>>>> Hashmi Fatima
> >>>>>
> >>>>> On Wed, Sep 29, 2010 at 8:49 PM, RAL HelpDesk {for Paul
Oldenburg} <
> >>>>> met_help at ucar.edu> wrote:
> >>>>>
> >>>>>> Hashmi,
> >>>>>>
> >>>>>> Here is the website, sorry:
> >>>>>>
> >>>>>> http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >>>>>>
> >>>>>> If you cannot get this to work, please let us know and we
will have
> to
> >>>> make
> >>>>>> other arrangements for you to send us your
> >>>>>> data. Sorry for all the trouble!
> >>>>>>
> >>>>>> Thanks,
> >>>>>>
> >>>>>> Paul
> >>>>>>
> >>>>>>
> >>>>>> RAL HelpDesk {for Paul Oldenburg} wrote:
> >>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >>>>>>>
> >>>>>>> Hashmi,
> >>>>>>>
> >>>>>>> We apparently have a new FTP server configuration that we
are still
> >>>>>> getting used to. Can you please refer to the
> >>>>>>> instructions on the following website for uploading data and
let us
> >>>> know
> >>>>>> if you can get it to work:
> >>>>>>>
> >>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
> >>>>>>>>
> >>>>>>>> Dear Paul
> >>>>>>>>
> >>>>>>>> ftp -p ftp.rap.ucar.edu is not working. It is giving me
error
> >>>> message.
> >>>>>>>> ftp -p ftp.rap.ucar.edu
> >>>>>>>> ftp: p: unknown option
> >>>>>>>> Usage: ftp [-v] [-d] [-i] [-n] [-g] [-k realm] [-f] [-x] [-
u] [-t]
> >>>>>> [host]
> >>>>>>>> Regards
> >>>>>>>> Hashmi Fatima
> >>>>>>>>
> >>>>>>>>
> >>>>>>>>
> >>>>>>>> On Wed, Sep 29, 2010 at 12:22 AM, RAL HelpDesk {for Paul
> Oldenburg}
> >> <
> >>>>>>>> met_help at ucar.edu> wrote:
> >>>>>>>>
> >>>>>>>>> Hashmi,
> >>>>>>>>>
> >>>>>>>>> Please try to invoke your ftp client using the -p command
line
> >>>>>> argument,
> >>>>>>>>> like so:
> >>>>>>>>>
> >>>>>>>>> ftp -p ftp.rap.ucar.edu
> >>>>>>>>> username = anonymous
> >>>>>>>>> password = "your email address"
> >>>>>>>>> cd incoming/irap/met_help/fatima_data
> >>>>>>>>> put "your files"
> >>>>>>>>>
> >>>>>>>>> Let us know how this approach goes.
> >>>>>>>>>
> >>>>>>>>> Thanks,
> >>>>>>>>>
> >>>>>>>>> Paul
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>> Paul Oldenburg wrote:
> >>>>>>>>>> Hashmi,
> >>>>>>>>>>
> >>>>>>>>>> We are investigating a potential problem with our FTP
site. I
> >> will
> >>>>>> let
> >>>>>>>>> you know when it is fixed, and then we can
> >>>>>>>>>> proceed with addressing your wavelet_stat question.
> >>>>>>>>>>
> >>>>>>>>>> Thanks,
> >>>>>>>>>>
> >>>>>>>>>> Paul
> >>>>>>>>>>
> >>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
> >>>>>>>>>>>
> >>>>>>>>>>> Dear Paul
> >>>>>>>>>>>
> >>>>>>>>>>> I tried to upload my files at the FTP server you have
told me
> in
> >>>>>>>>> previous
> >>>>>>>>>>> mails. But it says " Permission denied on server".
> >>>>>>>>>>>
> >>>>>>>>>>> Regards
> >>>>>>>>>>> Hashmi Fatima
> >>>>>>>>>>>
> >>>>>>>>>>> On Mon, Sep 27, 2010 at 10:25 PM, RAL HelpDesk {for Paul
> >> Oldenburg}
> >>>> <
> >>>>>>>>>>> met_help at ucar.edu> wrote:
> >>>>>>>>>>>
> >>>>>>>>>>>> Hashmi,
> >>>>>>>>>>>>
> >>>>>>>>>>>> I (Paul) have returned to my met_help responsibilities
and I
> >> have
> >>>>>>>>> followed
> >>>>>>>>>>>> your interactions with John. In order to
> >>>>>>>>>>>> help you run stat_analysis, I would like to have you
send me
> >> some
> >>>> of
> >>>>>>>>> your
> >>>>>>>>>>>> data via FTP. First, can you give us the WPP
> >>>>>>>>>>>> output that you fed into wavelet_stat? Also, can you
give us
> >> your
> >>>>>>>>>>>> wavelet_stat output? Finally, please copy the error
> >>>>>>>>>>>> message that you get from stat_analysis in your
response.
> That
> >>>>>>>>> information
> >>>>>>>>>>>> will be very helpful.
> >>>>>>>>>>>>
> >>>>>>>>>>>> Concerning the masking region selected by wavelet_stat,
the
> >> entire
> >>>>>> grid
> >>>>>>>>>>>> cannot be processed because wavelet_stat can
> >>>>>>>>>>>> process only dyadic regions - that is square regions
whose
> sides
> >>>> are
> >>>>>> a
> >>>>>>>>>>>> power of 2 grid points (like 256 or 2048). By
> >>>>>>>>>>>> default, wavelet_stat should try to guess the best
regions to
> >> use
> >>>> in
> >>>>>>>>> your
> >>>>>>>>>>>> domain. If you send us your WPP output, we
> >>>>>>>>>>>> will be able to help you configure the best masking
region for
> >>>> your
> >>>>>>>>> case.
> >>>>>>>>>>>> We look forward to receiving your data. If you have
any
> >>>> questions,
> >>>>>>>>> please
> >>>>>>>>>>>> let me know.
> >>>>>>>>>>>>
> >>>>>>>>>>>> Thanks,
> >>>>>>>>>>>>
> >>>>>>>>>>>> Paul
> >>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>>>>>>>> <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> Dear Sir,
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> Thanks for your help.Its working now.
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> I have run wavelet tool for computation of RH. Now i
want to
> >> use
> >>>>>> the
> >>>>>>>>>>>> output
> >>>>>>>>>>>>> of wavelet for StatAnalysis. but its not working. I am
> sending
> >>>> you
> >>>>>>>>>>>>> configuration file of wavelet and StatAnalysis. Could
you
> >> please
> >>>>>>>>> check.
> >>>>>>>>>>>>> One more thing I want to know that in the wavelet
output, why
> >>>> does
> >>>>>> it
> >>>>>>>>>>>> take
> >>>>>>>>>>>>> that much of area. I have not assigned any area for
the
> >>>>>> verification
> >>>>>>>>>>>> except
> >>>>>>>>>>>>> full. I have masked the whole grid. I am attaching
that
> output
> >>>> too
> >>>>>> as
> >>>>>>>>>>>>> Wavelet_output.png.
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> Thanks and Regards
> >>>>>>>>>>>>> Hashmi fatima
> >>>>>>>>>>>>>
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for
John
> Halley
> >>>>>> Gotway}
> >>>>>>>>> <
> >>>>>>>>>>>>> met_help at ucar.edu> wrote:
> >>>>>>>>>>>>>
> >>>>>>>>>>>>>> Hashmi,
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> Paul will be out of the office for the next couple of
weeks,
> >> but
> >>>>>> I'll
> >>>>>>>>> be
> >>>>>>>>>>>>>> happy to help you with this issue.
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> There are a couple of points I'd like to make:
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> First, the NetCDF files you sent us are the NetCDF
output of
> >>>> WRF.
> >>>>>>>>> MET
> >>>>>>>>>>>> does
> >>>>>>>>>>>>>> not directly read the NetCDF output of WRF. Instead,
for
> >>>> METv2.0,
> >>>>>>>>> your
> >>>>>>>>>>>> WRF
> >>>>>>>>>>>>>> output must first be run through the
> >>>>>>>>>>>>>> WRF-PostProcessor (WPP). Details of running WPP can
be
> found
> >> in
> >>>>>>>>> chapter
> >>>>>>>>>>>> 7
> >>>>>>>>>>>>>> of the WRF user's guide:
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>>
> >>>>
> >>
> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
> >>>>>>>>>>>>>> WPP writes out GRIB data, which MET can read.
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> Second, once you've run your data through the WRF
> >> PostProcessor,
> >>>>>> you
> >>>>>>>>> can
> >>>>>>>>>>>>>> select multiple thresholds for Grid-Stat to use. You
do so
> by
> >>>>>>>>> setting
> >>>>>>>>>>>> up
> >>>>>>>>>>>>>> the Grid-Stat configuration file something like this:
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
> >>>>>>>>>>>>>> obs_field[] = [];
> >>>>>>>>>>>>>> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240
ge260
> >>>> ge280
> >>>>>>>>> ge300"
> >>>>>>>>>>>> ];
> >>>>>>>>>>>>>> obs_thresh[] = [];
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> The settings above tell Grid-Stat to verify
temperature at
> >>>>>> 2-meters
> >>>>>>>>> as
> >>>>>>>>>>>> well
> >>>>>>>>>>>>>> as temperature at 500mb. For 2-meter temp, it'll use
the 2
> >>>>>>>>> thresholds
> >>>>>>>>>>>>>> listed (ge273 ge283) and for temp at 500mb, it'll use
> >>>>>>>>>>>>>> the 6 thresholds listed (ge200.0 ge220 ge240 ge260
ge280
> >> ge300).
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> So the real task is running your WRF output files
through
> WPP.
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> Lastly, I wanted to tell you that in METv3.0 (to be
released
> >>>>>>>>> shortly),
> >>>>>>>>>>>>>> we've added support for reading the NetCDF output of
the
> >> pinterp
> >>>>>> tool
> >>>>>>>>>>>> (which
> >>>>>>>>>>>>>> only processes ARW output). So once METv3.0 is
> >>>>>>>>>>>>>> available, you can use either the GRIB output of WPP
or the
> >>>> NetCDF
> >>>>>>>>>>>> output
> >>>>>>>>>>>>>> of the pinterp tool.
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> Hope that helps.
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> Thanks,
> >>>>>>>>>>>>>> John Halley Gotway
> >>>>>>>>>>>>>> met_help at ucar.edu
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>>>>>>>>>> <URL:
> >> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>> Dear Sir,
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>> I have put my files at the FTP site told by you.
Please
> >> check.
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>> One more thing, i want to know that do i need to
give equal
> >>>>>> number
> >>>>>>>>> of
> >>>>>>>>>>>>>>> threshold values as fcst_field/obs_field? Suppose I
there
> is
> >>>> only
> >>>>>>>>> one
> >>>>>>>>>>>>>>> fcst_field/ob_field, but I want to check it for
various
> >>>> threshold
> >>>>>>>>>>>> values.
> >>>>>>>>>>>>>>> In some other case, I have given 5
fcst_field/obs_field and
> >>>>>>>>>>>> corresponding
> >>>>>>>>>>>>>> 5
> >>>>>>>>>>>>>>> threshold (fcst_thresh/obs_thresh) values. But I
did not
> get
> >>>>>>>>> anything
> >>>>>>>>>>>> in
> >>>>>>>>>>>>>>> the output. I have got netcdf file in the output but
> nothing
> >> in
> >>>>>> the
> >>>>>>>>>>>> STAT
> >>>>>>>>>>>>>> and
> >>>>>>>>>>>>>>> .txt files.
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>> Thanks and Regards
> >>>>>>>>>>>>>>> Hashmi Fatima
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for
Paul
> >>>>>> Oldenburg} <
> >>>>>>>>>>>>>>> met_help at ucar.edu> wrote:
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>> Hashmi,
> >>>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>> Can you please send us an example of your model and
obs
> data
> >>>>>> along
> >>>>>>>>>>>> with
> >>>>>>>>>>>>>>>> your grid_stat config file? If the data is too
> >>>>>>>>>>>>>>>> large to attach to an email, please upload it to
our FTP
> >> site:
> >>>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>> ftp ftp.rap.ucar.edu
> >>>>>>>>>>>>>>>> username = anonymous
> >>>>>>>>>>>>>>>> password = "your email address"
> >>>>>>>>>>>>>>>> cd incoming/irap/met_help/fatima_data
> >>>>>>>>>>>>>>>> put "your files"
> >>>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>> Thanks,
> >>>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>> Paul
> >>>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
> >>>>>>>>>>>>>>>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted
upon.
> >>>>>>>>>>>>>>>>> Transaction: Ticket created by
hashmifatima at gmail.com
> >>>>>>>>>>>>>>>>> Queue: met_help
> >>>>>>>>>>>>>>>>> Subject: netcdf file for grid_stat (grid code
error)
> >>>>>>>>>>>>>>>>> Owner: Nobody
> >>>>>>>>>>>>>>>>> Requestors: hashmifatima at gmail.com
> >>>>>>>>>>>>>>>>> Status: new
> >>>>>>>>>>>>>>>>> Ticket <URL:
> >>>>>>>>>>>>
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
> >>>>>>>>>>>>>>>>> Dear Sir,
> >>>>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>>> I want to use grid_stat tool to verify 2m
temperature.
> Both
> >>>>>>>>>>>>>> (observation
> >>>>>>>>>>>>>>>> and
> >>>>>>>>>>>>>>>>> forecast) files are in netcdf format. when i am
defining
> >>>>>>>>> fcst_field,
> >>>>>>>>>>>> it
> >>>>>>>>>>>>>>>> is
> >>>>>>>>>>>>>>>>> asking for grid code. I am not using any grib file
for
> this
> >>>>>> test.
> >>>>>>>>>>>> What
> >>>>>>>>>>>>>>>>> should I do ?
> >>>>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>>> Thanks and Regards,
> >>>>>>
--------------------------------------------------------------------
> >>
> >
>
>
>
--
Hashmi Fatima
Research Scholar
Centre for Atmospheric Sciences
Indian Institute of Technology Delhi
Hauz Khas
New Delhi - 110016
INDIA
------------------------------------------------
Subject: netcdf file for grid_stat (grid code error)
From: Paul Oldenburg
Time: Mon Oct 11 15:27:55 2010
Hashmi,
I attached the data file that you sent, which includes SL1L2 data.
The stat_analysis command that I sent you in the
previous email gives the following output when run on this file.
Please study this example carefully, since it contains
all the information you need to run stat_analysis on your data. Make
sure to copy the data file into the directory
where you run stat analysis, so that it can be found.
$ stat_analysis -lookin . -job aggregate_stat -line_type SL1L2
-out_line_type CNT -vx_mask FULL -fcst_lead 240000
JOB_LIST: -job aggregate_stat -fcst_lead 240000 -vx_mask FULL
-line_type SL1L2 -out_line_type CNT -out_alpha
0.050000 -rank_corr_flag 1
COL_NAME: TOTAL FBAR FBAR_NCL FBAR_NCU FBAR_BCL FBAR_BCU
FSTDEV FSTDEV_NCL FSTDEV_NCU FSTDEV_BCL ...
CNT: 259920 5519.15261 5517.77172 5520.53350 NA NA
359.19459 358.22082 360.17371 NA ...
Thanks,
Paul
On 10/11/2010 01:12 PM, RAL HelpDesk {for Hashmi Fatima} wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>
> Thanks
>
> I am able to get the output by reading few lines but not SL1L2. It
says "no
> valid data".
>
> Regards
>
>
> On Tue, Oct 12, 2010 at 12:32 AM, RAL HelpDesk {for Paul Oldenburg}
<
> met_help at ucar.edu> wrote:
>
>> Hashmi,
>>
>> We are puzzled by your question. Which verification statistics are
you
>> interested in? If you would like stat_analysis
>> to aggregate continuous statistics using certain criteria, try a
command
>> line job like this:
>>
>> $ stat_analysis -lookin . -job aggregate_stat -line_type SL1L2
>> -out_line_type CNT -vx_mask FULL -fcst_lead 240000
>>
>> This job reads SL1L2 lines from all stat files in the current
directory and
>> computes continuous statistics for the FULL
>> domain data with lead time 240000. SL1L2 lines contain partial sum
data
>> information. The output of this command is
>> continuous statistic values like mean error (ME).
>>
>> If you want to aggregate categorical or thresholded statistics, try
a
>> command line job like this:
>>
>> $ stat_analysis -lookin . -job aggregate_stat -line_type CTC
-out_line_type
>> CTS -vx_mask FULL -fcst_lead 240000
>> -fcst_thresh ">=6000.000"
>>
>> The output of this command is categorical statistic values like
probability
>> of detection (PODY) for the threshold
>>> =6000.000.
>>
>> Note that neither of the above two commands uses a configuration
file. It
>> is often simplest and easiest to start
>> building stat_analysis jobs on the command line, so that you can
test the
>> criteria. If you have any questions, please
>> let us know.
>>
>> Thanks,
>>
>> Paul
>>
>>
>>
>>
>> On 10/06/2010 03:34 AM, RAL HelpDesk {for Hashmi Fatima} wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>
>>> Dear Sir,
>>>
>>> Thanks for this information.
>>>
>>> Could you please tell me that what line_type should i use to get
output
>> in
>>> stat_analysis?
>>>
>>> Regards
>>> Hashmi Fatima
>>>
>>> On Wed, Oct 6, 2010 at 2:14 AM, RAL HelpDesk {for Paul Oldenburg}
<
>>> met_help at ucar.edu> wrote:
>>>
>>>> Hashmi,
>>>>
>>>> We inspected your stat_analysis jobs and saw that your criteria
resulted
>> in
>>>> no MET output data being found. The reason
>>>> for this is that you are trying to select data in line types
SL1L2 and
>> CNT,
>>>> yet you are also specifying thresholds in
>>>> the same job. Line types SL1L2 and CNT contain statistics that
are not
>>>> divided among thresholds, because they are
>>>> continuous (as opposed to categorical or thresholded). Please
look
>>>> carefully at your MET output and you will see that
>>>> these lines have a value of NA in the fcst_thresh line. If you
have any
>>>> questions, please let us know.
>>>>
>>>> Thanks,
>>>>
>>>> Paul
>>>>
>>>>
>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>
>>>>> Dear Sir
>>>>>
>>>>> I have configured STATconfig file for the previous case.
>>>>>
>>>>> But now I am trying STATanalysis with grid_stat output for HGT
at P500.
>>>> But
>>>>> its again giving some error message. I have uploaded the
grid_stat
>> output
>>>>> files to your FTP server, which are input for STATanalysis case
(all
>>>>> grid_stat_HGT_P500-240000L_20101001_00000V_* files). STATconfig
file is
>>>>> STATAnalysisConfig_HGT. I have configurered the STATconfig file
>> according
>>>> to
>>>>> the input. but still unable to get the output.
>>>>> The error message which I am getting is srored in the file
>>>>> grid_output_stat_analysis_error.txt.
>>>>>
>>>>> Please check and let me know about the error.
>>>>>
>>>>> Thanks and Regards
>>>>> Hashmi Fatima
>>>>>
>>>>>
>>>>>
>>>>> On Sat, Oct 2, 2010 at 2:49 AM, RAL HelpDesk {for John Halley
Gotway} <
>>>>> met_help at ucar.edu> wrote:
>>>>>
>>>>>> Hashmi,
>>>>>>
>>>>>> Paul's out of the office, so hopefully I can provide some
advice on
>> this
>>>>>> issue.
>>>>>>
>>>>>> It looks like you're using the STATAnalysis config file that's
>> included
>>>>>> with the MET test scripts. This is really only meant as an
example to
>>>> work
>>>>>> with the test data. You'll need to define the
>>>>>> analysis jobs explicitly to match the output data you have.
>>>>>>
>>>>>> Since you're running Wavelet-Stat, I assume you may want to
aggregate
>>>> the
>>>>>> output of Wavelet-Stat across multiple cases. Please keep in
mind
>> that
>>>>>> Stat-Analysis can be run on the command line OR using
>>>>>> a config file. When you're trying to set up your analysis job
at the
>>>>>> beginning, I'd suggest using the command line option. For
example,
>>>> here's
>>>>>> the output I get for the following command:
>>>>>>
>>>>>> +++++++++++++++++++++++++++++++++++++++
>>>>>> bin/stat_analysis -lookin
wavelet_stat_000000L_20100920_000000V.stat
>>>> -job
>>>>>> aggregate -line_type ISC -vx_mask TILE_TOT -fcst_thresh
">=50.000"
>>>> -dump_row
>>>>>> aggr_ISC_50.stat -v 2
>>>>>> STAT Lines read = 330
>>>>>> STAT Lines retained = 330
>>>>>>
>>>>>> Processing Job 1: -job aggregate -vx_mask TILE_TOT -fcst_thresh
>>> =50.000
>>>>>> -line_type ISC -dump_row aggr_ISC_50.stat
>>>>>> Creating STAT output file "aggr_ISC_50.stat"
>>>>>> JOB_LIST: -job aggregate -vx_mask TILE_TOT -fcst_thresh
>=50.000
>>>>>> -line_type ISC -dump_row aggr_ISC_50.stat
>>>>>> COL_NAME: TOTAL TILE_DIM TILE_XLL TILE_YLL NSCALE ISCALE MSE
ISC
>>>>>> FENERGY2 OENERGY2 BASER FBIAS
>>>>>> ISC: 131072 NA NA NA 8 0
0.09870
>> 0.63219
>>>>>> 0.16414 0.15513 0.15513 1.05808
>>>>>> 131072 NA NA NA 8 1
0.01761
>> 0.40938
>>>>>> 0.00935 0.00921 0.15513 1.05808
>>>>>> 131072 NA NA NA 8 2
0.01881
>> 0.36913
>>>>>> 0.01064 0.01070 0.15513 1.05808
>>>>>> 131072 NA NA NA 8 3
0.02372
>> 0.20445
>>>>>> 0.01519 0.01605 0.15513 1.05808
>>>>>> 131072 NA NA NA 8 4
0.01836
>> 0.38422
>>>>>> 0.01856 0.01899 0.15513 1.05808
>>>>>> 131072 NA NA NA 8 5
0.01390
>> 0.53381
>>>>>> 0.01933 0.01968 0.15513 1.05808
>>>>>> 131072 NA NA NA 8 6
0.00524
>> 0.82426
>>>>>> 0.02334 0.02096 0.15513 1.05808
>>>>>> 131072 NA NA NA 8 7
0.00082
>> 0.97250
>>>>>> 0.02947 0.02577 0.15513 1.05808
>>>>>> 131072 NA NA NA 8 8
0.00011
>> 0.99631
>>>>>> 0.00854 0.00764 0.15513 1.05808
>>>>>> 131072 NA NA NA 8 9
0.00013
>> 0.99564
>>>>>> 0.02972 0.02614 0.15513 1.05808
>>>>>>
>>>>>> Job 1 used 10 out of 330 STAT lines.
>>>>>> +++++++++++++++++++++++++++++++++++++++
>>>>>>
>>>>>> Since you just sent us one Wavelet-Stat output file, that's all
we
>> could
>>>>>> test with. But in this example, I'm aggregating together
output for
>> the
>>>>>> ">=50.000" threshold. I'm guessing you'll want to
>>>>>> rerun this type of job for each of the thresholds you used in
>>>> Wavelet-Stat.
>>>>>> When you're first setting up jobs, it's a good idea to use the
>>>> "-dump_row"
>>>>>> option. That'll write all of the STAT lines the
>>>>>> job actually used to the file name you specify. That way, you
can
>> make
>>>>>> sure that STAT-Analysis is operating on the exact subset of
data that
>>>> you
>>>>>> expect.
>>>>>>
>>>>>> Hope that helps clarify.
>>>>>>
>>>>>> John Halley Gotway
>>>>>> met_help at ucar.edu
>>>>>>
>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804
>
>>>>>>>
>>>>>>> Dear Paul,
>>>>>>>
>>>>>>> Sorry for the delay. I have put my files at ftp site in my
folder. I
>>>> have
>>>>>>> done the verufication through wavelet for 1 month data. Here I
am
>>>> sending
>>>>>>> few files.
>>>>>>>
>>>>>>> *input files (WPP output):*
>>>>>>> gfs.t00z.pgrbf00 (analysis file)
>>>>>>> gfs.t00z.pgrbf24 (forecast file)
>>>>>>>
>>>>>>> *wavelet output file (stat analysis input) :*
>>>>>>> wavelet_stat_000000L_20100920_000000V.nc
>>>>>>> wavelet_stat_000000L_20100920_000000V.ps
>>>>>>> wavelet_stat_000000L_20100920_000000V.stat
>>>>>>> wavelet_stat_000000L_20100920_000000V_isc.txt
>>>>>>>
>>>>>>> The error message which I am getting after running
stat_analysis tool
>>>>>> after
>>>>>>> using wavelet output is in file *error_msg.txt*.
>>>>>>>
>>>>>>> Thanks and Regards
>>>>>>> Hashmi Fatima
>>>>>>>
>>>>>>> On Wed, Sep 29, 2010 at 8:49 PM, RAL HelpDesk {for Paul
Oldenburg} <
>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>
>>>>>>>> Hashmi,
>>>>>>>>
>>>>>>>> Here is the website, sorry:
>>>>>>>>
>>>>>>>> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>>>>>>>>
>>>>>>>> If you cannot get this to work, please let us know and we
will have
>> to
>>>>>> make
>>>>>>>> other arrangements for you to send us your
>>>>>>>> data. Sorry for all the trouble!
>>>>>>>>
>>>>>>>> Thanks,
>>>>>>>>
>>>>>>>> Paul
>>>>>>>>
>>>>>>>>
>>>>>>>> RAL HelpDesk {for Paul Oldenburg} wrote:
>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>>>>>
>>>>>>>>> Hashmi,
>>>>>>>>>
>>>>>>>>> We apparently have a new FTP server configuration that we
are still
>>>>>>>> getting used to. Can you please refer to the
>>>>>>>>> instructions on the following website for uploading data and
let us
>>>>>> know
>>>>>>>> if you can get it to work:
>>>>>>>>>
>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>>>>>>>>>>
>>>>>>>>>> Dear Paul
>>>>>>>>>>
>>>>>>>>>> ftp -p ftp.rap.ucar.edu is not working. It is giving me
error
>>>>>> message.
>>>>>>>>>> ftp -p ftp.rap.ucar.edu
>>>>>>>>>> ftp: p: unknown option
>>>>>>>>>> Usage: ftp [-v] [-d] [-i] [-n] [-g] [-k realm] [-f] [-x] [-
u] [-t]
>>>>>>>> [host]
>>>>>>>>>> Regards
>>>>>>>>>> Hashmi Fatima
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> On Wed, Sep 29, 2010 at 12:22 AM, RAL HelpDesk {for Paul
>> Oldenburg}
>>>> <
>>>>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>>>>
>>>>>>>>>>> Hashmi,
>>>>>>>>>>>
>>>>>>>>>>> Please try to invoke your ftp client using the -p command
line
>>>>>>>> argument,
>>>>>>>>>>> like so:
>>>>>>>>>>>
>>>>>>>>>>> ftp -p ftp.rap.ucar.edu
>>>>>>>>>>> username = anonymous
>>>>>>>>>>> password = "your email address"
>>>>>>>>>>> cd incoming/irap/met_help/fatima_data
>>>>>>>>>>> put "your files"
>>>>>>>>>>>
>>>>>>>>>>> Let us know how this approach goes.
>>>>>>>>>>>
>>>>>>>>>>> Thanks,
>>>>>>>>>>>
>>>>>>>>>>> Paul
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> Paul Oldenburg wrote:
>>>>>>>>>>>> Hashmi,
>>>>>>>>>>>>
>>>>>>>>>>>> We are investigating a potential problem with our FTP
site. I
>>>> will
>>>>>>>> let
>>>>>>>>>>> you know when it is fixed, and then we can
>>>>>>>>>>>> proceed with addressing your wavelet_stat question.
>>>>>>>>>>>>
>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>
>>>>>>>>>>>> Paul
>>>>>>>>>>>>
>>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
>>>>>>>>>>>>>
>>>>>>>>>>>>> Dear Paul
>>>>>>>>>>>>>
>>>>>>>>>>>>> I tried to upload my files at the FTP server you have
told me
>> in
>>>>>>>>>>> previous
>>>>>>>>>>>>> mails. But it says " Permission denied on server".
>>>>>>>>>>>>>
>>>>>>>>>>>>> Regards
>>>>>>>>>>>>> Hashmi Fatima
>>>>>>>>>>>>>
>>>>>>>>>>>>> On Mon, Sep 27, 2010 at 10:25 PM, RAL HelpDesk {for Paul
>>>> Oldenburg}
>>>>>> <
>>>>>>>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>>>>>>>
>>>>>>>>>>>>>> Hashmi,
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> I (Paul) have returned to my met_help responsibilities
and I
>>>> have
>>>>>>>>>>> followed
>>>>>>>>>>>>>> your interactions with John. In order to
>>>>>>>>>>>>>> help you run stat_analysis, I would like to have you
send me
>>>> some
>>>>>> of
>>>>>>>>>>> your
>>>>>>>>>>>>>> data via FTP. First, can you give us the WPP
>>>>>>>>>>>>>> output that you fed into wavelet_stat? Also, can you
give us
>>>> your
>>>>>>>>>>>>>> wavelet_stat output? Finally, please copy the error
>>>>>>>>>>>>>> message that you get from stat_analysis in your
response.
>> That
>>>>>>>>>>> information
>>>>>>>>>>>>>> will be very helpful.
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Concerning the masking region selected by wavelet_stat,
the
>>>> entire
>>>>>>>> grid
>>>>>>>>>>>>>> cannot be processed because wavelet_stat can
>>>>>>>>>>>>>> process only dyadic regions - that is square regions
whose
>> sides
>>>>>> are
>>>>>>>> a
>>>>>>>>>>>>>> power of 2 grid points (like 256 or 2048). By
>>>>>>>>>>>>>> default, wavelet_stat should try to guess the best
regions to
>>>> use
>>>>>> in
>>>>>>>>>>> your
>>>>>>>>>>>>>> domain. If you send us your WPP output, we
>>>>>>>>>>>>>> will be able to help you configure the best masking
region for
>>>>>> your
>>>>>>>>>>> case.
>>>>>>>>>>>>>> We look forward to receiving your data. If you have
any
>>>>>> questions,
>>>>>>>>>>> please
>>>>>>>>>>>>>> let me know.
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Paul
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>>>>>>> <URL:
>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Thanks for your help.Its working now.
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> I have run wavelet tool for computation of RH. Now i
want to
>>>> use
>>>>>>>> the
>>>>>>>>>>>>>> output
>>>>>>>>>>>>>>> of wavelet for StatAnalysis. but its not working. I am
>> sending
>>>>>> you
>>>>>>>>>>>>>>> configuration file of wavelet and StatAnalysis. Could
you
>>>> please
>>>>>>>>>>> check.
>>>>>>>>>>>>>>> One more thing I want to know that in the wavelet
output, why
>>>>>> does
>>>>>>>> it
>>>>>>>>>>>>>> take
>>>>>>>>>>>>>>> that much of area. I have not assigned any area for
the
>>>>>>>> verification
>>>>>>>>>>>>>> except
>>>>>>>>>>>>>>> full. I have masked the whole grid. I am attaching
that
>> output
>>>>>> too
>>>>>>>> as
>>>>>>>>>>>>>>> Wavelet_output.png.
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Thanks and Regards
>>>>>>>>>>>>>>> Hashmi fatima
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> On Tue, Sep 21, 2010 at 2:57 AM, RAL HelpDesk {for
John
>> Halley
>>>>>>>> Gotway}
>>>>>>>>>>> <
>>>>>>>>>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Hashmi,
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Paul will be out of the office for the next couple of
weeks,
>>>> but
>>>>>>>> I'll
>>>>>>>>>>> be
>>>>>>>>>>>>>>>> happy to help you with this issue.
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> There are a couple of points I'd like to make:
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> First, the NetCDF files you sent us are the NetCDF
output of
>>>>>> WRF.
>>>>>>>>>>> MET
>>>>>>>>>>>>>> does
>>>>>>>>>>>>>>>> not directly read the NetCDF output of WRF. Instead,
for
>>>>>> METv2.0,
>>>>>>>>>>> your
>>>>>>>>>>>>>> WRF
>>>>>>>>>>>>>>>> output must first be run through the
>>>>>>>>>>>>>>>> WRF-PostProcessor (WPP). Details of running WPP can
be
>> found
>>>> in
>>>>>>>>>>> chapter
>>>>>>>>>>>>>> 7
>>>>>>>>>>>>>>>> of the WRF user's guide:
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>
>>>>
>> http://www.dtcenter.org/wrf-
nmm/users/docs/user_guide/V3/users_guide_nmm_chap1-7.pdf
>>>>>>>>>>>>>>>> WPP writes out GRIB data, which MET can read.
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Second, once you've run your data through the WRF
>>>> PostProcessor,
>>>>>>>> you
>>>>>>>>>>> can
>>>>>>>>>>>>>>>> select multiple thresholds for Grid-Stat to use. You
do so
>> by
>>>>>>>>>>> setting
>>>>>>>>>>>>>> up
>>>>>>>>>>>>>>>> the Grid-Stat configuration file something like this:
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> fcst_field[] = [ "TMP/Z2", "TMP/P500" ];
>>>>>>>>>>>>>>>> obs_field[] = [];
>>>>>>>>>>>>>>>> fcst_thresh[] = [ "ge273 ge283", "ge200.0 ge220 ge240
ge260
>>>>>> ge280
>>>>>>>>>>> ge300"
>>>>>>>>>>>>>> ];
>>>>>>>>>>>>>>>> obs_thresh[] = [];
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> The settings above tell Grid-Stat to verify
temperature at
>>>>>>>> 2-meters
>>>>>>>>>>> as
>>>>>>>>>>>>>> well
>>>>>>>>>>>>>>>> as temperature at 500mb. For 2-meter temp, it'll use
the 2
>>>>>>>>>>> thresholds
>>>>>>>>>>>>>>>> listed (ge273 ge283) and for temp at 500mb, it'll use
>>>>>>>>>>>>>>>> the 6 thresholds listed (ge200.0 ge220 ge240 ge260
ge280
>>>> ge300).
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> So the real task is running your WRF output files
through
>> WPP.
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Lastly, I wanted to tell you that in METv3.0 (to be
released
>>>>>>>>>>> shortly),
>>>>>>>>>>>>>>>> we've added support for reading the NetCDF output of
the
>>>> pinterp
>>>>>>>> tool
>>>>>>>>>>>>>> (which
>>>>>>>>>>>>>>>> only processes ARW output). So once METv3.0 is
>>>>>>>>>>>>>>>> available, you can use either the GRIB output of WPP
or the
>>>>>> NetCDF
>>>>>>>>>>>>>> output
>>>>>>>>>>>>>>>> of the pinterp tool.
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Hope that helps.
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>>>>> John Halley Gotway
>>>>>>>>>>>>>>>> met_help at ucar.edu
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>>>>>>>>> <URL:
>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> I have put my files at the FTP site told by you.
Please
>>>> check.
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> One more thing, i want to know that do i need to
give equal
>>>>>>>> number
>>>>>>>>>>> of
>>>>>>>>>>>>>>>>> threshold values as fcst_field/obs_field? Suppose I
there
>> is
>>>>>> only
>>>>>>>>>>> one
>>>>>>>>>>>>>>>>> fcst_field/ob_field, but I want to check it for
various
>>>>>> threshold
>>>>>>>>>>>>>> values.
>>>>>>>>>>>>>>>>> In some other case, I have given 5
fcst_field/obs_field and
>>>>>>>>>>>>>> corresponding
>>>>>>>>>>>>>>>> 5
>>>>>>>>>>>>>>>>> threshold (fcst_thresh/obs_thresh) values. But I
did not
>> get
>>>>>>>>>>> anything
>>>>>>>>>>>>>> in
>>>>>>>>>>>>>>>>> the output. I have got netcdf file in the output but
>> nothing
>>>> in
>>>>>>>> the
>>>>>>>>>>>>>> STAT
>>>>>>>>>>>>>>>> and
>>>>>>>>>>>>>>>>> .txt files.
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> Thanks and Regards
>>>>>>>>>>>>>>>>> Hashmi Fatima
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> On Thu, Sep 16, 2010 at 8:27 PM, RAL HelpDesk {for
Paul
>>>>>>>> Oldenburg} <
>>>>>>>>>>>>>>>>> met_help at ucar.edu> wrote:
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> Hashmi,
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> Can you please send us an example of your model and
obs
>> data
>>>>>>>> along
>>>>>>>>>>>>>> with
>>>>>>>>>>>>>>>>>> your grid_stat config file? If the data is too
>>>>>>>>>>>>>>>>>> large to attach to an email, please upload it to
our FTP
>>>> site:
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> ftp ftp.rap.ucar.edu
>>>>>>>>>>>>>>>>>> username = anonymous
>>>>>>>>>>>>>>>>>> password = "your email address"
>>>>>>>>>>>>>>>>>> cd incoming/irap/met_help/fatima_data
>>>>>>>>>>>>>>>>>> put "your files"
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> Paul
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> RAL HelpDesk {for Hashmi Fatima} wrote:
>>>>>>>>>>>>>>>>>>> Thu Sep 16 01:44:32 2010: Request 40804 was acted
upon.
>>>>>>>>>>>>>>>>>>> Transaction: Ticket created by
hashmifatima at gmail.com
>>>>>>>>>>>>>>>>>>> Queue: met_help
>>>>>>>>>>>>>>>>>>> Subject: netcdf file for grid_stat (grid code
error)
>>>>>>>>>>>>>>>>>>> Owner: Nobody
>>>>>>>>>>>>>>>>>>> Requestors: hashmifatima at gmail.com
>>>>>>>>>>>>>>>>>>> Status: new
>>>>>>>>>>>>>>>>>>> Ticket <URL:
>>>>>>>>>>>>>>
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804>
>>>>>>>>>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> I want to use grid_stat tool to verify 2m
temperature.
>> Both
>>>>>>>>>>>>>>>> (observation
>>>>>>>>>>>>>>>>>> and
>>>>>>>>>>>>>>>>>>> forecast) files are in netcdf format. when i am
defining
>>>>>>>>>>> fcst_field,
>>>>>>>>>>>>>> it
>>>>>>>>>>>>>>>>>> is
>>>>>>>>>>>>>>>>>>> asking for grid code. I am not using any grib file
for
>> this
>>>>>>>> test.
>>>>>>>>>>>>>> What
>>>>>>>>>>>>>>>>>>> should I do ?
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> Thanks and Regards,
>>>>>>>>
--------------------------------------------------------------------
>>>>
>>>
>>
>>
>>
>
>
------------------------------------------------
Subject: netcdf file for grid_stat (grid code error)
From: Paul Oldenburg
Time: Mon Oct 11 15:27:55 2010
VERSION MODEL FCST_LEAD FCST_VALID_BEG FCST_VALID_END OBS_LEAD
OBS_VALID_BEG OBS_VALID_END FCST_VAR FCST_LEV OBS_VAR OBS_LEV
OBTYPE VX_MASK INTERP_MTHD INTERP_PNTS FCST_THRESH OBS_THRESH
COV_THRESH ALPHA LINE_TYPE
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL UW_MEAN 1 >=4000.000 >=4000.000 NA
NA FHO 259920 1.00000 1.00000 1.00000
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL UW_MEAN 1 >=4000.000 >=4000.000 NA
NA CTC 259920 259920 0 0 0
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL UW_MEAN 1 >=4000.000 >=4000.000 NA
0.10000 CTS 259920 1.00000 NA NA 1.00000
1.00000 1.00000 NA NA 1.00000 1.00000
1.00000 NA NA 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 NA NA 1.00000 1.00000 NA
NA NA NA NA NA NA NA NA NA
0.00000 NA NA 0.00000 0.00000 1.00000 NA NA
1.00000 1.00000 NA NA NA NA NA NA
NA NA NA NA NA NA NA
NA NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL UW_MEAN 1 >=5000.000 >=5000.000 NA
NA FHO 259920 0.85282 0.84121 0.85661
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL UW_MEAN 1 >=5000.000 >=5000.000 NA
NA CTC 259920 218648 3018 4003 34251
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL UW_MEAN 1 >=5000.000 >=5000.000 NA
0.10000 CTS 259920 0.85661 0.85549 0.85773 0.85552
0.85769 0.85282 0.85169 0.85395 0.85167 0.85391
0.97299 0.97248 0.97349 0.97245 0.97351 0.99558
0.99496 0.99618 0.98202 0.98161 0.98242 0.98156 0.98251
0.91902 0.91815 0.91988 0.91675 0.92137 0.08098 0.08012 0.08185
0.07863 0.08325 0.01362 0.01328 0.01396 0.01320 0.01401 0.96889
0.96834 0.96942 0.96826 0.96949 0.80381 0.80021 0.80718 0.90104
0.89982 0.90227 0.89873 0.90341 0.89124 0.88902
0.89330 619.88901 595.11330 645.69619 594.43222 646.17586
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL UW_MEAN 1 >=6000.000 >=6000.000 NA
NA FHO 259920 0.00000 0.00000 0.00000
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL UW_MEAN 1 >=6000.000 >=6000.000 NA
NA CTC 259920 0 0 0 259920
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL UW_MEAN 1 >=6000.000 >=6000.000 NA
0.10000 CTS 259920 0.00000 NA NA 0.00000
0.00000 0.00000 NA NA 0.00000 0.00000
1.00000 NA NA 1.00000 1.00000 NA NA
NA NA NA NA NA NA 1.00000 NA
NA 1.00000 1.00000 0.00000 NA NA 0.00000 0.00000
NA NA NA NA NA NA NA NA
NA NA NA NA NA NA NA NA
NA NA NA NA NA NA NA
NA NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL UW_MEAN 1 NA NA NA
0.10000 CNT 259920 5519.15261 5517.99373 5520.31148
5518.02204 5520.31532 359.19464 358.37720 360.01612 358.57697
359.79625 5520.20370 5519.04661 5521.36079 5519.10994 5521.32710
358.63971 357.82354 359.45992 358.05075 359.19672 0.98569 0.98560
0.98579 0.98559 0.98579 0.98377 -0.01483 259920 247508 247769
-1.05109 -1.24697 -0.85521 -1.22887 -0.86262 60.71305 60.57488
60.85190 60.52111 60.90417 0.99981 0.99978 0.99984 41.24238 41.10163
41.38311 3687.16463 3663.59913 3710.15439 3686.05983 3662.79055
3709.30395 60.72203 60.52767 60.91104 -70.89876 -71.49146 -70.30994
-24.50288 -24.80222 -24.28280 -1.50682 -1.60556 -1.40632 22.49984
22.18806 22.89759 75.01137 74.49398 75.68927
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL UW_MEAN 1 NA NA NA
NA SL1L2 259920 5519.15261 5520.20370 30593824.73777
30590065.79306 30601270.84710
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS UW_MEAN 1 >=4000.000 >=4000.000 NA
NA FHO 57520 1.00000 1.00000 1.00000
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS UW_MEAN 1 >=4000.000 >=4000.000 NA
NA CTC 57520 57520 0 0 0
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS UW_MEAN 1 >=4000.000 >=4000.000 NA
0.10000 CTS 57520 1.00000 0.99983 1.00012 1.00000
1.00000 1.00000 0.99983 1.00012 1.00000 1.00000
1.00000 0.99983 1.00012 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 0.99983 1.00012 1.00000 1.00000 NA
NA NA NA NA NA NA NA NA NA
0.00000 -0.00012 0.00017 0.00000 0.00000 1.00000 0.99983 1.00012
1.00000 1.00000 NA NA NA NA NA NA
NA NA NA NA NA NA NA
NA NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS UW_MEAN 1 >=5000.000 >=5000.000 NA
NA FHO 57520 1.00000 1.00000 1.00000
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS UW_MEAN 1 >=5000.000 >=5000.000 NA
NA CTC 57520 57520 0 0 0
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS UW_MEAN 1 >=5000.000 >=5000.000 NA
0.10000 CTS 57520 1.00000 0.99983 1.00012 1.00000
1.00000 1.00000 0.99983 1.00012 1.00000 1.00000
1.00000 0.99983 1.00012 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 0.99983 1.00012 1.00000 1.00000 NA
NA NA NA NA NA NA NA NA NA
0.00000 -0.00012 0.00017 0.00000 0.00000 1.00000 0.99983 1.00012
1.00000 1.00000 NA NA NA NA NA NA
NA NA NA NA NA NA NA
NA NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS UW_MEAN 1 >=6000.000 >=6000.000 NA
NA FHO 57520 0.00000 0.00000 0.00000
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS UW_MEAN 1 >=6000.000 >=6000.000 NA
NA CTC 57520 0 0 0 57520
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS UW_MEAN 1 >=6000.000 >=6000.000 NA
0.10000 CTS 57520 0.00000 -0.00012 0.00017 0.00000
0.00000 0.00000 -0.00012 0.00017 0.00000 0.00000
1.00000 0.99983 1.00012 1.00000 1.00000 NA NA
NA NA NA NA NA NA 1.00000 0.99983
1.00012 1.00000 1.00000 0.00000 -0.00012 0.00017 0.00000 0.00000
NA NA NA NA NA NA NA NA
NA NA NA NA NA NA NA NA
NA NA NA NA NA NA NA
NA NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS UW_MEAN 1 NA NA NA
0.10000 CNT 57520 5856.25195 5856.13512 5856.36878
5856.12889 5856.36357 17.03438 16.95220 17.11743 16.96187
17.10808 5857.00297 5856.90019 5857.10574 5856.89463 5857.10487
14.98538 14.91309 15.05844 14.90314 15.06071 0.84442 0.84244
0.84637 0.84225 0.84674 0.85920 0.67499 57520 56504 56498
-0.75102 -0.81374 -0.68830 -0.81208 -0.68763 9.14486 9.10074 9.18944
9.09081 9.19593 0.99987 0.99986 0.99988 7.18049 7.14210 7.22013
84.19101 83.21696 85.12348 83.62698 82.64135 84.56369
9.17557 9.12233 9.22624 -12.78700 -12.89244 -12.60786 -6.59916
-6.69820 -6.49865 0.00418 -0.08896 0.10124 5.58566 5.49863
5.60265 9.59227 9.49960 9.61027
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS UW_MEAN 1 NA NA NA
NA SL1L2 57520 5856.25195 5857.00297 34300300.61836
34295977.08733 34304708.34040
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 9 >=4000.000 >=4000.000 >=0.500
NA NBRCTC 257762 257762 0 0 0
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 9 >=4000.000 >=4000.000 >=0.500
0.10000 NBRCTS 257762 1.00000 NA NA 1.00000
1.00000 1.00000 NA NA 1.00000 1.00000
1.00000 NA NA 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 NA NA 1.00000 1.00000 NA
NA NA NA NA NA NA NA NA NA
0.00000 NA NA 0.00000 0.00000 1.00000 NA NA
1.00000 1.00000 NA NA NA NA NA NA
NA NA NA NA NA NA NA
NA NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 9 >=4000.000 >=4000.000 NA
0.10000 NBRCNT 257762 -0.00000 -0.00000 -0.00000 1.00000
1.00000 1.00000
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 9 >=4000.000 >=4000.000 >=0.500
NA NBRCTC 57360 57360 0 0 0
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 9 >=4000.000 >=4000.000 >=0.500
0.10000 NBRCTS 57360 1.00000 0.99981 1.00014 1.00000
1.00000 1.00000 0.99981 1.00014 1.00000 1.00000
1.00000 0.99981 1.00014 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 0.99981 1.00014 1.00000 1.00000 NA
NA NA NA NA NA NA NA NA NA
0.00000 -0.00014 0.00019 0.00000 0.00000 1.00000 0.99981 1.00014
1.00000 1.00000 NA NA NA NA NA NA
NA NA NA NA NA NA NA
NA NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 9 >=4000.000 >=4000.000 NA
0.10000 NBRCNT 57360 -0.00000 -0.00000 -0.00000 1.00000
1.00000 1.00000
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 9 >=5000.000 >=5000.000 >=0.500
NA NBRCTC 257762 217338 2996 4000 33428
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 9 >=5000.000 >=5000.000 >=0.500
0.10000 NBRCTS 257762 0.85869 0.85757 0.85981 0.85755
0.85982 0.85480 0.85366 0.85593 0.85365 0.85593
0.97286 0.97234 0.97337 0.97231 0.97339 0.99546
0.99482 0.99611 0.98193 0.98151 0.98234 0.98145 0.98239
0.91775 0.91686 0.91862 0.91528 0.91995 0.08225 0.08138 0.08314
0.08005 0.08472 0.01360 0.01325 0.01396 0.01320 0.01402 0.96881
0.96826 0.96936 0.96819 0.96941 0.80088 0.79719 0.80426 0.89967
0.89844 0.90091 0.89715 0.90191 0.88943 0.88715
0.89151 606.23954 581.94447 631.54888 581.03099 631.14947
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 9 >=5000.000 >=5000.000 NA
0.10000 NBRCNT 257762 0.02400 0.02352 0.02445 0.98596
0.98570 0.98625
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 9 >=5000.000 >=5000.000 >=0.500
NA NBRCTC 57360 57360 0 0 0
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 9 >=5000.000 >=5000.000 >=0.500
0.10000 NBRCTS 57360 1.00000 0.99981 1.00014 1.00000
1.00000 1.00000 0.99981 1.00014 1.00000 1.00000
1.00000 0.99981 1.00014 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 0.99981 1.00014 1.00000 1.00000 NA
NA NA NA NA NA NA NA NA NA
0.00000 -0.00014 0.00019 0.00000 0.00000 1.00000 0.99981 1.00014
1.00000 1.00000 NA NA NA NA NA NA
NA NA NA NA NA NA NA
NA NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 9 >=5000.000 >=5000.000 NA
0.10000 NBRCNT 57360 -0.00000 -0.00000 -0.00000 1.00000
1.00000 1.00000
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 9 >=6000.000 >=6000.000 >=0.500
NA NBRCTC 257762 0 0 0 257762
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 9 >=6000.000 >=6000.000 >=0.500
0.10000 NBRCTS 257762 0.00000 NA NA 0.00000
0.00000 0.00000 NA NA 0.00000 0.00000
1.00000 NA NA 1.00000 1.00000 NA NA
NA NA NA NA NA NA 1.00000 NA
NA 1.00000 1.00000 0.00000 NA NA 0.00000 0.00000
NA NA NA NA NA NA NA NA
NA NA NA NA NA NA NA NA
NA NA NA NA NA NA NA
NA NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 9 >=6000.000 >=6000.000 NA
0.10000 NBRCNT 257762 -0.00000 -0.00000 -0.00000 NA
NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 9 >=6000.000 >=6000.000 >=0.500
NA NBRCTC 57360 0 0 0 57360
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 9 >=6000.000 >=6000.000 >=0.500
0.10000 NBRCTS 57360 0.00000 -0.00014 0.00019 0.00000
0.00000 0.00000 -0.00014 0.00019 0.00000 0.00000
1.00000 0.99981 1.00014 1.00000 1.00000 NA NA
NA NA NA NA NA NA 1.00000 0.99981
1.00014 1.00000 1.00000 0.00000 -0.00014 0.00019 0.00000 0.00000
NA NA NA NA NA NA NA NA
NA NA NA NA NA NA NA NA
NA NA NA NA NA NA NA
NA NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 9 >=6000.000 >=6000.000 NA
0.10000 NBRCNT 57360 -0.00000 -0.00000 -0.00000 NA
NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 25 >=4000.000 >=4000.000 >=0.500
NA NBRCTC 255612 255612 0 0 0
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 25 >=4000.000 >=4000.000 >=0.500
0.10000 NBRCTS 255612 1.00000 NA NA 1.00000
1.00000 1.00000 NA NA 1.00000 1.00000
1.00000 NA NA 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 NA NA 1.00000 1.00000 NA
NA NA NA NA NA NA NA NA NA
0.00000 NA NA 0.00000 0.00000 1.00000 NA NA
1.00000 1.00000 NA NA NA NA NA NA
NA NA NA NA NA NA NA
NA NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 25 >=4000.000 >=4000.000 NA
0.10000 NBRCNT 255612 -0.00000 -0.00000 -0.00000 1.00000
1.00000 1.00000
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 25 >=4000.000 >=4000.000 >=0.500
NA NBRCTC 57200 57200 0 0 0
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 25 >=4000.000 >=4000.000 >=0.500
0.10000 NBRCTS 57200 1.00000 0.99979 1.00017 1.00000
1.00000 1.00000 0.99979 1.00017 1.00000 1.00000
1.00000 0.99979 1.00017 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 0.99979 1.00017 1.00000 1.00000 NA
NA NA NA NA NA NA NA NA NA
0.00000 -0.00017 0.00021 0.00000 0.00000 1.00000 0.99979 1.00017
1.00000 1.00000 NA NA NA NA NA NA
NA NA NA NA NA NA NA
NA NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 25 >=4000.000 >=4000.000 NA
0.10000 NBRCNT 57200 -0.00000 -0.00000 -0.00000 1.00000
1.00000 1.00000
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 25 >=5000.000 >=5000.000 >=0.500
NA NBRCTC 255612 216026 2963 3995 32628
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 25 >=5000.000 >=5000.000 >=0.500
0.10000 NBRCTS 255612 0.86076 0.85964 0.86188 0.85962
0.86191 0.85672 0.85559 0.85786 0.85558 0.85792
0.97278 0.97226 0.97329 0.97226 0.97327 0.99531
0.99471 0.99594 0.98184 0.98142 0.98226 0.98135 0.98231
0.91675 0.91585 0.91764 0.91436 0.91897 0.08325 0.08236 0.08415
0.08103 0.08564 0.01353 0.01318 0.01390 0.01314 0.01394 0.96880
0.96824 0.96935 0.96821 0.96936 0.79824 0.79473 0.80152 0.89859
0.89734 0.89985 0.89614 0.90090 0.88780 0.88563
0.88983 595.45376 571.50053 620.41095 571.80328 619.00176
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 25 >=5000.000 >=5000.000 NA
0.10000 NBRCNT 255612 0.02198 0.02155 0.02240 0.98716
0.98690 0.98742
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 25 >=5000.000 >=5000.000 >=0.500
NA NBRCTC 57200 57200 0 0 0
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 25 >=5000.000 >=5000.000 >=0.500
0.10000 NBRCTS 57200 1.00000 0.99979 1.00017 1.00000
1.00000 1.00000 0.99979 1.00017 1.00000 1.00000
1.00000 0.99979 1.00017 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 0.99979 1.00017 1.00000 1.00000 NA
NA NA NA NA NA NA NA NA NA
0.00000 -0.00017 0.00021 0.00000 0.00000 1.00000 0.99979 1.00017
1.00000 1.00000 NA NA NA NA NA NA
NA NA NA NA NA NA NA
NA NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 25 >=5000.000 >=5000.000 NA
0.10000 NBRCNT 57200 -0.00000 -0.00000 -0.00000 1.00000
1.00000 1.00000
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 25 >=6000.000 >=6000.000 >=0.500
NA NBRCTC 255612 0 0 0 255612
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 25 >=6000.000 >=6000.000 >=0.500
0.10000 NBRCTS 255612 0.00000 NA NA 0.00000
0.00000 0.00000 NA NA 0.00000 0.00000
1.00000 NA NA 1.00000 1.00000 NA NA
NA NA NA NA NA NA 1.00000 NA
NA 1.00000 1.00000 0.00000 NA NA 0.00000 0.00000
NA NA NA NA NA NA NA NA
NA NA NA NA NA NA NA NA
NA NA NA NA NA NA NA
NA NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS FULL NBRHD 25 >=6000.000 >=6000.000 NA
0.10000 NBRCNT 255612 -0.00000 -0.00000 -0.00000 NA
NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 25 >=6000.000 >=6000.000 >=0.500
NA NBRCTC 57200 0 0 0 57200
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 25 >=6000.000 >=6000.000 >=0.500
0.10000 NBRCTS 57200 0.00000 -0.00017 0.00021 0.00000
0.00000 0.00000 -0.00017 0.00021 0.00000 0.00000
1.00000 0.99979 1.00017 1.00000 1.00000 NA NA
NA NA NA NA NA NA 1.00000 0.99979
1.00017 1.00000 1.00000 0.00000 -0.00017 0.00021 0.00000 0.00000
NA NA NA NA NA NA NA NA
NA NA NA NA NA NA NA NA
NA NA NA NA NA NA NA
NA NA NA
V2.0 GFS 240000 20101001_000000 20101001_000000 000000
20100930_000000 20100930_000000 HGT P500 HGT P500
ANALYS TROPICS NBRHD 25 >=6000.000 >=6000.000 NA
0.10000 NBRCNT 57200 -0.00000 -0.00000 -0.00000 NA
NA NA
------------------------------------------------
More information about the Met_help
mailing list