#!/bin/sh
set -e
# subtract dhf_fslo from dhf_olrc

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

pfx=$1
inp1=$2
out=$3

access olrc ${pfx}_${inp1}
access fslo sa_extract_fslo1_1979_1979_dhf_fslo

# Convert from DATA TO GRID
echo "C*RELABL   DATA
C*         GRID" | ccc relabl olrc dg_olrc
echo "C*RELABL   DATA
C*         GRID" | ccc relabl fslo dg_fslo

#/home/esm/PRODUCTION/cmip5_scripts/gsubl dg_olrc dg_fslo ${pfx}_${out}
gsub dg_olrc dg_fslo ${pfx}_${out}

release olrc fslo dg_olrc dg_fslo x