#!/bin/ksh -e
#ryj: for qboi1* runs, select 30 pressure levels requested by QBOi project
#ryj: for qboix* runs, don't select, i.e output all 50 levels

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

plv="   30";
p01="-1400"; p02="-0100"; p03="-0150"; p04="-0200"; p05="-0300";
p06="-0500"; p07="-0700"; p08="   10"; p09="   15"; p10="   20";
p11="   30"; p12="   40"; p13="   50"; p14="   60"; p15="   70";
p16="   85"; p17="  100"; p18="  120"; p19="  150"; p20="  175";
p21="  200"; p22="  250"; p23="  300"; p24="  400"; p25="  500";
p26="  600"; p27="  700"; p28="  850"; p29="  925"; p30=" 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
echo "SELLEV    $p28$p29$p30                                            ">> ic.sellev

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

release inp ${pfx}_${out}
access  inp ${pfx}_${inp}

if [[ $pfx == sa_qboix* ]]; then
  cp inp ${pfx}_${out} 
  release inp
else
  sellev  inp ${pfx}_${out} input=ic.sellev
  release inp
fi