echo "------------------------------------------------"
echo "  Verification - ANLVER Processing"
echo "------------------------------------------------"

###################################################
# COMPUTES MONTHLY SUMMARY OF S1 SCORES FROM      #
# VARIOUS MODELS                                  #
###################################################

msg="VERFM_ANLVER_RPTS HAS BEGUN"
postmsg "$jlogfile" "$msg"

set -x

cd $DATA

export XLFRTEOPTS="unit_vars=yes"

########################################################
# COMPUTES MONTHLY SUMMARY OF S1 SCORES FOR MODELS:    #
#    NAM - 12HR FORECASTS                              #
#    GFS - 12HR FORECASTS                              #
#    GFL - 24HR FORECASTS                              #
########################################################

model_list=" GFL NAM GFS "

for mfile in $model_list
 do

MODEL=$DATA/${mfile}.${mm}.in
MAX=$DATA/s1max.${mfile}.${mm}
TRU=$DATA/s1tru.${mfile}.${mm}

#################
#  input files  #
#################

export XLFUNIT_10="arkv_verdate.in"
export XLFUNIT_12="$FIXhrly/anlver/arkv_sum${mfile}.in"
export XLFUNIT_50="$MODEL"

#################
# output files  #
#################

export XLFUNIT_30="$MAX"
export XLFUNIT_40="$TRU"

$EXEChrly/verf_s1sum  1> s1sum${mfile}.${mm} 2>s2sum${mfile}.err
 export err=$?; err_chk

done

##############################################################
# COMPUTES MONTHLY SUMMARY OF WMO MSLP,HGT,TEMP SCORES AND   #
# COMPUTES MONTHLY SUMMARY OF WMO WIND SCORES FOR BOTH THE   #
# GFL AND GFL MODELS                                         #
##############################################################

cyc_list=" 00 12 "

for ct in $cyc_list
 do

XMSLP=${mm}.xmslp.in${ct}
X850Z=${mm}.x850z.in${ct}
X850T=${mm}.x850t.in${ct}
X500Z=${mm}.x500z.in${ct}
X500T=${mm}.x500t.in${ct}
X250Z=${mm}.x250z.in${ct}
X250T=${mm}.x250t.in${ct}
V850=${mm}.w850.in${ct}
V500=${mm}.w500.in${ct}
V250=${mm}.w250.in${ct}
XYZ=xyzsum.${mm}
WND=windsum.${mm}

################
# input files  #
################

export XLFUNIT_10="arkv_verdate.in"
export XLFUNIT_21="$ANLVER_MONTHLYdata/$XMSLP"
export XLFUNIT_22="$ANLVER_MONTHLYdata/$X850Z"
export XLFUNIT_23="$ANLVER_MONTHLYdata/$X850T"
export XLFUNIT_24="$ANLVER_MONTHLYdata/$X500Z"
export XLFUNIT_25="$ANLVER_MONTHLYdata/$X500T"
export XLFUNIT_26="$ANLVER_MONTHLYdata/$X250Z"
export XLFUNIT_27="$ANLVER_MONTHLYdata/$X250T"

export XLFUNIT_28="$ANLVER_MONTHLYdata/$V850"
export XLFUNIT_29="$ANLVER_MONTHLYdata/$V500"
export XLFUNIT_30="$ANLVER_MONTHLYdata/$V250"

#################
# output files  #
#################

export XLFUNIT_70="$XYZ"
export XLFUNIT_80="$WND"

$EXEChrly/verf_xsum  1>xsum${ct}.${mm} 2>xsum${ct}.${mm}.err
 export err=$?; err_chk

$EXEChrly/verf_wsum  1>wsum${ct}.${mm} 2>wsum${ct}.${mm}.err
 export err=$?; err_chk

cat $XYZ $WND > vwmomon${ct}.${mm}

done

################
# input files  #
################

export XLFUNIT_20="vwmomon00.${mm}"
export XLFUNIT_30="vwmomon12.${mm}"

#################
# output files  #
#################

export XLFUNIT_40="vwmomon.${mm}"

$EXEChrly/verf_combintbls 
export err=$?; err_chk

##########################################################
# GENERATES MONTHLY VERIFICATION IN NON-STANDARD TABLES  #
# FOR NORTHERN AND SOUTHERN HEMISHPERE AND THE TROPICS   # 
##########################################################

VW=$DATA/vwmomon.${mm}
NHX=$DATA/wmonhx${mm}
SHX=$DATA/wmoshx${mm}
TRX=$DATA/wmotrx${mm}
NHW=$DATA/wmonhw${mm}
SHW=$DATA/wmoshw${mm}
TRW=$DATA/wmotrw${mm}
NHM=$DATA/wmonhm${mm}
SHM=$DATA/wmoshm${mm}
TRP=$DATA/wmotrp${mm}

################
# input file   #
################

export XLFUNIT_40="$VW"

###########################################################
#  output files (temp -- will be appended to arch files   #
#  and deleted, see end of script)                        #
###########################################################

export XLFUNIT_41="$NHX"
export XLFUNIT_42="$SHX"
export XLFUNIT_43="$TRX"
export XLFUNIT_44="$NHW"
export XLFUNIT_45="$SHW"
export XLFUNIT_46="$TRW"
export XLFUNIT_47="$NHM"
export XLFUNIT_48="$SHM"
export XLFUNIT_49="$TRP"
 
########################
# EXECUTING wmoptbls   #
########################

$EXEChrly/verf_wmoptbls 1>wmoptbls.${mm} 2>wmoptbls.err
 export err=$?; err_chk

########################
# EXECUTING wmoptblw   #
########################

$EXEChrly/verf_wmoptblw 1>wmoptblw.${mm} 2>wmoptblw.err
 export err=$?; err_chk

cp wmoptblw.${mm} $SUMAC4_MONTHLYreports/ryhwmx${mm}

cat $NHM >> $DATA/wmonhm$YYYY
cat $SHM >> $DATA/wmoshm$YYYY
cat $TRP >> $DATA/wmotrp$YYYY
cat $NHX >> $DATA/wmonhx$YYYY
cat $SHX >> $DATA/wmoshx$YYYY
cat $TRX >> $DATA/wmotrx$YYYY
cat $NHW >> $DATA/wmonhw$YYYY
cat $SHW >> $DATA/wmoshw$YYYY
cat $TRW >> $DATA/wmotrw$YYYY

####################################################################
# Appends current month's s1/s2 scores to archive for each model   #
####################################################################

  if [ ${mm} -lt '10' ]; then
     export msd=`echo ${mm} | cut -c2`
     sed -e 's/MSL    ${msd}/MSL   0${msd}/g' s1max.GFS.${mm} | sed -e 's/500MB  ${msd}/500MB 0${msd}/g' > h1
     sed -e 's/MSL    ${msd}/MSL   0${msd}/g' s1tru.GFS.${mm} | sed -e 's/500MB  ${msd}/500MB 0${msd}/g' > j1
     sed -e 's/MSL    ${msd}/MSL   0${msd}/g' s1max.NAM.${mm} | sed -e 's/500MB  ${msd}/500MB 0${msd}/g' > h2
     sed -e 's/MSL    ${msd}/MSL   0${msd}/g' s1tru.NAM.${mm} | sed -e 's/500MB  ${msd}/500MB 0${msd}/g' > j2
     sed -e 's/MSL    ${msd}/MSL   0${msd}/g' s1max.GFL.${mm} | sed -e 's/500MB  ${msd}/500MB 0${msd}/g' > h4
     sed -e 's/MSL    ${msd}/MSL   0${msd}/g' s1tru.GFL.${mm} | sed -e 's/500MB  ${msd}/500MB 0${msd}/g' > j4
  else
     export msd=`echo ${mm} | cut -c2`
     sed -e 's/MSL    ${msd}/MSL   ${msd}/g' s1max.GFS.${mm} | sed -e 's/500MB  ${msd}/500MB ${msd}/g' > h1
     sed -e 's/MSL    ${msd}/MSL   ${msd}/g' s1tru.GFS.${mm} | sed -e 's/500MB  ${msd}/500MB ${msd}/g' > j1
     sed -e 's/MSL    ${msd}/MSL   ${msd}/g' s1max.NAM.${mm} | sed -e 's/500MB  ${msd}/500MB ${msd}/g' > h2
     sed -e 's/MSL    ${msd}/MSL   ${msd}/g' s1tru.NAM.${mm} | sed -e 's/500MB  ${msd}/500MB ${msd}/g' > j2
     sed -e 's/MSL    ${msd}/MSL   ${msd}/g' s1max.GFL.${mm} | sed -e 's/500MB  ${msd}/500MB ${msd}/g' > h4
     sed -e 's/MSL    ${msd}/MSL   ${msd}/g' s1tru.GFL.${mm} | sed -e 's/500MB  ${msd}/500MB ${msd}/g' > j4
  fi

    mv h1 $ANLVER_MONTHLYdata/s1max.GFS.${mm}
    mv j1 $ANLVER_MONTHLYdata/s1tru.GFS.${mm}
    mv h2 $ANLVER_MONTHLYdata/s1max.NAM.${mm}
    mv j2 $ANLVER_MONTHLYdata/s1tru.NAM.${mm}
    mv h4 $ANLVER_MONTHLYdata/s1max.GFL.${mm}
    mv j4 $ANLVER_MONTHLYdata/s1tru.GFL.${mm}


  if [ ${mm} -eq '01' ]; then
   echo 'MONYR NN GFSMAX MSL 12-72HR  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> $ANLVER_MONTHLYreports/s1GFSmsl
   echo 'MONYR NN GFSMAX 500 12-72HR  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> $ANLVER_MONTHLYreports/s1GFS500 
   echo 'MONYR NN NAMMAX MSL 12-48HR  59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> $ANLVER_MONTHLYreports/s1NAMmsl
   echo 'MONYR NN NAMMAX 500 12-48HR  59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> $ANLVER_MONTHLYreports/s1NAM500
   echo 'MONYR NN GFLMAX MSL 24-240HR (X24)  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33'|cat >> $ANLVER_MONTHLYreports/s1GFLmsl
   echo 'MONYR NN GFLMAX 500 24-240HR (X24)  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' |cat >> $ANLVER_MONTHLYreports/s1GFL500
   echo 'MONYR NN GFSTRU MSL 12-72HR  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> $ANLVER_MONTHLYreports/s2GFSmsl
   echo 'MONYR NN GFSTRU 500 12-72HR  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> $ANLVER_MONTHLYreports/s2GFS500
   echo 'MONYR NN NAMTRU MSL 12-48HR  59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> $ANLVER_MONTHLYreports/s2NAMmsl
   echo 'MONYR NN NAMTRU 500 12-48HR  59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> $ANLVER_MONTHLYreports/s2NAM500
   echo 'MONYR NN GFLTRU MSL 24-240HR (X24)  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' |cat >> $ANLVER_MONTHLYreports/s2GFLmsl
   echo 'MONYR NN GFLTRU 500 24-240HR (X24)  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' |cat >> $ANLVER_MONTHLYreports/s2GFL500
  fi

  grep -v "MONYR" $ANLVER_MONTHLYreports/s1GFSmsl | cat >> a1
  head -6 $ANLVER_MONTHLYdata/s1max.GFS.${mm} >> a1
  rm -f $ANLVER_MONTHLYreports/s1GFSmsl
  echo 'MONYR NN GFSMAX MSL 12-72HR  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> a1
  cat a1 | cat >> $ANLVER_MONTHLYreports/s1GFSmsl

  grep -v "MONYR" $ANLVER_MONTHLYreports/s1GFS500 |cat >> b1
  tail -6 $ANLVER_MONTHLYdata/s1max.GFS.${mm} >> b1
  rm -f $ANLVER_MONTHLYreports/s1GFS500
  echo 'MONYR NN GFSMAX 500 12-72HR  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> b1
  cat b1 | cat >> $ANLVER_MONTHLYreports/s1GFS500

  grep -v "MONYR" $ANLVER_MONTHLYreports/s1NAMmsl |cat >> c1
  head -5 $ANLVER_MONTHLYdata/s1max.NAM.${mm} >> c1
  rm -f $ANLVER_MONTHLYreports/s1NAMmsl
  echo 'MONYR NN NAMMAX MSL 12-48HR  59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> c1
  cat c1 | cat >> $ANLVER_MONTHLYreports/s1NAMmsl

  grep -v "MONYR" $ANLVER_MONTHLYreports/s1NAM500 |cat >> d1
  tail -5 $ANLVER_MONTHLYdata/s1max.NAM.${mm} >> d1
  rm -f $ANLVER_MONTHLYreports/s1NAM500
  echo 'MONYR NN NAMMAX 500 12-48HR  59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> d1
  cat d1 | cat >> $ANLVER_MONTHLYreports/s1NAM500

  grep -v "MONYR" $ANLVER_MONTHLYreports/s1GFLmsl |cat >> e1
  head -6 $ANLVER_MONTHLYdata/s1max.GFL.${mm} >> e1
  rm -f $ANLVER_MONTHLYreports/s1GFLmsl
  echo 'MONYR NN GFLMAX MSL 24-240HR (X24)  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33'|cat >> e1
  cat e1 | cat >> $ANLVER_MONTHLYreports/s1GFLmsl

  grep -v "MONYR" $ANLVER_MONTHLYreports/s1GFL500 |cat >> f1
  tail -6 $ANLVER_MONTHLYdata/s1max.GFL.${mm} >> f1
  rm -f $ANLVER_MONTHLYreports/s1GFL500
  echo 'MONYR NN GFLMAX 500 24-240HR (X24)  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' |cat >> f1
  cat f1 | cat >> $ANLVER_MONTHLYreports/s1GFL500

  grep -v "MONYR" $ANLVER_MONTHLYreports/s2GFSmsl |cat >> a2
  head -6 $ANLVER_MONTHLYdata/s1tru.GFS.${mm} >> a2
  rm -f $ANLVER_MONTHLYreports/s2GFSmsl
  echo 'MONYR NN GFSTRU MSL 12-72HR  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> a2
  cat a2 | cat >> $ANLVER_MONTHLYreports/s2GFSmsl

  grep -v "MONYR" $ANLVER_MONTHLYreports/s2GFS500 |cat >> b2
  tail -6 $ANLVER_MONTHLYdata/s1tru.GFS.${mm} >> b2
  rm -f $ANLVER_MONTHLYreports/s2GFS500
  echo 'MONYR NN GFSTRU 500 12-72HR  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> b2
  cat b2 | cat >> $ANLVER_MONTHLYreports/s2GFS500

  grep -v "MONYR" $ANLVER_MONTHLYreports/s2NAMmsl |cat >> c2
  head -5 $ANLVER_MONTHLYdata/s1tru.NAM.${mm} >> c2
  rm -f $ANLVER_MONTHLYreports/s2NAMmsl
  echo 'MONYR NN NAMTRU MSL 12-48HR  59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> c2
  cat c2 | cat >> $ANLVER_MONTHLYreports/s2NAMmsl

  grep -v "MONYR" $ANLVER_MONTHLYreports/s2NAM500 |cat >> d2
  tail -5 $ANLVER_MONTHLYdata/s1tru.NAM.${mm} >> d2
  rm -f $ANLVER_MONTHLYreports/s2NAM500
  echo 'MONYR NN NAMTRU 500 12-48HR  59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' | cat >> d2
  cat d2 | cat >> $ANLVER_MONTHLYreports/s2NAM500

  grep -v "MONYR" $ANLVER_MONTHLYreports/s2GFLmsl |cat >> e2
  head -6 $ANLVER_MONTHLYdata/s1tru.GFL.${mm} >> e2
  rm -f $ANLVER_MONTHLYreports/s2GFLmsl
  echo 'MONYR NN GFLTRU MSL 24-240HR (X24)  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' |cat >> e2
  cat  e2 | cat >> $ANLVER_MONTHLYreports/s2GFLmsl

  grep -v "MONYR" $ANLVER_MONTHLYreports/s2GFL500 |cat >> f2
  tail -6 $ANLVER_MONTHLYdata/s1tru.GFL.${mm} >> f2
  rm -f $ANLVER_MONTHLYreports/s2GFL500
  echo 'MONYR NN GFLTRU 500 24-240HR (X24)  GBL,59,W33,E33,49,RGNAR1,LL59,LLW33,LLE33' |cat >> f2
  cat f2 | cat >> $ANLVER_MONTHLYreports/s2GFL500

##########################################################################################
# Saving the monthly s1max.GFS and s1max.NAM 36hr and 72hr S1 scores for skill's report  #
##########################################################################################

cd $DATA

  echo "NAM" | cat >> s1nam
  tail -4 $ANLVER_MONTHLYdata/s1max.NAM.${mm} >> s1nam

  echo "GFS" | cat >> s1gfs
  tail -4 $ANLVER_MONTHLYdata/s1max.GFS.${mm} >> s1gfs

model_list=" s1gfs s1nam "

for mfile in $model_list
 do

MODEL=${mfile}
S136=${mfile}"36"
S172=${mfile}"72"


#################
#  input files  #
#################

export XLFUNIT_10="gfs"
export XLFUNIT_50="$MODEL"


#################
# output files  #
#################

export XLFUNIT_30="$S136"
export XLFUNIT_40="$S172"

$EXEChrly/verf_gets1scores 1>output 2>gets1.err
 export err=$?; err_chk

cat $DATA/${mfile}"36" | cat >> $S1_SKILLscores/${mfile}"36500"

cat $DATA/${mfile}"72" | cat >> $S1_SKILLscores/${mfile}"72500"

done

msg="VERFM_ANLVER_RPTS completed normally"
postmsg "$jlogfile" "$msg"

###########################################################################
# GOOD RUN
set +x
echo "**************Script verfm_anlver_rpts COMPLETED NORMALLY ON THE IBM"
echo "**************Script verfm_anlver_rpts COMPLETED NORMALLY ON THE IBM"
echo "**************Script verfm_anlver_rpts COMPLETED NORMALLY ON THE IBM"
set -x
###########################################################################

#------------- end script verfm_anlver_rpts.sh ---------------------