#!/bin/sh
set -e
# soil water capacity

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

pfx=$1
desc=$2
out=$3

release desc
access  desc ${pfx}_${desc}

if [[ $(ggstat desc|grep MASK) ]]; then
  # for CanESM2; land fraction is mask
  echo "C*SELECT.  STEP         0 999999999    1         1    3 NAME MASK  DZG PORG" | ccc select desc mask dzg porg
  mlt dzg porg dzgporg
  echo "          3" | ccc bins dzgporg dzgporgbin
  echo "                  3." | ccc xlin dzgporgbin dzgporgbin1
  echo "C*FMSKPLT         -1 NEXT   EQ     -1.00    1     1.E38    1" | ccc fmskplt dzgporgbin1 ${pfx}_${out} mask
  
  release desc dzg porg dzgporg dzgporgbin dzgporgbin1 mask
else
  # for CanESM5; land fraction is flnd (dry land only)
  echo "C*SELECT.  STEP         0 999999999    1         1    3 NAME FLND  DZG PORG" | ccc select desc flnd dzg porg
  mlt dzg porg dzgporg
  echo "          3" | ccc bins dzgporg dzgporgbin
  echo "                  3." | ccc xlin dzgporgbin dzgporgbin1

  echo "C*FMSKPLT         -1 NEXT   GT      0.00    1     1.E38    1" | ccc fmskplt dzgporgbin1 ${pfx}_${out} flnd

  release desc dzg porg dzgporg dzgporgbin dzgporgbin1 flnd 
fi