WebSVN – lagranto.um – Diff – /trunk/bin/lagranto


#!/bin/csh
 
# ---------------------------------------------------------------------
# Usage
# ---------------------------------------------------------------------
 
# Write usage information
if ( (${#argv} == 0) ) then
  echo 
  echo "-----------------------------------------------------------------------------------"
  echo "Usage: "
  echo 
  echo "  1. lagranto {caseid[.label]} startdate enddate {startf} {select} [flags] "
  echo
  echo "         caseid     : case     (name of case directories; label appears in output directoy)"
  echo "                      local    (everything is run in local directory)"
  echo "                      interim  (take wind fields from Era-Interim)"
  echo "                      analysis (take wind fields from ECMWF operational analysis)"
  echo "                      forecast (take wind fields from ECMWF operational forecast)" 
  echo
  echo "         startdate  : start date in format yyyymmdd_hh(mm)"
  echo 
  echo "         enddate    : end date in format yyyymmdd_hh(mm)"
  echo  
  echo "         startf     : [ file | criterion ] with starting coordinates (lon,lat,p)"
  echo
  echo "         select     : [ file | criterion ] for trajectory selection (set 'nil' for none)"
  echo
  echo "         flags      : -o     outfile   (name of output trajectory file) "
  echo "                      -j               (jumping flag for crossing lower boundary)"
  echo "                      -v     filename  (name of tracing file)"
  echo "                      -r     filename  (name of region file)"
  echo "                      -noclean         (do not clean trajectory directory after run)"
  echo "                      -prep            (prepare run script, but do not run lagranto)"
  echo "                      -log             (logfile is written to screen instead of logfile)"
  echo "                      -changet         (change times on netCDF files realtive to startdate)"
  echo
  echo "  2. lagranto [ -open | -remove ] caseid[.label] "
  echo
  echo "        -open   : change to the run directory [open a new Xterm window]"
  echo "        -remove : remove the whole run directory"
  echo
  echo "Examples:"
  echo 
  echo "-----------------------------------------------------------------------------------"
  exit 0
endif
 
# Check whether LAGRANTO was set with setenv
set test = `setenv | grep LAGRANTO` 
if ( "${test}" == "" ) then
   echo "---------------------------------------------------------------------"
   echo "Please set environmental variable <LAGRANTO> to your source directory"
   echo 'csh > setenv LAGRANTO ${TO_YOUR_LAGRANTO_INSTALLATION}'
   echo "---------------------------------------------------------------------"
   exit 1
endif
  
# ---------------------------------------------------------------------
# Set some general parameters
# ---------------------------------------------------------------------
 
# Set the prefix of the primary and secondary data files
set charp = 'P'
set chars = 'S'
 
# Set the name of the run script
set runscript = "runscript"
 
# Remember the calling directory
set calldir = $PWD
 
# ---------------------------------------------------------------------
# Parameter handling
# ---------------------------------------------------------------------
 
# Write title
echo
echo '========================================================='
echo '             *** START OF LAGRANTO ***             '
echo
 
# Save the complete argument list
set lagranto_call = "$*"
 
echo '---- HANDLE ARGUMENTS -----------------------------------'
echo
 
# Get fixed arguments
if (${#argv} == 1) then
  set mode = '-open'
  set caseid    = "$1"
  shift
 
else if (${#argv} == 2) then
  set mode      = "$1"
  set caseid    = "$2"
  shift
  shift
 
else
  set caseid    = "$1"
  set startdate = "$2"
  set enddate   = "$3"
  set startf    = "$4"
  set select    = "$5"
  set mode      = "runlagranto"
  shift
  shift
  shift
  shift
  shift
 
endif
 
# Split <caseid> into base and appendix 
set casebase = `(echo $caseid | awk -F "." '{print $1}')`   
set appendix = `(echo $caseid | awk -F "." '{print $2}')`   
if ( "${appendix}" != "" ) then
  set appendix = ".${appendix}"
endif
 
# Set the directories (depending on <caseid>)
if ( "${casebase}" == "local" ) then
  set cdfdir = $PWD
  set tradir = $PWD
 
else if ( "${casebase}" == "interim" ) then
  set cdfdir = /net/dansgaard/atmosdyn/erainterim/cdf/
  set tradir = $PWD
 
else if ( "${casebase}" == "analysis" ) then
  set cdfdir = 'TO_BE_SPECIFIED'
  set tradir = $PWD
  echo "ERROR: caseid (analysis) not yet implemented"
  exit 1  
 
else if ( "${casebase}" == "forecast" ) then
  set cdfdir = 'TO_BE_SPECIFIED'
  set tradir = $PWD
  echo "ERROR: caseid (forecast) not yet implemented"
  exit 1  
 
else
  set cdfdir = ${HOME}/cdf/${casebase}
  set tradir = ${HOME}/tra/${casebase}
endif
 
# Change to the run directory if requested
if ( "${mode}" == "-open" ) then
   set count = 0   
   foreach folder ( ${tradir}/ntr_*_${caseid}_* )
     @ count = ${count} + 1
     echo "[${count}] ${folder}"
   end
   echo
   echo "Please select directory [0...${count}]"
   set id = $<
   set count = 0   
   foreach folder ( ${tradir}/ntr_*_${caseid}_* )
     @ count = ${count} + 1
     if ( "${count}" == "${id}" ) then 
        cd ${folder}
        xterm &
     endif
   end
   exit 0
endif
 
# Remove the whole run directory if requested
if ( "${mode}" == "-remove" ) then
   set count = 0   
   foreach folder ( ${tradir}/ntr_*_${caseid}_* )
     @ count = ${count} + 1
     echo "[${count}] ${folder}"
   end
   echo
   echo "Please select directory [0...${count}, all]"
   set id = $<
   set count = 0   
   foreach folder ( ${tradir}/ntr_*_${caseid}_* )
     @ count = ${count} + 1
     if ( ( "${count}" == "${id}" ) | ( "${id}" == "all" ) ) then 
        \rm -r ${folder}
     endif
   end
   exit 0
endif
 
# Show the trajectory file
if ( "${mode}" == "-show" ) then
   set count = 0   
   foreach folder ( ${tradir}/ntr_*_${caseid}_* )
     @ count = ${count} + 1
     echo "[${count}] ${folder}"
   end
   echo
   echo "Please select directory [0...${count}]"
   set id = $<
   set count = 0   
   foreach folder ( ${tradir}/ntr_*_${caseid}_* )
     @ count = ${count} + 1
     if ( "${count}" == "${id}" ) then 
        cd ${folder}
        set filename = ` sed -ne '9,9p' runscript.sh | cut -c 2-` 
	${LAGRANTO}/goodies/trainfo.sh ${filename} list | less
     endif
   end
   exit 0
endif
 
# Decide whether <select> is a file or an explicit criterion
set flag_select = 'file'
set test = `echo ${select} | grep ':' | wc -c`
if ( "${test}" != "0" ) then 
  set flag_select = 'criterion'
endif
 
# Decide whether <startf> is a file or an explicit criterion
set flag_startf = 'file'
 
set test = `echo ${startf} | grep ' ' | wc -c`
if ( "${test}" != "0" ) then 
  set flag_startf = 'point'
endif
 
set test = `echo ${startf} | grep '@' | wc -c`
if ( "${test}" != "0" ) then 
  set flag_startf = 'criterion'
endif
 
 
# Write some status information
echo " Case ID       : ${caseid}"
echo " Start date    : ${startdate}"
echo " End date      : ${enddate}"
 
if ( "${flag_startf}" == "criterion" ) then
    echo " Start file    : ${startf} [criterion] -> startf.criterion [file]"
else if ( "${flag_startf}" == "file" ) then
    echo " Start points  : ${startf} [file]"
else if ( "${flag_startf}" == "point" ) then
    echo " Start points  : ${startf} [point] -> startf [file]"
endif
 
if ( "${flag_select}" == "criterion" ) then
    echo " Selection     : ${select} [criterion] -> selectf.criterion [file]"
else
    echo " Selection     : ${select} [file]"
endif    
 
echo
 
# Set defaults for flags
set jumpflag   = ""
set tracefile  = "tracevars"
set outfile    = "lsl_${startdate}"
set keepflag   = ""
set prepflag   = ""
set logfile    = "${runscript}.logfile"
set regionf    = "regionf"
set changet    = "false"
 
# Handle flags
while ( $#argv > 0 )
 
  switch ( $argv[1] )
 
  echo $argv[1]
 
   case -j
     set jumpflag = "-j"
     echo " Flag '-j [jump]'             -> True "
     echo
   breaksw
 
   case -prep
     set prepflag = "prep"
     echo " Flag '-prep [prepare only]'  -> ${prepflag} "
     echo
   breaksw
 
   case -noclean
    set keepflag = "keep"
     echo " Flag '-noclean [no cleaning]' -> ${keepflag} "
     echo
   breaksw
 
   case -o
     set outfile =  $argv[2]
     echo " Flag '-o [output]'            -> ${outfile} "
     echo
   breaksw
 
   case -v
     set tracefile = $argv[2]
     echo " Flag '-v [tracevars]'         -> ${tracefile} "
     echo
     shift
   breaksw
 
   case -r
     set regionfile = $argv[2]
     echo " Flag '-r [regionfile]'        -> ${regionfile} "
     echo
     shift
   breaksw
 
   case -log
     set logfile = "/dev/stdout"
     echo " Flag '-v [log]'               -> ${logfile} "
   breaksw
 
   case -changet
     set changet = "true"
     echo " Flag '-changet [change times]'-> True "
     echo
   breaksw
 
  endsw
 
  shift;
 
end
 
# Decide whether a tracevars file is given or not
if ( -f ${tradir}/${tracefile} ) then
   set flag_tracevars = 'file'
else
   set flag_tracevars = 'nil'
endif
 
echo
if ( "${flag_tracevars}" == "file" ) then
    echo " Trace Vars    : ${tracefile} [file]"
else
    echo " Trace Vars    : [no fields are traced]"
endif    
echo
 
# Set the name of the run directory (forward/backward distinction)
set dir   = 'f'
set range = `${LAGRANTO}/bin/gettidiff.sh ${enddate} ${startdate}`
if ( ${range} < 0 ) then
  set dir   = 'b'
  set range = `${LAGRANTO}/bin/gettidiff.sh ${startdate} ${enddate}`
endif
set rundir = ${tradir}/ntr_${startdate}_${dir}${range}_${caseid}
 
if ( "${flag_startf}" == "criterion" ) then
  set rundir = "${rundir}_startf"
else if ( "${flag_startf}" == "point" ) then
  set rundir = "${rundir}_startf"
else if ( "${flag_startf}" == "file" ) then
  set rundir = "${rundir}_${startf}"
endif
 
if (  "${flag_select}" != "file" ) then
  set rundir = "${rundir}_selectf"
else
 set rundir = "${rundir}_${select}"
endif
 
# Write directory names
echo '---- DIRECTORIES ----------------------------------------'
echo
echo " CDF = ${cdfdir}"
echo " TRA = ${tradir}"
echo " RUN = ${rundir}"
 
# Check whether CDF and TRA directories exist
if ( ! -d ${cdfdir} ) then
    echo 'Cannot find CDF directory ... Stop'
    exit 1
endif
if ( ! -d ${tradir} ) then
    echo 'Cannot find TRA directory ... Stop'
    exit 1
endif
 
# Create the run directory
if ( ! -d ${rundir} ) then
   mkdir ${rundir}
endif 
if ( ! -d ${rundir} ) then
    echo 'Cannot create RUN directory ... Stop'
    exit 1
endif 
 
# ---------------------------------------------------------------------
# Handle the time specifier - startdate, enddate
# ---------------------------------------------------------------------
 
echo
echo '---- TIME RANGE -----------------------------------------'
echo
 
# Check format of start and end date - must be the same
set ns=`echo $startdate | sed -e 's/_[0-9]*//' | wc -c`
set ne=`echo $enddate   | sed -e 's/_[0-9]*//' | wc -c`
if ( $ns != $ne ) then
  echo " ERROR: start and end date must be in the same format ***"
  exit 1
endif
if ( $ns != 9 ) then
  echo " ERROR: Date format must be yyyymmdd ***"
  exit 1
endif
set ns=`echo $startdate | sed -e 's/[0-9]*_//' | wc -c`
set ne=`echo $enddate   | sed -e 's/[0-9]*_//' | wc -c`
if ( $ns != $ne ) then
  echo " ERROR: start and end date must be in the same format ***"
  exit 1
endif
if ( ( $ns != 5 ) & ( $ns != 3 ) ) then
  echo " ERROR: Time format must be hh(mm) ***"
  exit 1
endif
 
# Split the start and end date into <yymmdd_hh and mm>
set startdate_ymdh = `echo $startdate | cut -c 1-11`
set startdate_min  = `echo $startdate | cut -c 12-13`
if ( $startdate_min == "" ) set startdate_min = 00
 
set enddate_ymdh = `echo $enddate | cut -c 1-11`
set enddate_min  = `echo $enddate | cut -c 12-13`
if ( $enddate_min == "" ) set enddate_min = 00
 
# Get the time difference between <start_ymdh> and <end_ymdh> date
# Decide whether trajectoriesare forward or backward
set timediff_hh = `${LAGRANTO}/goodies/gettidiff ${enddate_ymdh} ${startdate_ymdh}`
 
if ( $timediff_hh == 0 ) then
  if ( $enddate_min > $startdate_min ) then
    set direction = f
    set idir      = 1
  else
    set direction = b
    set idir      = -1
  endif
else if ( $timediff_hh > 0 ) then
  set direction = f
  set idir      = 1
else
  set direction = b
  set idir      = -1
  @ timediff_hh = $idir * $timediff_hh
endif
 
# Get also minutes for time difference, if <start_min> or <end_min> != 0
set timediff_mm=
 
if ( $startdate_min != 00 || $enddate_min != 00 ) then
  @ min = ( $enddate_min - $startdate_min )
  if ( $min == 0 ) then
    set timediff_mm=
  else if ( $min > 0 ) then
    if ( $idir == 1 ) then
      set timediff_mm=$min
    else
      @ timediff_hh --
      @ timediff_mm = 60 - $min
    endif
  else
    if ( $idir == 1 ) then
      @ timediff_hh --
      @ timediff_mm = 60 + $min
    else
      @ timediff_mm = 0 - $min
    endif
  endif
endif
 
# Set the reference date equal to the startdate
set refdate=${startdate}
 
# Write status information
echo "Time range      : ${startdate} -> ${enddate}"
if ( ${timediff_mm} !=  "" ) then
   echo "Time difference : ${timediff_hh} h ${timediff_mm} min"
else
   echo "Time difference : ${timediff_hh} h"
endif
echo "Direction       : ${direction} (${idir})"
echo "Reference date  : ${refdate}"
 
# ---------------------------------------------------------------------
# Decide which files are needed and check availability
# ---------------------------------------------------------------------
 
echo
echo '---- INPUT FILES ----------------------------------------'
echo
 
# Change to input data (cdf) directory
cd ${cdfdir}
echo
 
# Take the time increment from flag list ('nil', if not defined)
set timeinc = 'nil'
 
# Find a first data file (if possible corresponding to start/end date
# If starttime is not a data time, take the first file in the direectory
if ( $direction == "f" ) then
  set file=${charp}${startdate_ymdh}
else
  set file=${charp}${enddate_ymdh}
endif
if ( ! -f $file ) then
  set file=`ls ${charp}[0-9_]*[0-9] | head -1 | sed -e 's/@//'`
endif
 
# Determine timeinc (the time difference in hours between two data file)
# if not already defined with option -i
if ( ${timeinc} == "nil" ) then
  set date1=`echo $file | cut -c 2-12`
  set n=`ls ${charp}[0-9_]*[0-9] | grep -n $date1 | awk -F: '{print $1}'`
  @ n ++
  set date2=`ls ${charp}[0-9_]*[0-9] | head -$n | tail -1 | cut -c 2-12`
  set timeinc=`${LAGRANTO}/goodies/gettidiff $date2 $date1`
endif
if ( $timeinc == 0 ) then
    echo " ERROR: cannot set the time increment between input files ***"
    exit 1
endif
 
# Search the first file to use
set flag=0
set td=
foreach i ( ${charp}????????_?? )
 
  set date = `echo $i | cut -c 2-12`
  set td1  = `${LAGRANTO}/goodies/gettidiff ${startdate_ymdh} ${date}`
  set td2  = `${LAGRANTO}/goodies/gettidiff ${enddate_ymdh}   ${date}`
 
  if (( $td1 < $timeinc || $td2 < $timeinc ) && ( $td1 >= 0 || $td2 >= 0 )) then
      set datfiles=$date
      if ( $td1 < $timeinc    ) set td=$td1
      if ( $td2 < $timeinc    ) set td=$td2
      if ( ( $startdate_min > 0 ) || ( $enddate_min > 0 ) ) @ td ++ 
      goto label2
  endif
 
end
 
# if no P/T-files are available for the specified time period, then $td is
# still undefined
if ( $td == "" ) then
  echo " ERROR: no data files available for the specified time period"
  exit 1
endif
 
# Everything is fine so far: proceed
label2:
 
# Calculate the number of further files
@ num = ( $timediff_hh + $td ) / $timeinc + 1
@ dum1 = ( $num - 1 ) * $timeinc
@ dum2 = $timediff_hh + $td
if ( $dum1 != $dum2 ) @ num ++
 
# Get a list of all needed files
set numfiles=$num
set sfiles=1
while ( $num > 1 )
  set date=`${LAGRANTO}/goodies/newtime $date $timeinc`
  if ( ! -f ${charp}${date} ) then
    echo " ERROR: file with primary data is missing for $date"
    exit 1
  else if ( ! -f ${chars}${date} ) then
    set sfiles=0
    set datfiles=`echo $datfiles $date`
  else
    set datfiles=`echo $datfiles $date`
  endif
  @ num --
end
 
# Write some status information
echo "Primary file prefix            : ${charp}"
echo "Secondary file prefix          : ${chars}"
echo "Time increment for input files : ${timeinc}"
echo "# input files                  : ${numfiles}"
echo "${charp} files availability           : 1"  
echo "${chars} files availability           : ${sfiles}"     
echo "First input file               : $datfiles[1] " 
echo "Last input file                : $datfiles[$numfiles] " 
 
 
# ---------------------------------------------------------------------
# Prepare run directory - copy all needed files
# ---------------------------------------------------------------------
 
echo
echo '---- PREPARE RUN DIRECTORY ------------------------------'
echo
 
# Change to run directory
cd ${rundir}
echo
 
 
# Link cdf input files
foreach date ( $datfiles )
  if ( -f ${cdfdir}/${charp}${date} ) then
     if ( ! -f  ${charp}${date} ) then
        ln -s  ${cdfdir}/${charp}${date}  ${charp}${date}
        echo "Linking  ${date} "
     else
        echo "Existing ${date} "
     endif
  endif
  if ( -f ${cdfdir}/${chars}${date} ) then
     if ( ! -f  ${chars}${date} ) then
        ln -sf  ${cdfdir}/${chars}${date}  ${chars}${date}
     endif
  endif
end
 
# Change times on netCDF files relative to reference date
if ( "${changet}" == "true" ) then
 echo
 foreach i ( $datfiles )
   ${LAGRANTO}/goodies/changet.sh ${refdate} ${charp}${i}
 end
 if ( ${sfiles} == 1 ) then
   echo
   foreach i ( $datfiles )
     ${LAGRANTO}/goodies/changet.sh ${refdate} ${chars}${i}
   end
 endif
endif
 
# Get the constants file for each file
set first = 1
foreach date ( $datfiles )
  set cstname = `ncdump -h ${charp}${date} | grep constants_file_name`
  set cstname = `echo \'${cstname}\' |   awk -F\" '{print $2}'`
  if ( -f ${cdfdir}/${cstname} ) then
    if ( ! -f ${cstname} ) then
       ln -s  ${cdfdir}/${cstname}  ${cstname}
       echo "Linking  constants file ${cstname}"
       set first = 0
    else if ( "${first}" == "1" ) then
       echo "Existing constants file ${cstname}"
       set first = 0
    endif
  endif
end
 
# Copy the tracevars file
if ( "${flag_tracevars}" == "file" ) then
   if ( -f ${tradir}/${tracefile} ) then
      \cp ${tradir}/${tracefile} .
   endif
endif
 
# Copy the start file/criterion
if ( "${flag_startf}" == "criterion" ) then
   \rm -f startf.criterion
   echo ${startf} > startf.criterion
   set startf =  "startf.criterion"
else if ( "${flag_startf}" == "file" ) then
   if ( -f ${tradir}/${startf} ) then
     \cp ${tradir}/${startf} .
   endif
else if ( "${flag_startf}" == "point" ) then
   \rm -f startf
   echo ${startf}  > startf
   set startf      =  "startf"
   set flag_startf = "file"
endif 
 
# Copy the select file/criterion
if ( "${flag_select}" == "criterion" ) then
   \rm -f select.criterion
   echo ${select} > select.criterion
   set select = "select.criterion"
else
   if ( -f ${tradir}/${select} ) then
     \cp ${tradir}/${select} .
   endif
endif 
 
 
# ---------------------------------------------------------------------
# Prepare the run script
# ---------------------------------------------------------------------
 
echo
echo '---- PREPARE RUN SCRIPT ---------------------------------'
echo
 
# Remove runscript if it already exists
if ( -f ${runscript}.sh ) then
   \rm -f ${runscript}.sh
endif 
 
# Write header
echo '#\!/bin/csh'                    >! ${runscript}.sh
echo "#"                              >> ${runscript}.sh
echo "#----- Calling command"         >> ${runscript}.sh
echo "#"                              >> ${runscript}.sh
echo "# lagranto ${lagranto_call} "   >> ${runscript}.sh
echo "#"                              >> ${runscript}.sh 
echo "#----- Output file    "         >> ${runscript}.sh
echo "#"                              >> ${runscript}.sh
echo "# ${outfile}                "   >> ${runscript}.sh
echo "#"                              >> ${runscript}.sh 
 
# Decide whether <create_startf> has to be run
if ( "${flag_startf}" == "criterion" ) then
  echo "#------ Run <startf>"                                                                                       >> ${runscript}.sh
  echo "#"                                                                                                                 >> ${runscript}.sh
  echo ${LAGRANTO}/bin/create.sh ${startdate} startf \"`cat startf.criterion`\" -t ${tracefile} -regionf ${regionf} >> ${runscript}.sh 
  echo "#"                                                                                                                 >> ${runscript}.sh
  set startf = "startf"
endif
 
# Check whether startfile is ready
cat >> ${runscript}.sh <<EOF 
#------ Abort if no startf is available
#
if ( ! -f ${startf} ) then
  echo " ERROR: no start file available .... Stop"
  exit 1
endif
EOF
 
# Decide wich intermediate format is used (needed for efficient communication between different programs)
set iform1 = ""
set iform2 = ""
if ( "${tracefile}" != ""    ) set iform1 = ".4"
if ( "${select}"    != "nil" ) set iform2 = ".4"
 
# Remove all trajectory files which might already be there
cat >> ${runscript}.sh <<EOF 
#
#------ Remove existing trajectory files
#
if ( -f ${outfile}${iform1} ) then
  \rm -f ${outfile}${iform1}
endif
if ( -f ${outfile}${iform2} ) then
  \rm -f ${outfile}${iform2}
endif
if ( -f ${outfile} ) then
  \rm -f ${outfile}
endif
 
EOF
 
# Write call to <caltra> into runscript
echo "#"                                                                                      >> ${runscript}.sh
echo "#------ Run <caltra>"                                                                   >> ${runscript}.sh
echo "#"                                                                                      >> ${runscript}.sh
echo ${LAGRANTO}/bin/caltra.sh ${startdate} ${enddate} startf ${outfile}${iform1} ${jumpflag} >> ${runscript}.sh  
 
# Check whether caltra was successful
cat >> ${runscript}.sh <<EOF 
#
#------ Abort if caltra was not successful
#
if ( ! -f ${outfile}${iform1} ) then
  echo " ERROR: caltra failed .... Stop"
  exit 1
endif
EOF
 
#  Decide whether <trace> has to be run
if ( "${tracefile}" != "" ) then
  echo "#"                                                                                    >> ${runscript}.sh
  echo "#------ Run <trace>"                                                                  >> ${runscript}.sh
  echo "#"                                                                                    >> ${runscript}.sh
  echo ${LAGRANTO}/bin/trace.sh ${outfile}${iform1} ${outfile}${iform2} -v ${tracefile}       >> ${runscript}.sh 
 
endif
 
# Check whether trace was successful
cat >> ${runscript}.sh <<EOF 
#
#------ Abort if trace was not successful
#
if ( ! -f ${outfile}${iform2} ) then
  echo " ERROR: trace failed .... Stop"
  exit 1
endif
EOF
 
#  Decide whether <select> has to be run
if ( "${select}" != "nil" ) then                        
  echo "#"                                                                                               >> ${runscript}.sh
  echo "#------ Run <select>"                                                                            >> ${runscript}.sh
  echo "#"                                                                                               >> ${runscript}.sh
  echo ${LAGRANTO}/bin/select.sh ${outfile}${iform2} ${outfile}  \"`cat ${select}`\" -regionf ${regionf} >> ${runscript}.sh 
  echo "#"                                                                                               >> ${runscript}.sh
endif
 
# Show runscript on screen
more ${runscript}.sh
 
# Change runscript to executable
chmod u+x ${runscript}.sh
 
# Decide whether to stop at this place (only preparation for Lagranto run)
if ( "${prepflag}" == "prep" ) then
  echo
  echo '---- READY TO RUN ----- ---------------------------------'
  echo
  echo " Would you like to change to the rundirectory [y/n] "
 
  set id = $<
  if ( ( "${id}" == "y" ) | ( "${id}" == "Y" ) ) then
     cd ${rundir}
     xterm &
  endif
 
  exit 0
 
endif  
 
# Define environment variable for the case
echo
echo '---- REMEMBER MOST RECENT CALL --------------------------'
echo
cd ${calldir}
if ( -f lagranto.param ) \rm lagranto.param
echo "${rundir}"  >! lagranto.param
echo "${outfile}" >> lagranto.param
cd ${rundir}
echo
echo " run directory + output file -> lagranto.param"
 
# ---------------------------------------------------------------------
# Run lagranto
# ---------------------------------------------------------------------
 
echo
echo '---- RUN LAGRANTO ---------------------------------------'
echo
 
if ( "${logfile}" != "/dev/stdout" ) then
    echo "log goes to ${logfile}"
endif
 
# Remove logfile
if ( "${logfile}" != "/dev/stdout" ) then
    \rm -f    ${logfile}
endif
 
# Run it
${rundir}/${runscript}.sh > ${logfile}
 
# Check for errors in logfile
if ( "${logfile}" != "/dev/stdout" ) then
   grep ERROR ${logfile} 
endif
 
# ---------------------------------------------------------------------
# Make clean and finish lagranto run
# ---------------------------------------------------------------------
 
# Remove unneeded files
\rm -f fort.9
 
# Remove files if requested
if ( "${keepflag}" != "keep" ) then
 
   echo 
   echo '---- MAKING CLEAN ---------------------------------------'
   echo
 
   cd ${rundir}
 
       
   set runscript = "${runscript}.sh"
 
   foreach file ( `\ls -L1` )
    
      set id = "0"
 
      if ( "${file}" == "${outfile}"   ) set id="1" 
      if ( "${file}" == "${runscript}" ) set id="1" 
      if ( "${file}" == "${logfile}"   ) set id="1" 
  
      if ( "${id}" == "0" ) \rm -f ${file}      
 
   end
  
endif
 
# Everything is fine!
echo
echo '             *** END OF LAGRANTO ***             '
echo '========================================================='
echo
 
exit 0 
   
   

Rev 3	Rev 7
1	`#!/bin/csh`	1	`#!/bin/csh`
2		2
3	`# ---------------------------------------------------------------------`	3	`# ---------------------------------------------------------------------`
4	`# Usage`	4	`# Usage`
5	`# ---------------------------------------------------------------------`	5	`# ---------------------------------------------------------------------`
6		6
7	`# Write usage information`	7	`# Write usage information`
8	`if ( (${#argv} == 0) ) then`	8	`if ( (${#argv} == 0) ) then`
9	`echo`	9	`echo`
10	`echo "-----------------------------------------------------------------------------------"`	10	`echo "-----------------------------------------------------------------------------------"`
11	`echo "Usage: "`	11	`echo "Usage: "`
12	`echo`	12	`echo`
13	`echo " 1. lagranto {caseid[.label]} startdate enddate {startf} {select} [flags] "`	13	`echo " 1. lagranto {caseid[.label]} startdate enddate {startf} {select} [flags] "`
14	`echo`	14	`echo`
15	`echo " caseid : case (name of case directories; label appears in output directoy)"`	15	`echo " caseid : case (name of case directories; label appears in output directoy)"`
16	`echo " local (everything is run in local directory)"`	16	`echo " local (everything is run in local directory)"`
17	`echo " interim (take wind fields from Era-Interim)"`	17	`echo " interim (take wind fields from Era-Interim)"`
18	`echo " analysis (take wind fields from ECMWF operational analysis)"`	18	`echo " analysis (take wind fields from ECMWF operational analysis)"`
19	`echo " forecast (take wind fields from ECMWF operational forecast)"`	19	`echo " forecast (take wind fields from ECMWF operational forecast)"`
20	`echo`	20	`echo`
21	`echo " startdate : start date in format yyyymmdd_hh(mm)"`	21	`echo " startdate : start date in format yyyymmdd_hh(mm)"`
22	`echo`	22	`echo`
23	`echo " enddate : end date in format yyyymmdd_hh(mm)"`	23	`echo " enddate : end date in format yyyymmdd_hh(mm)"`
24	`echo`	24	`echo`
25	`echo " startf : [ file \| criterion ] with starting coordinates (lon,lat,p)"`	25	`echo " startf : [ file \| criterion ] with starting coordinates (lon,lat,p)"`
26	`echo`	26	`echo`
27	`echo " select : [ file \| criterion ] for trajectory selection (set 'nil' for none)"`	27	`echo " select : [ file \| criterion ] for trajectory selection (set 'nil' for none)"`
28	`echo`	28	`echo`
29	`echo " flags : -o outfile (name of output trajectory file) "`	29	`echo " flags : -o outfile (name of output trajectory file) "`
30	`echo " -j (jumping flag for crossing lower boundary)"`	30	`echo " -j (jumping flag for crossing lower boundary)"`
31	`echo " -v filename (name of tracing file)"`	31	`echo " -v filename (name of tracing file)"`
32	`echo " -r filename (name of region file)"`	32	`echo " -r filename (name of region file)"`
33	`echo " -noclean (do not clean trajectory directory after run)"`	33	`echo " -noclean (do not clean trajectory directory after run)"`
34	`echo " -prep (prepare run script, but do not run lagranto)"`	34	`echo " -prep (prepare run script, but do not run lagranto)"`
35	`echo " -log (logfile is written to screen instead of logfile)"`	35	`echo " -log (logfile is written to screen instead of logfile)"`
36	`echo " -changet (change times on netCDF files realtive to startdate)"`	36	`echo " -changet (change times on netCDF files realtive to startdate)"`
37	`echo`	37	`echo`
38	`echo " 2. lagranto [ -open \| -remove ] caseid[.label] "`	38	`echo " 2. lagranto [ -open \| -remove ] caseid[.label] "`
39	`echo`	39	`echo`
40	`echo " -open : change to the run directory [open a new Xterm window]"`	40	`echo " -open : change to the run directory [open a new Xterm window]"`
41	`echo " -remove : remove the whole run directory"`	41	`echo " -remove : remove the whole run directory"`
42	`echo`	42	`echo`
43	`echo "Examples:"`	43	`echo "Examples:"`
44	`echo`	44	`echo`
45	`echo "-----------------------------------------------------------------------------------"`	45	`echo "-----------------------------------------------------------------------------------"`
46	`exit 0`	46	`exit 0`
47	`endif`	47	`endif`
48		48
49	`# Check whether LAGRANTO was set with setenv`	49	`# Check whether LAGRANTO was set with setenv`
50	set test = `setenv \| grep LAGRANTO`	50	set test = `setenv \| grep LAGRANTO`
51	`if ( "${test}" == "" ) then`	51	`if ( "${test}" == "" ) then`
52	`echo "---------------------------------------------------------------------"`	52	`echo "---------------------------------------------------------------------"`
53	`echo "Please set environmental variable <LAGRANTO> to your source directory"`	53	`echo "Please set environmental variable <LAGRANTO> to your source directory"`
54	`echo 'csh > setenv LAGRANTO ${TO_YOUR_LAGRANTO_INSTALLATION}'`	54	`echo 'csh > setenv LAGRANTO ${TO_YOUR_LAGRANTO_INSTALLATION}'`
55	`echo "---------------------------------------------------------------------"`	55	`echo "---------------------------------------------------------------------"`
56	`exit 1`	56	`exit 1`
57	`endif`	57	`endif`
58		58
59	`# ---------------------------------------------------------------------`	59	`# ---------------------------------------------------------------------`
60	`# Set some general parameters`	60	`# Set some general parameters`
61	`# ---------------------------------------------------------------------`	61	`# ---------------------------------------------------------------------`
62		62
63	`# Set the prefix of the primary and secondary data files`	63	`# Set the prefix of the primary and secondary data files`
64	`set charp = 'P'`	64	`set charp = 'P'`
65	`set chars = 'S'`	65	`set chars = 'S'`
66		66
67	`# Set the name of the run script`	67	`# Set the name of the run script`
68	`set runscript = "runscript"`	68	`set runscript = "runscript"`
69		69
70	`# Remember the calling directory`	70	`# Remember the calling directory`
71	`set calldir = $PWD`	71	`set calldir = $PWD`
72		72
73	`# ---------------------------------------------------------------------`	73	`# ---------------------------------------------------------------------`
74	`# Parameter handling`	74	`# Parameter handling`
75	`# ---------------------------------------------------------------------`	75	`# ---------------------------------------------------------------------`
76		76
77	`# Write title`	77	`# Write title`
78	`echo`	78	`echo`
79	`echo '========================================================='`	79	`echo '========================================================='`
80	`echo ' * START OF LAGRANTO * '`	80	`echo ' * START OF LAGRANTO * '`
81	`echo`	81	`echo`
82		82
83	`# Save the complete argument list`	83	`# Save the complete argument list`
84	`set lagranto_call = "$*"`	84	`set lagranto_call = "$*"`
85		85
86	`echo '---- HANDLE ARGUMENTS -----------------------------------'`	86	`echo '---- HANDLE ARGUMENTS -----------------------------------'`
87	`echo`	87	`echo`
88		88
89	`# Get fixed arguments`	89	`# Get fixed arguments`
90	`if (${#argv} == 1) then`	90	`if (${#argv} == 1) then`
91	`set mode = '-open'`	91	`set mode = '-open'`
92	`set caseid = "$1"`	92	`set caseid = "$1"`
93	`shift`	93	`shift`
94		94
95	`else if (${#argv} == 2) then`	95	`else if (${#argv} == 2) then`
96	`set mode = "$1"`	96	`set mode = "$1"`
97	`set caseid = "$2"`	97	`set caseid = "$2"`
98	`shift`	98	`shift`
99	`shift`	99	`shift`
100		100
101	`else`	101	`else`
102	`set caseid = "$1"`	102	`set caseid = "$1"`
103	`set startdate = "$2"`	103	`set startdate = "$2"`
104	`set enddate = "$3"`	104	`set enddate = "$3"`
105	`set startf = "$4"`	105	`set startf = "$4"`
106	`set select = "$5"`	106	`set select = "$5"`
107	`set mode = "runlagranto"`	107	`set mode = "runlagranto"`
108	`shift`	108	`shift`
109	`shift`	109	`shift`
110	`shift`	110	`shift`
111	`shift`	111	`shift`
112	`shift`	112	`shift`
113		113
114	`endif`	114	`endif`
115		115
116	`# Split <caseid> into base and appendix`	116	`# Split <caseid> into base and appendix`
117	set casebase = `(echo $caseid \| awk -F "." '{print $1}')`	117	set casebase = `(echo $caseid \| awk -F "." '{print $1}')`
118	set appendix = `(echo $caseid \| awk -F "." '{print $2}')`	118	set appendix = `(echo $caseid \| awk -F "." '{print $2}')`
119	`if ( "${appendix}" != "" ) then`	119	`if ( "${appendix}" != "" ) then`
120	`set appendix = ".${appendix}"`	120	`set appendix = ".${appendix}"`
121	`endif`	121	`endif`
122		122
123	`# Set the directories (depending on <caseid>)`	123	`# Set the directories (depending on <caseid>)`
124	`if ( "${casebase}" == "local" ) then`	124	`if ( "${casebase}" == "local" ) then`
125	`set cdfdir = $PWD`	125	`set cdfdir = $PWD`
126	`set tradir = $PWD`	126	`set tradir = $PWD`
127		127
128	`else if ( "${casebase}" == "interim" ) then`	128	`else if ( "${casebase}" == "interim" ) then`
129	`set cdfdir = /net/dansgaard/atmosdyn/erainterim/cdf/`	129	`set cdfdir = /net/dansgaard/atmosdyn/erainterim/cdf/`
130	`set tradir = $PWD`	130	`set tradir = $PWD`
131		131
132	`else if ( "${casebase}" == "analysis" ) then`	132	`else if ( "${casebase}" == "analysis" ) then`
133	`set cdfdir = 'TO_BE_SPECIFIED'`	133	`set cdfdir = 'TO_BE_SPECIFIED'`
134	`set tradir = $PWD`	134	`set tradir = $PWD`
135	`echo "ERROR: caseid (analysis) not yet implemented"`	135	`echo "ERROR: caseid (analysis) not yet implemented"`
136	`exit 1`	136	`exit 1`
137		137
138	`else if ( "${casebase}" == "forecast" ) then`	138	`else if ( "${casebase}" == "forecast" ) then`
139	`set cdfdir = 'TO_BE_SPECIFIED'`	139	`set cdfdir = 'TO_BE_SPECIFIED'`
140	`set tradir = $PWD`	140	`set tradir = $PWD`
141	`echo "ERROR: caseid (forecast) not yet implemented"`	141	`echo "ERROR: caseid (forecast) not yet implemented"`
142	`exit 1`	142	`exit 1`
143		143
144	`else`	144	`else`
145	`set cdfdir = ${HOME}/cdf/${casebase}`	145	`set cdfdir = ${HOME}/cdf/${casebase}`
146	`set tradir = ${HOME}/tra/${casebase}`	146	`set tradir = ${HOME}/tra/${casebase}`
147	`endif`	147	`endif`
148		148
149	`# Change to the run directory if requested`	149	`# Change to the run directory if requested`
150	`if ( "${mode}" == "-open" ) then`	150	`if ( "${mode}" == "-open" ) then`
151	`set count = 0`	151	`set count = 0`
152	`foreach folder ( ${tradir}/ntr__${caseid}_ )`	152	`foreach folder ( ${tradir}/ntr__${caseid}_ )`
153	`@ count = ${count} + 1`	153	`@ count = ${count} + 1`
154	`echo "[${count}] ${folder}"`	154	`echo "[${count}] ${folder}"`
155	`end`	155	`end`
156	`echo`	156	`echo`
157	`echo "Please select directory [0...${count}]"`	157	`echo "Please select directory [0...${count}]"`
158	`set id = $<`	158	`set id = $<`
159	`set count = 0`	159	`set count = 0`
160	`foreach folder ( ${tradir}/ntr__${caseid}_ )`	160	`foreach folder ( ${tradir}/ntr__${caseid}_ )`
161	`@ count = ${count} + 1`	161	`@ count = ${count} + 1`
162	`if ( "${count}" == "${id}" ) then`	162	`if ( "${count}" == "${id}" ) then`
163	`cd ${folder}`	163	`cd ${folder}`
164	`xterm &`	164	`xterm &`
165	`endif`	165	`endif`
166	`end`	166	`end`
167	`exit 0`	167	`exit 0`
168	`endif`	168	`endif`
169		169
170	`# Remove the whole run directory if requested`	170	`# Remove the whole run directory if requested`
171	`if ( "${mode}" == "-remove" ) then`	171	`if ( "${mode}" == "-remove" ) then`
172	`set count = 0`	172	`set count = 0`
173	`foreach folder ( ${tradir}/ntr__${caseid}_ )`	173	`foreach folder ( ${tradir}/ntr__${caseid}_ )`
174	`@ count = ${count} + 1`	174	`@ count = ${count} + 1`
175	`echo "[${count}] ${folder}"`	175	`echo "[${count}] ${folder}"`
176	`end`	176	`end`
177	`echo`	177	`echo`
178	`echo "Please select directory [0...${count}, all]"`	178	`echo "Please select directory [0...${count}, all]"`
179	`set id = $<`	179	`set id = $<`
180	`set count = 0`	180	`set count = 0`
181	`foreach folder ( ${tradir}/ntr__${caseid}_ )`	181	`foreach folder ( ${tradir}/ntr__${caseid}_ )`
182	`@ count = ${count} + 1`	182	`@ count = ${count} + 1`
183	`if ( ( "${count}" == "${id}" ) \| ( "${id}" == "all" ) ) then`	183	`if ( ( "${count}" == "${id}" ) \| ( "${id}" == "all" ) ) then`
184	`\rm -r ${folder}`	184	`\rm -r ${folder}`
185	`endif`	185	`endif`
186	`end`	186	`end`
187	`exit 0`	187	`exit 0`
188	`endif`	188	`endif`
189		189
190	`# Show the trajectory file`	190	`# Show the trajectory file`
191	`if ( "${mode}" == "-show" ) then`	191	`if ( "${mode}" == "-show" ) then`
192	`set count = 0`	192	`set count = 0`
193	`foreach folder ( ${tradir}/ntr__${caseid}_ )`	193	`foreach folder ( ${tradir}/ntr__${caseid}_ )`
194	`@ count = ${count} + 1`	194	`@ count = ${count} + 1`
195	`echo "[${count}] ${folder}"`	195	`echo "[${count}] ${folder}"`
196	`end`	196	`end`
197	`echo`	197	`echo`
198	`echo "Please select directory [0...${count}]"`	198	`echo "Please select directory [0...${count}]"`
199	`set id = $<`	199	`set id = $<`
200	`set count = 0`	200	`set count = 0`
201	`foreach folder ( ${tradir}/ntr__${caseid}_ )`	201	`foreach folder ( ${tradir}/ntr__${caseid}_ )`
202	`@ count = ${count} + 1`	202	`@ count = ${count} + 1`
203	`if ( "${count}" == "${id}" ) then`	203	`if ( "${count}" == "${id}" ) then`
204	`cd ${folder}`	204	`cd ${folder}`
205	set filename = ` sed -ne '9,9p' runscript.sh \| cut -c 2-`	205	set filename = ` sed -ne '9,9p' runscript.sh \| cut -c 2-`
206	`${LAGRANTO}/goodies/trainfo.sh ${filename} list \| less`	206	`${LAGRANTO}/goodies/trainfo.sh ${filename} list \| less`
207	`endif`	207	`endif`
208	`end`	208	`end`
209	`exit 0`	209	`exit 0`
210	`endif`	210	`endif`
211		211
212	`# Decide whether <select> is a file or an explicit criterion`	212	`# Decide whether <select> is a file or an explicit criterion`
213	`set flag_select = 'file'`	213	`set flag_select = 'file'`
214	set test = `echo ${select} \| grep ':' \| wc -c`	214	set test = `echo ${select} \| grep ':' \| wc -c`
215	`if ( "${test}" != "0" ) then`	215	`if ( "${test}" != "0" ) then`
216	`set flag_select = 'criterion'`	216	`set flag_select = 'criterion'`
217	`endif`	217	`endif`
218		218
219	`# Decide whether <startf> is a file or an explicit criterion`	219	`# Decide whether <startf> is a file or an explicit criterion`
220	`set flag_startf = 'file'`	220	`set flag_startf = 'file'`
221		221
222	set test = `echo ${startf} \| grep ' ' \| wc -c`	222	set test = `echo ${startf} \| grep ' ' \| wc -c`
223	`if ( "${test}" != "0" ) then`	223	`if ( "${test}" != "0" ) then`
224	`set flag_startf = 'point'`	224	`set flag_startf = 'point'`
225	`endif`	225	`endif`
226		226
227	set test = `echo ${startf} \| grep '@' \| wc -c`	227	set test = `echo ${startf} \| grep '@' \| wc -c`
228	`if ( "${test}" != "0" ) then`	228	`if ( "${test}" != "0" ) then`
229	`set flag_startf = 'criterion'`	229	`set flag_startf = 'criterion'`
230	`endif`	230	`endif`
231		231
232		232
233	`# Write some status information`	233	`# Write some status information`
234	`echo " Case ID : ${caseid}"`	234	`echo " Case ID : ${caseid}"`
235	`echo " Start date : ${startdate}"`	235	`echo " Start date : ${startdate}"`
236	`echo " End date : ${enddate}"`	236	`echo " End date : ${enddate}"`
237		237
238	`if ( "${flag_startf}" == "criterion" ) then`	238	`if ( "${flag_startf}" == "criterion" ) then`
239	`echo " Start file : ${startf} [criterion] -> startf.criterion [file]"`	239	`echo " Start file : ${startf} [criterion] -> startf.criterion [file]"`
240	`else if ( "${flag_startf}" == "file" ) then`	240	`else if ( "${flag_startf}" == "file" ) then`
241	`echo " Start points : ${startf} [file]"`	241	`echo " Start points : ${startf} [file]"`
242	`else if ( "${flag_startf}" == "point" ) then`	242	`else if ( "${flag_startf}" == "point" ) then`
243	`echo " Start points : ${startf} [point] -> startf [file]"`	243	`echo " Start points : ${startf} [point] -> startf [file]"`
244	`endif`	244	`endif`
245		245
246	`if ( "${flag_select}" == "criterion" ) then`	246	`if ( "${flag_select}" == "criterion" ) then`
247	`echo " Selection : ${select} [criterion] -> selectf.criterion [file]"`	247	`echo " Selection : ${select} [criterion] -> selectf.criterion [file]"`
248	`else`	248	`else`
249	`echo " Selection : ${select} [file]"`	249	`echo " Selection : ${select} [file]"`
250	`endif`	250	`endif`
251		251
252	`echo`	252	`echo`
253		253
254	`# Set defaults for flags`	254	`# Set defaults for flags`
255	`set jumpflag = ""`	255	`set jumpflag = ""`
256	`set tracefile = "tracevars"`	256	`set tracefile = "tracevars"`
257	`set outfile = "lsl_${startdate}"`	257	`set outfile = "lsl_${startdate}"`
258	`set keepflag = ""`	258	`set keepflag = ""`
259	`set prepflag = ""`	259	`set prepflag = ""`
260	`set logfile = "${runscript}.logfile"`	260	`set logfile = "${runscript}.logfile"`
261	`set regionf = "regionf"`	261	`set regionf = "regionf"`
262	`set changet = "false"`	262	`set changet = "false"`
263		263
264	`# Handle flags`	264	`# Handle flags`
265	`while ( $#argv > 0 )`	265	`while ( $#argv > 0 )`
266		266
267	`switch ( $argv[1] )`	267	`switch ( $argv[1] )`
268		268
269	`echo $argv[1]`	269	`echo $argv[1]`
270		270
271	`case -j`	271	`case -j`
272	`set jumpflag = "-j"`	272	`set jumpflag = "-j"`
273	`echo " Flag '-j [jump]' -> True "`	273	`echo " Flag '-j [jump]' -> True "`
274	`echo`	274	`echo`
275	`breaksw`	275	`breaksw`
276		276
277	`case -prep`	277	`case -prep`
278	`set prepflag = "prep"`	278	`set prepflag = "prep"`
279	`echo " Flag '-prep [prepare only]' -> ${prepflag} "`	279	`echo " Flag '-prep [prepare only]' -> ${prepflag} "`
280	`echo`	280	`echo`
281	`breaksw`	281	`breaksw`
282		282
283	`case -noclean`	283	`case -noclean`
284	`set keepflag = "keep"`	284	`set keepflag = "keep"`
285	`echo " Flag '-noclean [no cleaning]' -> ${keepflag} "`	285	`echo " Flag '-noclean [no cleaning]' -> ${keepflag} "`
286	`echo`	286	`echo`
287	`breaksw`	287	`breaksw`
288		288
289	`case -o`	289	`case -o`
290	`set outfile = $argv[2]`	290	`set outfile = $argv[2]`
291	`echo " Flag '-o [output]' -> ${outfile} "`	291	`echo " Flag '-o [output]' -> ${outfile} "`
292	`echo`	292	`echo`
293	`breaksw`	293	`breaksw`
294		294
295	`case -v`	295	`case -v`
296	`set tracefile = $argv[2]`	296	`set tracefile = $argv[2]`
297	`echo " Flag '-v [tracevars]' -> ${tracefile} "`	297	`echo " Flag '-v [tracevars]' -> ${tracefile} "`
298	`echo`	298	`echo`
299	`shift`	299	`shift`
300	`breaksw`	300	`breaksw`
301		301
302	`case -r`	302	`case -r`
303	`set regionfile = $argv[2]`	303	`set regionfile = $argv[2]`
304	`echo " Flag '-r [regionfile]' -> ${regionfile} "`	304	`echo " Flag '-r [regionfile]' -> ${regionfile} "`
305	`echo`	305	`echo`
306	`shift`	306	`shift`
307	`breaksw`	307	`breaksw`
308		308
309	`case -log`	309	`case -log`
310	`set logfile = "/dev/stdout"`	310	`set logfile = "/dev/stdout"`
311	`echo " Flag '-v [log]' -> ${logfile} "`	311	`echo " Flag '-v [log]' -> ${logfile} "`
312	`breaksw`	312	`breaksw`
313		313
314	`case -changet`	314	`case -changet`
315	`set changet = "true"`	315	`set changet = "true"`
316	`echo " Flag '-changet [change times]'-> True "`	316	`echo " Flag '-changet [change times]'-> True "`
317	`echo`	317	`echo`
318	`breaksw`	318	`breaksw`
319		319
320	`endsw`	320	`endsw`
321		321
322	`shift;`	322	`shift;`
323		323
324	`end`	324	`end`
325		325
326	`# Decide whether a tracevars file is given or not`	326	`# Decide whether a tracevars file is given or not`
327	`if ( -f ${tradir}/${tracefile} ) then`	327	`if ( -f ${tradir}/${tracefile} ) then`
328	`set flag_tracevars = 'file'`	328	`set flag_tracevars = 'file'`
329	`else`	329	`else`
330	`set flag_tracevars = 'nil'`	330	`set flag_tracevars = 'nil'`
331	`endif`	331	`endif`
332		332
333	`echo`	333	`echo`
334	`if ( "${flag_tracevars}" == "file" ) then`	334	`if ( "${flag_tracevars}" == "file" ) then`
335	`echo " Trace Vars : ${tracefile} [file]"`	335	`echo " Trace Vars : ${tracefile} [file]"`
336	`else`	336	`else`
337	`echo " Trace Vars : [no fields are traced]"`	337	`echo " Trace Vars : [no fields are traced]"`
338	`endif`	338	`endif`
339	`echo`	339	`echo`
340		340
341	`# Set the name of the run directory (forward/backward distinction)`	341	`# Set the name of the run directory (forward/backward distinction)`
342	`set dir = 'f'`	342	`set dir = 'f'`
343	set range = `${LAGRANTO}/bin/gettidiff.sh ${enddate} ${startdate}`	343	set range = `${LAGRANTO}/bin/gettidiff.sh ${enddate} ${startdate}`
344	`if ( ${range} < 0 ) then`	344	`if ( ${range} < 0 ) then`
345	`set dir = 'b'`	345	`set dir = 'b'`
346	set range = `${LAGRANTO}/bin/gettidiff.sh ${startdate} ${enddate}`	346	set range = `${LAGRANTO}/bin/gettidiff.sh ${startdate} ${enddate}`
347	`endif`	347	`endif`
348	`set rundir = ${tradir}/ntr_${startdate}_${dir}${range}_${caseid}`	348	`set rundir = ${tradir}/ntr_${startdate}_${dir}${range}_${caseid}`
349		349
350	`if ( "${flag_startf}" == "criterion" ) then`	350	`if ( "${flag_startf}" == "criterion" ) then`
351	`set rundir = "${rundir}_startf"`	351	`set rundir = "${rundir}_startf"`
352	`else if ( "${flag_startf}" == "point" ) then`	352	`else if ( "${flag_startf}" == "point" ) then`
353	`set rundir = "${rundir}_startf"`	353	`set rundir = "${rundir}_startf"`
354	`else if ( "${flag_startf}" == "file" ) then`	354	`else if ( "${flag_startf}" == "file" ) then`
355	`set rundir = "${rundir}_${startf}"`	355	`set rundir = "${rundir}_${startf}"`
356	`endif`	356	`endif`
357		357
358	`if ( "${flag_select}" != "file" ) then`	358	`if ( "${flag_select}" != "file" ) then`
359	`set rundir = "${rundir}_selectf"`	359	`set rundir = "${rundir}_selectf"`
360	`else`	360	`else`
361	`set rundir = "${rundir}_${select}"`	361	`set rundir = "${rundir}_${select}"`
362	`endif`	362	`endif`
363		363
364	`# Write directory names`	364	`# Write directory names`
365	`echo '---- DIRECTORIES ----------------------------------------'`	365	`echo '---- DIRECTORIES ----------------------------------------'`
366	`echo`	366	`echo`
367	`echo " CDF = ${cdfdir}"`	367	`echo " CDF = ${cdfdir}"`
368	`echo " TRA = ${tradir}"`	368	`echo " TRA = ${tradir}"`
369	`echo " RUN = ${rundir}"`	369	`echo " RUN = ${rundir}"`
370		370
371	`# Check whether CDF and TRA directories exist`	371	`# Check whether CDF and TRA directories exist`
372	`if ( ! -d ${cdfdir} ) then`	372	`if ( ! -d ${cdfdir} ) then`
373	`echo 'Cannot find CDF directory ... Stop'`	373	`echo 'Cannot find CDF directory ... Stop'`
374	`exit 1`	374	`exit 1`
375	`endif`	375	`endif`
376	`if ( ! -d ${tradir} ) then`	376	`if ( ! -d ${tradir} ) then`
377	`echo 'Cannot find TRA directory ... Stop'`	377	`echo 'Cannot find TRA directory ... Stop'`
378	`exit 1`	378	`exit 1`
379	`endif`	379	`endif`
380		380
381	`# Create the run directory`	381	`# Create the run directory`
382	`if ( ! -d ${rundir} ) then`	382	`if ( ! -d ${rundir} ) then`
383	`mkdir ${rundir}`	383	`mkdir ${rundir}`
384	`endif`	384	`endif`
385	`if ( ! -d ${rundir} ) then`	385	`if ( ! -d ${rundir} ) then`
386	`echo 'Cannot create RUN directory ... Stop'`	386	`echo 'Cannot create RUN directory ... Stop'`
387	`exit 1`	387	`exit 1`
388	`endif`	388	`endif`
389		389
390	`# ---------------------------------------------------------------------`	390	`# ---------------------------------------------------------------------`
391	`# Handle the time specifier - startdate, enddate`	391	`# Handle the time specifier - startdate, enddate`
392	`# ---------------------------------------------------------------------`	392	`# ---------------------------------------------------------------------`
393		393
394	`echo`	394	`echo`
395	`echo '---- TIME RANGE -----------------------------------------'`	395	`echo '---- TIME RANGE -----------------------------------------'`
396	`echo`	396	`echo`
397		397
398	`# Check format of start and end date - must be the same`	398	`# Check format of start and end date - must be the same`
399	set ns=`echo $startdate \| sed -e 's/_[0-9]*//' \| wc -c`	399	set ns=`echo $startdate \| sed -e 's/_[0-9]*//' \| wc -c`
400	set ne=`echo $enddate \| sed -e 's/_[0-9]*//' \| wc -c`	400	set ne=`echo $enddate \| sed -e 's/_[0-9]*//' \| wc -c`
401	`if ( $ns != $ne ) then`	401	`if ( $ns != $ne ) then`
402	`echo " ERROR: start and end date must be in the same format ***"`	402	`echo " ERROR: start and end date must be in the same format ***"`
403	`exit 1`	403	`exit 1`
404	`endif`	404	`endif`
405	`if ( $ns != 9 ) then`	405	`if ( $ns != 9 ) then`
406	`echo " ERROR: Date format must be yyyymmdd ***"`	406	`echo " ERROR: Date format must be yyyymmdd ***"`
407	`exit 1`	407	`exit 1`
408	`endif`	408	`endif`
409	set ns=`echo $startdate \| sed -e 's/[0-9]*_//' \| wc -c`	409	set ns=`echo $startdate \| sed -e 's/[0-9]*_//' \| wc -c`
410	set ne=`echo $enddate \| sed -e 's/[0-9]*_//' \| wc -c`	410	set ne=`echo $enddate \| sed -e 's/[0-9]*_//' \| wc -c`
411	`if ( $ns != $ne ) then`	411	`if ( $ns != $ne ) then`
412	`echo " ERROR: start and end date must be in the same format ***"`	412	`echo " ERROR: start and end date must be in the same format ***"`
413	`exit 1`	413	`exit 1`
414	`endif`	414	`endif`
415	`if ( ( $ns != 5 ) & ( $ns != 3 ) ) then`	415	`if ( ( $ns != 5 ) & ( $ns != 3 ) ) then`
416	`echo " ERROR: Time format must be hh(mm) ***"`	416	`echo " ERROR: Time format must be hh(mm) ***"`
417	`exit 1`	417	`exit 1`
418	`endif`	418	`endif`
419		419
420	`# Split the start and end date into <yymmdd_hh and mm>`	420	`# Split the start and end date into <yymmdd_hh and mm>`
421	set startdate_ymdh = `echo $startdate \| cut -c 1-11`	421	set startdate_ymdh = `echo $startdate \| cut -c 1-11`
422	set startdate_min = `echo $startdate \| cut -c 12-13`	422	set startdate_min = `echo $startdate \| cut -c 12-13`
423	`if ( $startdate_min == "" ) set startdate_min = 00`	423	`if ( $startdate_min == "" ) set startdate_min = 00`
424		424
425	set enddate_ymdh = `echo $enddate \| cut -c 1-11`	425	set enddate_ymdh = `echo $enddate \| cut -c 1-11`
426	set enddate_min = `echo $enddate \| cut -c 12-13`	426	set enddate_min = `echo $enddate \| cut -c 12-13`
427	`if ( $enddate_min == "" ) set enddate_min = 00`	427	`if ( $enddate_min == "" ) set enddate_min = 00`
428		428
429	`# Get the time difference between <start_ymdh> and <end_ymdh> date`	429	`# Get the time difference between <start_ymdh> and <end_ymdh> date`
430	`# Decide whether trajectoriesare forward or backward`	430	`# Decide whether trajectoriesare forward or backward`
431	set timediff_hh = `${LAGRANTO}/goodies/gettidiff ${enddate_ymdh} ${startdate_ymdh}`	431	set timediff_hh = `${LAGRANTO}/goodies/gettidiff ${enddate_ymdh} ${startdate_ymdh}`
432		432
433	`if ( $timediff_hh == 0 ) then`	433	`if ( $timediff_hh == 0 ) then`
434	`if ( $enddate_min > $startdate_min ) then`	434	`if ( $enddate_min > $startdate_min ) then`
435	`set direction = f`	435	`set direction = f`
436	`set idir = 1`	436	`set idir = 1`
437	`else`	437	`else`
438	`set direction = b`	438	`set direction = b`
439	`set idir = -1`	439	`set idir = -1`
440	`endif`	440	`endif`
441	`else if ( $timediff_hh > 0 ) then`	441	`else if ( $timediff_hh > 0 ) then`
442	`set direction = f`	442	`set direction = f`
443	`set idir = 1`	443	`set idir = 1`
444	`else`	444	`else`
445	`set direction = b`	445	`set direction = b`
446	`set idir = -1`	446	`set idir = -1`
447	`@ timediff_hh = $idir * $timediff_hh`	447	`@ timediff_hh = $idir * $timediff_hh`
448	`endif`	448	`endif`
449		449
450	`# Get also minutes for time difference, if <start_min> or <end_min> != 0`	450	`# Get also minutes for time difference, if <start_min> or <end_min> != 0`
451	`set timediff_mm=`	451	`set timediff_mm=`
452		452
453	`if ( $startdate_min != 00 \|\| $enddate_min != 00 ) then`	453	`if ( $startdate_min != 00 \|\| $enddate_min != 00 ) then`
454	`@ min = ( $enddate_min - $startdate_min )`	454	`@ min = ( $enddate_min - $startdate_min )`
455	`if ( $min == 0 ) then`	455	`if ( $min == 0 ) then`
456	`set timediff_mm=`	456	`set timediff_mm=`
457	`else if ( $min > 0 ) then`	457	`else if ( $min > 0 ) then`
458	`if ( $idir == 1 ) then`	458	`if ( $idir == 1 ) then`
459	`set timediff_mm=$min`	459	`set timediff_mm=$min`
460	`else`	460	`else`
461	`@ timediff_hh --`	461	`@ timediff_hh --`
462	`@ timediff_mm = 60 - $min`	462	`@ timediff_mm = 60 - $min`
463	`endif`	463	`endif`
464	`else`	464	`else`
465	`if ( $idir == 1 ) then`	465	`if ( $idir == 1 ) then`
466	`@ timediff_hh --`	466	`@ timediff_hh --`
467	`@ timediff_mm = 60 + $min`	467	`@ timediff_mm = 60 + $min`
468	`else`	468	`else`
469	`@ timediff_mm = 0 - $min`	469	`@ timediff_mm = 0 - $min`
470	`endif`	470	`endif`
471	`endif`	471	`endif`
472	`endif`	472	`endif`
473		473
474	`# Set the reference date equal to the startdate`	474	`# Set the reference date equal to the startdate`
475	`set refdate=${startdate}`	475	`set refdate=${startdate}`
476		476
477	`# Write status information`	477	`# Write status information`
478	`echo "Time range : ${startdate} -> ${enddate}"`	478	`echo "Time range : ${startdate} -> ${enddate}"`
479	`if ( ${timediff_mm} != "" ) then`	479	`if ( ${timediff_mm} != "" ) then`
480	`echo "Time difference : ${timediff_hh} h ${timediff_mm} min"`	480	`echo "Time difference : ${timediff_hh} h ${timediff_mm} min"`
481	`else`	481	`else`
482	`echo "Time difference : ${timediff_hh} h"`	482	`echo "Time difference : ${timediff_hh} h"`
483	`endif`	483	`endif`
484	`echo "Direction : ${direction} (${idir})"`	484	`echo "Direction : ${direction} (${idir})"`
485	`echo "Reference date : ${refdate}"`	485	`echo "Reference date : ${refdate}"`
486		486
487	`# ---------------------------------------------------------------------`	487	`# ---------------------------------------------------------------------`
488	`# Decide which files are needed and check availability`	488	`# Decide which files are needed and check availability`
489	`# ---------------------------------------------------------------------`	489	`# ---------------------------------------------------------------------`
490		490
491	`echo`	491	`echo`
492	`echo '---- INPUT FILES ----------------------------------------'`	492	`echo '---- INPUT FILES ----------------------------------------'`
493	`echo`	493	`echo`
494		494
495	`# Change to input data (cdf) directory`	495	`# Change to input data (cdf) directory`
496	`cd ${cdfdir}`	496	`cd ${cdfdir}`
497	`echo`	497	`echo`
498		498
499	`# Take the time increment from flag list ('nil', if not defined)`	499	`# Take the time increment from flag list ('nil', if not defined)`
500	`set timeinc = 'nil'`	500	`set timeinc = 'nil'`

Subversion Repositories lagranto.um

(root)/trunk/bin/lagranto – Rev 3 → 7