#! /bin/csh -f

set start_time = `date +%s`
source ~/.cshrc
cd /network/rit/lab/mesonet/mesonet/radar
setenv RAD /network/daes/data7/nsat/radar/noaaport

set mi = `date -u +%M`
set mm = `date -u +%m`
set dd = `date -u +%d`
set hh = `date -u +%H`
set yy = `date -u +%y`
echo The Current Date Is $yy$mm$dd$hh$mi
set date0 = $yy$mm$dd$hh$mi

set mi = `date -u -d -5minutes +%M`; set mm = `date -u -d -5minutes +%m`; set dd = `date -u -d -5minutes +%d`
set hh = `date -u -d -5minutes +%H`; set yy = `date -u -d -5minutes +%y`; set date1 = $yy$mm$dd$hh$mi

set mi = `date -u -d -10minutes +%M`; set mm = `date -u -d -10minutes +%m`; set dd = `date -u -d -10minutes +%d`
set hh = `date -u -d -10minutes +%H`; set yy = `date -u -d -10minutes +%y`; set date2 = $yy$mm$dd$hh$mi

set mi = `date -u -d -15minutes +%M`; set mm = `date -u -d -15minutes +%m`; set dd = `date -u -d -15minutes +%d`
set hh = `date -u -d -15minutes +%H`; set yy = `date -u -d -15minutes +%y`; set date3 = $yy$mm$dd$hh$mi

set mi = `date -u -d -20minutes +%M`; set mm = `date -u -d -20minutes +%m`; set dd = `date -u -d -20minutes +%d`
set hh = `date -u -d -20minutes +%H`; set yy = `date -u -d -20minutes +%y`; set date4 = $yy$mm$dd$hh$mi

set mi = `date -u -d -25minutes +%M`; set mm = `date -u -d -25minutes +%m`; set dd = `date -u -d -25minutes +%d`
set hh = `date -u -d -25minutes +%H`; set yy = `date -u -d -25minutes +%y`; set date5 = $yy$mm$dd$hh$mi

set mi = `date -u -d -30minutes +%M`; set mm = `date -u -d -30minutes +%m`; set dd = `date -u -d -30minutes +%d`
set hh = `date -u -d -30minutes +%H`; set yy = `date -u -d -30minutes +%y`; set date6 = $yy$mm$dd$hh$mi

set mi = `date -u -d -35minutes +%M`; set mm = `date -u -d -35minutes +%m`; set dd = `date -u -d -35minutes +%d`
set hh = `date -u -d -35minutes +%H`; set yy = `date -u -d -35minutes +%y`; set date7 = $yy$mm$dd$hh$mi

set mi = `date -u -d -40minutes +%M`; set mm = `date -u -d -40minutes +%m`; set dd = `date -u -d -40minutes +%d`
set hh = `date -u -d -40minutes +%H`; set yy = `date -u -d -40minutes +%y`; set date8 = $yy$mm$dd$hh$mi

set mi = `date -u -d -45minutes +%M`; set mm = `date -u -d -45minutes +%m`; set dd = `date -u -d -45minutes +%d`
set hh = `date -u -d -45minutes +%H`; set yy = `date -u -d -45minutes +%y`; set date9 = $yy$mm$dd$hh$mi

set mi = `date -u -d -50minutes +%M`; set mm = `date -u -d -50minutes +%m`; set dd = `date -u -d -50minutes +%d`
set hh = `date -u -d -50minutes +%H`; set yy = `date -u -d -50minutes +%y`; set date10 = $yy$mm$dd$hh$mi

set mi = `date -u -d -55minutes +%M`; set mm = `date -u -d -55minutes +%m`; set dd = `date -u -d -55minutes +%d`
set hh = `date -u -d -55minutes +%H`; set yy = `date -u -d -55minutes +%y`; set date11 = $yy$mm$dd$hh$mi

set mi = `date -u -d -60minutes +%M`; set mm = `date -u -d -60minutes +%m`; set dd = `date -u -d -60minutes +%d`
set hh = `date -u -d -60minutes +%H`; set yy = `date -u -d -60minutes +%y`; set date12 = $yy$mm$dd$hh$mi

set mi = `date -u -d -65minutes +%M`; set mm = `date -u -d -65minutes +%m`; set dd = `date -u -d -65minutes +%d`
set hh = `date -u -d -65minutes +%H`; set yy = `date -u -d -65minutes +%y`; set date13 = $yy$mm$dd$hh$mi

set mi = `date -u -d -70minutes +%M`; set mm = `date -u -d -70minutes +%m`; set dd = `date -u -d -70minutes +%d`
set hh = `date -u -d -70minutes +%H`; set yy = `date -u -d -70minutes +%y`; set date14 = $yy$mm$dd$hh$mi

set mi = `date -u -d -75minutes +%M`; set mm = `date -u -d -75minutes +%m`; set dd = `date -u -d -75minutes +%d`
set hh = `date -u -d -75minutes +%H`; set yy = `date -u -d -75minutes +%y`; set date15 = $yy$mm$dd$hh$mi

set mi = `date -u -d -80minutes +%M`; set mm = `date -u -d -80minutes +%m`; set dd = `date -u -d -80minutes +%d`
set hh = `date -u -d -80minutes +%H`; set yy = `date -u -d -70minutes +%y`; set date16 = $yy$mm$dd$hh$mi

set mi = `date -u -d -85minutes +%M`; set mm = `date -u -d -85minutes +%m`; set dd = `date -u -d -85minutes +%d`
set hh = `date -u -d -85minutes +%H`; set yy = `date -u -d -85minutes +%y`; set date17 = $yy$mm$dd$hh$mi

set mi = `date -u -d -90minutes +%M`; set mm = `date -u -d -90minutes +%m`; set dd = `date -u -d -90minutes +%d`
set hh = `date -u -d -90minutes +%H`; set yy = `date -u -d -90minutes +%y`; set date18 = $yy$mm$dd$hh$mi

set mi = `date -u -d -95minutes +%M`; set mm = `date -u -d -95minutes +%m`; set dd = `date -u -d -95minutes +%d`
set hh = `date -u -d -95minutes +%H`; set yy = `date -u -d -95minutes +%y`; set date19 = $yy$mm$dd$hh$mi

set mi = `date -u -d -100minutes +%M`; set mm = `date -u -d -100minutes +%m`; set dd = `date -u -d -100minutes +%d`
set hh = `date -u -d -100minutes +%H`; set yy = `date -u -d -100minutes +%y`; set date20 = $yy$mm$dd$hh$mi

#mv radar.gem $date1.gem
#mv meso.gem meso$date1.gem
rm radar.gem

gdradr<< EOF


 GRDAREA  = 39.8;-80.5;46;-70.5
 PROJ     = LCC/30.5;-75.5;55.5
 KXKY     = 790;680
 GDPFUN   = DHR
 GDFILE   = radar.gem
 RADTIM   = current
 RADDUR   = 60
 RADFRQ   = 0
 CPYFIL   =  
 STNFIL   = nexrad.tbl
 MAXGRD   = 2000
 RADMODE  = PC
 NDVAL    =  

r

e

EOF

gdcsv<< EOF
 GDATTIM  = last
 GDFILE   = radar.gem
 GLEVEL   = 0
 GAREA    = grid
 PROJ     = LCC/30.5;-75.5;55.5
 GVCORD   = none
 GFUNC    = DHR
 SCALE    = 0
 OUTPUT   = f/radar.csv
r

e

EOF

rm meso.gem

sfcfil<< EOF

 SFOUTF   = meso.gem
 SFPRMF   = TMPF;DWPF;SMPH;DRCT;PRES;WSYM;GUST
 STNFIL   = mesonet.tbl
 SHIPFL   = no
 TIMSTN   = 1000/130
 SFFSRC   = AIRW|TEXT

r

e

EOF

sfedit<< EOF
 SFEFIL   = surface.dat
 SFFILE   = meso.gem
r

e

EOF

gdplot<< EOF
 GDFILE   = radar.gem
 GDATTIM  = last
 GLEVEL   = 0
 GVCORD   = none
 PANEL    = 0
 SKIP     = 0
 SCALE    = 0
 GFUNC    = DHR
 CTYPE    = f
 CONTUR   = 3/1
 CINT     = 1;2;5;10;20;30
 LINE     = 31/-2/2
 FINT     = 9;12;15;18;21;24;27;30;33;36;39;42;45;48;51;54;57;60;63;66;69;72;75
 FLINE    = 0;30-19;10;17-11;7;6;1
 HILO     =  
 HLSYM    =  
 CLRBAR   = 31/v///3;.08/1
 CLRBAR   = 31/h///.98/-1
 GVECT    =  
 WIND     = BM1
 REFVEC   =  
 TITLE    = 31/-2/^ Base Radar Composite & Current Observations
 TEXT     = 1.25/22/HW
 CLEAR    = n
 GAREA    = grid
 PROJ     = LCC/30.5;-75.5;55.5
 MAP      = 8//2+31/-2/2+31/-2/4
 MSCALE   = 0
 LATLON   =  
 DEVICE   = gif|radar.gif|1500;1500|C
 STNPLT   =  
 \$MAPFIL  = hicnus.nws + hiisus.nws +hipowo.cia

r

e

EOF

sfmap<< EOF
 AREA     = 39.8;-80.5;46;-70.5
 GAREA    =39.8;-80.5;46;-70.5
 SATFIL   =  
 RADFIL   = $RAD/NEXRCOMP/1km/n0r/n0r_151223_1835
 IMCBAR   =  
 SFPARM   = skyc:.75;tmpf:1.25:3;wsym:.75;rmsl:1.25:3;ptnd;dwpf:1.25:3;stid:1.25:3;;vsby*10;brbk::3;mark:4
 SFPARM   = skyc:.75;tmpf:1.2:3;wsym:1.1:2;smsl:1.2:3;ptnd;dwpf:1.2:3;;;vsby*10;brbk::3;mark:4
 DATTIM   = last
 SFFILE   = meso.gem
 COLORS   = 26;2;(19-90-10/17;6;8;17;3;3;1;22;2);25;20;22;7;24;6
 COLORS   = 26;2;31;31;20;22;7;31;31
# Colors   = 31
 MSCALE   = 0
 LATLON   =  
 FILTER   = no
 TEXT     = 1/22/HW
 LUTFIL   =  
 STNPLT   =  
 CLRBAR   =  
 LSTPRM   =  
 CLEAR    = n
 MAP      = 0
 TITLE    = 
r

e

EOF

gpend

#ncl radar.ncl
#convert -shave 150x150 radar.png temp.png
#mv temp.png radar.png

scp  *.gif *.png nbassill@operations.nysmesonet.org:/home/nbassill/public_html/plots/live

set end_time   = `date  +%s`
@ length_time  = $end_time - $start_time
set mm = `date -u -d -${length_time}seconds +%m`
set dd = `date -u -d -${length_time}seconds +%d`
set yy = `date -u -d -${length_time}seconds +%y`
scp radar.gif nbassill@operations.nysmesonet.org:/home/nbassill/public_html/archive/20$yy/$mm/$dd/radar/20$date0-radar.gif
scp *csh *ncl *tbl nbassill@operations.nysmesonet.org:/home/nbassill/public_html/archive/20$yy/$mm/$dd/radar/scripts

mv radar.gif gifs/20$date0-radar.gif
mv radar.png gifs/20$date0-radar.png

ln -s gifs/20$date0-radar.gif 20$date0-radar0.gif
ln -s gifs/20$date0-radar.gif 20$date0-radar1.gif
ln -s gifs/20$date0-radar.gif 20$date0-radar2.gif
ln -s gifs/20$date0-radar.gif 20$date0-radar3.gif
ln -s gifs/20$date0-radar.gif 20$date0-radar4.gif
ln -s gifs/20$date1-radar.gif .
ln -s gifs/20$date2-radar.gif .
ln -s gifs/20$date3-radar.gif .
ln -s gifs/20$date4-radar.gif .
ln -s gifs/20$date5-radar.gif .
ln -s gifs/20$date6-radar.gif .
ln -s gifs/20$date7-radar.gif .
ln -s gifs/20$date8-radar.gif .
ln -s gifs/20$date9-radar.gif .
ln -s gifs/20$date10-radar.gif .
ln -s gifs/20$date11-radar.gif .
ln -s gifs/20$date10-radar.gif .
ln -s gifs/20$date11-radar.gif .
ln -s gifs/20$date12-radar.gif .
ln -s gifs/20$date13-radar.gif .
ln -s gifs/20$date14-radar.gif .
ln -s gifs/20$date15-radar.gif .
ln -s gifs/20$date16-radar.gif .
ln -s gifs/20$date17-radar.gif .
ln -s gifs/20$date18-radar.gif .
ln -s gifs/20$date19-radar.gif .
ln -s gifs/20$date20-radar.gif .

set d = 20 # delay length

convert -delay 20 -loop 0 *gif animatedradar.gif 
#rm /home/nb965428/public_html/live/animatedradar.gif
scp animatedradar.gif nbassill@operations.nysmesonet.org:/home/nbassill/public_html/plots/live
#cp radar.gem ../analysis

rm *.png *.gif

set end_time   = `date  +%s`
@ length_time  = $end_time - $start_time
echo "It took this long to run radar.csh: $length_time seconds"