Subversion Repositories lagranto.wrf

#!/bin/csh

# Set Lagranto 
set LAGRANTO = ${LAGRANTOBASE}.${MODEL}/

# Set some parameters
set crifile=${PWD}/select.parsed

# Set base directories (run+prog)
set cdfdir=${PWD}
set tradir=${PWD}

# Write usage information
if ( ${#argv} ==  0 ) then
  echo 
  ${LAGRANTO}/bin/lagrantohelp select short
  echo  
  exit 0
endif

# Write list of special selection criteria
if ( "$1" == "-special" ) then
    grep '%)' ${LAGRANTO}/select/special.f
    exit 0
endif

# Write title
echo 
echo '========================================================='
echo '       *** START OF PREPROCESSOR SELECT ***              '
echo

# Save input arguments
set inpfile=$1
set outfile=$2
set crit="$3"
shift
shift
shift

# Handle optional arguments
echo
echo '---- OPTIONAL FLAGS -------------------------------------'
echo  

set noclean = 'false'
set format  = 'trajectory'
set regionf = 'regionf'
set trigger = 'nil'

while ( $#argv > 0 )

  switch ( $argv[1] )

   case -noclean
     set noclean = 'true'
     echo "noclean       -> true (user defined)"
   breaksw

   case -trigger
     set trigger = '-trigger'
     echo "trigger      ->  true (user defined)"
   breaksw

   case -boolean
     set format = 'boolean'
     echo "format       ->  boolean (user defined)"
   breaksw

   case -index
     set format = 'index'
     echo "format       ->  index (user defined)"
   breaksw

   case -count
     set format = 'count'
     echo "format       ->  count (user defined)"
   breaksw

   case -startf
     set format = 'startf'
     echo "format       ->  startf (user defined)"
   breaksw

   case -regionf
     set regionf = $argv[2]
     echo "regionf                -> ${regionf} (user defined)"
     shift;
   breaksw

  endsw
 
  shift;

end

# Decide whether <select> is a file or an explicit criterion
set flag_select = 'criterion'
set test = `echo ${crit} | grep ':' | wc -c`
if ( "${test}" == "0" ) then 
  set flag_select     = 'file'
  set flag_selectfile = $crit
  if ( -f ${flag_selectfile} ) then 
     set crit             = `cat ${flag_selectfile}`
  else
     echo " ERROR: criterion file ${flag_selectfile} is missing... Stop"
     exit 1
  endif
endif

# Get the start, end and reference date for the tracing
set ntra      =  `${LAGRANTO}/goodies/trainfo.sh ${inpfile} ntra`
set ntim      =  `${LAGRANTO}/goodies/trainfo.sh ${inpfile} ntim`
set ncol      =  `${LAGRANTO}/goodies/trainfo.sh ${inpfile} ncol`
set times     =  `${LAGRANTO}/goodies/trainfo.sh ${inpfile} times`
set vars      =  `${LAGRANTO}/goodies/trainfo.sh ${inpfile} vars`

# Split the criterion into subunits (with Perl)
${LAGRANTO}/select/select.perl "${crit}" >! ${crifile} 
if ( "${status}" != "0" ) then
     echo "ERROR:  Parser <select> failed"
     exit 1
endif

# Write some status information
echo
echo '---- DIRECTORIES AND PROGRAMS ---------------------------'
echo    
echo "CDF directory         : ${cdfdir}"
echo "TRA directory         : ${tradir}"
echo "PROGRAM SELECT        : ${LAGRANTO}/select/select"
echo "PARSER                : ${LAGRANTO}/select/select.perl"
echo "CRITERION FILE        : ${crifile}"
echo
echo '---- INPUT PARAMETERS -----------------------------------'
echo    
echo "Input file            : ${inpfile}"
echo "Output file           : ${outfile}"
echo "Output format         : ${format}"
if ( "${flag_select}" == "criterion" ) then
    echo "Criterion             : ${crit}"
else
    echo "Criterion             : ${crit} (from file ${flag_selectfile})"
endif
echo
echo '---- INPUT FILE -----------------------------------------'
echo    
echo "# TRA                 : ${ntra}"
echo "# TIMES               : ${ntim}"
echo "# COLUMNS             : ${ncol}"
echo "Times                 : ${times}"
echo "Variables             : ${vars}"
echo
echo '---- PARSED CRITERION ------------------------------------'
echo 
cat ${crifile}

# Finish the preprocessor
echo 
echo '       *** END OF PREPROCESSOR SELECT ***              '
echo '========================================================='
echo

# Run the selection programme
echo \"${inpfile}\" >! select.param
echo \"${outfile}\" >> select.param
echo \"${format}\"  >> select.param
echo \"${crifile}\" >> select.param
echo ${ntra}        >> select.param
echo ${ntim}        >> select.param
echo ${ncol}        >> select.param
echo \"${regionf}\" >> select.param
echo \"${trigger}\" >> select.param

${LAGRANTO}/select/select

if ( "${status}" != "0" ) then
  echo "ERROR:  Program <select> failed"
  exit 1
endif

# Make clean
if ( "${noclean}" == "false" ) then
  \rm -f select.param
  \rm -f ${crifile}
endif

exit 0