[Met_help] [rt.rap.ucar.edu #40804] History for netcdf file for grid_stat (grid code error)

RAL HelpDesk {for John Halley Gotway} met_help at ucar.edu
Fri Oct 29 13:30:29 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

------------------------------------------------
Subject: netcdf file for grid_stat (grid code error)
From: Hashmi Fatima
Time: Mon Oct 18 06:29:05 2010

Dear Sir,

Thanks for all you help. I wanted to use NetCDF output of WRF for
grid_stat.
So I have install METv3.0. I have configuered config file accoding to
NetCDF
file format. But it is giving a message to me " NetCDF: Attribute not
found". I think there is something lacking in defining the
"fcst_field" in
config file, which I am not able to figure out.  I am sending you the
files
which I am using i.e script, config file, and NetCDF file. I have put
these
file to your ftp server in my folder
"incoming/irap/met_help/fatima_data".

Please check and let me know, how should I write fcst_field for NetCDF
files.

Thanks and Regards
Hashmi Fatima

On Tue, Oct 12, 2010 at 2:58 AM, RAL HelpDesk {for Paul Oldenburg} <
met_help at ucar.edu> wrote:

> According to our records, your request has been resolved. If you
have any
> further questions or concerns, please respond to this message.
>



--
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] Resolved: netcdf file for grid_stat (grid code error)
From: John Halley Gotway
Time: Mon Oct 18 09:05:17 2010

Hashimi,

METv3.0 does NOT read the NetCDF output directly from WRF.  You must
post-process the data first.  We recommend post-processing it using
the WRF Post Processor (WPP) which writes it's output in GRIB
format.  However, as of METv3.0, you may also post process it using
the pinterp tool (only works for WRF-ARW):
   http://www.mmm.ucar.edu/wrf/OnLineTutorial/Tools/p_interp.htm

METv3.0 can read the GRIB output of WPP or the NetCDF output of
pinterp.

Thanks,
John

On 10/18/2010 06:29 AM, RAL HelpDesk {for Hashmi Fatima} wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=40804 >
>
> Dear Sir,
>
> Thanks for all you help. I wanted to use NetCDF output of WRF for
grid_stat.
> So I have install METv3.0. I have configuered config file accoding
to NetCDF
> file format. But it is giving a message to me " NetCDF: Attribute
not
> found". I think there is something lacking in defining the
"fcst_field" in
> config file, which I am not able to figure out.  I am sending you
the files
> which I am using i.e script, config file, and NetCDF file. I have
put these
> file to your ftp server in my folder
"incoming/irap/met_help/fatima_data".
>
> Please check and let me know, how should I write fcst_field for
NetCDF
> files.
>
> Thanks and Regards
> Hashmi Fatima
>
> On Tue, Oct 12, 2010 at 2:58 AM, RAL HelpDesk {for Paul Oldenburg} <
> met_help at ucar.edu> wrote:
>
>> According to our records, your request has been resolved. If you
have any
>> further questions or concerns, please respond to this message.
>>
>
>
>

------------------------------------------------


More information about the Met_help mailing list