Rev 3 | Blame | Compare with Previous | Last modification | View Log | Download | RSS feed
#!/bin/csh
# Set Lagranto
set LAGRANTO = ${LAGRANTOBASE}.${MODEL}/
# Write usage information
if ( ${#argv} == 0) then
echo
echo "----------------------------------------------------------------"
echo "Usage: difference.sh trainp1 trainp2 traout field [-single|-max]"
echo
echo " trainp1,2 : Input trajectory files"
echo " traout : Output trajectory file"
echo " field : Field for which difference is taken;"
echo " 'latlon' if spherical distance required"
echo
echo " -single : difference for all times"
echo " -max : only time for maximal difference"
echo "----------------------------------------------------------------"
exit 0
endif
# Get input and output trajectory file
set inpfile1 = $1
set inpfile2 = $2
set outfile = $3
set fieldname = $4
# Handle optional arguments
set mode = "single"
while ( $#argv > 0 )
switch ( $argv[1] )
case -single
set mode = "single"
breaksw
case -max
set mode = "max"
breaksw
endsw
shift;
end
# Get the dimensions of the trajectory files
set dim1=`${LAGRANTO}/goodies/trainfo.sh ${inpfile1} dim`
set dim2=`${LAGRANTO}/goodies/trainfo.sh ${inpfile2} dim`
# Prepare parameter file and run program
\rm -f difference.param
echo \"${inpfile1}\" >! difference.param
echo \"${inpfile2}\" >> difference.param
echo \"${outfile}\" >> difference.param
echo ${dim1} >> difference.param
echo ${dim2} >> difference.param
echo \"${mode}\" >> difference.param
echo \"${fieldname}\" >> difference.param
${LAGRANTO}/goodies/difference
# Make clean
#\rm -f difference.param
exit 0