#!/bin/sh
set -x

############################################
# Script to setup and compute monthly means
# Originally written by Shrinivas Moorthi
# Updated by Patrick Tripp - Sept 2010
############################################


export OUTDIR=$MONTHDIR

cd $RUNDIR

errs=0

export dogribchk=${dogribchk:-YES}
echo "dogribchk is $dogribchk"

$USHcfs/cfs_cdas_timeavg_anl.sh $type $cdate1 $cdate2 $nhours $prefix1 $prefix2 $fhini $fhmax $fhout
if [[ $? -ne 0 ]] ; then
  echo "ERROR: cfs_cdas_timeavg_anl.sh returned non-zero exit"
  ((errs+=1))
fi

# Create low-res versions
if [[ $makelow -eq 1 ]] ; then

  if [ $type = "fcst" ] ; then

    for fh in $fhlist
    do
      ifile=$OUTDIR/${prefix1}${prefix2}${fh}$SUFOUT
      if [[ ! -s $ifile ]] ; then
        echo "ERROR: $ifile does not exist"
        ((errs+=1))
      fi
      ofile=$OUTDIR/${prefix1}${prefixo}${fh}$SUFOUT
      /bin/rm -f $ofile

      # pgb
      if [[ $prefix1 == "pgb" ]] ; then
        /bin/rm -f $RUNDIR/pgb25a $RUNDIR/pgb25b

        $wgrib $ifile|grep -i :kpds6=100:|$copygb -xkw -i4,0,35 -g$ogrid $ifile $RUNDIR/pgb25a
        if [[ $? -ne 0 ]] ; then
          echo "ERROR: copygb returned non-zero exit for $ifile $ofile"
          ((errs+=1))
        fi

        $wgrib $ifile|grep -vi :kpds6=100:|$copygb -xkw -i0 -g$ogrid $ifile $RUNDIR/pgb25b
        if [[ $? -ne 0 ]] ; then
          echo "ERROR: copygb returned non-zero exit for $ifile $ofile"
          ((errs+=1))
        fi

        cat $RUNDIR/pgb25a $RUNDIR/pgb25b > $ofile
        if [[ $? -ne 0 ]] ; then
          echo "ERROR: cat returned non-zero exit for $ofile"
          ((errs+=1))
        fi

      else
        # non pgb
        $copygb -g$ogrid -i0 -x $ifile $ofile
        if [[ $? -ne 0 ]] ; then
          echo "ERROR: copygb returned non-zero exit for $ifile $ofile"
          ((errs+=1))
        fi
      fi
      if [[ ! -s $ofile ]] ; then
        echo "ERROR: $ofile does not exist"
        ((errs+=1))
      fi
    done
  else   # non fcst, no fhlist

    ifile=$OUTDIR/${prefix1}${prefix2}$SUFOUT
    if [[ ! -s $ifile ]] ; then
      echo "ERROR: $ifile does not exist"
      ((errs+=1))
    fi

    ofile=$OUTDIR/${prefixo}$SUFOUT
    /bin/rm -f $ofile

    # pgb
    if [[ $prefix1 == "pgb" ]] ; then

      /bin/rm -f $RUNDIR/pgb25a $RUNDIR/pgb25b

      $wgrib $ifile|grep -i :kpds6=100:|$copygb -xkw -i4,0,35 -g$ogrid $ifile $RUNDIR/pgb25a
      if [[ $? -ne 0 ]] ; then
        echo "ERROR: copygb returned non-zero exit for $ifile $ofile"
        ((errs+=1))
      fi

      $wgrib $ifile|grep -vi :kpds6=100:|$copygb -xkw -i0 -g$ogrid $ifile $RUNDIR/pgb25b
      if [[ $? -ne 0 ]] ; then
        echo "ERROR: copygb returned non-zero exit for $ifile $ofile"
        ((errs+=1))
      fi

      cat $RUNDIR/pgb25a $RUNDIR/pgb25b > $ofile
      if [[ $? -ne 0 ]] ; then
        echo "ERROR: cat returned non-zero exit for $ofile"
        ((errs+=1))
      fi

    else
      # non pgb
      $copygb -g$ogrid -i0 -x $ifile $ofile
      if [[ $? -ne 0 ]] ; then
        echo "ERROR: copygb returned non-zero exit for $ifile $ofile"
        ((errs+=1))
      fi
    fi

    if [[ ! -s $ofile ]] ; then
      echo "ERROR: $ofile does not exist"
      ((errs+=1))
    fi
  fi
fi

# Exit with error count

exit $errs