#!/bin/sh
####################################################
#
# script: clim_climate_avrg.sh
# purpose: pick up the file for ensemble mean
#
#  $1 ... year month day hour
#  $2 ... start year, say 1979
#  $3 ... end year, say 1999
#  $4 ... ensemble mean directory
#
if [ $# -ne 4 ] ; then
        $DATA/err_exit "$0: argument error"
fi
yyyymmddhh=$1
yyyystart=$2
yyyyend=$3
ensmdir=$4
year=` echo $yyyymmddhh | cut -c1-4`
month=`echo $yyyymmddhh | cut -c5-6`
day=`  echo $yyyymmddhh | cut -c7-8`
hour=` echo $yyyymmddhh | cut -c9-10`
ys=` echo $yyyystart | cut -c3-4`
ye=` echo $yyyyend   | cut -c3-4`
nmon=1

vmonth=` echo $VYYYYMM | cut -c5-6`
vyear=` echo $VYYYYMM | cut -c1-4`
mkdir -m 775 -p $ensmdir/$VYYYYMM

while [ $nmon -le 6 ]
do

   yyyy=$yyyystart
   mm=`expr $vmonth + $nmon`
   hhyy=$yyyy
   vy=$vyear
   if [ $mm -gt 12 ]
   then
        mm=`expr $mm - 12`
        hhyy=`expr $hhyy + 1`
        vy=`expr $vy + 1`
   fi
   if [ $mm -lt 10 ]; then mm=0$mm; fi
   rm -f flx.list$$ pgb.list$$ var.list$$

   while [ $yyyy -le $yyyyend ]
   do

     for file in flx pgb
     do
       ensmfile=ensmean.$file.$hhyy$mm.avrg.grib
       if [ -s $ensmdir/$yyyy$month/$ensmfile ]
       then
           cp $ensmdir/$yyyy$month/$ensmfile $ensmfile
           echo $ensmfile >>$file.list$$
       fi
     done
     
     yyyy=`expr $yyyy + 1`
     hhyy=`expr $hhyy + 1`

   done

   for file in flx pgb
   do
     ensmfile=climmean.$file.$vy$mm.avrg.grib
     $USHclim/clim_grib_average.sh $file.list$$ $ensmfile
     if [ -s $ensmfile ]
     then
         cp $ensmfile $ensmdir/$vyear$vmonth/$ensmfile
         if [ $SENDDBN = YES ]
         then
             $DBNROOT/bin/dbn_alert MODEL CFS_FCST $job $ensmdir/$vyear$vmonth/$ensmfile
         fi
   fi
   done

   nmon=`expr $nmon + 1`

done
~