#!/bin/csh

# transform the coast line data
sed 's/nan nan/371 371/' < coastll.dat | trans_pole 1 0 285.6 50.1 276 54 > coastll_trans.dat
awk '{print $6, $7}' coastll_trans.dat | grep -v 0.000000 |  grep -v nan > coast.por.dat

# transform the los data
trans_pole 1 1 285.6 50.1 276 54 < los > alos_asc.trans.llv
awk '{print $6,$7,$3}' alos_asc.trans.llv > alos_asc.por.yxv

trans_pole 1 1 285.6 50.1 276 54 < u > alos_asc.trans.llu
awk '{print $6,$7,$3}' alos_asc.trans.llu > alos_asc.por.yxu

trans_pole 1 1 285.6 50.1 276 54 < std > alos_asc.trans.lls
awk '{print $6,$7,$3}' alos_asc.trans.lls > alos_asc.por.yxs

paste e n | awk '{print $1,$2,$3,$6,0,0,0}' > alos_asc.llen
trans_pole 1 2 285.6 50.1 276 54 < alos_asc.llen > alos_asc.trans.llen

awk '{print $15, $16, $10}' alos_asc.trans.llen > alos_asc.por.yxyy
awk '{print $15, $16, $11}' alos_asc.trans.llen > alos_asc.por.yxxx

#
awk '{print $3}' alos_asc.por.yxs > tmps
awk '{print $3}' alos_asc.por.yxv > tmpv
awk '{print $3}' alos_asc.por.yxu > tmpu

# shift the X Y coordinates
awk '{print $1+850,$2+400,$3}' alos_asc.por.yxxx > xx
awk '{print $3}' alos_asc.por.yxyy > yy
paste xx yy tmpu tmpv tmps > alos_asc.por.dat

# cleanup
rm alos_asc.trans.*
rm tmps tmpv tmpu xx yy
rm alos_asc.llen
rm coastll_trans.dat