#!/bin/sh
# runs ssi

# input files: $prepqm $sges $oldsig
# output files: $sanl

cd $DATA
set -x
date

divterrs=$FIXcdas/cdas_divterrs2812645
satv28=$FIXcdas/cdas_v28newx
eofs12628=$FIXcdas/cdas_eofs28126

# convert sges, bges and oldsig to double precision

export pgm=cdas_sig2sig; . prep_step
startmsg
$EXECcdas/cdas_sig2sig -toR8 $bges bges.dbl
export err=$?; err_chk
$EXECcdas/cdas_sig2sig -toR8 $sges sges.dbl
export err=$?; err_chk
$EXECcdas/cdas_sig2sig -toR8 $oldsig oldsig.dbl
export err=$?; err_chk

export pgm=cdas_ssi; . prep_step
# input
export XLFUNIT_30="$prepqm"
export XLFUNIT_35="sges.dbl"
export XLFUNIT_36="oldsig.dbl"
export XLFUNIT_37="bges.dbl"
export XLFUNIT_47="$divterrs"
export XLFUNIT_48="$satv28"
export XLFUNIT_49="$eofs12628"

# output
export XLFUNIT_51="sanl.dbl"

IGEN=2
startmsg
$EXECcdas/cdas_ssi <<EOF >> $pgmout 2> errfile
&NAMANAL igen=$IGEN,JCAP=62,NLATH=48,NLON=192,NSIG=28,niter=100,miter=1,
   a=.25,.33,.42,.45,ampdivt=.7,dampdivt=.8,grosst=10.,grossst=10.,
    grossw=10.,grossp=10.,grossq=5.,grosspw=10.,
/
EOF
export err=$?; err_chk

# convert sigma analysis to single precision
export pgm=cdas_sig2sig; . prep_step
startmsg
$EXECcdas/cdas_sig2sig -toR4 sanl.dbl $sanl
export err=$?; err_chk

# clean up
rm bges.dbl sges.dbl oldsig.dbl sanl.dbl