#!/bin/csh -f

source /star/etc/login  >& /dev/null
source /star/etc/cshrc  >& /dev/null
kappa >& /dev/null

if (`alias kappa` == "") then
  echo ""
  echo "Fatal: Initialization error of kappa environment"
  echo ""
  exit
endif

set kap = $KAPPA_DIR

if ($#argv < 2) set argv = "help"

if ($argv[1] =~ *"help"* || $argv[1] == "-h") then
  echo ""
  echo "Use: sp_col  file  v1:v2"
  echo "     e.g. 'sp_col myfile -300:-100'"
  echo ""
  echo "Note: the boundaries must be in increasing order!"
  echo ""
  exit
endif

set infile = $argv[1]
set infile = `echo "$infile" | cut -d\. -f1`
set infile = "${infile}.sdf"
shift argv
if !(-e "$infile") then
  echo ""
  echo "Fatal: file $infile not found"
  echo ""
  exit
endif

echo ""
echo -n "Collapsing ${infile} over range "

# ----------------------------------------------------------------------
# Collapse cube over range

    set parts = `echo "$argv" | sed s/":"/" "/g`
    set regions = ""
    @ ii = 1
    while ($ii <= 2) 
      set part = ${parts[$ii]}
      if ("$part" !~ *'.'*) then
        set part = "${part}."
      endif
      if ($ii == 1) then
        set regions = "low=${part}"
      else
        set regions = "${regions} high=${part}"
      endif
      @ ii ++
    end
    unset parts

    echo "$regions"
    set outfile = `echo "${infile}" | cut -d\. -f1`"_col"
    ${kap}/collapse in=${infile} out=${outfile} axis=3 ${regions}

    echo "The collapse file has been saved as:"
    echo "   '${outfile}.sdf'"
    echo ""