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

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

plv="   17";

p01="   10"; 
p02="   20";
p03="   30"; 
p04="   50"; 
p05="   70"; 
p06="  100"; 
p07="  150"; 
p08="  200";
p09="  250"; 
p10="  300"; 
p11="  400"; 
p12="  500"; 
p13="  600"; 
p14="  700";
p15="  850"; 
p16="  925"; 
p17=" 1000"; 

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