#!/bin/ksh -e
#ryj: select 20 pressure levels (geopotential height) requested by CMAMdyn project

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

plv="   20";
p01="-0100"; p02="-0200"; p03="-0300"; p04="-0500"; p05="-0700";
p06="   10"; p07="   20"; p08="   30"; p09="   40"; p10="   50";
p11="   70"; p12="  100"; p13="  150"; p14="  200"; p15="  300";
p16="  400"; p17="  500"; p18="  700"; p19="  850"; p20=" 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                            ">> ic.sellev

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

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

#change unit into metre
echo "             9.80616        0.         1" | ccc xlin inp tmp

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