#!/bin/sh
set -x
#
#########################################################################
# Usage: ofs_correct_forcing.sh                                         #
#                                                                       #
# Description: Correct GFS near-surface air temperature                 #
#              Originally, air temperature is in the file               #
#              forcing.airtmp.[ab]It is moved to forcing.airtm1.[ab]    #
#              GFS surface temperature is read from forcing.surtem.[    #
#              ab] and air temperature is corrected where surface       #
#              temperature is below freezing. The result is written to  #
#              forcing.airtmp.[ab]                                      #
#                                                                       #
# History:                                                              #
#    02-02-2007 Ilya Rivin                                              #
#########################################################################

echo "*** Started script $0"

# It is assumed that forcing.imput files
# forcing.airtmp.[ab] forcing.surtem.[ab] and regional.mask.[ab]
# are already in the working directory $DATA

cd $DATA

# prepare input file 
blines=`cat ${DATA}/forcing.surtem.b | wc -l`
num_frames=`grep surtem ${DATA}/forcing.surtem.b | wc -l`
ihead=`expr ${blines} \- ${num_frames} `
sname=`grep span ${DATA}/forcing.surtem.b |head -1 |cut -c1-10`
aname=`grep span ${DATA}/forcing.airtmp.b |head -1 |cut -c1-10`

cat > correct_forcing.in <<EOF
$ihead
$num_frames
$aname
$sname
EOF

mv ${DATA}/forcing.airtmp.a ${DATA}/forcing.airtm1.a 
mv ${DATA}/forcing.airtmp.b ${DATA}/forcing.airtm1.b 

$EXECofs/ofs_correct_forcing < correct_forcing.in >>$pgmout 2>errfile
export err=$?; err_chk

echo "*** Finished script $0"