[Met_help] [rt.rap.ucar.edu #55241] History for Using point_stat tool to verify the rainfall data

Paul Oldenburg via RT met_help at ucar.edu
Tue Apr 3 08:44:39 MDT 2012


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

Hi

I want to verified the rainfall data by using the point_stat.sh.
I have several questions, hope you can help me.

My observation data includes hourly rainfall data, 12hr accumulated
rainfall data and the daily accumulate rainfall data.
My forecast data is the WRF-ARW output data


Below is my processing steps:

   1. Prepare the observation rainfall data in 10 columns ascii format.

   For 24hr accumulate rainfall data.  The data is excerpted as following.
   ---------------
   ADPSFC    C1I230   20110829_000000   23.9639  120.8372 917   61 12
   -9999.00     35.5
   ADPSFC    466940   20110829_000000   25.1348  121.7321 27 61 12
   -9999.00     80
   ADPSFC    466900   20110829_000000   25.1656  121.44   19 61 12
   -9999.00    14.5
   ......
   ---------------
   This date set is for 24-hr accumulated rainfall from 20110828 000000UTC
   to 20110829 000000UTC.


   For 12hr accumulate rainfall data.  The data is excerpted as following.
   ---------------
   ADPSFC   C1I230   20110828_120000   23.9639  120.8372 917   61 12
   -9999.00    21.5
   ADPSFC   466940   20110828_120000   25.1348  121.7321 27 61 12 -9999.00
      29
   ADPSFC   466900   20110828_120000   25.1656  121.44   19 61 12 -9999.00
      7
   ......
   ADPSFC   C1I230   20110829_000000   23.9639  120.8372 917   61 12
   -9999.00    14
   ADPSFC   466940   20110829_000000   25.1348  121.7321 27 61 12 -9999.00
      51
   ADPSFC   466900   20110829_000000   25.1656  121.44   19 61 12 -9999.00
      7.5
   ---------------


   2. Using the ascii2nc to produced the observation data in netcdf format.
   For 24-hr accumulated rainfall data, the netcdf file is *
   precip24_2011082900.nc*
   For 12-hr accumulated rainfall data, the netcdf file is *
   precip24_20110828_0829_12hr.nc*

   3. Transfer the WRF-ARW output by using UPPV1.0
   My output data is saved as one file per hour.
   So, it produced several file from UPP program.

   for example:* wrfprs_d03.000*, *wrfprs_d03.001*,   *wrfprs_d03.002* .......
   *wrfprs_d03.024*

   4. If I want to verify the first 24hr accumulated rainfall, by the 24-hr
   observational accumulated rainfall

   The config file was excerpted as following
   ----------------
   fcst_field[] = [ "APCP/A24" , "APCP/A24" , "APCP/A24" , "APCP/A24" ];
   obs_field[]  = [ "61/A24" , "61/A24" , "61/A24" , "61/A24" ];
   ----------------

   ../bin/point_stat \
   ../../POST_UPP/M15_2011082800/*wrfprs_d03.024 *\
   ../out/ascii2nc/*precip24_2011082800.nc *\
   config/PointStat_Rain \
   -outdir ../out/point_stat -v 2

   5. If I want to verify the first 24hr accumulated rainfall, by the 12-hr
   observational accumulated rainfall

   The config file was excerpted as following
   ----------------
   fcst_field[] = [ "APCP/A24" , "APCP/A24" , "APCP/A24" , "APCP/A24" ];
   obs_field[]  = [ "61/A12" , "61/A12" , "61/A12" , "61/A12" ];
   ----------------

   ../bin/point_stat \
   ../../POST_UPP/M15_2011082800/*wrfprs_d03.024* \
   ../out/ascii2nc/*precip24_20110828_0829_12hr.nc* \
   config/PointStat_Rain \
   -outdir ../out/point_stat -v 2


My question is:

   - Does there have some procedure I am doing wrong?
   - The step 4 and step 5 should be identical, but the results shows that
   the step 5 only compare with 12hr accumulate.
   How can I using the 12-hr observational data to verify the 24-hr
   accumulated rainfall?
   - If I want to verify the continues hourly rainfall data with the
   observational hourly rainfall data.
   How do I put the fcst_file in the execute line?  (Since I have several
   output file. )


Thank you very much.

-- 
Delia, Yen-Chu Chen
Department of Atmospheric Sciences
National Taiwan University
Email : deliachen at ntu.edu.tw


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

Subject: Re: [rt.rap.ucar.edu #55241] Using point_stat tool to verify the rainfall data
From: Paul Oldenburg
Time: Mon Mar 19 08:49:30 2012

Delia,

I see that your email address is not currently a registered user of
the MET software.  Prior to providing direct support
in the use of MET, we require that users fill out the free online
registration form.  That enables us to ensure that we
remain in  compliance with the international export laws of the United
States.

To fill out the registration, just go here:
    http://www.dtcenter.org/met/users/downloads/index.php
Enter you email address and hit enter to be taken to the form.

Thanks,

Paul


On 03/19/2012 02:33 AM, DeliaChen via RT wrote:
>
> Mon Mar 19 02:33:28 2012: Request 55241 was acted upon.
> Transaction: Ticket created by deliachen at ntu.edu.tw
>         Queue: met_help
>       Subject: Using point_stat tool to verify the rainfall data
>         Owner: Nobody
>    Requestors: deliachen at ntu.edu.tw
>        Status: new
>   Ticket<URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=55241>
>
>
> Hi
>
> I want to verified the rainfall data by using the point_stat.sh.
> I have several questions, hope you can help me.
>
> My observation data includes hourly rainfall data, 12hr accumulated
> rainfall data and the daily accumulate rainfall data.
> My forecast data is the WRF-ARW output data
>
>
> Below is my processing steps:
>
>     1. Prepare the observation rainfall data in 10 columns ascii
format.
>
>     For 24hr accumulate rainfall data.  The data is excerpted as
following.
>     ---------------
>     ADPSFC    C1I230   20110829_000000   23.9639  120.8372 917   61
12
>     -9999.00     35.5
>     ADPSFC    466940   20110829_000000   25.1348  121.7321 27 61 12
>     -9999.00     80
>     ADPSFC    466900   20110829_000000   25.1656  121.44   19 61 12
>     -9999.00    14.5
>     ......
>     ---------------
>     This date set is for 24-hr accumulated rainfall from 20110828
000000UTC
>     to 20110829 000000UTC.
>
>
>     For 12hr accumulate rainfall data.  The data is excerpted as
following.
>     ---------------
>     ADPSFC   C1I230   20110828_120000   23.9639  120.8372 917   61
12
>     -9999.00    21.5
>     ADPSFC   466940   20110828_120000   25.1348  121.7321 27 61 12
-9999.00
>        29
>     ADPSFC   466900   20110828_120000   25.1656  121.44   19 61 12
-9999.00
>        7
>     ......
>     ADPSFC   C1I230   20110829_000000   23.9639  120.8372 917   61
12
>     -9999.00    14
>     ADPSFC   466940   20110829_000000   25.1348  121.7321 27 61 12
-9999.00
>        51
>     ADPSFC   466900   20110829_000000   25.1656  121.44   19 61 12
-9999.00
>        7.5
>     ---------------
>
>
>     2. Using the ascii2nc to produced the observation data in netcdf
format.
>     For 24-hr accumulated rainfall data, the netcdf file is *
>     precip24_2011082900.nc*
>     For 12-hr accumulated rainfall data, the netcdf file is *
>     precip24_20110828_0829_12hr.nc*
>
>     3. Transfer the WRF-ARW output by using UPPV1.0
>     My output data is saved as one file per hour.
>     So, it produced several file from UPP program.
>
>     for example:* wrfprs_d03.000*, *wrfprs_d03.001*,
*wrfprs_d03.002* .......
>     *wrfprs_d03.024*
>
>     4. If I want to verify the first 24hr accumulated rainfall, by
the 24-hr
>     observational accumulated rainfall
>
>     The config file was excerpted as following
>     ----------------
>     fcst_field[] = [ "APCP/A24" , "APCP/A24" , "APCP/A24" ,
"APCP/A24" ];
>     obs_field[]  = [ "61/A24" , "61/A24" , "61/A24" , "61/A24" ];
>     ----------------
>
>     ../bin/point_stat \
>     ../../POST_UPP/M15_2011082800/*wrfprs_d03.024 *\
>     ../out/ascii2nc/*precip24_2011082800.nc *\
>     config/PointStat_Rain \
>     -outdir ../out/point_stat -v 2
>
>     5. If I want to verify the first 24hr accumulated rainfall, by
the 12-hr
>     observational accumulated rainfall
>
>     The config file was excerpted as following
>     ----------------
>     fcst_field[] = [ "APCP/A24" , "APCP/A24" , "APCP/A24" ,
"APCP/A24" ];
>     obs_field[]  = [ "61/A12" , "61/A12" , "61/A12" , "61/A12" ];
>     ----------------
>
>     ../bin/point_stat \
>     ../../POST_UPP/M15_2011082800/*wrfprs_d03.024* \
>     ../out/ascii2nc/*precip24_20110828_0829_12hr.nc* \
>     config/PointStat_Rain \
>     -outdir ../out/point_stat -v 2
>
>
> My question is:
>
>     - Does there have some procedure I am doing wrong?
>     - The step 4 and step 5 should be identical, but the results
shows that
>     the step 5 only compare with 12hr accumulate.
>     How can I using the 12-hr observational data to verify the 24-hr
>     accumulated rainfall?
>     - If I want to verify the continues hourly rainfall data with
the
>     observational hourly rainfall data.
>     How do I put the fcst_file in the execute line?  (Since I have
several
>     output file. )
>
>
> Thank you very much.
>


------------------------------------------------
Subject: Using point_stat tool to verify the rainfall data
From: Delia Chen
Time: Mon Mar 19 18:17:53 2012

Paul,

Sorry for the misunderstanding,
I have several email accounts, I have already registered the other
one,
which is
 "delia at nat.as.ntu.edu.tw".

Thank you very much for your help.

Sincerely,
Delia

2012/3/19 Paul Oldenburg via RT <met_help at ucar.edu>

> Delia,
>
> I see that your email address is not currently a registered user of
the
> MET software.  Prior to providing direct support
> in the use of MET, we require that users fill out the free online
> registration form.  That enables us to ensure that we
> remain in  compliance with the international export laws of the
United
> States.
>
> To fill out the registration, just go here:
>    http://www.dtcenter.org/met/users/downloads/index.php
> Enter you email address and hit enter to be taken to the form.
>
> Thanks,
>
> Paul
>
>
> On 03/19/2012 02:33 AM, DeliaChen via RT wrote:
> >
> > Mon Mar 19 02:33:28 2012: Request 55241 was acted upon.
> > Transaction: Ticket created by deliachen at ntu.edu.tw
> >         Queue: met_help
> >       Subject: Using point_stat tool to verify the rainfall data
> >         Owner: Nobody
> >    Requestors: deliachen at ntu.edu.tw
> >        Status: new
> >   Ticket<URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=55241>
> >
> >
> > Hi
> >
> > I want to verified the rainfall data by using the point_stat.sh.
> > I have several questions, hope you can help me.
> >
> > My observation data includes hourly rainfall data, 12hr
accumulated
> > rainfall data and the daily accumulate rainfall data.
> > My forecast data is the WRF-ARW output data
> >
> >
> > Below is my processing steps:
> >
> >     1. Prepare the observation rainfall data in 10 columns ascii
format.
> >
> >     For 24hr accumulate rainfall data.  The data is excerpted as
> following.
> >     ---------------
> >     ADPSFC    C1I230   20110829_000000   23.9639  120.8372 917
61 12
> >     -9999.00     35.5
> >     ADPSFC    466940   20110829_000000   25.1348  121.7321 27 61
12
> >     -9999.00     80
> >     ADPSFC    466900   20110829_000000   25.1656  121.44   19 61
12
> >     -9999.00    14.5
> >     ......
> >     ---------------
> >     This date set is for 24-hr accumulated rainfall from 20110828
> 000000UTC
> >     to 20110829 000000UTC.
> >
> >
> >     For 12hr accumulate rainfall data.  The data is excerpted as
> following.
> >     ---------------
> >     ADPSFC   C1I230   20110828_120000   23.9639  120.8372 917   61
12
> >     -9999.00    21.5
> >     ADPSFC   466940   20110828_120000   25.1348  121.7321 27 61 12
> -9999.00
> >        29
> >     ADPSFC   466900   20110828_120000   25.1656  121.44   19 61 12
> -9999.00
> >        7
> >     ......
> >     ADPSFC   C1I230   20110829_000000   23.9639  120.8372 917   61
12
> >     -9999.00    14
> >     ADPSFC   466940   20110829_000000   25.1348  121.7321 27 61 12
> -9999.00
> >        51
> >     ADPSFC   466900   20110829_000000   25.1656  121.44   19 61 12
> -9999.00
> >        7.5
> >     ---------------
> >
> >
> >     2. Using the ascii2nc to produced the observation data in
netcdf
> format.
> >     For 24-hr accumulated rainfall data, the netcdf file is *
> >     precip24_2011082900.nc*
> >     For 12-hr accumulated rainfall data, the netcdf file is *
> >     precip24_20110828_0829_12hr.nc*
> >
> >     3. Transfer the WRF-ARW output by using UPPV1.0
> >     My output data is saved as one file per hour.
> >     So, it produced several file from UPP program.
> >
> >     for example:* wrfprs_d03.000*, *wrfprs_d03.001*,
*wrfprs_d03.002*
> .......
> >     *wrfprs_d03.024*
> >
> >     4. If I want to verify the first 24hr accumulated rainfall, by
the
> 24-hr
> >     observational accumulated rainfall
> >
> >     The config file was excerpted as following
> >     ----------------
> >     fcst_field[] = [ "APCP/A24" , "APCP/A24" , "APCP/A24" ,
"APCP/A24" ];
> >     obs_field[]  = [ "61/A24" , "61/A24" , "61/A24" , "61/A24" ];
> >     ----------------
> >
> >     ../bin/point_stat \
> >     ../../POST_UPP/M15_2011082800/*wrfprs_d03.024 *\
> >     ../out/ascii2nc/*precip24_2011082800.nc *\
> >     config/PointStat_Rain \
> >     -outdir ../out/point_stat -v 2
> >
> >     5. If I want to verify the first 24hr accumulated rainfall, by
the
> 12-hr
> >     observational accumulated rainfall
> >
> >     The config file was excerpted as following
> >     ----------------
> >     fcst_field[] = [ "APCP/A24" , "APCP/A24" , "APCP/A24" ,
"APCP/A24" ];
> >     obs_field[]  = [ "61/A12" , "61/A12" , "61/A12" , "61/A12" ];
> >     ----------------
> >
> >     ../bin/point_stat \
> >     ../../POST_UPP/M15_2011082800/*wrfprs_d03.024* \
> >     ../out/ascii2nc/*precip24_20110828_0829_12hr.nc* \
> >     config/PointStat_Rain \
> >     -outdir ../out/point_stat -v 2
> >
> >
> > My question is:
> >
> >     - Does there have some procedure I am doing wrong?
> >     - The step 4 and step 5 should be identical, but the results
shows
> that
> >     the step 5 only compare with 12hr accumulate.
> >     How can I using the 12-hr observational data to verify the 24-
hr
> >     accumulated rainfall?
> >     - If I want to verify the continues hourly rainfall data with
the
> >     observational hourly rainfall data.
> >     How do I put the fcst_file in the execute line?  (Since I have
> several
> >     output file. )
> >
> >
> > Thank you very much.
> >
>
>
>


--
Delia, Yen-Chu Chen
Email : deliachen at ntu.edu.tw

------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #55241] Using point_stat tool to verify the  rainfall data
From: Paul Oldenburg
Time: Mon Mar 19 20:38:15 2012

Delia,

I am sorry to ask you this, but I cannot provide MET help to an
unregistered email address for
legal reasons.  Can you either register your current email address or
re-send me your MET help
inquiry from your other email account?  Sorry for the inconvenience.
The reason for it is because
the US federal government considers MET help an "export".

Thanks,

Paul


>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=55241 >
>
> Paul,
>
> Sorry for the misunderstanding,
> I have several email accounts, I have already registered the other
one,
> which is
>  "delia at nat.as.ntu.edu.tw".
>
> Thank you very much for your help.
>
> Sincerely,
> Delia
>
> 2012/3/19 Paul Oldenburg via RT <met_help at ucar.edu>
>
>> Delia,
>>
>> I see that your email address is not currently a registered user of
the
>> MET software.  Prior to providing direct support
>> in the use of MET, we require that users fill out the free online
>> registration form.  That enables us to ensure that we
>> remain in  compliance with the international export laws of the
United
>> States.
>>
>> To fill out the registration, just go here:
>>    http://www.dtcenter.org/met/users/downloads/index.php
>> Enter you email address and hit enter to be taken to the form.
>>
>> Thanks,
>>
>> Paul
>>
>>
>> On 03/19/2012 02:33 AM, DeliaChen via RT wrote:
>> >
>> > Mon Mar 19 02:33:28 2012: Request 55241 was acted upon.
>> > Transaction: Ticket created by deliachen at ntu.edu.tw
>> >         Queue: met_help
>> >       Subject: Using point_stat tool to verify the rainfall data
>> >         Owner: Nobody
>> >    Requestors: deliachen at ntu.edu.tw
>> >        Status: new
>> >   Ticket<URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=55241>
>> >
>> >
>> > Hi
>> >
>> > I want to verified the rainfall data by using the point_stat.sh.
>> > I have several questions, hope you can help me.
>> >
>> > My observation data includes hourly rainfall data, 12hr
accumulated
>> > rainfall data and the daily accumulate rainfall data.
>> > My forecast data is the WRF-ARW output data
>> >
>> >
>> > Below is my processing steps:
>> >
>> >     1. Prepare the observation rainfall data in 10 columns ascii
format.
>> >
>> >     For 24hr accumulate rainfall data.  The data is excerpted as
>> following.
>> >     ---------------
>> >     ADPSFC    C1I230   20110829_000000   23.9639  120.8372 917
61 12
>> >     -9999.00     35.5
>> >     ADPSFC    466940   20110829_000000   25.1348  121.7321 27 61
12
>> >     -9999.00     80
>> >     ADPSFC    466900   20110829_000000   25.1656  121.44   19 61
12
>> >     -9999.00    14.5
>> >     ......
>> >     ---------------
>> >     This date set is for 24-hr accumulated rainfall from 20110828
>> 000000UTC
>> >     to 20110829 000000UTC.
>> >
>> >
>> >     For 12hr accumulate rainfall data.  The data is excerpted as
>> following.
>> >     ---------------
>> >     ADPSFC   C1I230   20110828_120000   23.9639  120.8372 917
61 12
>> >     -9999.00    21.5
>> >     ADPSFC   466940   20110828_120000   25.1348  121.7321 27 61
12
>> -9999.00
>> >        29
>> >     ADPSFC   466900   20110828_120000   25.1656  121.44   19 61
12
>> -9999.00
>> >        7
>> >     ......
>> >     ADPSFC   C1I230   20110829_000000   23.9639  120.8372 917
61 12
>> >     -9999.00    14
>> >     ADPSFC   466940   20110829_000000   25.1348  121.7321 27 61
12
>> -9999.00
>> >        51
>> >     ADPSFC   466900   20110829_000000   25.1656  121.44   19 61
12
>> -9999.00
>> >        7.5
>> >     ---------------
>> >
>> >
>> >     2. Using the ascii2nc to produced the observation data in
netcdf
>> format.
>> >     For 24-hr accumulated rainfall data, the netcdf file is *
>> >     precip24_2011082900.nc*
>> >     For 12-hr accumulated rainfall data, the netcdf file is *
>> >     precip24_20110828_0829_12hr.nc*
>> >
>> >     3. Transfer the WRF-ARW output by using UPPV1.0
>> >     My output data is saved as one file per hour.
>> >     So, it produced several file from UPP program.
>> >
>> >     for example:* wrfprs_d03.000*, *wrfprs_d03.001*,
*wrfprs_d03.002*
>> .......
>> >     *wrfprs_d03.024*
>> >
>> >     4. If I want to verify the first 24hr accumulated rainfall,
by the
>> 24-hr
>> >     observational accumulated rainfall
>> >
>> >     The config file was excerpted as following
>> >     ----------------
>> >     fcst_field[] = [ "APCP/A24" , "APCP/A24" , "APCP/A24" ,
"APCP/A24" ];
>> >     obs_field[]  = [ "61/A24" , "61/A24" , "61/A24" , "61/A24" ];
>> >     ----------------
>> >
>> >     ../bin/point_stat \
>> >     ../../POST_UPP/M15_2011082800/*wrfprs_d03.024 *\
>> >     ../out/ascii2nc/*precip24_2011082800.nc *\
>> >     config/PointStat_Rain \
>> >     -outdir ../out/point_stat -v 2
>> >
>> >     5. If I want to verify the first 24hr accumulated rainfall,
by the
>> 12-hr
>> >     observational accumulated rainfall
>> >
>> >     The config file was excerpted as following
>> >     ----------------
>> >     fcst_field[] = [ "APCP/A24" , "APCP/A24" , "APCP/A24" ,
"APCP/A24" ];
>> >     obs_field[]  = [ "61/A12" , "61/A12" , "61/A12" , "61/A12" ];
>> >     ----------------
>> >
>> >     ../bin/point_stat \
>> >     ../../POST_UPP/M15_2011082800/*wrfprs_d03.024* \
>> >     ../out/ascii2nc/*precip24_20110828_0829_12hr.nc* \
>> >     config/PointStat_Rain \
>> >     -outdir ../out/point_stat -v 2
>> >
>> >
>> > My question is:
>> >
>> >     - Does there have some procedure I am doing wrong?
>> >     - The step 4 and step 5 should be identical, but the results
shows
>> that
>> >     the step 5 only compare with 12hr accumulate.
>> >     How can I using the 12-hr observational data to verify the
24-hr
>> >     accumulated rainfall?
>> >     - If I want to verify the continues hourly rainfall data with
the
>> >     observational hourly rainfall data.
>> >     How do I put the fcst_file in the execute line?  (Since I
have
>> several
>> >     output file. )
>> >
>> >
>> > Thank you very much.
>> >
>>
>>
>>
>
>
> --
> Delia, Yen-Chu Chen
> Email : deliachen at ntu.edu.tw
>



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


More information about the Met_help mailing list