#!/bin/sh
set -e
# glacier fraction, %

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; glacier is dpth(-4) over mask
  echo "C*SELECT.  STEP         0         0    1         1    1 NAME DPTH" | ccc select desc dpth
  echo "C*FMASK           -1 NEXT   EQ       0.0" | ccc fmask dpth fmsk
  echo "C*XLIN          100.        0." | ccc xlin fmsk ${pfx}_${out}
  
  release desc dpth fmsk
else
  # for CanESM5; glacier is dpth(-4) over flnd
  echo "C*SELECT.  STEP         0         0    1         1    1 NAME DPTH FLND" | ccc select desc dpth flnd
  echo "C*FMASK           -1 NEXT   LT       0.0" | ccc fmask dpth fmsk

  mlt flnd fmsk glac 
  echo "C*XLIN          100.        0." | ccc xlin glac ${pfx}_${out}

  release desc flnd dpth glac fmsk
fi