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


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

Subversion Repositories lagranto.um

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