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
