#!/bin/ksh -e
#ryj: select 31 standard pressure levels for CCMI project

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

plv="   31";
p01="-1100"; p02="-1200"; p03="-1300"; p04="-1500"; p05="-0100"; p06="-0150";
p07="-0200"; p08="-0300"; p09="-0500"; p10="-0700"; p11="   10"; p12="   15";
p13="   20"; p14="   30"; p15="   50"; p16="   70"; p17="   80"; p18="   90";
p19="  100"; p20="  115"; p21="  130"; p22="  150"; p23="  170"; p24="  200";
p25="  250"; p26="  300"; p27="  400"; p28="  500"; p29="  700"; p30="  850";
p31=" 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$p31                                        ">> ic.sellev

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

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

#commented out by ryj
#several variables need extra processes before selecting 31 pressure levels
#var=`ggstat inp| egrep "GRID|ZONL"|head -1|cut -c25-28`
#if [ $var = "GZ" -o $var = "PHI" ]; then
#  #for zg in day table and ZAmon table
#  echo "C*XLIN       9.80616        0.         1" | ccc xlin inp tmp
#  sellev tmp ${pfx}_${out} input=ic.sellev
#elif [ $var = "TDNX" ]; then
#  #for accelgw in ZAmon table
#  access inp2 ${pfx}_tdox_p
#  add inp inp2 tmp
#  sellev tmp ${pfx}_${out} input=ic.sellev
#  release inp2
#elif [ $var = "VPTP" ]; then
#  #for vt100 in ZAmon table
#  echo "SELECT.   STEPS         1 999999999    1 LEVS  100  100 NAME  ALL"|ccc select inp ${pfx}_${out} 
#else
#  sellev inp ${pfx}_${out} input=ic.sellev
#fi
#
#release inp