#!/bin/sh
set -e
# 1. heat flux into sea water due to sea ice thermodynamics (melting/growing of sea ice + conductive heat flux)  = OBEG - BEG.
# 2. Interpolate from agcm to ogcm grid.
# inp1 is obeg; inp2 is beg; desc is gz_data_description.

if [ $# -ne 5 ] ; then
  echo "Error in $0: invalid number of arguments $#"
  exit 1
fi

pfx=$1
inp1=$2
inp2=$3
desc=$4
out=$5

access inp1 ${pfx}_${inp1}
access inp2 ${pfx}_${inp2}
access desc ${pfx}_${desc}

sub inp1 inp2 inpd
echo "C*SELECT   STEP         0 999999999    1     -900199999 NAME  LAT" | ccc select desc olat
echo "C*SELECT   STEP         0 999999999    1     -900199999 NAME  LON" | ccc select desc olon
grido desc betao dao dx dy dz
echo "C* A2OGRID     BILIN" | ccc a2ogrid inpd olon olat dao betao inp3
echo "C*NEWNAM   HFSI" | ccc newnam inp3 ${pfx}_${out}

release inp1 inp2 desc
rm -f olat olon betao dao dx dy dz inpd inp3