[Met_help] [rt.rap.ucar.edu #65926] History for Need help to convert PREPBUFR format data using MET tool pb2nc

John Halley Gotway via RT met_help at ucar.edu
Mon Jun 2 14:16:03 MDT 2014


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

Dear Met Help Desk,

I am trying to convert NCEP ADP Global Upper Air and Surface Weather  
Observations (PREPBUFR format) data from CISL research data archive (   
http://rda.ucar.edu/datasets/ds337.0/index.html#!description) using  
DTC's MET tool pb2nc but having difficulties to success although I am  
able to finish converting the example on-line tutorial data set  
(http://www.dtcenter.org/met/users/support/online_tutorial/METv3.0/pb2nc/run.php).

I have downloaded and extracted the above data set and get a file with  
name  "prepbufr.gdas.2002072306.wo40". Could you please tell me how  
can I convert the above data set into NetCDF format using the pb2nc  
tool. Any help would be greatly appreciated!


Thanks
Kamal
PhD Candidate
Earth and Environmental Sciences
University of Waterloo
Canada








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

Subject: Re: [rt.rap.ucar.edu #65926] Need help to convert PREPBUFR format data using MET tool pb2nc
From: John Halley Gotway
Time: Fri Mar 21 09:58:01 2014

Kamal,

I went to the CISL website and retrieved this tar file:
    prepbufr.20020723.wo40.tar.gz

I extracted the file using this command:
    tar -xvzf prepbufr.20020723.wo40.tar.gz

And then I ran the PB2NC utility using the following command:
    METv4.1/bin/pb2nc \
      20020723.wo40/prepbufr.gdas.2002072306.wo40 \
      20020723.wo40/prepbufr.gdas.2002072306.wo40.nc \
      METv4.1/data/config/PB2NCConfig_default

It all worked fine, but did take about 3.5 minutes to run.  Here's the
output I see printed to the screen:

DEBUG 1: Default Config File: /d1/johnhg/MET/MET_development/svn-met-
dev.cgd.ucar.edu/trunk/met/data/config/PB2NCConfig_default
DEBUG 1: User Config File: PB2NCConfig_default
DEBUG 1: Creating NetCDF File:
prepbufr.gdas.2002072306.wo40.nc
DEBUG 1: Processing PrepBufr File:      prepbufr.gdas.2002072306.wo40
DEBUG 1: Blocking PrepBufr file to:     /tmp/tmp_pb2nc_blk_10827_0
DEBUG 2: PrepBufr Time Center:          20020723_060000
DEBUG 2: Searching Time Window:         20020723_043000 to
20020723_073000
DEBUG 2: Processing 267962 PrepBufr messages...
5% 10% 15% 20% 25% 30% 35% 40% 45% 50% 55% 60% 65% 70% 75% 80% 85% 90%
95% 100%
DEBUG 2: Total PrepBufr Messages processed      = 267962
DEBUG 2: Rejected based on message type         = 0
DEBUG 2: Rejected based on station id           = 0
DEBUG 2: Rejected based on valid time           = 113337
DEBUG 2: Rejected based on masking grid         = 0
DEBUG 2: Rejected based on masking polygon      = 0
DEBUG 2: Rejected based on elevation            = 0
DEBUG 2: Rejected based on pb report type       = 0
DEBUG 2: Rejected based on input report type    = 0
DEBUG 2: Rejected based on instrument type      = 0
DEBUG 2: Rejected based on zero observations    = 73203
DEBUG 2: Total PrepBufr Messages retained       = 81422
DEBUG 2: Total observations retained or derived = 211818

However, I see that you're running the METv3.0 tutorial - so I tried
doing the same thing with the METv3.0 version of pb2nc.  It also ran
fine and took about the same amount of time to run.

Two things...
(1) If possible, I'd suggest using the latest available version of
MET, METv4.1.  There have been several significant bugfixes and
updates since METv3.0.
     You can download it here:
http://www.dtcenter.org/met/users/downloads/index.php
     And here's the corresponding tutorial:
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/index.php

(2) If you're still having trouble running PB2NC, please send a
detailed description of the problem, and I'll do my best to help.

Thanks,
John Halley Gotway
met_help at ucar.edu





On 03/20/2014 04:28 PM, mmkamal via RT wrote:
>
> Thu Mar 20 16:28:05 2014: Request 65926 was acted upon.
> Transaction: Ticket created by mmkamal at uwaterloo.ca
>         Queue: met_help
>       Subject: Need help to convert PREPBUFR format data using MET
tool pb2nc
>         Owner: Nobody
>    Requestors: mmkamal at uwaterloo.ca
>        Status: new
>   Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65926 >
>
>
> Dear Met Help Desk,
>
> I am trying to convert NCEP ADP Global Upper Air and Surface Weather
> Observations (PREPBUFR format) data from CISL research data archive
(
> http://rda.ucar.edu/datasets/ds337.0/index.html#!description) using
> DTC's MET tool pb2nc but having difficulties to success although I
am
> able to finish converting the example on-line tutorial data set
>
(http://www.dtcenter.org/met/users/support/online_tutorial/METv3.0/pb2nc/run.php).
>
> I have downloaded and extracted the above data set and get a file
with
> name  "prepbufr.gdas.2002072306.wo40". Could you please tell me how
> can I convert the above data set into NetCDF format using the pb2nc
> tool. Any help would be greatly appreciated!
>
>
> Thanks
> Kamal
> PhD Candidate
> Earth and Environmental Sciences
> University of Waterloo
> Canada
>
>
>
>
>
>

------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #65926] Need help to convert PREPBUFR format data using MET tool pb2nc
From: mmkamal
Time: Fri Mar 21 18:37:13 2014

Hi John,


Thanks for your quick reply and suggestion. Unfortunately I have less
expertise in software installation in High Performance Computing
Cluster (SciNet, operated by university of Toronto). But, definitely I
will give a try next month.  Fortunately I am able to get an output
using my existing one and when check I got the following header
information from my output. During compilation I get the following
information:


===================================================
gpc-f101n084-$ ./pb2nc
/scratch/j/jcl/mkamal/DATA/Observation/Upper_Air_Data/July/20020723.wo40/prepbufr.gdas.2002072312.wo40
/scratch/j/jcl/mkamal/METv3.0.1/tutorial/out/pb2nc/gdas.2002072312.nc
/scratch/j/jcl/mkamal/METv3.0.1/tutorial/config/PB2NCConfig_default
Reading Config File:
/scratch/j/jcl/mkamal/METv3.0.1/tutorial/config/PB2NCConfig_default
Creating NetCDF File:
/scratch/j/jcl/mkamal/METv3.0.1/tutorial/out/pb2nc/gdas.2002072312.nc
Processing PrepBufr File:
/scratch/j/jcl/mkamal/DATA/Observation/Upper_Air_Data/July/20020723.wo40/prepbufr.gdas.2002072312.wo40

===============================================================

And my output file contain the following header information
===============================================================

gpc-f101n084-$
/scinet/gpc/Libraries/netcdf-4.1.3/serial_gcc/bin/ncdump -h
netcdf.gdas.2002072312.nc

netcdf netcdf.gdas.2002072312 {
dimensions:
          mxstr = 16 ;
          hdr_arr_len = 3 ;
          obs_arr_len = 5 ;
          nobs = UNLIMITED ; // (56403 currently)
          nhdr = 9253 ;
variables:
          float obs_arr(nobs, obs_arr_len) ;
                  obs_arr:long_name = "array of observation values" ;
                  obs_arr:_fill_value = -9999.f ;
                  obs_arr:columns = "hdr_id gc lvl hgt ob" ;
                  obs_arr:hdr_id_long_name = "index of matching header
data" ;
                  obs_arr:gc_long_name = "grib code corresponding to
the observation type" ;
                  obs_arr:lvl_long_name = "pressure level (hPa) or
accumulation interval (sec)" ;
                  obs_arr:hgt_long_name = "height in meters above sea
level (msl)" ;
                  obs_arr:ob_long_name = "observation value" ;
          char hdr_typ(nhdr, mxstr) ;
                  hdr_typ:long_name = "message type" ;
          char hdr_sid(nhdr, mxstr) ;
                  hdr_sid:long_name = "station identification" ;
          char hdr_vld(nhdr, mxstr) ;
                  hdr_vld:long_name = "valid time" ;
                  hdr_vld:units = "YYYYMMDD_HHMMSS" ;
          float hdr_arr(nhdr, hdr_arr_len) ;
                  hdr_arr:long_name = "array of observation station
header values" ;
                  hdr_arr:_fill_value = -9999.f ;
                  hdr_arr:columns = "lat lon elv" ;
                  hdr_arr:lat_long_name = "latitude" ;
                  hdr_arr:lat_units = "degrees_north" ;
                  hdr_arr:lon_long_name = "longitude" ;
                  hdr_arr:lon_units = "degrees_east" ;
                  hdr_arr:elv_long_name = "elevation" ;
                  hdr_arr:elv_units = "meters above sea level (msl)" ;

// global attributes:
                  :FileOrigins = "File
/scratch/j/jcl/mkamal/METv3.0.1/tutorial/out/pb2nc/gdas.2002072312.nc
generated 20140321_233334 UTC on host gpc-f101n084 by the MET pb2nc
tool" ;
                  :MET_version = "V3.0.1" ;
                  :MET_tool = "pb2nc" ;
}

===========================================================================

I was wondering whether I need to change any configuration inside the
"PB2NCConfig_default" file to fit my data set or not? Because, when I
tried to use the output NetCDF data set to plot  skew-T log-p diagram
for Buffalo (72528), NY, USA
station using the following NCL script
(http://www.ncl.ucar.edu/Support/talk_archives/2012/att-
2060/skewt_station.ncl) I get the following error
message:

==================================================
==============================================


gpc-f101n084-$ ncl 'ncfile="netcdf.gdas.2002072312.nc"'
'station_icao="KBUF"' 'station_synop="72528"' skewt_station.ncl

   Copyright (C) 1995-2012 - All Rights Reserved
   University Corporation for Atmospheric Research
   NCAR Command Language Version 6.1.0
   The use of this software is governed by a License Agreement.
   See http://www.ncl.ucar.edu/ for more details.


Variable: where_station
Type: integer
Total Size: 4 bytes
              1 values
Number of Dimensions: 1
Dimensions and sizes:   [1]
Coordinates:
Number Of Attributes: 1
    _FillValue :  -2147483647
(0)     -2147483647


Variable: count_station
Type: integer
Total Size: 4 bytes
              1 values
Number of Dimensions: 1
Dimensions and sizes:   [1]
Coordinates:
(0)     1

Variable: obs_arr
Type: float
Total Size: 1128060 bytes
              282015 values
Number of Dimensions: 2
Dimensions and sizes:   [nobs | 56403] x [obs_arr_len | 5]
Coordinates:
Number Of Attributes: 8
    long_name :   array of observation values
    _fill_value : -9999
    columns :     hdr_id gc lvl hgt ob
    hdr_id_long_name :    index of matching header data
    gc_long_name :        grib code corresponding to the observation
type
    lvl_long_name :       pressure level (hPa) or accumulation
interval (sec)
    hgt_long_name :       height in meters above sea level (msl)
    ob_long_name :        observation value

Variable: obs
Type: integer
Total Size: 4 bytes
              1 values
Number of Dimensions: 1
Dimensions and sizes:   [1]
Coordinates:
Number Of Attributes: 1
    _FillValue :  -2147483647
fatal:Subscript out of range, error in subscript #0
fatal:An error occurred reading hdr_vld
fatal:["Execute.c":7743]:Execute: Error occurred at or near line 78 in
file skewt_station.ncl
==============================================================================


I look forward to hearing from you.


Thanks
Kamal





Quoting John Halley Gotway via RT <met_help at ucar.edu>:

> Kamal,
>
> I went to the CISL website and retrieved this tar file:
>     prepbufr.20020723.wo40.tar.gz
>
> I extracted the file using this command:
>     tar -xvzf prepbufr.20020723.wo40.tar.gz
>
> And then I ran the PB2NC utility using the following command:
>     METv4.1/bin/pb2nc \
>       20020723.wo40/prepbufr.gdas.2002072306.wo40 \
>       20020723.wo40/prepbufr.gdas.2002072306.wo40.nc \
>       METv4.1/data/config/PB2NCConfig_default
>
> It all worked fine, but did take about 3.5 minutes to run.  Here's
> the output I see printed to the screen:
>
> DEBUG 1: Default Config File:
> /d1/johnhg/MET/MET_development/svn-met-
dev.cgd.ucar.edu/trunk/met/data/config/PB2NCConfig_default
> DEBUG 1: User Config File: PB2NCConfig_default
> DEBUG 1: Creating NetCDF File:
prepbufr.gdas.2002072306.wo40.nc
> DEBUG 1: Processing PrepBufr File:
prepbufr.gdas.2002072306.wo40
> DEBUG 1: Blocking PrepBufr file to:     /tmp/tmp_pb2nc_blk_10827_0
> DEBUG 2: PrepBufr Time Center:          20020723_060000
> DEBUG 2: Searching Time Window:         20020723_043000 to
20020723_073000
> DEBUG 2: Processing 267962 PrepBufr messages...
> 5% 10% 15% 20% 25% 30% 35% 40% 45% 50% 55% 60% 65% 70% 75% 80% 85%
> 90% 95% 100%
> DEBUG 2: Total PrepBufr Messages processed      = 267962
> DEBUG 2: Rejected based on message type         = 0
> DEBUG 2: Rejected based on station id           = 0
> DEBUG 2: Rejected based on valid time           = 113337
> DEBUG 2: Rejected based on masking grid         = 0
> DEBUG 2: Rejected based on masking polygon      = 0
> DEBUG 2: Rejected based on elevation            = 0
> DEBUG 2: Rejected based on pb report type       = 0
> DEBUG 2: Rejected based on input report type    = 0
> DEBUG 2: Rejected based on instrument type      = 0
> DEBUG 2: Rejected based on zero observations    = 73203
> DEBUG 2: Total PrepBufr Messages retained       = 81422
> DEBUG 2: Total observations retained or derived = 211818
>
> However, I see that you're running the METv3.0 tutorial - so I tried
> doing the same thing with the METv3.0 version of pb2nc.  It also ran
> fine and took about the same amount of time to run.
>
> Two things...
> (1) If possible, I'd suggest using the latest available version of
> MET, METv4.1.  There have been several significant bugfixes and
> updates since METv3.0.
>      You can download it here:
> http://www.dtcenter.org/met/users/downloads/index.php
>      And here's the corresponding tutorial:
>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/index.php
>
> (2) If you're still having trouble running PB2NC, please send a
> detailed description of the problem, and I'll do my best to help.
>
> Thanks,
> John Halley Gotway
> met_help at ucar.edu
>
>
>
>
>
> On 03/20/2014 04:28 PM, mmkamal via RT wrote:
>>
>> Thu Mar 20 16:28:05 2014: Request 65926 was acted upon.
>> Transaction: Ticket created by mmkamal at uwaterloo.ca
>>         Queue: met_help
>>       Subject: Need help to convert PREPBUFR format data using MET
>> tool pb2nc
>>         Owner: Nobody
>>    Requestors: mmkamal at uwaterloo.ca
>>        Status: new
>>   Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65926 >
>>
>>
>> Dear Met Help Desk,
>>
>> I am trying to convert NCEP ADP Global Upper Air and Surface
Weather
>> Observations (PREPBUFR format) data from CISL research data archive
(
>> http://rda.ucar.edu/datasets/ds337.0/index.html#!description) using
>> DTC's MET tool pb2nc but having difficulties to success although I
am
>> able to finish converting the example on-line tutorial data set
>>
(http://www.dtcenter.org/met/users/support/online_tutorial/METv3.0/pb2nc/run.php).
>>
>> I have downloaded and extracted the above data set and get a file
with
>> name  "prepbufr.gdas.2002072306.wo40". Could you please tell me how
>> can I convert the above data set into NetCDF format using the pb2nc
>> tool. Any help would be greatly appreciated!
>>
>>
>> Thanks
>> Kamal
>> PhD Candidate
>> Earth and Environmental Sciences
>> University of Waterloo
>> Canada
>>
>>
>>
>>
>>
>>
>
>






------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #65926] Need help to convert PREPBUFR format data using MET tool pb2nc
From: John Halley Gotway
Time: Mon Mar 24 10:16:20 2014

Kamal,

I'm sorry but I'm really not familiar with that NCL script you
referenced.  I did try running it using the output of the PB2NC tool
for MET versions 3.0, 3.0.1, 3.1, 4.0, and 4.1.  All of them
resulted in an error.  So there may be a problem in that script.

If you're just trying to see if PB2NC successfully ran, you could run
the output through the plot_point_obs utility.  That tool creates a
postscript plot with a red dot at every lat/lon where there
was an observation.  It'll show you the locations of the observations.
You can run it like this:

   METv3.0.1/bin/plot_point_obs netcdf.gdas.2002072312.nc
netcdf.gdas.2002072312.ps

As for changing the default configuration settings, it's really up to
you to decide based on the type of verification you're trying to do.
I'd suggest reading through the comments in the PB2NC config
file for an explanation of each option.  The GDAS dataset contains
global observations, but you're probably not doing verification over
the whole globe!  I'd suggest deciding the area over which you
want to perform verification.  Then set the "mask_grid" or "mask_poly"
option in the config file so as to only retain observations that fall
within your area of interest.  Also look at the
"obs_grib_code" setting and check to see if that's the list of
observation types you want to retain.

Hope that helps.

Thanks,
John

On 03/21/2014 06:37 PM, mmkamal via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65926 >
>
> Hi John,
>
>
> Thanks for your quick reply and suggestion. Unfortunately I have
less
> expertise in software installation in High Performance Computing
> Cluster (SciNet, operated by university of Toronto). But, definitely
I
> will give a try next month.  Fortunately I am able to get an output
> using my existing one and when check I got the following header
> information from my output. During compilation I get the following
> information:
>
>
> ===================================================
> gpc-f101n084-$ ./pb2nc
>
/scratch/j/jcl/mkamal/DATA/Observation/Upper_Air_Data/July/20020723.wo40/prepbufr.gdas.2002072312.wo40
/scratch/j/jcl/mkamal/METv3.0.1/tutorial/out/pb2nc/gdas.2002072312.nc
> /scratch/j/jcl/mkamal/METv3.0.1/tutorial/config/PB2NCConfig_default
> Reading Config File:
> /scratch/j/jcl/mkamal/METv3.0.1/tutorial/config/PB2NCConfig_default
> Creating NetCDF File:
>
/scratch/j/jcl/mkamal/METv3.0.1/tutorial/out/pb2nc/gdas.2002072312.nc
> Processing PrepBufr File:
>
/scratch/j/jcl/mkamal/DATA/Observation/Upper_Air_Data/July/20020723.wo40/prepbufr.gdas.2002072312.wo40
>
> ===============================================================
>
> And my output file contain the following header information
> ===============================================================
>
> gpc-f101n084-$
> /scinet/gpc/Libraries/netcdf-4.1.3/serial_gcc/bin/ncdump -h
> netcdf.gdas.2002072312.nc
>
> netcdf netcdf.gdas.2002072312 {
> dimensions:
>            mxstr = 16 ;
>            hdr_arr_len = 3 ;
>            obs_arr_len = 5 ;
>            nobs = UNLIMITED ; // (56403 currently)
>            nhdr = 9253 ;
> variables:
>            float obs_arr(nobs, obs_arr_len) ;
>                    obs_arr:long_name = "array of observation values"
;
>                    obs_arr:_fill_value = -9999.f ;
>                    obs_arr:columns = "hdr_id gc lvl hgt ob" ;
>                    obs_arr:hdr_id_long_name = "index of matching
header data" ;
>                    obs_arr:gc_long_name = "grib code corresponding
to
> the observation type" ;
>                    obs_arr:lvl_long_name = "pressure level (hPa) or
> accumulation interval (sec)" ;
>                    obs_arr:hgt_long_name = "height in meters above
sea
> level (msl)" ;
>                    obs_arr:ob_long_name = "observation value" ;
>            char hdr_typ(nhdr, mxstr) ;
>                    hdr_typ:long_name = "message type" ;
>            char hdr_sid(nhdr, mxstr) ;
>                    hdr_sid:long_name = "station identification" ;
>            char hdr_vld(nhdr, mxstr) ;
>                    hdr_vld:long_name = "valid time" ;
>                    hdr_vld:units = "YYYYMMDD_HHMMSS" ;
>            float hdr_arr(nhdr, hdr_arr_len) ;
>                    hdr_arr:long_name = "array of observation station
> header values" ;
>                    hdr_arr:_fill_value = -9999.f ;
>                    hdr_arr:columns = "lat lon elv" ;
>                    hdr_arr:lat_long_name = "latitude" ;
>                    hdr_arr:lat_units = "degrees_north" ;
>                    hdr_arr:lon_long_name = "longitude" ;
>                    hdr_arr:lon_units = "degrees_east" ;
>                    hdr_arr:elv_long_name = "elevation" ;
>                    hdr_arr:elv_units = "meters above sea level
(msl)" ;
>
> // global attributes:
>                    :FileOrigins = "File
>
/scratch/j/jcl/mkamal/METv3.0.1/tutorial/out/pb2nc/gdas.2002072312.nc
> generated 20140321_233334 UTC on host gpc-f101n084 by the MET pb2nc
> tool" ;
>                    :MET_version = "V3.0.1" ;
>                    :MET_tool = "pb2nc" ;
> }
>
>
===========================================================================
>
> I was wondering whether I need to change any configuration inside
the
> "PB2NCConfig_default" file to fit my data set or not? Because, when
I
> tried to use the output NetCDF data set to plot  skew-T log-p
diagram
> for Buffalo (72528), NY, USA
> station using the following NCL script
> (http://www.ncl.ucar.edu/Support/talk_archives/2012/att-
2060/skewt_station.ncl) I get the following error
> message:
>
> ==================================================
> ==============================================
>
>
> gpc-f101n084-$ ncl 'ncfile="netcdf.gdas.2002072312.nc"'
> 'station_icao="KBUF"' 'station_synop="72528"' skewt_station.ncl
>
>     Copyright (C) 1995-2012 - All Rights Reserved
>     University Corporation for Atmospheric Research
>     NCAR Command Language Version 6.1.0
>     The use of this software is governed by a License Agreement.
>     See http://www.ncl.ucar.edu/ for more details.
>
>
> Variable: where_station
> Type: integer
> Total Size: 4 bytes
>                1 values
> Number of Dimensions: 1
> Dimensions and sizes:   [1]
> Coordinates:
> Number Of Attributes: 1
>      _FillValue :  -2147483647
> (0)     -2147483647
>
>
> Variable: count_station
> Type: integer
> Total Size: 4 bytes
>                1 values
> Number of Dimensions: 1
> Dimensions and sizes:   [1]
> Coordinates:
> (0)     1
>
> Variable: obs_arr
> Type: float
> Total Size: 1128060 bytes
>                282015 values
> Number of Dimensions: 2
> Dimensions and sizes:   [nobs | 56403] x [obs_arr_len | 5]
> Coordinates:
> Number Of Attributes: 8
>      long_name :   array of observation values
>      _fill_value : -9999
>      columns :     hdr_id gc lvl hgt ob
>      hdr_id_long_name :    index of matching header data
>      gc_long_name :        grib code corresponding to the
observation type
>      lvl_long_name :       pressure level (hPa) or accumulation
interval (sec)
>      hgt_long_name :       height in meters above sea level (msl)
>      ob_long_name :        observation value
>
> Variable: obs
> Type: integer
> Total Size: 4 bytes
>                1 values
> Number of Dimensions: 1
> Dimensions and sizes:   [1]
> Coordinates:
> Number Of Attributes: 1
>      _FillValue :  -2147483647
> fatal:Subscript out of range, error in subscript #0
> fatal:An error occurred reading hdr_vld
> fatal:["Execute.c":7743]:Execute: Error occurred at or near line 78
in
> file skewt_station.ncl
>
==============================================================================
>
>
> I look forward to hearing from you.
>
>
> Thanks
> Kamal
>
>
>
>
>
> Quoting John Halley Gotway via RT <met_help at ucar.edu>:
>
>> Kamal,
>>
>> I went to the CISL website and retrieved this tar file:
>>      prepbufr.20020723.wo40.tar.gz
>>
>> I extracted the file using this command:
>>      tar -xvzf prepbufr.20020723.wo40.tar.gz
>>
>> And then I ran the PB2NC utility using the following command:
>>      METv4.1/bin/pb2nc \
>>        20020723.wo40/prepbufr.gdas.2002072306.wo40 \
>>        20020723.wo40/prepbufr.gdas.2002072306.wo40.nc \
>>        METv4.1/data/config/PB2NCConfig_default
>>
>> It all worked fine, but did take about 3.5 minutes to run.  Here's
>> the output I see printed to the screen:
>>
>> DEBUG 1: Default Config File:
>> /d1/johnhg/MET/MET_development/svn-met-
dev.cgd.ucar.edu/trunk/met/data/config/PB2NCConfig_default
>> DEBUG 1: User Config File: PB2NCConfig_default
>> DEBUG 1: Creating NetCDF File:
prepbufr.gdas.2002072306.wo40.nc
>> DEBUG 1: Processing PrepBufr File:
prepbufr.gdas.2002072306.wo40
>> DEBUG 1: Blocking PrepBufr file to:     /tmp/tmp_pb2nc_blk_10827_0
>> DEBUG 2: PrepBufr Time Center:          20020723_060000
>> DEBUG 2: Searching Time Window:         20020723_043000 to
20020723_073000
>> DEBUG 2: Processing 267962 PrepBufr messages...
>> 5% 10% 15% 20% 25% 30% 35% 40% 45% 50% 55% 60% 65% 70% 75% 80% 85%
>> 90% 95% 100%
>> DEBUG 2: Total PrepBufr Messages processed      = 267962
>> DEBUG 2: Rejected based on message type         = 0
>> DEBUG 2: Rejected based on station id           = 0
>> DEBUG 2: Rejected based on valid time           = 113337
>> DEBUG 2: Rejected based on masking grid         = 0
>> DEBUG 2: Rejected based on masking polygon      = 0
>> DEBUG 2: Rejected based on elevation            = 0
>> DEBUG 2: Rejected based on pb report type       = 0
>> DEBUG 2: Rejected based on input report type    = 0
>> DEBUG 2: Rejected based on instrument type      = 0
>> DEBUG 2: Rejected based on zero observations    = 73203
>> DEBUG 2: Total PrepBufr Messages retained       = 81422
>> DEBUG 2: Total observations retained or derived = 211818
>>
>> However, I see that you're running the METv3.0 tutorial - so I
tried
>> doing the same thing with the METv3.0 version of pb2nc.  It also
ran
>> fine and took about the same amount of time to run.
>>
>> Two things...
>> (1) If possible, I'd suggest using the latest available version of
>> MET, METv4.1.  There have been several significant bugfixes and
>> updates since METv3.0.
>>       You can download it here:
>> http://www.dtcenter.org/met/users/downloads/index.php
>>       And here's the corresponding tutorial:
>>
http://www.dtcenter.org/met/users/support/online_tutorial/METv4.1/index.php
>>
>> (2) If you're still having trouble running PB2NC, please send a
>> detailed description of the problem, and I'll do my best to help.
>>
>> Thanks,
>> John Halley Gotway
>> met_help at ucar.edu
>>
>>
>>
>>
>>
>> On 03/20/2014 04:28 PM, mmkamal via RT wrote:
>>>
>>> Thu Mar 20 16:28:05 2014: Request 65926 was acted upon.
>>> Transaction: Ticket created by mmkamal at uwaterloo.ca
>>>          Queue: met_help
>>>        Subject: Need help to convert PREPBUFR format data using
MET
>>> tool pb2nc
>>>          Owner: Nobody
>>>     Requestors: mmkamal at uwaterloo.ca
>>>         Status: new
>>>    Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=65926 >
>>>
>>>
>>> Dear Met Help Desk,
>>>
>>> I am trying to convert NCEP ADP Global Upper Air and Surface
Weather
>>> Observations (PREPBUFR format) data from CISL research data
archive (
>>> http://rda.ucar.edu/datasets/ds337.0/index.html#!description)
using
>>> DTC's MET tool pb2nc but having difficulties to success although I
am
>>> able to finish converting the example on-line tutorial data set
>>>
(http://www.dtcenter.org/met/users/support/online_tutorial/METv3.0/pb2nc/run.php).
>>>
>>> I have downloaded and extracted the above data set and get a file
with
>>> name  "prepbufr.gdas.2002072306.wo40". Could you please tell me
how
>>> can I convert the above data set into NetCDF format using the
pb2nc
>>> tool. Any help would be greatly appreciated!
>>>
>>>
>>> Thanks
>>> Kamal
>>> PhD Candidate
>>> Earth and Environmental Sciences
>>> University of Waterloo
>>> Canada
>>>
>>>
>>>
>>>
>>>
>>>
>>
>>
>
>
>
>
>

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


More information about the Met_help mailing list