#!/bin/sh
set -e
# Compute dry static energy s=cp*t+gz at the lowest model level

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

pfx=$1
inp1=$2 # temp
inp2=$3 # phi
out=$4 # dses

release t gz
access  t  ${pfx}_${inp1}
access  gz ${pfx}_${inp2}

# select lowest model level
echo "C*SELLEV      1  995" | ccc sellev t  t1
echo "C*SELLEV      1  995" | ccc sellev gz gz1

# s=cp*t+gz
echo "C*XLIN        1004.5" | ccc xlin t1 cpt
add cpt gz1 ${pfx}_${out}

release t t1 gz gz1 cpt .ccc_cards