[Met_help] pb2nc: -pbfile option not working right

Case, Jonathan (MSFC)[] Jonathan.Case-1 at nasa.gov
Wed Feb 18 09:40:38 MST 2009


Thanks.  I must have copied the parameter over from the old config file,
overlooking the fact that the name changed. 
Jonathan

> -----Original Message-----
> From: John Halley Gotway [mailto:johnhg at rap.ucar.edu]
> Sent: Wednesday, February 18, 2009 10:32 AM
> To: Case, Jonathan (MSFC)[]; met_help at ucar.edu
> Subject: Re: pb2nc: -pbfile option not working right
> 
> Jonathan,
> 
> It would appear that we've renamed the "mask_stations" parameter to
> "mask_sid".  And the config file parser is complaining that it can't
> find an entry for "mask_sid".  So in your PointStat config
> file, change:
> 
> mask_stations = "MET_BASE/data/stations/KHSV.stations";
> to
> mask_sid = "MET_BASE/data/stations/KHSV.stations";
> 
> I suppose we wanted to emphasize that this should contain station
ID's.
> If you look in the comments in
> METv2.0beta6/data/config/PointStatConfig_default, it should list
> "mask_sid".
> 
> Sorry for the confusion.
> 
> John
> 
> Case, Jonathan (MSFC)[] wrote:
> > Thanks John,
> >
> > I have pb2nc working now for v2Beta.
> >
> > Now I'm having a problem in using mask_stations, which worked in
v1.1
> >
> > In the config file, I have set:
> > mask_stations = "MET_BASE/data/stations/KHSV.stations";
> >
> > This file simply contains "KHSV" in it to process a single station.
> >
> > When I ran point_stat, I received the following error:
> > lc2:/raid1/casejl/MET/POINTSTAT > point_stat WRFPRS15.tm00.d01
> > gdas1.20080601.SFONLY.WRF.nc PointStatConfig_test_v2 -valid_beg
> > 20080601_174500 -valid_end 20080601_181500 -v 2
> >
> >
> >   point_stat_Conf::read(const char *) -> can't get symbol table
entry
> > for variable "mask_sid"
> >
> > Any idea what this means?
> > Jonathan
> >
> >> -----Original Message-----
> >> From: John Halley Gotway [mailto:johnhg at rap.ucar.edu]
> >> Sent: Tuesday, February 17, 2009 10:15 PM
> >> To: Case, Jonathan (MSFC)[]
> >> Cc: met_help at ucar.edu
> >> Subject: RE: pb2nc: -pbfile option not working right
> >>
> >> Jonathan,
> >>
> >> Looks like I have too many versions of code running around!  I've
> >> posted this fix but only for METv1.1.  So I applied it to the
> >> METv2.0beta6 version and have attached it to this message.  Also,
> > it'll
> >> be included in the next beta release.
> >>
> >> Please use the attached version of pb2nc.cc to replace the file
> >> METv2.0beta6/src/pb2nc/pb2nc.cc.
> >>
> >> Hopefully that'll do the trick.
> >>
> >> Thanks,
> >> John
> >>
> >>> Hi again John,
> >>>
> >>> Could you possibly direct me to the bug fix in pb2ncc.cc because I
> >>> encountered the same error in the v2 Beta pb2nc when processing
> >>> multiple
> >>> pbfiles:
> >>>
> >>> ERROR: process_pbfile() -> the observation time should remain the
> >> same
> >>> for all PrepBufr messages: -62169984000 != 1212300000
> >>>
> >>> If it's a relatively simply fix, I can make it myself.
> >>> Thanks!
> >>> Jon
> >>>
> >>>
> >>>> -----Original Message-----
> >>>> From: John Halley Gotway [mailto:johnhg at rap.ucar.edu]
> >>>> Sent: Thursday, February 12, 2009 1:00 PM
> >>>> To: Case, Jonathan (MSFC)[]; met_help at ucar.edu
> >>>> Subject: Re: pb2nc: -pbfile option not working right
> >>>>
> >>>> Jonathan,
> >>>>
> >>>> Good news.  I found a bug in PB2NC which I think is causing this
> >>>> problem.
> >>>>
> >>>> When PB2NC reads the first PrepBufr file, it figures out how many
> >>>> messages that file contains.  Then it attempts to read that
number
> >> of
> >>>> messages from all subsequent PrepBufr files, rather than figuring
> >> out
> >>>> how many messages are actually in those files.  I think the error
> >>>> you're seeing is PB2NC attempting to read more messages than are
> >>>> actually in the file!
> >>>>
> >>>> I'll work up a fix for this, post it to the MET website, and send
> > it
> >>> to
> >>>> you.
> >>>>
> >>>> Thanks for finding this issue, and sorry for the inconvenience.
> >> I'll
> >>>> write you when I have the fix ready.
> >>>>
> >>>> Thanks,
> >>>> John
> >>>>
> >>>> Case, Jonathan (MSFC)[] wrote:
> >>>>> John,
> >>>>>
> >>>>>
> >>>>>
> >>>>> I'm having issues trying to process multiple PREPBUFR files into
> > a
> >>>>> single daily netcdf file using pb2nc.
> >>>>>
> >>>>>
> >>>>>
> >>>>> So far, I've only had one successfully attempt when I first
> >>> processed
> >>>>> the 12z gdas1 file, then the 00z, 06z, 18z, and next day's 00z,
> >>>>> each using the -pbfile option.  [I had to specify the beg_ds and
> >>>>> end_ds
> >>> to
> >>>> be
> >>>>> +/- 12 hours (43200 seconds) in order for this to work.  Smaller
> >>>> values
> >>>>> assigned to beg_ds and end_ds resulted in a failure.]  Even so,
> > it
> >>>> does
> >>>>> not seem as if the -valid_begin and -valid_end options have any
> >>>> effect.
> >>>>> In my data, I had a few OBS times from before the specified -
> >>>> valid_begin
> >>>>> time of 20080531_234500:
> >>>>>
> >>>>>
> >>>>>
> >>>>> lc2:/raid1/casejl/MET/PB2NC > pb2nc
> >>>> gdas1.20080601.t12z.prepbufr.blk.nr
> >>>>> test_daily.nc PB2NCConfig_SFONLY_WRF -valid_begin
20080531_234500
> >>>>> -valid_end 20080601_234500 -pbfile
> >>>> gdas1.20080601.t00z.prepbufr.blk.nr
> >>>>> -pbfile gdas1.20080601.t06z.prepbufr.blk.nr -pbfile
> >>>>> gdas1.20080601.t18z.prepbufr.blk.nr -pbfile
> >>>>> gdas1.20080602.t00z.prepbufr.blk.nr
> >>>>>
> >>>>> Reading Config File:    PB2NCConfig_SFONLY_WRF
> >>>>>
> >>>>> Creating NetCDF File:   test_daily.nc
> >>>>>
> >>>>> Reading PrepBufr File:  gdas1.20080601.t12z.prepbufr.blk.nr
> >>>>>
> >>>>> Reading PrepBufr File:  gdas1.20080601.t00z.prepbufr.blk.nr
> >>>>>
> >>>>> Reading PrepBufr File:  gdas1.20080601.t06z.prepbufr.blk.nr
> >>>>>
> >>>>> Reading PrepBufr File:  gdas1.20080601.t18z.prepbufr.blk.nr
> >>>>>
> >>>>> Reading PrepBufr File:  gdas1.20080602.t00z.prepbufr.blk.nr
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>> Then I tried processing the files in temporal order as follows
> >>>> (keeping
> >>>>> everything else the same):
> >>>>>
> >>>>>
> >>>>>
> >>>>> lc2:/raid1/casejl/MET/PB2NC > pb2nc
> >>>> gdas1.20080601.t00z.prepbufr.blk.nr
> >>>>> test_daily.nc PB2NCConfig_SFONLY_WRF -valid_begin
20080531_234500
> >>>>> -valid_end 20080601_234500 -pbfile
> >>>> gdas1.20080601.t06z.prepbufr.blk.nr
> >>>>> -pbfile gdas1.20080601.t12z.prepbufr.blk.nr -pbfile
> >>>>> gdas1.20080601.t18z.prepbufr.blk.nr -pbfile
> >>>>> gdas1.20080602.t00z.prepbufr.blk.nr
> >>>>>
> >>>>> Reading Config File:    PB2NCConfig_SFONLY_WRF
> >>>>>
> >>>>> Creating NetCDF File:   test_daily.nc
> >>>>>
> >>>>> Reading PrepBufr File:  gdas1.20080601.t00z.prepbufr.blk.nr
> >>>>>
> >>>>> Reading PrepBufr File:  gdas1.20080601.t06z.prepbufr.blk.nr
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>> ERROR: process_pbfile() -> the observation time should remain
the
> >>>> same
> >>>>> for all PrepBufr messages: -62169984000 != 1212300000
> >>>>>
> >>>>> Your advice is much appreciated,
> >>>>>
> >>>>> Jon
> >>>>>
> >>>>>
> >>>>>
> >>>>> ***********************************************************
> >>>>> Jonathan Case, ENSCO, Inc.
> >>>>> Aerospace Sciences & Engineering Division Short-term Prediction
> >>>>> Research and Transition Center 320 Sparkman Drive, Room 3062
> >>>>> Huntsville, AL 35805-1912
> >>>>> Voice: (256) 961-7504   Fax: (256) 961-7788
> >>>>> Emails: Jonathan.Case-1 at nasa.gov
> >>>>>
> >>>>>              case.jonathan at ensco.com
> >>>>>
> >>>>> ***********************************************************
> >>>>>
> >>>>>
> >>>>>
> >>>>>


More information about the Met_help mailing list