#!/bin/csh 

# generate a file of XYZ differences for GPPS
# example output
# 2  300   30   30000.1193    .1221   0   0  ARO0  CCC0   311

# get X, in columns 5 and 7, after chopping 2-part line name 
cat patr.x.xyz.bsl | colrm 120 | sed 's/_/ /' | awk '{OFS="\t";print 2,"N"$1,"N"$2,$5,3*$7,0,0,$1,$2,$3"Kclen"}' | sed -f codes.sed  >! tmp.x

# get Y, in columns 9 and 11, after chopping 2-part line name 
cat patr.x.xyz.bsl | colrm 120 | sed 's/_/ /' | awk '{OFS="\t";print 3,"N"$1,"N"$2,$9,3*$11,0,0,$1,$2,$3"Kclen"}' | sed -f codes.sed  >! tmp.y

# get Z, in columns 13 and 15, after chopping 2-part line name 
cat patr.x.xyz.bsl | colrm 120 | sed 's/_/ /' | awk '{OFS="\t";print 4,"N"$1,"N"$2,$13,3*$15,0,0,$1,$2,$3"Kclen"}' | sed -f codes.sed  >! tmp.z

cat tmp.x tmp.y tmp.z | sort +7 +8 +9 +0 | expand -4 | sed 's/1991.//' >! bsl2gpps.out

head bsl2gpps.out
echo created bsl2gpps.out