[Met_help] [rt.rap.ucar.edu #71092] History for Sir,I got a problem using Wavelet_stat

John Halley Gotway via RT met_help at ucar.edu
Thu Mar 26 09:37:11 MDT 2015


----------------------------------------------------------------
  Initial Request
----------------------------------------------------------------

Good morning sir:
          I'm a new user of MET. Recently I'd like to use it to perform a scale decomposition. My question is that my forecast data is from ARW-WRF(wrfout) which contain many variable and levels, I'd like to do it with RAINC+RAINNC. It seems that MET cannot directly run with that.
         How should I do to transform my wrfout file or is there any other way?
         Thanks for your consideration !




 

----------------------------------------------------------------
  Complete Ticket History
----------------------------------------------------------------

Subject: Sir,I got a problem using Wavelet_stat
From: John Halley Gotway
Time: Fri Mar 13 16:09:20 2015

Hello,

MET cannot read raw WRFOUT files directly.  Instead, they must first
be
post-processed using either the Unified PostProcessor (UPP) or the
pinterp
utility for WRF-ARW.  We strongly recommend using UPP which writes
GRIB
output files that MET was designed to handle.  While MET should be
able to
read pinterp output, it can only handle variables defined on a
non-staggered dimension.  That means that MET won't read winds from a
pinterp file since they're staggered.

Information about UPP can be found here:
   http://www.dtcenter.org/wrf-nmm/users/overview/upp_overview.php

It can read output from WRF-NMM and WRF-ARW.

I'd suggest giving that a shot and let us know what other questions
arise.
Questions specifically about UPP should be directed to
wrfhelp at ucar.edu.

Thanks,
John Halley Gotway
met_help at ucar.edu

On Thu, Mar 12, 2015 at 9:05 PM, Zhuang from Nuist via RT
<met_help at ucar.edu
> wrote:

>
> Thu Mar 12 21:05:41 2015: Request 71092 was acted upon.
> Transaction: Ticket created by zxrxz3212009 at 163.com
>        Queue: met_help
>      Subject: Sir,I got a problem using Wavelet_stat
>        Owner: Nobody
>   Requestors: zxrxz3212009 at 163.com
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71092 >
>
>
> Good morning sir:
>           I'm a new user of MET. Recently I'd like to use it to
perform a
> scale decomposition. My question is that my forecast data is from
> ARW-WRF(wrfout) which contain many variable and levels, I'd like to
do it
> with RAINC+RAINNC. It seems that MET cannot directly run with that.
>          How should I do to transform my wrfout file or is there any
other
> way?
>          Thanks for your consideration !
>
>
>
>
>
>

------------------------------------------------
Subject: Sir,I got a problem using Wavelet_stat
From: Zhuang from Nuist
Time: Mon Mar 16 08:28:43 2015

Hello,
      In my test I transform wrfout using UPP and the total
precipitation variable in wrfprs_d02.00 is A_PCP_GDS3_SFC_acc.
      But when I use met to perform a wavelet_stat, the error is ERROR
:
ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation 'R_H_GDS3_HTGL' for table version 2
ERROR  :
 I have no idea what does that mean.

 Thanks for your consideration !






At 2015-03-14 06:09:20, "John Halley Gotway via RT"
<met_help at ucar.edu> wrote:
>Hello,
>
>MET cannot read raw WRFOUT files directly.  Instead, they must first
be
>post-processed using either the Unified PostProcessor (UPP) or the
pinterp
>utility for WRF-ARW.  We strongly recommend using UPP which writes
GRIB
>output files that MET was designed to handle.  While MET should be
able to
>read pinterp output, it can only handle variables defined on a
>non-staggered dimension.  That means that MET won't read winds from a
>pinterp file since they're staggered.
>
>Information about UPP can be found here:
>   http://www.dtcenter.org/wrf-nmm/users/overview/upp_overview.php
>
>It can read output from WRF-NMM and WRF-ARW.
>
>I'd suggest giving that a shot and let us know what other questions
arise.
>Questions specifically about UPP should be directed to
wrfhelp at ucar.edu.
>
>Thanks,
>John Halley Gotway
>met_help at ucar.edu
>
>On Thu, Mar 12, 2015 at 9:05 PM, Zhuang from Nuist via RT
<met_help at ucar.edu
>> wrote:
>
>>
>> Thu Mar 12 21:05:41 2015: Request 71092 was acted upon.
>> Transaction: Ticket created by zxrxz3212009 at 163.com
>>        Queue: met_help
>>      Subject: Sir,I got a problem using Wavelet_stat
>>        Owner: Nobody
>>   Requestors: zxrxz3212009 at 163.com
>>       Status: new
>>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71092 >
>>
>>
>> Good morning sir:
>>           I'm a new user of MET. Recently I'd like to use it to
perform a
>> scale decomposition. My question is that my forecast data is from
>> ARW-WRF(wrfout) which contain many variable and levels, I'd like to
do it
>> with RAINC+RAINNC. It seems that MET cannot directly run with that.
>>          How should I do to transform my wrfout file or is there
any other
>> way?
>>          Thanks for your consideration !
>>
>>
>>
>>
>>
>>
>

------------------------------------------------
Subject: Sir,I got a problem using Wavelet_stat
From: John Halley Gotway
Time: Mon Mar 16 10:42:42 2015

Do you have the wgrib utility available on your system?  It can be
used to
dump the contents of a GRIB file.  If so, please run that file through
wgrib and send me the output:
   wgrib wrfprs_d02.00 > wrfprs_d02.00.dump

Then send me that file "wrfprs_d02.00.dump" along with a description
of
what field you'd like to run through Wavelet-Stat.  That'll help me
advise
how you should set up the configuration file.

Thanks,
John Halley Gotway
met_help at ucar.edu

On Monday, March 16, 2015, Zhuang from Nuist via RT
<met_help at ucar.edu>
wrote:

>
> Mon Mar 16 08:39:05 2015: Request 71115 was acted upon.
> Transaction: Ticket created by zxrxz3212009 at 163.com
>        Queue: met_help
>      Subject: Sir,I got a problem using Wavelet_stat
>        Owner: Nobody
>   Requestors: zxrxz3212009 at 163.com
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71115 >
>
>
> Hello,
>       In my test I transform wrfout using UPP and the total
precipitation
> variable in wrfprs_d02.00 is A_PCP_GDS3_SFC_acc.
>       But when I use met to perform a wavelet_stat, and set       "
name
>      = "A_PCP_GDS3_SFC_acc";"   in WaveletStatConfig,
> the error is ERROR  :
> ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
> 'R_H_GDS3_HTGL' for table version 2
> ERROR  :
>       It can be seen from file WaveletStatConfig_default that
> "name="APCP"" and "level=[ "A03" ]" ,but the sample data
> "data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212" doesn't have
a var
> named "APCP".
>      I have no idea what does that mean and I do like to know how to
fill
> "name" and "level" in WaveletStatConfig.
>
>
>     Thanks for your consideration!
>
>
>

------------------------------------------------
Subject: Sir,I got a problem using Wavelet_stat
From: Zhuang from Nuist
Time: Tue Mar 17 21:15:01 2015

Hello sir:
         Thanks very much for your help and I have successfully run
wavelrt_stat by its default config.
         Now I have some questions:
1. is there any way for me to only decompose the fcst field without a
matched obs_field like USER_GUIDE figure9.2.
As we know, the picture out from default config perform a scale
decomposition of the difference between fcst_field and obs_field.


2.how can I perform a forecast energy squared picture like figure 9.6
in user guide.


Thanks for your kind consideration!






At 2015-03-17 00:42:42, "John Halley Gotway via RT"
<met_help at ucar.edu> wrote:
>Do you have the wgrib utility available on your system?  It can be
used to
>dump the contents of a GRIB file.  If so, please run that file
through
>wgrib and send me the output:
>   wgrib wrfprs_d02.00 > wrfprs_d02.00.dump
>
>Then send me that file "wrfprs_d02.00.dump" along with a description
of
>what field you'd like to run through Wavelet-Stat.  That'll help me
advise
>how you should set up the configuration file.
>
>Thanks,
>John Halley Gotway
>met_help at ucar.edu
>
>On Monday, March 16, 2015, Zhuang from Nuist via RT
<met_help at ucar.edu>
>wrote:
>
>>
>> Mon Mar 16 08:39:05 2015: Request 71115 was acted upon.
>> Transaction: Ticket created by zxrxz3212009 at 163.com
>>        Queue: met_help
>>      Subject: Sir,I got a problem using Wavelet_stat
>>        Owner: Nobody
>>   Requestors: zxrxz3212009 at 163.com
>>       Status: new
>>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71115 >
>>
>>
>> Hello,
>>       In my test I transform wrfout using UPP and the total
precipitation
>> variable in wrfprs_d02.00 is A_PCP_GDS3_SFC_acc.
>>       But when I use met to perform a wavelet_stat, and set       "
name
>>      = "A_PCP_GDS3_SFC_acc";"   in WaveletStatConfig,
>> the error is ERROR  :
>> ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
>> 'R_H_GDS3_HTGL' for table version 2
>> ERROR  :
>>       It can be seen from file WaveletStatConfig_default that
>> "name="APCP"" and "level=[ "A03" ]" ,but the sample data
>> "data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212" doesn't
have a var
>> named "APCP".
>>      I have no idea what does that mean and I do like to know how
to fill
>> "name" and "level" in WaveletStatConfig.
>>
>>
>>     Thanks for your consideration!
>>
>>
>>
>

------------------------------------------------
Subject: Sir,I got a problem using Wavelet_stat
From: John Halley Gotway
Time: Wed Mar 18 11:17:49 2015

Zhuang,

I see that you'd like to use Wavelet-Stat to decompose a single field.
Here's one thing you could do... When you run Wavelet-Stat pass in the
forecast file as both the forecast and observation file.  In the
configuration file, set the categorical threshold (cat_thresh) for the
forecast field in the way you'd like.  Then set cat_thresh for the
observation field very, very high.  That'll result in the forecast
field
being what you want and the observation field being all 0's.  So when
you
subtract the forecast minus observation, you'll end up with the
forecast
field.

For example...

fcst = {
   field = [
      {
        name       = "APCP";
        level      = [ "A12" ];
        cat_thresh = [ >0.0 ];
      }
   ];
};
obs = {
   field = [
      {
        name       = "APCP";
        level      = [ "A12" ];
        cat_thresh = [ >99999 ];
      }
   ];
};

This will define events in the forecast field but no event in the
observation field.  So fcst - obs = fcst.

Make sense?

Thanks,
John

On Tue, Mar 17, 2015 at 9:29 PM, Zhuang from Nuist via RT
<met_help at ucar.edu
> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71092 >
>
> Hello sir:
>          Thanks very much for your help and I have successfully run
> wavelrt_stat by its default config.
>          Now I have some questions:
> 1. is there any way for me to only decompose the fcst field without
a
> matched obs_field like USER_GUIDE figure9.2.
> As we know, the picture out from default config perform a scale
> decomposition of the difference between fcst_field and obs_field.
>
>
> 2.how can I perform a forecast energy squared picture like figure
9.6  in
> user guide.
>
>
> Thanks for your kind consideration!
>
>
>
>
>
>
> At 2015-03-17 00:42:42, "John Halley Gotway via RT"
<met_help at ucar.edu>
> wrote:
> >Do you have the wgrib utility available on your system?  It can be
used to
> >dump the contents of a GRIB file.  If so, please run that file
through
> >wgrib and send me the output:
> >   wgrib wrfprs_d02.00 > wrfprs_d02.00.dump
> >
> >Then send me that file "wrfprs_d02.00.dump" along with a
description of
> >what field you'd like to run through Wavelet-Stat.  That'll help me
advise
> >how you should set up the configuration file.
> >
> >Thanks,
> >John Halley Gotway
> >met_help at ucar.edu
> >
> >On Monday, March 16, 2015, Zhuang from Nuist via RT
<met_help at ucar.edu>
> >wrote:
> >
> >>
> >> Mon Mar 16 08:39:05 2015: Request 71115 was acted upon.
> >> Transaction: Ticket created by zxrxz3212009 at 163.com
> >>        Queue: met_help
> >>      Subject: Sir,I got a problem using Wavelet_stat
> >>        Owner: Nobody
> >>   Requestors: zxrxz3212009 at 163.com
> >>       Status: new
> >>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71115 >
> >>
> >>
> >> Hello,
> >>       In my test I transform wrfout using UPP and the total
> precipitation
> >> variable in wrfprs_d02.00 is A_PCP_GDS3_SFC_acc.
> >>       But when I use met to perform a wavelet_stat, and set
" name
> >>      = "A_PCP_GDS3_SFC_acc";"   in WaveletStatConfig,
> >> the error is ERROR  :
> >> ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
> >> 'R_H_GDS3_HTGL' for table version 2
> >> ERROR  :
> >>       It can be seen from file WaveletStatConfig_default that
> >> "name="APCP"" and "level=[ "A03" ]" ,but the sample data
> >> "data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212" doesn't
have a
> var
> >> named "APCP".
> >>      I have no idea what does that mean and I do like to know how
to
> fill
> >> "name" and "level" in WaveletStatConfig.
> >>
> >>
> >>     Thanks for your consideration!
> >>
> >>
> >>
> >
>
>

------------------------------------------------
Subject: Sir,I got a problem using Wavelet_stat
From: John Halley Gotway
Time: Wed Mar 18 11:22:48 2015

Also, the images shown in Figure 9.6 were created by plotting the
gridded
NetCDF output from the Wavelet-Stat tool.  It was probably done using
NCL
or R.  I don't have the script that was used, but it should be pretty
straight-forward for you to plot the contents of those NetCDF output
files.

You could start by just displaying that data using the ncview utility.

Thanks,
John

On Wed, Mar 18, 2015 at 11:17 AM, John Halley Gotway <johnhg at ucar.edu>
wrote:

> Zhuang,
>
> I see that you'd like to use Wavelet-Stat to decompose a single
field.
> Here's one thing you could do... When you run Wavelet-Stat pass in
the
> forecast file as both the forecast and observation file.  In the
> configuration file, set the categorical threshold (cat_thresh) for
the
> forecast field in the way you'd like.  Then set cat_thresh for the
> observation field very, very high.  That'll result in the forecast
field
> being what you want and the observation field being all 0's.  So
when you
> subtract the forecast minus observation, you'll end up with the
forecast
> field.
>
> For example...
>
> fcst = {
>    field = [
>       {
>         name       = "APCP";
>         level      = [ "A12" ];
>         cat_thresh = [ >0.0 ];
>       }
>    ];
> };
> obs = {
>    field = [
>       {
>         name       = "APCP";
>         level      = [ "A12" ];
>         cat_thresh = [ >99999 ];
>       }
>    ];
> };
>
> This will define events in the forecast field but no event in the
> observation field.  So fcst - obs = fcst.
>
> Make sense?
>
> Thanks,
> John
>
> On Tue, Mar 17, 2015 at 9:29 PM, Zhuang from Nuist via RT <
> met_help at ucar.edu> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71092 >
>>
>> Hello sir:
>>          Thanks very much for your help and I have successfully run
>> wavelrt_stat by its default config.
>>          Now I have some questions:
>> 1. is there any way for me to only decompose the fcst field without
a
>> matched obs_field like USER_GUIDE figure9.2.
>> As we know, the picture out from default config perform a scale
>> decomposition of the difference between fcst_field and obs_field.
>>
>>
>> 2.how can I perform a forecast energy squared picture like figure
9.6  in
>> user guide.
>>
>>
>> Thanks for your kind consideration!
>>
>>
>>
>>
>>
>>
>> At 2015-03-17 00:42:42, "John Halley Gotway via RT"
<met_help at ucar.edu>
>> wrote:
>> >Do you have the wgrib utility available on your system?  It can be
used
>> to
>> >dump the contents of a GRIB file.  If so, please run that file
through
>> >wgrib and send me the output:
>> >   wgrib wrfprs_d02.00 > wrfprs_d02.00.dump
>> >
>> >Then send me that file "wrfprs_d02.00.dump" along with a
description of
>> >what field you'd like to run through Wavelet-Stat.  That'll help
me
>> advise
>> >how you should set up the configuration file.
>> >
>> >Thanks,
>> >John Halley Gotway
>> >met_help at ucar.edu
>> >
>> >On Monday, March 16, 2015, Zhuang from Nuist via RT
<met_help at ucar.edu>
>> >wrote:
>> >
>> >>
>> >> Mon Mar 16 08:39:05 2015: Request 71115 was acted upon.
>> >> Transaction: Ticket created by zxrxz3212009 at 163.com
>> >>        Queue: met_help
>> >>      Subject: Sir,I got a problem using Wavelet_stat
>> >>        Owner: Nobody
>> >>   Requestors: zxrxz3212009 at 163.com
>> >>       Status: new
>> >>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71115
>> >
>> >>
>> >>
>> >> Hello,
>> >>       In my test I transform wrfout using UPP and the total
>> precipitation
>> >> variable in wrfprs_d02.00 is A_PCP_GDS3_SFC_acc.
>> >>       But when I use met to perform a wavelet_stat, and set
"
>> name
>> >>      = "A_PCP_GDS3_SFC_acc";"   in WaveletStatConfig,
>> >> the error is ERROR  :
>> >> ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field
>> abbreviation
>> >> 'R_H_GDS3_HTGL' for table version 2
>> >> ERROR  :
>> >>       It can be seen from file WaveletStatConfig_default that
>> >> "name="APCP"" and "level=[ "A03" ]" ,but the sample data
>> >> "data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212" doesn't
have a
>> var
>> >> named "APCP".
>> >>      I have no idea what does that mean and I do like to know
how to
>> fill
>> >> "name" and "level" in WaveletStatConfig.
>> >>
>> >>
>> >>     Thanks for your consideration!
>> >>
>> >>
>> >>
>> >
>>
>>
>

------------------------------------------------
Subject: Sir,I got a problem using Wavelet_stat
From: Zhuang from Nuist
Time: Fri Mar 20 02:39:15 2015

Hello sir::
        First thanks a lot for your kind consideration and it's really
helped me a lot.
        Now I'd like to do wavelet_Stat while using trmm pcp data as
the obs_data. I interpolate trmm data to grid that is matched with my
wrfout data(post-possec by UPP with type of grb),
as your konw, the wrfout data post-pocessed by UPP cannot be wrote, so
I transformed it to nc use ncl so that I can wrfite the intepolated
trmm data to it. Then I did the wavelet_stat step, but it comes out
with the error ""NetCDF: Attribute not found"". My config file is :
fcst = {
   field = [
      {
        name       = "APCP";
        level      = [  "A24" ];
        cat_thresh = [ >100 ];
      }
   ];
};
obs = {
   field = [
      {
        name       = "APCP";
        level      = [ "A24" ];
        cat_thresh = [ >99999 ];
      }
   ];
};


I have no idea why it comes out with this error when I only changed my
file type grb to nc(for write intepolated trmm data), the original
file named wrfprs_d02.24 and the chaged file named wrfprs_d02_24.nc


I'm looking forward to your early reply ! Thanks very much!









At 2015-03-19 01:17:50, "John Halley Gotway via RT"
<met_help at ucar.edu> wrote:
>Zhuang,
>
>I see that you'd like to use Wavelet-Stat to decompose a single
field.
>Here's one thing you could do... When you run Wavelet-Stat pass in
the
>forecast file as both the forecast and observation file.  In the
>configuration file, set the categorical threshold (cat_thresh) for
the
>forecast field in the way you'd like.  Then set cat_thresh for the
>observation field very, very high.  That'll result in the forecast
field
>being what you want and the observation field being all 0's.  So when
you
>subtract the forecast minus observation, you'll end up with the
forecast
>field.
>
>For example...
>
>fcst = {
>   field = [
>      {
>        name       = "APCP";
>        level      = [ "A12" ];
>        cat_thresh = [ >0.0 ];
>      }
>   ];
>};
>obs = {
>   field = [
>      {
>        name       = "APCP";
>        level      = [ "A12" ];
>        cat_thresh = [ >99999 ];
>      }
>   ];
>};
>
>This will define events in the forecast field but no event in the
>observation field.  So fcst - obs = fcst.
>
>Make sense?
>
>Thanks,
>John
>
>On Tue, Mar 17, 2015 at 9:29 PM, Zhuang from Nuist via RT
<met_help at ucar.edu
>> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71092 >
>>
>> Hello sir:
>>          Thanks very much for your help and I have successfully run
>> wavelrt_stat by its default config.
>>          Now I have some questions:
>> 1. is there any way for me to only decompose the fcst field without
a
>> matched obs_field like USER_GUIDE figure9.2.
>> As we know, the picture out from default config perform a scale
>> decomposition of the difference between fcst_field and obs_field.
>>
>>
>> 2.how can I perform a forecast energy squared picture like figure
9.6  in
>> user guide.
>>
>>
>> Thanks for your kind consideration!
>>
>>
>>
>>
>>
>>
>> At 2015-03-17 00:42:42, "John Halley Gotway via RT"
<met_help at ucar.edu>
>> wrote:
>> >Do you have the wgrib utility available on your system?  It can be
used to
>> >dump the contents of a GRIB file.  If so, please run that file
through
>> >wgrib and send me the output:
>> >   wgrib wrfprs_d02.00 > wrfprs_d02.00.dump
>> >
>> >Then send me that file "wrfprs_d02.00.dump" along with a
description of
>> >what field you'd like to run through Wavelet-Stat.  That'll help
me advise
>> >how you should set up the configuration file.
>> >
>> >Thanks,
>> >John Halley Gotway
>> >met_help at ucar.edu
>> >
>> >On Monday, March 16, 2015, Zhuang from Nuist via RT
<met_help at ucar.edu>
>> >wrote:
>> >
>> >>
>> >> Mon Mar 16 08:39:05 2015: Request 71115 was acted upon.
>> >> Transaction: Ticket created by zxrxz3212009 at 163.com
>> >>        Queue: met_help
>> >>      Subject: Sir,I got a problem using Wavelet_stat
>> >>        Owner: Nobody
>> >>   Requestors: zxrxz3212009 at 163.com
>> >>       Status: new
>> >>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71115 >
>> >>
>> >>
>> >> Hello,
>> >>       In my test I transform wrfout using UPP and the total
>> precipitation
>> >> variable in wrfprs_d02.00 is A_PCP_GDS3_SFC_acc.
>> >>       But when I use met to perform a wavelet_stat, and set
" name
>> >>      = "A_PCP_GDS3_SFC_acc";"   in WaveletStatConfig,
>> >> the error is ERROR  :
>> >> ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
>> >> 'R_H_GDS3_HTGL' for table version 2
>> >> ERROR  :
>> >>       It can be seen from file WaveletStatConfig_default that
>> >> "name="APCP"" and "level=[ "A03" ]" ,but the sample data
>> >> "data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212" doesn't
have a
>> var
>> >> named "APCP".
>> >>      I have no idea what does that mean and I do like to know
how to
>> fill
>> >> "name" and "level" in WaveletStatConfig.
>> >>
>> >>
>> >>     Thanks for your consideration!
>> >>
>> >>
>> >>
>> >
>>
>>
>

------------------------------------------------
Subject: Sir,I got a problem using Wavelet_stat
From: John Halley Gotway
Time: Fri Mar 20 13:37:25 2015

Hello,

I see that you're having trouble using the NetCDF output of NCL in
MET.
MET can easily read GRIB data and can read some flavors of NetCDF.
But
NetCDF is a very generic format, and MET is looking for certain
information
to process the data correctly.  That's why you're getting an attribute
not
found error.

Let's take a step back.  It sounds like you want to compare your WRF
output
to TRMM data.  If that's correct, you could do the following...

(1) Obtain 3-hourly or daily TRMM data in binary format.  Look for
those
links on this page:
   http://www.dtcenter.org/met/users/downloads/observation_data.php

(2) Run the "trmmbin2nc.R" Rscript mentioned on that page to convert
the
binary TRMM data to NetCDF.  Alos, you can edit the top-section of
that
script to define the range of lat/lon's you're studying.

(3) Run the GRIB output of UPP through the copygb utility to regrid
your
forecast data to the same grid as the TRMM data.  An examples or using
copygb to regrid to lat/lon can be found here:

http://www.dtcenter.org/met/users/support/online_tutorial/METv5.0/copygb/run2.php

Once you have the forecast and observation fields on the same grid,
you can
run them through any of the MET grid-to-grid tools.

Lastly, I want to mention that we're currently working on an automated
regridding feature for MET.  It'll be included in the next release,
version
5.1.  When that's available, you'll be able to pass 2 files to
wavelet_stat
that are not already on the same grid.  And you can tell it to regrid
the
forecast grid to the observations or vice-versa.

Hope that helps get you going.

Thanks,
John

On Fri, Mar 20, 2015 at 2:39 AM, Zhuang from Nuist via RT
<met_help at ucar.edu
> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71092 >
>
> Hello sir::
>         First thanks a lot for your kind consideration and it's
really
> helped me a lot.
>         Now I'd like to do wavelet_Stat while using trmm pcp data as
the
> obs_data. I interpolate trmm data to grid that is matched with my
wrfout
> data(post-possec by UPP with type of grb),
> as your konw, the wrfout data post-pocessed by UPP cannot be wrote,
so I
> transformed it to nc use ncl so that I can wrfite the intepolated
trmm data
> to it. Then I did the wavelet_stat step, but it comes out with the
error
> ""NetCDF: Attribute not found"". My config file is :
> fcst = {
>    field = [
>       {
>         name       = "APCP";
>         level      = [  "A24" ];
>         cat_thresh = [ >100 ];
>       }
>    ];
> };
> obs = {
>    field = [
>       {
>         name       = "APCP";
>         level      = [ "A24" ];
>         cat_thresh = [ >99999 ];
>       }
>    ];
> };
>
>
> I have no idea why it comes out with this error when I only changed
my
> file type grb to nc(for write intepolated trmm data), the original
file
> named wrfprs_d02.24 and the chaged file named wrfprs_d02_24.nc
>
>
> I'm looking forward to your early reply ! Thanks very much!
>
>
>
>
>
>
>
>
>
> At 2015-03-19 01:17:50, "John Halley Gotway via RT"
<met_help at ucar.edu>
> wrote:
> >Zhuang,
> >
> >I see that you'd like to use Wavelet-Stat to decompose a single
field.
> >Here's one thing you could do... When you run Wavelet-Stat pass in
the
> >forecast file as both the forecast and observation file.  In the
> >configuration file, set the categorical threshold (cat_thresh) for
the
> >forecast field in the way you'd like.  Then set cat_thresh for the
> >observation field very, very high.  That'll result in the forecast
field
> >being what you want and the observation field being all 0's.  So
when you
> >subtract the forecast minus observation, you'll end up with the
forecast
> >field.
> >
> >For example...
> >
> >fcst = {
> >   field = [
> >      {
> >        name       = "APCP";
> >        level      = [ "A12" ];
> >        cat_thresh = [ >0.0 ];
> >      }
> >   ];
> >};
> >obs = {
> >   field = [
> >      {
> >        name       = "APCP";
> >        level      = [ "A12" ];
> >        cat_thresh = [ >99999 ];
> >      }
> >   ];
> >};
> >
> >This will define events in the forecast field but no event in the
> >observation field.  So fcst - obs = fcst.
> >
> >Make sense?
> >
> >Thanks,
> >John
> >
> >On Tue, Mar 17, 2015 at 9:29 PM, Zhuang from Nuist via RT <
> met_help at ucar.edu
> >> wrote:
> >
> >>
> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71092 >
> >>
> >> Hello sir:
> >>          Thanks very much for your help and I have successfully
run
> >> wavelrt_stat by its default config.
> >>          Now I have some questions:
> >> 1. is there any way for me to only decompose the fcst field
without a
> >> matched obs_field like USER_GUIDE figure9.2.
> >> As we know, the picture out from default config perform a scale
> >> decomposition of the difference between fcst_field and obs_field.
> >>
> >>
> >> 2.how can I perform a forecast energy squared picture like figure
9.6
> in
> >> user guide.
> >>
> >>
> >> Thanks for your kind consideration!
> >>
> >>
> >>
> >>
> >>
> >>
> >> At 2015-03-17 00:42:42, "John Halley Gotway via RT"
<met_help at ucar.edu>
> >> wrote:
> >> >Do you have the wgrib utility available on your system?  It can
be
> used to
> >> >dump the contents of a GRIB file.  If so, please run that file
through
> >> >wgrib and send me the output:
> >> >   wgrib wrfprs_d02.00 > wrfprs_d02.00.dump
> >> >
> >> >Then send me that file "wrfprs_d02.00.dump" along with a
description of
> >> >what field you'd like to run through Wavelet-Stat.  That'll help
me
> advise
> >> >how you should set up the configuration file.
> >> >
> >> >Thanks,
> >> >John Halley Gotway
> >> >met_help at ucar.edu
> >> >
> >> >On Monday, March 16, 2015, Zhuang from Nuist via RT
<met_help at ucar.edu
> >
> >> >wrote:
> >> >
> >> >>
> >> >> Mon Mar 16 08:39:05 2015: Request 71115 was acted upon.
> >> >> Transaction: Ticket created by zxrxz3212009 at 163.com
> >> >>        Queue: met_help
> >> >>      Subject: Sir,I got a problem using Wavelet_stat
> >> >>        Owner: Nobody
> >> >>   Requestors: zxrxz3212009 at 163.com
> >> >>       Status: new
> >> >>  Ticket <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71115 >
> >> >>
> >> >>
> >> >> Hello,
> >> >>       In my test I transform wrfout using UPP and the total
> >> precipitation
> >> >> variable in wrfprs_d02.00 is A_PCP_GDS3_SFC_acc.
> >> >>       But when I use met to perform a wavelet_stat, and set
"
> name
> >> >>      = "A_PCP_GDS3_SFC_acc";"   in WaveletStatConfig,
> >> >> the error is ERROR  :
> >> >> ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field
> abbreviation
> >> >> 'R_H_GDS3_HTGL' for table version 2
> >> >> ERROR  :
> >> >>       It can be seen from file WaveletStatConfig_default that
> >> >> "name="APCP"" and "level=[ "A03" ]" ,but the sample data
> >> >> "data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212"
doesn't have
> a
> >> var
> >> >> named "APCP".
> >> >>      I have no idea what does that mean and I do like to know
how to
> >> fill
> >> >> "name" and "level" in WaveletStatConfig.
> >> >>
> >> >>
> >> >>     Thanks for your consideration!
> >> >>
> >> >>
> >> >>
> >> >
> >>
> >>
> >
>
>

------------------------------------------------
Subject: Sir,I got a problem using Wavelet_stat
From: Zhuang from Nuist
Time: Sun Mar 22 02:30:30 2015

       Thanks very much for your kindly help and your method really
helped me a lot.
       There is one more question for wavelet_stat. As 'Multiscale
characteristics and evalution of perturbations for warm season
convection-allowing precipitation forecasts: dependence on background
flow and method of perturbation', the author used MET Wavelet_stat to
perform a scale decompositon of total precipitation field insdead of
divide the pcp field into binary field(ie. 1 or 0) primarily due to a
giving threhold.
       My question is that if I'd like to do a scale decompostion only
to the original field(within or without a giving threhold) instead of
its binary field. Then what should I do?
      Thanks for your considerition.!








At 2015-03-21 03:37:26, "John Halley Gotway via RT"
<met_help at ucar.edu> wrote:
>Hello,
>
>I see that you're having trouble using the NetCDF output of NCL in
MET.
>MET can easily read GRIB data and can read some flavors of NetCDF.
But
>NetCDF is a very generic format, and MET is looking for certain
information
>to process the data correctly.  That's why you're getting an
attribute not
>found error.
>
>Let's take a step back.  It sounds like you want to compare your WRF
output
>to TRMM data.  If that's correct, you could do the following...
>
>(1) Obtain 3-hourly or daily TRMM data in binary format.  Look for
those
>links on this page:
>   http://www.dtcenter.org/met/users/downloads/observation_data.php
>
>(2) Run the "trmmbin2nc.R" Rscript mentioned on that page to convert
the
>binary TRMM data to NetCDF.  Alos, you can edit the top-section of
that
>script to define the range of lat/lon's you're studying.
>
>(3) Run the GRIB output of UPP through the copygb utility to regrid
your
>forecast data to the same grid as the TRMM data.  An examples or
using
>copygb to regrid to lat/lon can be found here:
>
>http://www.dtcenter.org/met/users/support/online_tutorial/METv5.0/copygb/run2.php
>
>Once you have the forecast and observation fields on the same grid,
you can
>run them through any of the MET grid-to-grid tools.
>
>Lastly, I want to mention that we're currently working on an
automated
>regridding feature for MET.  It'll be included in the next release,
version
>5.1.  When that's available, you'll be able to pass 2 files to
wavelet_stat
>that are not already on the same grid.  And you can tell it to regrid
the
>forecast grid to the observations or vice-versa.
>
>Hope that helps get you going.
>
>Thanks,
>John
>
>On Fri, Mar 20, 2015 at 2:39 AM, Zhuang from Nuist via RT
<met_help at ucar.edu
>> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71092 >
>>
>> Hello sir::
>>         First thanks a lot for your kind consideration and it's
really
>> helped me a lot.
>>         Now I'd like to do wavelet_Stat while using trmm pcp data
as the
>> obs_data. I interpolate trmm data to grid that is matched with my
wrfout
>> data(post-possec by UPP with type of grb),
>> as your konw, the wrfout data post-pocessed by UPP cannot be wrote,
so I
>> transformed it to nc use ncl so that I can wrfite the intepolated
trmm data
>> to it. Then I did the wavelet_stat step, but it comes out with the
error
>> ""NetCDF: Attribute not found"". My config file is :
>> fcst = {
>>    field = [
>>       {
>>         name       = "APCP";
>>         level      = [  "A24" ];
>>         cat_thresh = [ >100 ];
>>       }
>>    ];
>> };
>> obs = {
>>    field = [
>>       {
>>         name       = "APCP";
>>         level      = [ "A24" ];
>>         cat_thresh = [ >99999 ];
>>       }
>>    ];
>> };
>>
>>
>> I have no idea why it comes out with this error when I only changed
my
>> file type grb to nc(for write intepolated trmm data), the original
file
>> named wrfprs_d02.24 and the chaged file named wrfprs_d02_24.nc
>>
>>
>> I'm looking forward to your early reply ! Thanks very much!
>>
>>
>>
>>
>>
>>
>>
>>
>>
>> At 2015-03-19 01:17:50, "John Halley Gotway via RT"
<met_help at ucar.edu>
>> wrote:
>> >Zhuang,
>> >
>> >I see that you'd like to use Wavelet-Stat to decompose a single
field.
>> >Here's one thing you could do... When you run Wavelet-Stat pass in
the
>> >forecast file as both the forecast and observation file.  In the
>> >configuration file, set the categorical threshold (cat_thresh) for
the
>> >forecast field in the way you'd like.  Then set cat_thresh for the
>> >observation field very, very high.  That'll result in the forecast
field
>> >being what you want and the observation field being all 0's.  So
when you
>> >subtract the forecast minus observation, you'll end up with the
forecast
>> >field.
>> >
>> >For example...
>> >
>> >fcst = {
>> >   field = [
>> >      {
>> >        name       = "APCP";
>> >        level      = [ "A12" ];
>> >        cat_thresh = [ >0.0 ];
>> >      }
>> >   ];
>> >};
>> >obs = {
>> >   field = [
>> >      {
>> >        name       = "APCP";
>> >        level      = [ "A12" ];
>> >        cat_thresh = [ >99999 ];
>> >      }
>> >   ];
>> >};
>> >
>> >This will define events in the forecast field but no event in the
>> >observation field.  So fcst - obs = fcst.
>> >
>> >Make sense?
>> >
>> >Thanks,
>> >John
>> >
>> >On Tue, Mar 17, 2015 at 9:29 PM, Zhuang from Nuist via RT <
>> met_help at ucar.edu
>> >> wrote:
>> >
>> >>
>> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71092 >
>> >>
>> >> Hello sir:
>> >>          Thanks very much for your help and I have successfully
run
>> >> wavelrt_stat by its default config.
>> >>          Now I have some questions:
>> >> 1. is there any way for me to only decompose the fcst field
without a
>> >> matched obs_field like USER_GUIDE figure9.2.
>> >> As we know, the picture out from default config perform a scale
>> >> decomposition of the difference between fcst_field and
obs_field.
>> >>
>> >>
>> >> 2.how can I perform a forecast energy squared picture like
figure 9.6
>> in
>> >> user guide.
>> >>
>> >>
>> >> Thanks for your kind consideration!
>> >>
>> >>
>> >>
>> >>
>> >>
>> >>
>> >> At 2015-03-17 00:42:42, "John Halley Gotway via RT"
<met_help at ucar.edu>
>> >> wrote:
>> >> >Do you have the wgrib utility available on your system?  It can
be
>> used to
>> >> >dump the contents of a GRIB file.  If so, please run that file
through
>> >> >wgrib and send me the output:
>> >> >   wgrib wrfprs_d02.00 > wrfprs_d02.00.dump
>> >> >
>> >> >Then send me that file "wrfprs_d02.00.dump" along with a
description of
>> >> >what field you'd like to run through Wavelet-Stat.  That'll
help me
>> advise
>> >> >how you should set up the configuration file.
>> >> >
>> >> >Thanks,
>> >> >John Halley Gotway
>> >> >met_help at ucar.edu
>> >> >
>> >> >On Monday, March 16, 2015, Zhuang from Nuist via RT
<met_help at ucar.edu
>> >
>> >> >wrote:
>> >> >
>> >> >>
>> >> >> Mon Mar 16 08:39:05 2015: Request 71115 was acted upon.
>> >> >> Transaction: Ticket created by zxrxz3212009 at 163.com
>> >> >>        Queue: met_help
>> >> >>      Subject: Sir,I got a problem using Wavelet_stat
>> >> >>        Owner: Nobody
>> >> >>   Requestors: zxrxz3212009 at 163.com
>> >> >>       Status: new
>> >> >>  Ticket <URL:
>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71115 >
>> >> >>
>> >> >>
>> >> >> Hello,
>> >> >>       In my test I transform wrfout using UPP and the total
>> >> precipitation
>> >> >> variable in wrfprs_d02.00 is A_PCP_GDS3_SFC_acc.
>> >> >>       But when I use met to perform a wavelet_stat, and set
"
>> name
>> >> >>      = "A_PCP_GDS3_SFC_acc";"   in WaveletStatConfig,
>> >> >> the error is ERROR  :
>> >> >> ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field
>> abbreviation
>> >> >> 'R_H_GDS3_HTGL' for table version 2
>> >> >> ERROR  :
>> >> >>       It can be seen from file WaveletStatConfig_default that
>> >> >> "name="APCP"" and "level=[ "A03" ]" ,but the sample data
>> >> >> "data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212"
doesn't have
>> a
>> >> var
>> >> >> named "APCP".
>> >> >>      I have no idea what does that mean and I do like to know
how to
>> >> fill
>> >> >> "name" and "level" in WaveletStatConfig.
>> >> >>
>> >> >>
>> >> >>     Thanks for your consideration!
>> >> >>
>> >> >>
>> >> >>
>> >> >
>> >>
>> >>
>> >
>>
>>
>

------------------------------------------------
Subject: Sir,I got a problem using Wavelet_stat
From: John Halley Gotway
Time: Mon Mar 23 10:36:17 2015

Zhuang,

The wavelet-stat tool is set up to be run on binary fields by applying
a
threshold to them.  If you want to change that behavior, you'd need to
edit
the code and recompile.  But it'll be up to you to interpret the
output.
Here's how you could do that...

(1) Edit lines 976 and 977 of file wavelet_stat.cc to skip the
thresholding
step:
OLD...
    976          f_dat[j] =
isc_info[i].cts_fcst_thresh.check(f_na[j]);
    977          o_dat[j] = isc_info[i].cts_obs_thresh.check(o_na[j]);
NEW...
    976          f_dat[j] = f_na[j];
    977          o_dat[j] = o_na[j];

(2) Rebuild MET: make install

(3) Edit the wavelet_stat configuration file by changing the plotting
range:
wvlt_plot = {
   color_table = "MET_BASE/colortables/NCL_colortables/BlWhRe.ctable";
   plot_min = -25;
   plot_max = 25;
};

(4) Choose a different wavelet type for the decomposition.  The Haar
wavelet is a good choice for binary fields, but not for more
continuous
fields.  Search for "wavelet" in the file met-5.0/data/config/README
for a
listing of the available wavelet types.
wavelet = {
   type   = HAAR;
   member = 2;
};

Thanks,
John

On Sun, Mar 22, 2015 at 2:30 AM, Zhuang from Nuist via RT
<met_help at ucar.edu
> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71092 >
>
>        Thanks very much for your kindly help and your method really
helped
> me a lot.
>        There is one more question for wavelet_stat. As 'Multiscale
> characteristics and evalution of perturbations for warm season
> convection-allowing precipitation forecasts: dependence on
background flow
> and method of perturbation', the author used MET Wavelet_stat to
perform a
> scale decompositon of total precipitation field insdead of divide
the pcp
> field into binary field(ie. 1 or 0) primarily due to a giving
threhold.
>        My question is that if I'd like to do a scale decompostion
only to
> the original field(within or without a giving threhold) instead of
its
> binary field. Then what should I do?
>       Thanks for your considerition.!
>
>
>
>
>
>
>
>
> At 2015-03-21 03:37:26, "John Halley Gotway via RT"
<met_help at ucar.edu>
> wrote:
> >Hello,
> >
> >I see that you're having trouble using the NetCDF output of NCL in
MET.
> >MET can easily read GRIB data and can read some flavors of NetCDF.
But
> >NetCDF is a very generic format, and MET is looking for certain
> information
> >to process the data correctly.  That's why you're getting an
attribute not
> >found error.
> >
> >Let's take a step back.  It sounds like you want to compare your
WRF
> output
> >to TRMM data.  If that's correct, you could do the following...
> >
> >(1) Obtain 3-hourly or daily TRMM data in binary format.  Look for
those
> >links on this page:
> >   http://www.dtcenter.org/met/users/downloads/observation_data.php
> >
> >(2) Run the "trmmbin2nc.R" Rscript mentioned on that page to
convert the
> >binary TRMM data to NetCDF.  Alos, you can edit the top-section of
that
> >script to define the range of lat/lon's you're studying.
> >
> >(3) Run the GRIB output of UPP through the copygb utility to regrid
your
> >forecast data to the same grid as the TRMM data.  An examples or
using
> >copygb to regrid to lat/lon can be found here:
> >
> >
>
http://www.dtcenter.org/met/users/support/online_tutorial/METv5.0/copygb/run2.php
> >
> >Once you have the forecast and observation fields on the same grid,
you
> can
> >run them through any of the MET grid-to-grid tools.
> >
> >Lastly, I want to mention that we're currently working on an
automated
> >regridding feature for MET.  It'll be included in the next release,
> version
> >5.1.  When that's available, you'll be able to pass 2 files to
> wavelet_stat
> >that are not already on the same grid.  And you can tell it to
regrid the
> >forecast grid to the observations or vice-versa.
> >
> >Hope that helps get you going.
> >
> >Thanks,
> >John
> >
> >On Fri, Mar 20, 2015 at 2:39 AM, Zhuang from Nuist via RT <
> met_help at ucar.edu
> >> wrote:
> >
> >>
> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71092 >
> >>
> >> Hello sir::
> >>         First thanks a lot for your kind consideration and it's
really
> >> helped me a lot.
> >>         Now I'd like to do wavelet_Stat while using trmm pcp data
as the
> >> obs_data. I interpolate trmm data to grid that is matched with my
wrfout
> >> data(post-possec by UPP with type of grb),
> >> as your konw, the wrfout data post-pocessed by UPP cannot be
wrote, so I
> >> transformed it to nc use ncl so that I can wrfite the intepolated
trmm
> data
> >> to it. Then I did the wavelet_stat step, but it comes out with
the error
> >> ""NetCDF: Attribute not found"". My config file is :
> >> fcst = {
> >>    field = [
> >>       {
> >>         name       = "APCP";
> >>         level      = [  "A24" ];
> >>         cat_thresh = [ >100 ];
> >>       }
> >>    ];
> >> };
> >> obs = {
> >>    field = [
> >>       {
> >>         name       = "APCP";
> >>         level      = [ "A24" ];
> >>         cat_thresh = [ >99999 ];
> >>       }
> >>    ];
> >> };
> >>
> >>
> >> I have no idea why it comes out with this error when I only
changed my
> >> file type grb to nc(for write intepolated trmm data), the
original file
> >> named wrfprs_d02.24 and the chaged file named wrfprs_d02_24.nc
> >>
> >>
> >> I'm looking forward to your early reply ! Thanks very much!
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >> At 2015-03-19 01:17:50, "John Halley Gotway via RT"
<met_help at ucar.edu>
> >> wrote:
> >> >Zhuang,
> >> >
> >> >I see that you'd like to use Wavelet-Stat to decompose a single
field.
> >> >Here's one thing you could do... When you run Wavelet-Stat pass
in the
> >> >forecast file as both the forecast and observation file.  In the
> >> >configuration file, set the categorical threshold (cat_thresh)
for the
> >> >forecast field in the way you'd like.  Then set cat_thresh for
the
> >> >observation field very, very high.  That'll result in the
forecast
> field
> >> >being what you want and the observation field being all 0's.  So
when
> you
> >> >subtract the forecast minus observation, you'll end up with the
> forecast
> >> >field.
> >> >
> >> >For example...
> >> >
> >> >fcst = {
> >> >   field = [
> >> >      {
> >> >        name       = "APCP";
> >> >        level      = [ "A12" ];
> >> >        cat_thresh = [ >0.0 ];
> >> >      }
> >> >   ];
> >> >};
> >> >obs = {
> >> >   field = [
> >> >      {
> >> >        name       = "APCP";
> >> >        level      = [ "A12" ];
> >> >        cat_thresh = [ >99999 ];
> >> >      }
> >> >   ];
> >> >};
> >> >
> >> >This will define events in the forecast field but no event in
the
> >> >observation field.  So fcst - obs = fcst.
> >> >
> >> >Make sense?
> >> >
> >> >Thanks,
> >> >John
> >> >
> >> >On Tue, Mar 17, 2015 at 9:29 PM, Zhuang from Nuist via RT <
> >> met_help at ucar.edu
> >> >> wrote:
> >> >
> >> >>
> >> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71092
>
> >> >>
> >> >> Hello sir:
> >> >>          Thanks very much for your help and I have
successfully run
> >> >> wavelrt_stat by its default config.
> >> >>          Now I have some questions:
> >> >> 1. is there any way for me to only decompose the fcst field
without a
> >> >> matched obs_field like USER_GUIDE figure9.2.
> >> >> As we know, the picture out from default config perform a
scale
> >> >> decomposition of the difference between fcst_field and
obs_field.
> >> >>
> >> >>
> >> >> 2.how can I perform a forecast energy squared picture like
figure 9.6
> >> in
> >> >> user guide.
> >> >>
> >> >>
> >> >> Thanks for your kind consideration!
> >> >>
> >> >>
> >> >>
> >> >>
> >> >>
> >> >>
> >> >> At 2015-03-17 00:42:42, "John Halley Gotway via RT" <
> met_help at ucar.edu>
> >> >> wrote:
> >> >> >Do you have the wgrib utility available on your system?  It
can be
> >> used to
> >> >> >dump the contents of a GRIB file.  If so, please run that
file
> through
> >> >> >wgrib and send me the output:
> >> >> >   wgrib wrfprs_d02.00 > wrfprs_d02.00.dump
> >> >> >
> >> >> >Then send me that file "wrfprs_d02.00.dump" along with a
> description of
> >> >> >what field you'd like to run through Wavelet-Stat.  That'll
help me
> >> advise
> >> >> >how you should set up the configuration file.
> >> >> >
> >> >> >Thanks,
> >> >> >John Halley Gotway
> >> >> >met_help at ucar.edu
> >> >> >
> >> >> >On Monday, March 16, 2015, Zhuang from Nuist via RT <
> met_help at ucar.edu
> >> >
> >> >> >wrote:
> >> >> >
> >> >> >>
> >> >> >> Mon Mar 16 08:39:05 2015: Request 71115 was acted upon.
> >> >> >> Transaction: Ticket created by zxrxz3212009 at 163.com
> >> >> >>        Queue: met_help
> >> >> >>      Subject: Sir,I got a problem using Wavelet_stat
> >> >> >>        Owner: Nobody
> >> >> >>   Requestors: zxrxz3212009 at 163.com
> >> >> >>       Status: new
> >> >> >>  Ticket <URL:
> >> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71115 >
> >> >> >>
> >> >> >>
> >> >> >> Hello,
> >> >> >>       In my test I transform wrfout using UPP and the total
> >> >> precipitation
> >> >> >> variable in wrfprs_d02.00 is A_PCP_GDS3_SFC_acc.
> >> >> >>       But when I use met to perform a wavelet_stat, and set
>  "
> >> name
> >> >> >>      = "A_PCP_GDS3_SFC_acc";"   in WaveletStatConfig,
> >> >> >> the error is ERROR  :
> >> >> >> ERROR  : VarInfoGrib::set_dict() - unrecognized GRIB1 field
> >> abbreviation
> >> >> >> 'R_H_GDS3_HTGL' for table version 2
> >> >> >> ERROR  :
> >> >> >>       It can be seen from file WaveletStatConfig_default
that
> >> >> >> "name="APCP"" and "level=[ "A03" ]" ,but the sample data
> >> >> >> "data/sample_fcst/2005080700/wrfprs_ruc13_12.tm00_G212"
doesn't
> have
> >> a
> >> >> var
> >> >> >> named "APCP".
> >> >> >>      I have no idea what does that mean and I do like to
know how
> to
> >> >> fill
> >> >> >> "name" and "level" in WaveletStatConfig.
> >> >> >>
> >> >> >>
> >> >> >>     Thanks for your consideration!
> >> >> >>
> >> >> >>
> >> >> >>
> >> >> >
> >> >>
> >> >>
> >> >
> >>
> >>
> >
>
>

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


More information about the Met_help mailing list