NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.

To learn about what's going on, see About the Archive Site.

[gembud] wrf2gem and problem with gempak file

hi,
I converted a WRF output file into a gempak file choosing only 2 m-temperature variable. Now I am having problems with creating a grid file. Can anyone out there who has had experience with wrf2gem and gempak plotting provide some plotting scripts or examples of how to / what to use for the gdcfil, gdedit and gdplot2 etc.? Once you read the u,v, temperature and pressure and write them out to a gempak file (by making changes in the file convert.nl ), could you please explain step by step what to do, how to do to plot say temperature advection at 850 mb pressure level for a given time. 

I have attached a script which supposedly makes a grid file.
After it is ran, I have got the following messages
*** Generating Grid file 043005_0000_grd **
GDOUTF    Output grid file                  043005_0000_grd
PROJ      Map projection/angles/margins|dr  mer
GRDAREA   Area covered by grid              25.0;-102.0;36.0;-88.0
KXKY      Number of grid points in x;y      50;50
MAXGRD    Maximum number of grids           10000
CPYFIL Grid file whose navigation is to ANLYSS Grid analysis block Parameters requested: GDOUTF,PROJ,GRDAREA,KXKY,MAXGRD,CPYFIL,ANLYSS. GEMPAK-GDCFIL> GEMPAK-GDCFIL> GEMPAK-GDCFIL> GEMPAK-GDCFIL> GEMPAK-GDCFIL> GEMPAK-GDCFIL> GEMPAK-GDCFIL> GEMPAK-GDCFIL>Creating process: gplt for queue 18055174 

GDCFIL PARAMETERS:
New grid file: 043005_0000_grd Grid projection: mer Grid area: 25.0;-102.0;36.0;-88.0 Grid dimension: 50;50 Grid anlysis: Maximum number of grids that can be added to file 10000 
grid file has just been prepared
fill the grid file with the data from gempak file
GDEFIL    Grid edit file                    wrf_k05_sfc.gem
GDFILE    Grid file                         043005_0000_grd
GPACK     Packing type/number of bits       NONE
Parameters requested: GDEFIL,GDFILE,GPACK.
GEMPAK-GDEDIT> GEMPAK-GDEDIT> GEMPAK-GDEDIT> GEMPAK-GDEDIT> [GDEDIT -3] Valid grid identifier not found. 
grid file has just been filled*

thanks in advance
Haldun


#!/bin/csh

 set yr=05 
 set mon=04
 set day=30
 set hr=00 
 set lat1=25.0
 set lon1=-102.0
 set lat2=36.0
 set lon2=-88.0
 set gemfile=wrf_k05_sfc.gem

 echo " " 
 echo "** Generating Grid file $mon${day}${yr}_${hr}00_grd **"
 echo "GDOUTF   = $mon${day}${yr}_${hr}00_grd" >> MKGRD
 echo "PROJ     = mer" >> MKGRD
 echo "GRDAREA  = $lat1;$lon1;$lat2;$lon2" >> MKGRD
 echo "KXKY     = 50;50" >> MKGRD
 echo "MAXGRD   = 10000" >> MKGRD
 echo "CPYFIL   = " >> MKGRD
 echo "ANLYSS   = " >> MKGRD
 echo "r" >> MKGRD
 echo " " >> MKGRD
 echo " " >> MKGRD
 echo "e" >> MKGRD
 $GEMEXE/gdcfil < MKGRD
 rm -f MKGRD
 echo "grid file has just been prepared"

 echo " fill the grid file with the data from gempak file "
 echo "GDEFIL = $gemfile" >> GRDFIL
 echo "GDFILE = $mon${day}${yr}_${hr}00_grd" >> GRDFIL
 echo "GPACK = NONE" >> GRDFIL
 echo "r" >> GRDFIL
 echo " " >> GRDFIL
 echo " " >> GRDFIL
 echo "e" >> GRDFIL
 $GEMEXE/gdedit < GRDFIL
 rm -f GRDFIL
 echo "grid file has just been filled"

 $GEMEXE/gpend
 exit
  • 2007 messages navigation, sorted by:
    1. Thread
    2. Subject
    3. Author
    4. Date
    5. ↑ Table Of Contents
  • Search the gembud archives: