#!/bin/sh

if [ $# -lt 3 ]
then
  echo "Usage: $0 [sso_freeze ocsol output file] [# of regions] [outfile name]"
  exit
fi

ocsol=$1
nreg=$2
outf=$3

#...check that input OC solution file exists
if [ ! -f $ocsol ]
then
  echo "$0: File not found: $ocsol"
  exit -1
fi

echo "ocsol, nreg are: $ocsol $nreg"

nrows=`dmlist ${ocsol} counts`

echo "nrows is $nrows"

dt=`echo "scale=4;${nrows-1}/(${nreg}-1)" | bc -l`

echo "dt is $dt"

echo "# Region file format: DS9 version 4.0" > $outf
echo "global color=white font=\"helvetica 10 normal\" select=1 highlite=1 edit=1 move=1 delete=1 include=1 fixed=0 source" >> $outf
echo "fk5" >> $outf

row=1
irow=1

while [ $irow -le ${nrows} ]
do
  ra=`dmlist $ocsol"[col ssora]" opt=data,clean rows=$row:$row | tail -1 | awk '{print $1}'`
  dec=`dmlist $ocsol"[col ssodec]" opt=data,clean rows=$row:$row | tail -1 | awk '{print $1}'`
  echo "ssora,dec are: $ra $dec"
  row=`echo "$row+$dt" | bc -l`
  irow=`echo "scale=0;$row/1" | bc -l`
  if [ x$prior_ra != "x" ]
  then
    rad=`echo "scale=8;(sqrt(($ra-$prior_ra)*($ra-$prior_ra))+sqrt(($dec-$prior_dec)*($dec-$prior_dec)))/2.0" | bc -l`
    radasec=`echo "scale=2;$rad*3600.0" | bc -l`
    echo "circle($prior_ra,$prior_dec,$radasec\")" >> $outf
  fi
  prior_ra=$ra
  prior_dec=$dec
done

echo "circle($prior_ra,$prior_dec,$radasec\") # color=green" >> $outf