[Met_help] PB2NC config file

Ruifang Li lir at ucar.edu
Tue Jul 14 23:15:21 MDT 2009


John,
I get this error when I run point_stat tool:

*** Running POINT_STAT on WRF OUTPUT ***
***          UGRD                ***
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=1666342218
Forecast File:
/ptmp/lir/test/data/2009_t8_15km/verify_met2.0/gsi/wpp/2007081512/postprd/wrfprs_d01.000
Climatology File: none
Configuration File: ./config/PointStatConfig.UGRD
Observation File:
/ptmp/lir/test/data/2009_t8_15km/verify_met2.0/ob_nc/2007081512/ob.nc

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

Reading records for UGRD/P1000.
For UGRD/P1000 found 1 forecast levels and 0 climatology levels.


ERROR: is_grid_relative() -> Unsupported grid type value: ^A


But point_stat works well with TMP. I looked at code read_grib.cc and found
out

// FUTURE WORK: ADD SUPPORT FOR THE MERCATOR PROJECTION

   //
   // Mercator Projection Grid
   //

else if(r.gds->type == 1) {

      cerr << "\n\nERROR: read_gds() -> "
           << "Mercator grid type " << (int) r.gds->type
           << " not currently supported.\n\n"
           << flush;
      exit(1);
   }

My wrf forecast use Mercator projection, It looks like that I can not do U,
V, WIND using MET2.0, right?
Are there other ways that I can try?

Thanks,
Ruifang


On Tue, Jul 14, 2009 at 10:18 AM, John Halley Gotway <johnhg at rap.ucar.edu>wrote:

> Ruifang,
>
> Sorry about that, I forgot the attachment.  It's attached now.
>
> You're right that you could use a polyline to define the masking region for
> the observations you'd like to retain.  But to answer you question, does
> G181 = CAR.poly?  I'd say the answer is NO.
>
> Look in the file "METv2.0/data/poly/ncep_vx_regions.pdf", at the region
> marked "CAR".  I'm sure that doesn't line up exactly with NCEP Grid 181,
> although I can't find a plot of Grid 181 right now.
>
> You can do the masking with a grid or a polyline, whichever you prefer.
>
> Thanks,
> John
>
> Ruifang Li wrote:
> > John,
> > Thanks for quick reply.
> > One more question, If I define these variables in the config file
> >
> > mask_grid=" "
> > maks_poly="MET_BASE/data/poly/CAR.poly"
> >
> > Is it as same as
> >
> > mask_grid="G181 "
> > maks_poly=" "
> >
> > By the way, I could not see your attached file.
> >
> > Thanks again,
> > Ruifang
> >
> >
> > On Tue, Jul 14, 2009 at 6:37 AM, John Halley Gotway <johnhg at rap.ucar.edu
> >wrote:
> >
> >> Hello,
> >>
> >> You raise a good point.  Here's the webpage from which we originally
> >> extracted the grid definitions:
> >> http://www.nco.ncep.noaa.gov/pmb/docs/on388/tableb.html
> >>
> >> I see that the information on there was updated on 9/15/2008, and that
> the
> >> changes are listed in red, including the addition of grid #181.  This
> will
> >> require some code changes to grab all of these
> >> updated grid definitions.  In the meantime, I've modified the file
> >> "METv2.0/lib/vx_data_grids/grid_defs.f" to add NCEP Grid 181 to the list
> of
> >> definitions.  Please try replacing your original version
> >> of the file with the attached one, and then recompile MET being sure to
> do
> >> a "make clean" first.  Please let me know if you have any troubles using
> >> grid #181 after recompiling.
> >>
> >> When I get a chance sometime over the next month, I'll update that file
> >> "METv2.0/lib/vx_data_grids/grid_defs.f" with all of the changes and post
> it
> >> to the MET website.
> >>
> >> Thanks for raising this issue!
> >>
> >> John Halley Gotway
> >> johnhg at ucar.edu
> >>
> >> Ruifang Li wrote:
> >>> Hi, John,
> >>> I have a question about  variable mask_grid = "G212" in the default
> >> config
> >>> file of PB2NC. I tried to use mask_grid = "G181"
> >>> for Caribean domain and it says:
> >>> ERROR: process_config_file() -> the mask_grid requested "G181" is not
> >>> defined.
> >>>
> >>> But G212 does work.
> >>>
> >>> How to use mask_grid = "G181"
> >>>
> >>> Thanks,
> >>>
> >>>
> >>>
> ------------------------------------------------------------------------
> >>>
> >>> _______________________________________________
> >>> Met_help mailing list
> >>> Met_help at mailman.ucar.edu
> >>> http://mailman.ucar.edu/mailman/listinfo/met_help
> >
> >
> >
>



-- 
Ruifang

Mesoscale & Microscale Meteorology Division
Phone: 303-497-8938
Office:  FL3-3085
Email:  lir at ucar.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/met_help/attachments/20090714/b2992df2/attachment.html 


More information about the Met_help mailing list