#!/bin/sh
set -e
# convert from volume (Sv) to mass (kg/s) transport - times 1.025E+9;
# combine basins global, arc-atl and ind-pac into one file,
# indicated by ibuf(7)=0(global), 1(arc-atl),2(ind-pac).

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

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

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

echo "C*FMASK           -1 NEXT   EQ     1.E38    1     1.E38    1    1" |ccc fmask inp1 inpmask1
echo "C*FMASK           -1 NEXT   EQ     1.E38    1     1.E38    1    1" |ccc fmask inp2 inpmask2
echo "C*FMASK           -1 NEXT   EQ     1.E38    1     1.E38    1    1" |ccc fmask inp3 inpmask3

echo "            1.025E+9        0." | ccc xlin inp1 inp10
echo "            1.025E+9        0." | ccc xlin inp2 inp20
echo "            1.025E+9        0." | ccc xlin inp3 inp30

echo "C* FMSKPLT        -1 NEXT   LT     1.E38                   1 " |ccc fmskplt inp10 inp11 inpmask1
echo "C* FMSKPLT        -1 NEXT   LT     1.E38                   1 " |ccc fmskplt inp20 inp21 inpmask2
echo "C* FMSKPLT        -1 NEXT   LT     1.E38                   1 " |ccc fmskplt inp30 inp31 inpmask3

echo "                                                 0
C*                                               0" | ccc relabl inp11 inp12
echo "                                                 0
C*                                               1" | ccc relabl inp21 inp22
echo "C*                                               0
C*                                               2" | ccc relabl inp31 inp32

#ryj on July 06, 2015: adjust the output order (from 123 into 231) so as to be consistent with cmor 
rmerge ${pfx}_${out} inp22 inp32 inp12 

release inp1 inp2 inp3
rm -f inpmask1 inpmask2 inpmask3 inp10 inp20 inp30 inp11 inp21 inp31 inp12 inp22 inp32