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
