#!/bin/sh
set -e
# Set the small non-zero values under sea ice that were caused by the 1:4 packing density as 0.
# inp1 is oflg/ohfs/ohfl; inp2 is gc; desc is gz_data_description.

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

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

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

nm=`ggstat inp1 | grep GRID | awk '{print $4}' | head -1`
echo "C*SELECT   STEP         0 999999999    1     -900199999 NAME $nm" | ccc select inp1 inp3
echo "C* FMSKPLT        -1 NEXT   GT      -1.0    1        1.    1" | ccc fmskplt inp2 mask1
echo "C*SELECT   STEP         0 999999999    1     -900199999 NAME  LAT" | ccc select desc olat
echo "C*SELECT   STEP         0 999999999    1     -900199999 NAME  LON" | ccc select desc olon
grido desc betao dao dx dy dz
echo "C* A2OGRID     BILIN" | ccc a2ogrid mask1 olon olat dao betao mask2

echo "C* FMSKPLT        -1 NEXT   LE       0.5    1        0.    1" | ccc fmskplt inp3 ${pfx}_${out} mask2

release inp1 inp2 desc
rm -f inp3 mask1 mask2 olat olon betao dao dx dy dz