#!/bin/sh
set -e
# 3D field multiplie 2D field -- DY, Feb/2019.

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

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

# Get CDO / TEMPORARY!
export PATH=/fs/ssm/hpco/exp/mib002/anaconda/anaconda-4.4.0/anaconda_4.4.0_ubuntu-14.04-amd64-64/envs/cdo-1.9.0/bin:$PATH

release inp1 inp2
rm -f tmp*

# access inp1 & inp2
access inp1 ${pfx}_${inp1}.nc
access inp2 ${pfx}_${inp2}.nc

# Append inp2 to inp1
cp inp1 tmp1
chmod ug+w tmp1
ncks -A inp2 tmp1

# inp1 (3D) * inp2 (2D)
ncap2 -O -s "${inp1}=${inp1}*${inp2}" tmp1 tmp2

# remove $inp2
ncks -v ${inp2} -x tmp2 ${pfx}_${out}

release inp1 inp2
rm -f tmp*