#!/bin/sh
set -e
# divide  input file by rho0*g/1.0E+4 (=1025*9.806/1.0E+4=1.005115)
if [ $# -ne 3 ] ; then
  echo "Error in $0: invalid number of arguments $#"
  exit 1
fi
pfx=$1
inp=$2
out=$3
access inp ${pfx}_${inp}
echo "C*FMASK           -1 NEXT   EQ     1.E38    1     1.E38    1    1" |ccc fmask inp inpmask
echo "            1.005115        0.         1" | ccc xlin inp inp2
echo "C* FMSKPLT        -1 NEXT   LT     1.E20                   1 " |ccc fmskplt inp2 ${pfx}_${out} inpmask
release inp
