#!/bin/ksh +x
#*********************************#
#                                 #
#     Hough_Functions_uvh_run     #
#                                 #
#*********************************#
#
if [ -z "${1}" ]
then
  echo " "
  echo "   Usage:"
  echo " "
  echo "   ./Hough_Functions_uvh_run cs [Mmax]"
  echo " "
  echo "    cs : bc for Baroclinic Mode or br for Barotropic Mode"
  echo "  Mmax : four digits number of zonal waves (Rhomb+1)"
  echo "         if not given, default Mmax = 126"
  echo " "
  echo " Logics corret the number of digits for Mmax"
  echo " "
  exit
fi
#
# Arguments for ./Hough_Functions_uvh
# ${1} = bc for Baroclinic Mode and br for Barotropic Mode
# ${2} = Wave : rb, kv, mx, gw, ge
# ${3} = Wave Number Plus 1
# ${4} = Meridional Index
# ${5} = Symmetry (1 - Symmetric ; 2- Antisymmetric)
# ${6} = Mmax (default = 126) Rhomboidal Truncations (Rhomb=Mmax-1)
#
# Suitable Cases to Analyse the Horizontal Structure of the Selected Modes
#
cs=${1}
if [ "${cs}" = "" ]
then
  cs="br"
fi
mmx=${2}
if [ "${mmx}" = "" ]
then
  mmx=126
else
  mmx=${2}
fi
#
let lrb1=${mmx}+1
let lrb2=${mmx}+4
let lkv=${mmx}/2+1
let lmx=${mmx}+1
let lge=${mmx}/2+2
#
set -x
#
./Hough_Functions_uvh ${cs} rb 2 ${lrb1} 1 ${mmx}
sleep 3s
./Hough_Functions_uvh ${cs} rb 2 ${lrb2} 2 ${mmx}
sleep 3s
./Hough_Functions_uvh ${cs} rb 7 ${lrb1} 1 ${mmx}
sleep 3s
./Hough_Functions_uvh ${cs} rb 7 ${lrb2} 2 ${mmx}
sleep 3s
./Hough_Functions_uvh ${cs} kv 2 ${lkv} 1 ${mmx}
sleep 3s
./Hough_Functions_uvh ${cs} kv 7 ${lkv} 1 ${mmx}
sleep 3s
./Hough_Functions_uvh ${cs} mx 2 ${lmx} 2 ${mmx}
sleep 3s
./Hough_Functions_uvh ${cs} mx 7 ${lmx} 2 ${mmx}
sleep 3s
./Hough_Functions_uvh ${cs} gw 2 1 1 ${mmx}
sleep 3s
./Hough_Functions_uvh ${cs} gw 7 1 1 ${mmx}
sleep 3s
./Hough_Functions_uvh ${cs} ge 2 ${lge} 1 ${mmx}
sleep 3s
./Hough_Functions_uvh ${cs} ge 7 ${lge} 1 ${mmx}
#
