#!/bin/ksh -e
#ryj: select 3 pressure levels for HAPPI project

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

plv="    3";

p01="   10"; 
p02="  250"; 
p03="  850"; 

echo "SELLEV    $plv$p01$p02$p03$p04$p05$p06$p07$p08$p09$p10$p11$p12$p13" > ic.sellev
echo "SELLEV    $p14$p15$p16$p17$p18$p19$p20$p21$p22$p23$p24$p25$p26$p27">> ic.sellev

pfx=$1
inp=$2
out=$3

release inp ${pfx}_${out}
access  inp ${pfx}_${inp}
sellev  inp outmp input=ic.sellev

if [[ $inp = gz ]]; then
  echo "C*XLIN       9.80616        0.         1" | ccc xlin outmp ${pfx}_${out}
elif [[ $inp = rh ]]; then
  echo "C*XLIN          100." | ccc xlin outmp ${pfx}_${out}
else
  mv outmp ${pfx}_${out}
fi

release inp outmp