[gembud] Model help

• To: <gembud@xxxxxxxxxxxxxxxx>
• Subject: [gembud] Model help
• From: "Jeff Lake - Admin" <admin@xxxxxxxxxxxxxxxxxxxx>
• Date: Mon, 2 Aug 2010 09:30:08 -0400

greetings all..
attached is my simple GFS precip/pmsl map script..
not really finding any examples to create maps from the models,
I wrote this..
Is there a better way?
Or is this the basic way I will need to do it for whatever I want out of the 
models?

here is my end result after my php scripting is done
http://images.michiganwxsystem.net/grads/models/gfs/

thanks..
Jeff Lake K8JSL
MichiganWxSystem.com
WeatherMichigan.net
TheWeatherCenter.net << Watches/Warnings for your inbox!!

#!/bin/sh


cd ~/nawips/imgcreate/gfs

. ~/nawips/Gemenviron.profile

rm -fr p06m
mkdir p06m
cd p06m
cp ~/nawips/imgcreate/gfs/precp.tbl coltbl.xwp

gdinfo << EOFE
GDFILE = gfs211
LSTALL = YES
OUTPUT = f/211.txt
GDATTIM = all
GLEVEL = all
GVCORD = all
GFUNC = all

r

e

EOFE

cp 211.txt /var/www/html/images/grads/gfs/



DATE[0]=F006
DATE[1]=F012
DATE[2]=F018
DATE[3]=F024
DATE[4]=F030
DATE[5]=F036
DATE[6]=F042
DATE[7]=F048
DATE[8]=F054
DATE[9]=F060
DATE[10]=F066
DATE[11]=F072
DATE[12]=F078
DATE[13]=F084
DATE[14]=F090
DATE[15]=F096
DATE[16]=F102
DATE[17]=F108
DATE[18]=F114
DATE[19]=F120
DATE[20]=F126
DATE[21]=F132
DATE[22]=F138
DATE[23]=F144
DATE[24]=F150

for ((i=0; i < ${#DATE[*]}; i++))
   do






gpend

gdplot3 << EOF2
device   = gif|gfsprecip_${DATE[$i]}.gif|800;600
gdfile   = gfs211
\$mapfil =
proj     = lcc/20;-95;60/0;0;0;0
garea    = uslcc
gdpfun  =  p06i
map      = 32//1
clear    = y
gdattim  =${DATE[$i]}
text     = 1/21
TITLE    = 0
GLEVEL   = 0
GVCORD   = none
TYPE    = F
PANEL    = 0
SKIP     = 0
SCALE    = 0
HILO     =
LATLON   = 0
STNPLT   =
WOU     =
WARN =
SATFIL   =
RADFIL   =
LUTFIL   =
STREAM   =
CLRBAR = 0
POSN     = 0
MARKER   = 0
GRDLBL   = 0
FILTER   = NO
CONTUR = 3/3
FINT =.01;.1;.25;.5;.75;1;1.25;1.5;1.75;2;2.5;3;4;5;6;7;8;9
FLINE =0;1;2;3;4;5;6;7;8;9;10;11;12;13;14;15;16;17;18;19;20;21

r

gdpfun = pmsl
cint = 2/880/1100
line = 30//1
type = c
text = .83/1
clear = n

r

e

EOF2

gpend

/usr/bin/convert gfsprecip_${DATE[$i]}.gif 
/var/www/html/images/grads/gfs/us_precipmslp_${DATE[$i]}.png
php ~/decoders/gfs/gfsprecip.php run=211 image=us_precipmslp_${DATE[$i]}


done