[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

20041207: 20040928: 20040927: converting WRF grib output to GEMPAK



Michael,

would it be possible for you to post your file:
000012406_4km_wrf_f06.gem

I can see what your data looks like.

Steve Chiswell
Unidata User Support


>From: Michael Brennan <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200412071525.iB7FPclI021103

>Hi Steve,
>
>I've finally gotten GEMPAK5.7.4 to compile and install correctly.  After 
>taking your suggestion and modifying MAXPRM and MAXLVL in:
>
>$GEMPAK/source/programs/gd/gdvint/viprm.prm
>
>and typing
>
>make clean
>
>make all
>
>make install
>
>I'm still unable to get GDVINT to interpolate anymore than the hardwired 
>number of parameters.  I know I'm pointing to the correct executable and 
>everything.  I've even tried to modify vglvpr.f to get MAXPRM to write 
>out to the screen and I've had no luck with that either.
>
>  GEMPAK-GDVINT>l
>  GDFILE   = 2000012406_4km_wrf_f06.gem
>  GDOUTF   = test.gem
>  GDATTIM  = f06
>  GVCORD   = elvl:100/pres
>  GLEVEL   = 1000-100-25
>  MAXGRD   = 5000
>  GAREA    = grid
>  VCOORD   = none
>
>  GEMPAK-GDVINT>r
>   Excessive number of parameters.  Parameter QCLD         skipped.
>   Excessive number of parameters.  Parameter QRAI         skipped.
>   Excessive number of parameters.  Parameter QICE         skipped.
>   Excessive number of parameters.  Parameter QSNO         skipped.
>   Excessive number of parameters.  Parameter QGRA         skipped.
>   Excessive number of parameters.  Parameter QCLD         skipped.
>   Excessive number of parameters.  Parameter QRAI         skipped.
>   Excessive number of parameters.  Parameter QICE         skipped.
>   Excessive number of parameters.  Parameter QSNO         skipped.
>   Excessive number of parameters.  Parameter QGRA         skipped.
>   Excessive number of parameters.  Parameter QCLD         skipped.
>   Excessive number of parameters.  Parameter QRAI         skipped.
>
>Here is my viprm.prm file:
>
>C************************************************************************
>C* VIPRM.PRM
>C*
>C* This include file contains parameters used in the vertical
>C* interpolation software.
>C************************************************************************
>C!
>C!      Parameter names.
>C!
>         CHARACTER*4     NTMPC, NTMPF, NTHTA, NPSYM, NHGHT, NPRES
>         CHARACTER*4     NSPFH, NMIXR, NRMIX, NRELH, NDWPT, NDWPC
>         CHARACTER*4     NDWPK, NTMPK
>C!
>         PARAMETER       ( NTMPC = 'TMPC' )
>C!                                              Temperature in deg C
>         PARAMETER       ( NTMPF = 'TMPF' )
>C!                                              Temperature in deg F
>         PARAMETER       ( NTHTA = 'THTA' )
>C!                                              Pot. Temperature in K
>         PARAMETER       ( NPSYM = 'PSYM' )
>C!                                              Mont. Str. Func in mm/ss
>         PARAMETER       ( NHGHT = 'HGHT' )
>C!                                              Height in m
>         PARAMETER       ( NPRES = 'PRES' )
>C!                                              Press. in hpa
>         PARAMETER       ( NSPFH = 'SPFH' )
>C!                                              Specific humidity
>         PARAMETER       ( NMIXR = 'MIXR' )
>C!                                              Mixing ratio * 1000
>         PARAMETER       ( NRMIX = 'RMIX' )
>C!                                              Mixing ratio
>         PARAMETER       ( NRELH = 'RELH' )
>C!                                              Relative humidity
>         PARAMETER       ( NDWPT = 'DWPT' )
>C!                                              Dewpt temp in C
>         PARAMETER       ( NDWPC = 'DWPC' )
>C!                                              Dewpt temp in C
>         PARAMETER       ( NDWPK = 'DWPK' )
>C!                                              Dewpt temp in K
>         PARAMETER       ( NTMPK = 'TMPK' )
>C!                                              Temperature in deg K
>C------------------------------------------------------------------------
>C!
>C!      Coordinate names.
>C!
>         CHARACTER*4     NESFC, NNONE, NSGMA, NETA, NZAGL
>C!
>         PARAMETER       ( NESFC = 'ESFC' )
>C!                                              Earth's surface
>         PARAMETER       ( NNONE = 'NONE' )
>C!                                              None (surface)
>         PARAMETER       ( NSGMA = 'SGMA' )
>C!                                              Sigma
>         PARAMETER       ( NETA  = 'ELVL'  )
>C!                                              Eta
>         PARAMETER       ( NZAGL = 'ZAGL' )
>C!                                              Height above grnd lvl
>C------------------------------------------------------------------------
>C!
>C!      Miscellaneous parameters and GEMPAK parameters.
>C!
>C************************************************************************
>         INCLUDE 'GEMPRM.PRM'
>         PARAMETER ( MXGD = LLMXGD )
>C!
>C!                                              Maximum grid size
>         PARAMETER ( MAXLVL = 100 )
>C!                                              Maximum number of lvls
>         PARAMETER ( MAXPRM = 50 )
>C!                                              Maximum number of parms
>         PARAMETER ( MXSFVC = 8 )
>C!                                              Maximum number of extra
>C!                                              vertical coordinates
>         PARAMETER ( LNAV = LLNNAV, LANL = LLNANL )
>C!                                              Max nav. & anal. blk sz
>
>Michael J. Brennan
>Graduate Research Assistant
>Department of Marine, Earth and Atmospheric Sciences
>North Carolina State University
>Email: address@hidden
>Phone: (919) 515-1447
>
>Unidata Support wrote:
>> Michale,
>> 
>> Yes, they are compiletime defined in 
>> $GEMPAK/source/programs/gd/gdvint/viprm.prm.
>> 
>> You should be able to redefine and compile just gdvint.
>> 
>> Steve Chiswell
>> 
>> 
>> 
>>>From: Michael Brennan <address@hidden>
>>>Organization: UCAR/Unidata
>>>Keywords: 200409282143.i8SLhPHP017603
>> 
>> 
>>>Hi Steve,
>>>
>>>Thanks, your suggestions were very helpful.
>>>
>>>The conversion to GEMPAK from grib was pretty straightforward from that 
>>>point.
>>>
>>>Now I'm trying to use gdvint to try and interpolate from this WRF eta 
>>>level data to pressure level data.  I've had some partial success, but 
>>>I'm running into some hard limits for:
>>>
>>>1) the number of eta levels that I can process
>>>2) the number of pressure levels I can write output to
>>>3) the number of parameters I can interpolate
>>>
>>>Are these hardwired somewhere in the source code for gdvint?
>>>
>>>Thanks,
>>>Mike
>>>
>>>Michael J. Brennan
>>>Graduate Research Assistant
>>>Department of Marine, Earth and Atmospheric Sciences
>>>North Carolina State University
>>>Email: address@hidden
>>>Phone: (919) 515-1447
>>>
>>>Unidata Support wrote:
>>>
>>>>Michael,
>>>>
>>>>The vcrdgrib1.tbl provides for the ETA level and Layer coordinate.
>>>>
>>>>Your output below: WMO parameter table to wmogrib200.tbl
>>>>shows that your grib is using parameters in table 200.
>>>>This is not a table that exists in the current NCEP grib
>>>>documentation so you will have to create that table for parameters 1-127, 
>>>>as well as the local (xxxxgrib200.tbl) file  for parameters 128-255.
>>>>What model center is being used? The model center number will have to
>>>>be defined in cntrgrib1.tbl (if you aren't using an existing number),
>>>>and that will define the string used for your local table name.
>>>>
>>>>Without those tables, you wont have any defined parameters to decode,
>>>>so you won't have anything in your output file..
>>>>
>>>>Steve Chiswell
>>>>
>>>>
>>>>
>>>>>From: Michael Brennan <address@hidden>
>>>>>Organization: UCAR/Unidata
>>>>>Keywords: 200409271712.i8RHC5nJ025022
>>>>
>>>>
>>>>>Steve,
>>>>>
>>>>>I've been able to get WRF model output in grib format by using an grib 
>>>>>output module from Todd Hutchinson at WSI.  I'd like to try and convert 
>>>>>the grib files to GEMPAK.  However, the WRF output is in the WRF's "Eta" 
>>>>>vertical coordinate, and I think that trips up dcgrib2.  I don't get an 
>>>>>error message when I run dcgrib2, but the GEMPAK file doesn't get created.
>>>>>
>>>>>[mjbrenna@compute-2-35 mjbrenna]$ dcgrib2 wrf_test.gem < wrfout.grb
>>>>>Changing center table to cntrgrib1.tbl
>>>>>Changing vertical coord table to vcrdgrib1.tbl
>>>>>Changing WMO parameter table to wmogrib200.tbl
>>>>>
>>>>>Do you have any thoughts or suggestions?
>>>>>
>>>>>Thanks,
>>>>>Mike
>>>>>-- 
>>>>>Michael J. Brennan
>>>>>Graduate Research Assistant
>>>>>Department of Marine, Earth and Atmospheric Sciences
>>>>>North Carolina State University
>>>>>Email: address@hidden
>>>>>Phone: (919) 515-1447
>>>>>
>>>>
>>>>--
>>>>***************************************************************************
> *
>>>
>>>
>>>>Unidata User Support                                    UCAR Unidata Progra
> m
>>>
>>>
>>>>(303)497-8643                                                  P.O. Box 300
> 0
>>>
>>>
>>>>address@hidden                                   Boulder, CO 8030
> 7
>>>
>>>
>>>>---------------------------------------------------------------------------
> -
>>>
>>>
>>>>Unidata WWW Service              http://my.unidata.ucar.edu/content/support
>  
>>>
>>>
>>>>---------------------------------------------------------------------------
> -
>>>
>>>
>>>>NOTE: All email exchanges with Unidata User Support are recorded in the
>>>>Unidata inquiry tracking system and then made publicly available
>>>>through the web.  If you do not want to have your interactions made
>>>>available in this way, you must let us know in each email you send to us.
>>>
>> --
>> ****************************************************************************
>> Unidata User Support                                    UCAR Unidata Program
>> (303)497-8643                                                  P.O. Box 3000
>> address@hidden                                   Boulder, CO 80307
>> ----------------------------------------------------------------------------
>> Unidata WWW Service              http://my.unidata.ucar.edu/content/support 
>> ----------------------------------------------------------------------------
>> NOTE: All email exchanges with Unidata User Support are recorded in the
>> Unidata inquiry tracking system and then made publicly available
>> through the web.  If you do not want to have your interactions made
>> available in this way, you must let us know in each email you send to us.
>
--
NOTE: All email exchanges with Unidata User Support are recorded in the
Unidata inquiry tracking system and then made publicly available
through the web.  If you do not want to have your interactions made
available in this way, you must let us know in each email you send to us.