#/bin/csh -f
# Script to plot track results
set f1 = TRAK356_19.NEU.crbt.LC
set f2 = TRAK356_19.NEU.lows.LC
set f3 = TRAK356_19.NEU.pomm.LC
#
set ps = TRAK356_19_eq.ps
echo "Creating $ps"

#
set f1r = `head -3 $f1 | tail -n -1 | awk '{print $7}'`
tail -n +3 $f1 | awk -v f1r=$f1r '{print $18,($7-f1r)*1000-100,$8*1000}'  >! t.1

set f2r = `head -3 $f2 | tail -n -1 | awk '{print $7}'`
tail -n +3 $f2 | awk -v f2r=$f2r '{print $18,($7-f2r)*1000+50,$8*1000}' >! t.2

set f3r = `head -3 $f3 | tail -n -1 | awk '{print $7}'`
tail -n +3 $f3 | awk -v f3r=$f3r '{print $18,($7-f3r)*1000+150,$8*1000}'  >! t.3

psxy t.1 -R900/1200/-300/250 -JX5/3 -Ba100f20/a100f50:"North (mm)":/SWen -Y6 -K -W1/255/0/0 -P >! $ps
awk '{if (NR-int(NR/10)*10 == 0 ) {print $0}}' t.1 | psxy  -R -JX -Sp -Ey/1/200/200/200 -G255/0/0 -K -O >> $ps
psxy t.2 -R -JX -K -O -W1/0/0/0 >> $ps
awk '{if (NR-int(NR/10)*10 == 0 ) {print $0}}' t.2 | psxy  -R -JX -Sp -Ey/1/200/200/200 -G128/0/0 -K -O >> $ps
psxy t.3 -R -JX -K -O -W1/0/0/255 >> $ps
awk '{if (NR-int(NR/10)*10 == 0 ) {print $0}}' t.3 | psxy  -R -JX -Sp -Ey/1/200/200/200 -G0/0/255 -K -O >> $ps
pstext -R -JX -K -O <<! >> $ps
910 -200 10 0 0 ML CRBT
910  -50 10 0 0 ML LOWS
910  200 10 0 0 ML POMM
!


# East
set f1r = `head -3 $f1 | tail -n -1 | awk '{print $9}'`
tail -n +3 $f1 | awk -v f1r=$f1r '{print $18,($9-f1r)*1000-100,$10*1000}'  >! t.1

set f2r = `head -3 $f2 | tail -n -1 | awk '{print $9}'`
tail -n +3 $f2 | awk -v f2r=$f2r '{print $18,($9-f2r)*1000+50,$10*1000}' >! t.2

set f3r = `head -3 $f3 | tail -n -1 | awk '{print $9}'`
tail -n +3 $f3 | awk -v f3r=$f3r '{print $18,($9-f3r)*1000+150,$10*1000}'  >! t.3

psxy t.1 -R900/1200/-300/250 -JX5/3 -Ba100f20:"Seconds from 19 00 Dec 21, 2003":/a100f50:"East (mm)":/SWen -Y-3.5 -K -O -W1/255/0/0 -P >> $ps
awk '{if (NR-int(NR/10)*10 == 0 ) {print $0}}' t.1 | psxy  -R -JX -Sp -Ey/1/200/200/200 -G255/0/0 -K -O >> $ps
psxy t.2 -R -JX -K -O -W1/0/0/0 >> $ps
awk '{if (NR-int(NR/10)*10 == 0 ) {print $0}}' t.2 | psxy  -R -JX -Sp -Ey/1/200/200/200 -G128/0/0 -K -O >> $ps
psxy t.3 -R -JX -K -O -W1/0/0/255 >> $ps
awk '{if (NR-int(NR/10)*10 == 0 ) {print $0}}' t.3 | psxy  -R -JX -Sp -Ey/1/200/200/200 -G0/0/255 -K -O >> $ps
pstext -R -JX -O <<! >> $ps
910 -200 10 0 0 ML CRBT
910    0 10 0 0 ML LOWS
910  200 10 0 0 ML POMM
!


#-----------------------------------------------------------------------------------------------------
# Now plot the whole 1-hour of data
#
set ps = TRAK356_19_all.ps
echo "Creating $ps"

#
set f1r = `head -3 $f1 | tail -n -1 | awk '{print $7}'`
tail -n +3 $f1 | awk -v f1r=$f1r '{print $18,($7-f1r)*1000-100,$8*1000}'  >! t.1

set f2r = `head -3 $f2 | tail -n -1 | awk '{print $7}'`
tail -n +3 $f2 | awk -v f2r=$f2r '{print $18,($7-f2r)*1000+50,$8*1000}' >! t.2

set f3r = `head -3 $f3 | tail -n -1 | awk '{print $7}'`
tail -n +3 $f3 | awk -v f3r=$f3r '{print $18,($7-f3r)*1000+150,$8*1000}'  >! t.3

psxy t.1 -R0/3600/-300/250 -JX5/3 -Ba600f300/a100f50:"North (mm)":/SWen -Y6 -K -W1/255/0/0 -P >! $ps
awk '{if (NR-int(NR/30)*30 == 0 ) {print $0}}' t.1 | psxy  -R -JX -Sp -Ey/1/200/200/200 -G255/0/0 -K -O >> $ps
psxy t.2 -R -JX -K -O -W1/0/0/0 >> $ps
awk '{if (NR-int(NR/30)*30 == 0 ) {print $0}}' t.2 | psxy  -R -JX -Sp -Ey/1/200/200/200 -G128/0/0 -K -O >> $ps
psxy t.3 -R -JX -K -O -W1/0/0/255 >> $ps
awk '{if (NR-int(NR/30)*30 == 0 ) {print $0}}' t.3 | psxy  -R -JX -Sp -Ey/1/200/200/200 -G0/0/255 -K -O >> $ps
pstext -R -JX -K -O <<! >> $ps
50 -200 10 0 0 ML CRBT
50  -50 10 0 0 ML LOWS
50  200 10 0 0 ML POMM
!


# East
set f1r = `head -3 $f1 | tail -n -1 | awk '{print $9}'`
tail -n +3 $f1 | awk -v f1r=$f1r '{print $18,($9-f1r)*1000-100,$10*1000}'  >! t.1

set f2r = `head -3 $f2 | tail -n -1 | awk '{print $9}'`
tail -n +3 $f2 | awk -v f2r=$f2r '{print $18,($9-f2r)*1000+50,$10*1000}' >! t.2

set f3r = `head -3 $f3 | tail -n -1 | awk '{print $9}'`
tail -n +3 $f3 | awk -v f3r=$f3r '{print $18,($9-f3r)*1000+150,$10*1000}'  >! t.3

psxy t.1 -R0/3600/-300/250 -JX5/3 -Ba600f300:"Seconds from 19 00 Dec 21, 2003":/a100f50:"East (mm)":/SWen -Y-3.5 -K -O -W1/255/0/0 -P >> $ps
awk '{if (NR-int(NR/30)*30 == 0 ) {print $0}}' t.1 | psxy  -R -JX -Sp -Ey/1/200/200/200 -G255/0/0 -K -O >> $ps
psxy t.2 -R -JX -K -O -W1/0/0/0 >> $ps
awk '{if (NR-int(NR/30)*30 == 0 ) {print $0}}' t.2 | psxy  -R -JX -Sp -Ey/1/200/200/200 -G128/0/0 -K -O >> $ps
psxy t.3 -R -JX -K -O -W1/0/0/255 >> $ps
awk '{if (NR-int(NR/30)*30 == 0 ) {print $0}}' t.3 | psxy  -R -JX -Sp -Ey/1/200/200/200 -G0/0/255 -K -O >> $ps
pstext -R -JX -O <<! >> $ps
50 -200 10 0 0 ML CRBT
50    0 10 0 0 ML LOWS
50  200 10 0 0 ML POMM
!

\rm t.?