#!/bin/ksh
set -aex

#######################################################################

set +e
/bin/rm -r *gr *dat
set -e
#
${NCP:-cp} $MASKDIR/cfs_gridsT62yg.dat ./gridsT62yg.dat
${NCP:-cp} $MASKDIR/cfs_gridsMOM3.dat ./gridsMOM3.dat
${NCP:-cp} $MASKDIR/cfs_maskMRFT62.gr ./maskMRFT62.gr
${NCP:-cp} $MASKDIR/cfs_maskMOM3.gr ./maskMOM3.gr
#
export NDATE=${NDATE:-/nwprod/util/exec/ndate}
export wgrib=${wgrib:-/nwprod/util/exec/wgrib}
#
# create fluxes for ocean component (MOM3)
#

if [ -f flxMOM3.gr ]; then
 rm flxMOM3.gr
fi

cpldate=$VDATE     ;# date of current coupling episode
fcpl=$FHCYC        ;# coupling frequency in hours

varlist=':UFLX:|:VFLX:|:DSWRF:sfc:|:USWRF:sfc:|:DLWRF:sfc:|:ULWRF:sfc:|:PRATE:|:LHTFL:|:SHTFL:|:LAND:'

yyyy0=`echo $cpldate | cut -c1-4`   ;#start year of this job
mm0=`echo $cpldate | cut -c5-6` 
dd0=`echo $cpldate | cut -c7-8` 

set +e
rm -r flxcplgrb*
set -e
hh=00
FH=$FHINI
until [[ $((FH+=$FHOUT)) -gt $FHMAX ]]; do
  [[ $FH -lt 10 ]]&&FH=0$FH
  cpldate=`$NDATE $FHOUT $cpldate`
  hh=$((FH-FHINI))
  [[ $hh -lt 10 ]]&&hh=0$hh
  flxcplgrb=flxcplgrb$hh
  if [ ! -f flxf$FH ] ; then
    echo "msg from coupA2O: flxf$FH does not exist."
    exit 8
  fi
# if [ ! -f sigf$cpldate ] ; then
#   echo "msg from coupA2O: sigf$cpldate does not exist."
#   exit 8
# fi
  $wgrib flxf$FH -s | egrep "($varlist)" | $wgrib flxf$FH -o $flxcplgrb -grib -append -s -i
  $force_grib_date_y2k $flxcplgrb $cpldate
  export err=$?;$DATA/err_chk
done

$coupA2Oexec >> $pgmout 2>errfile
export err=$?;$DATA/err_chk

${NCP:-cp} flxMOM3.gr flxMOM3$cpldate.gr