#!/bin/sh
set -x
#
REGION=$1
if [ "$REGION" = AK ] 
then
  $USHhrly/marine_stats_getrep_ak.sh $DATA/archive.sort
  cp $DATA/ak.ships $DATA/ships1
  EXEC2=$EXEChrly/marine_stats_genak
  prog=ak
  wc -l $DATA/ships1 > $DATA/numrep.$prog
elif [ "$REGION" = ISS ]
then
  $USHhrly/marine_stats_getrep_iss.sh $DATA/archive.sort
  cp $DATA/iss.ships $DATA/ships1
  EXEC2=$EXEChrly/marine_stats_geniss
  prog=iss
  wc -l $DATA/ships1 > $DATA/numrep.$prog
fi

#
export pgm=$EXEC2
. prep_step
startmsg
#
#input files
export XLFRTEOPTS="unit_vars=yes"     # Allow overriding default names.
#
#input files
export XLFUNIT_11="$DATA/ships1"   
export XLFUNIT_12="$DATA/numrep.$prog"   
#
#output files
export XLFUNIT_40="$DATA/slp.stats.$prog"
export XLFUNIT_41="$DATA/air.stats.$prog"
export XLFUNIT_42="$DATA/wdi.stats.$prog"
export XLFUNIT_43="$DATA/wsp.stats.$prog"
export XLFUNIT_44="$DATA/sst.stats.$prog"
#
#run compiled fortran
#
$EXEC2 > $DATA/out.$prog
#
#cat output files together
cat $DATA/slp.stats.$prog $DATA/air.stats.$prog > $DATA/tmp1                       
cat $DATA/wdi.stats.$prog $DATA/wsp.stats.$prog > $DATA/tmp2
cat $DATA/sst.stats.$prog > $DATA/tmp3       

cat $DATA/tmp1 $DATA/tmp2 $DATA/tmp3 > $DATA/$prog.stats

if [ "$REGION" = AK ] 
then
export pgm=newstats
. prep_step
startmsg
#
#input files
export XLFRTEOPTS="unit_vars=yes"     # Allow overriding default names.

export XLFUNIT_11="$FIXhrly/ship_names"
export XLFUNIT_12="$prog.stats"
export XLFUNIT_20="new$prog.stats"
$EXEChrly/marine_stats_newstats

fi
#
exit
#