#!/bin/sh
set -e
# maximum root depth, m

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
  release anfile
  access  anfile mc_igm_1850_m12_an
  
  echo "C*SELECT.  STEP         0 999999999    1         0    4 NAME MASK ROOT" | ccc select anfile mask root
  printf "C*RELABL\nC*RELABL                          1\n" | ccc relabl root root1
  
  timmax root1 root2
  printf "C*RELABL\nC*RELABL                0         1\n" | ccc relabl root2  root3
  
  echo "C*FMSKPLT         -1 NEXT   EQ     -1.00    1     1.E38    1" | ccc fmskplt root3 ${pfx}_${out} mask
  
  release anfile root root1 root2 root3 mask desc
else
  # for CanESM5; land fraction is flnd (dry land only)
  #release anfile
  #access  anfile mc_canesm5_1850_m01_an.001   #todo need access model an file
  ln -sf /fs/site[12]/dev/eccc/crd/ccrn/users/cmr203/cmip6nc/mc_canesm5_1850_m01_an.001 anfile
  
  echo "C*SELECT.  STEP         0 999999999    1         0 4001 NAME FLND ROOT" | ccc select anfile flnd root
  printf "C*RELABL\nC*RELABL                          1\n" | ccc relabl root root1
  
  timmax root1 root2
  printf "C*RELABL\nC*RELABL                0         1\n" | ccc relabl root2  root3
  
  echo "C*FMSKPLT         -1 NEXT   GT      0.00    1     1.E38    1" | ccc fmskplt root3 ${pfx}_${out} flnd
  
  release anfile root root1 root2 root3 flnd land desc
fi