#!/bin/sh -

#HH=`grep ' H= ' "$3" | awk '{print $2}'`
#H=`echo $HH | awk '{print $2}'`
#pV=$1
#eta=$2
#D=`echo $H $pV | /cygdrive/z/fitTM/h2d | awk '{print $3}'`
#G=0.15
#A=$pV*0.393
#echo $3  H=$H D=$D
#echo ------------wavelengths -----------------------------------------------------
#echo 3.6, 4.5, 5.8, 8.0, 16.0, 22.0  Date| 
#     awk '{printf("%8.3f %8.3f %8.3f %8.3f %8.3f %8.3f %s\n",$1,$2,$3,$4,$5,$6,$7)}'
#echo ' '
#echo ------------fluxes '(mJy)'-----------------------------------------------------

if [ $# = 0 ]; then
    echo "usage: $program file ..."
    exit 1
fi   

if [ $# = 1 ]; then
exec awk ' BEGIN { FS=","; OFMT="%14.7f"}
/^ *[0-9]/ && $2 >0  {
		      lambdaH=$7/180.*3.1415; betaH=$8/180*3.1415;
		      clH=cos(lambdaH); slH=sin(lambdaH);
		      cbH=cos(betaH); sbH=sin(betaH);
		      xH=clH*cbH*$9; yH=slH*cbH*$9; zH=sbH*$9;	    
		      
		      lambdaO=$14/180.*3.1415; betaO=$15/180*3.1415;
		      clO=cos(lambdaO); slO=sin(lambdaO);
		      cbO=cos(betaO); sbO=sin(betaO);
		      xO=clO*cbO*$11; yO=slO*cbO*$11; zO=sbO*$11
                      print $2,xH,yH,zH,xO,yO,zO;
                      next }
' "$1"
fi

if [ $# = 3 ]; then
exec awk -v jdi=$2 -v jde=$3 ' BEGIN { FS=","; OFMT="%14.7f"; jdini=$2}
/^ *[0-9]/ && $2 >= jdi && $2 <= jde {
		      lambdaH=$7/180.*3.1415; betaH=$8/180*3.1415;
		      clH=cos(lambdaH); slH=sin(lambdaH);
		      cbH=cos(betaH); sbH=sin(betaH);
		      xH=clH*cbH*$9; yH=slH*cbH*$9; zH=sbH*$9;	    
		      
		      lambdaO=$14/180.*3.1415; betaO=$15/180*3.1415;
		      clO=cos(lambdaO); slO=sin(lambdaO);
		      cbO=cos(betaO); sbO=sin(betaO);
		      xO=clO*cbO*$11; yO=slO*cbO*$11; zO=sbO*$11
                      print $2,xH,yH,zH,xO,yO,zO;
                      next }
' "$1"
fi
