#!/bin/sh
set -e

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 nemo_mesh_mask
rm -f tmp*

access nemo_mesh_mask nemo_mesh_mask_rc3.nc
access  inp1 ${pfx}_${inp1}.nc
access  inp2 ${pfx}_${inp2}.nc

# deptho = e3t * tmask
ncap2 -O -s "deptho=e3t*tmask" nemo_mesh_mask tmp0
cdo -select,name=deptho tmp0 y.nc

cdo mul inp1 y.nc z1.nc
cdo mul inp2 y.nc z2.nc
cdo add z1.nc z2.nc z.nc
rm z1.nc z2.nc y.nc

# integral of -e3t*tmask*data along depth
cdo vertsum z.nc ${pfx}_${out}

# Rename t to time_counter
#ncrename -d time,time_counter tmp2 ${pfx}_${out}

chmod uo+w ${pfx}_${out}

# Select nav_lon & nav_lat
#ncks -O -C -a -v nav_lon,nav_lat inp1 tmp3

# Append nav_lat & nav_lon
#ncks -A tmp3 ${pfx}_${out}

release inp1 inp2 nemo_mesh_mask
rm -f tmp* z.nc