--- MITgcm/verification/testreport	2008/04/02 23:41:22	1.114
+++ MITgcm/verification/testreport	2009/04/16 18:46:59	1.123
@@ -1,6 +1,6 @@
 #! /usr/bin/env bash
 #
-#  $Header: /home/ubuntu/mnt/e9_copy/MITgcm/verification/testreport,v 1.114 2008/04/02 23:41:22 jmc Exp $
+#  $Header: /home/ubuntu/mnt/e9_copy/MITgcm/verification/testreport,v 1.123 2009/04/16 18:46:59 jmc Exp $
 #  $Name:  $
 #
 
@@ -11,7 +11,7 @@
     echo
     echo "where possible OPTIONS are:"
     echo "  (-help|-h)               print usage"
-    echo "  (-mth)                   run multi threaded (using eedata.mth)"
+    echo "  (-mth)                   run multi-threaded (using eedata.mth)"
     echo "  (-mpi)                   compile and run using MPI"
     echo "  (-ieee|-noieee)          if possible, use IEEE compiler flags"
     echo "                             (DEF=\"-ieee\")"
@@ -30,6 +30,7 @@
     echo "                             Bourne-compatible \"sh\" shell"
     echo "                             (DEF=\"\" for \"bash\")"
     echo "  (-adm|-ad)               perform an adjoint run"
+    echo "  (-oad)                   perform an OpenAD adjoint run"
     echo "  (-command) STRING        command to run"
     echo "                             (DEF=\"make output.txt\")"
     echo "  (-m|-make) STRING        command to use for \"make\""
@@ -179,31 +180,6 @@
     return $digits_of_similarity
 }
 
-check_for_add_mon_output()
-{
-    # Check for additional types of monitor output
-    if test "x$1" = x ; then
-	return
-    fi
-
-    for ii in $PTRACERS_NUM ; do
-	eval "HAVE_PTR0"$ii"=f"
-    done
-
-    ptr_add="trcstat_ptracerXX_min trcstat_ptracerXX_max"
-    ptr_add="$ptr_add trcstat_ptracerXX_mean trcstat_ptracerXX_sd"
-    for ii in $PTRACERS_NUM ; do
-	for jj in $ptr_add ; do
-	    name=`eval "echo $jj | sed -e 's|XX|0"$ii"|g'"`
-	    tst=`grep $name $1 | wc -l | awk '{print $1}'`
-	    if test ! "x$tst" = x0 ; then
-		eval "HAVE_PTR0"$ii"=t"
-	    fi
-	done
-	#  eval 'echo "HAVE_PTR0'$ii' = $HAVE_PTR0'$ii'"'
-    done
-}
-
 testoutput_run()
 {
     # testoutput_run directory subdir reference_output
@@ -221,12 +197,30 @@
         listVar=`echo $listChk | sed 's/ [a-zA-Z0-9]*+/&mn &mx &av &sd/g' \
                                | sed 's/+//g' | sed "s/^$sVar//"`
 	if [ $debug -gt 0 ]; then echo "testoutput_run: listVar(I)='$listVar'" 1>&2 ; fi
- 	for ii in 1 2 3 4 5 6 7 8 9 ; do
-	  tst=`eval 'echo "$HAVE_PTR0'$ii'"'`
-	 #echo "-- ptr test=" $tst "number of var=" `echo $listVar | awk '{print NF}'` 1>&2
- 	  if test "x$tst" != xt ; then listVar=`echo "$listVar" | sed "s/ pt$ii..//g"` ; fi
-        done
-        tst=`echo $sVar $listVar | awk '{ for(i=2;i<=NF;i++){t+=($i==$1)}; print t }'`
+        # check for ptracer output in reference_output file :
+        outpref=$1/results/$3
+	ptr_mon="trcstat_ptracerXX_min trcstat_ptracerXX_max"
+	ptr_mon="$ptr_mon trcstat_ptracerXX_mean trcstat_ptracerXX_sd"
+	for ii in $PTRACERS_NUM ; do
+	    ptrfound=0
+	    for jj in $ptr_mon ; do
+		name=`eval "echo $jj | sed -e 's|XX|0"$ii"|g'"`
+		tst=`grep $name $outpref | wc -l | awk '{print $1}'`
+		if test ! "x$tst" = x0 ; then ptrfound=1 ; fi
+	    done
+	    if test $ptrfound = '1' ; then
+		eval "HAVE_PTR0"$ii"=t"
+	    else
+		eval "HAVE_PTR0"$ii"=f"
+	      if test "x$ADM" = x -a "x$OADM" = x; then
+	      # remove this ptr from the list of output variable to check
+	      # echo "-- ptr test=" $tst "number of var=" `echo $listVar | awk '{print NF}'` 1>&2
+		listVar=`echo "$listVar" | sed "s/ pt$ii..//g"`
+	      fi
+	    fi
+	#   eval 'echo "HAVE_PTR0'$ii' = $HAVE_PTR0'$ii'"' 1>&2
+	done
+        tst=`echo $sVar $listVar | awk '{ for(i=2;i<=NF;i++){if($i==$1)t+=1}; print t }'`
         if test $tst != 1 ; then
           if test $tst = 0 ; then echo "==> WARNING: selected var >$sVar< not found" 1>&2
                  else echo "==> WARNING: found selected var >$sVar< $tst times" 1>&2 ; fi
@@ -321,10 +315,12 @@
 	(
 	    cd $1;
 	    command="$GENMAKE2  -ds -m $MAKE"
-	    if test "x$ADM" = x ; then
-		command="$command --mods=../code"
-	    else
+	    if test "x$ADM" = xt ; then
 		command="$command --mods=../code_ad"
+	    elif test "x$OADM" = xt ; then
+		command="$command -adof ../../../tools/adjoint_options/adjoint_oad -mods '../code_oad ../../OpenAD/code_oad_all'"
+	    else
+		command="$command --mods=../code"
 	    fi
 	    if test "x$OPTFILE" != xNONE ; then
 		command="$command --optfile=$OPTFILE"
@@ -338,6 +334,10 @@
 	    if test "x$MPI" = xt ; then
 		command="$command -mpi"
 	    fi
+	    if test "x$MULTI_THREAD" = xt ; then
+	    #- run multi-threaded using OpenMP:
+		command="$command -omp"
+	    fi
 	    if test "x$TS" = xt ; then
 		command="$command -ts"
 	    fi
@@ -349,7 +349,8 @@
 	    fi
 	    fi
 	    printf 'genmake ... ' 1>&2
-	    $command > make.log 2>&1
+	    # $command > make.log 2>&1
+	    eval $command > make.log 2>&1
 	    RETVAL=$?
 	    #  Reduce the size of the testing emails!
 	    head -100 Makefile > $CDIR/Makefile_head
@@ -444,14 +445,16 @@
 	cd $1;
 	if test -r Makefile ; then
 	    printf 'make ... ' 1>&2
-	    if test "x$ADM" = x ; then
+	    if test "x$ADM" = xt ; then
+		$MAKE adall >> make.log 2>&1
+	    elif test "x$OADM" = xt ; then
+		$MAKE adAll >> make.log 2>&1
+            else
 		if test "x$JOBS" = x ; then
 		    $MAKE >> make.log 2>&1
 		else
 		    $MAKE -j $JOBS >> make.log 2>&1
 		fi
-	    else
-		$MAKE adall >> make.log 2>&1
 	    fi
 	    RETVAL=$?
 	    if test "x$RETVAL" != x0 ; then
@@ -573,52 +576,52 @@
 	echo
 	if test -L $EXECUTABLE -a -x "../"$builddir"/"$EXECUTABLE ; then
 	    diff -q $EXECUTABLE "../"$builddir"/"$EXECUTABLE > /dev/null 2>&1
-	    outD=$? ; if test $outD != 0 ; then rm -f $EXECUTABLE ; rm -f run.log ; fi
+	    outD=$? ; if test $outD != 0 ; then rm -f $EXECUTABLE ; fi
 	fi
         if test ! -x $EXECUTABLE -a -x "../"$builddir"/"$EXECUTABLE ; then
-	    echo " link" $EXECUTABLE "from dir ../"$builddir > run.log_00
+	    echo " link" $EXECUTABLE "from dir ../"$builddir > run.log_tmp
 	    ln -sf "../"$builddir"/"$EXECUTABLE .
 	fi
         if test ! -x $EXECUTABLE ; then
-		rm -f run.log ; touch run.log
-		if test -f run.log_00 ; then cat run.log_00 >> run.log ; fi
-		echo " no executable:" $EXECUTABLE >> run.log
+		rm -f $RUNLOG ; touch $RUNLOG
+		if test -f run.log_tmp ; then cat run.log_tmp >> $RUNLOG ; fi
+		echo " no executable:" $EXECUTABLE >> $RUNLOG
 		RETVAL=8
 		ENDVAL=-1
 	else
             if test ! -f $OUTPUTFILE -o $OUTPUTFILE -ot $EXECUTABLE ; then
               # output do not exist or is older than executable:
-		rm -f run.log ; touch run.log
-		if test -f run.log_00 ; then cat run.log_00 >> run.log ; fi
-		( eval $COMMAND ) >> run.log 2>&1
+		rm -f $RUNLOG ; touch $RUNLOG
+		if test -f run.log_tmp ; then cat run.log_tmp >> $RUNLOG ; fi
+		( eval $COMMAND ) >> $RUNLOG 2>&1
 		RETVAL=$?
 	    else
 		RETVAL=0
-		if test -f run.log ; then
-		    if test -f run.log_00 ; then cat run.log_00 >> run.log ; fi
-     		    echo "---------->> $OUTPUTFILE is up to date " >> run.log 2>&1
+		if test -f $RUNLOG ; then
+		    if test -f run.log_tmp ; then cat run.log_tmp >> $RUNLOG ; fi
+     		    echo "---------->> $OUTPUTFILE is up to date " >> $RUNLOG 2>&1
 		else 
-		    touch run.log
-		    if test -f run.log_00 ; then cat run.log_00 >> run.log ; fi
-     		    echo "---------->> $OUTPUTFILE is up to date " >> run.log 2>&1
-     		    echo " no previous run.log: assume NORMAL END" >> run.log 2>&1
+		    touch $RUNLOG
+		    if test -f run.log_tmp ; then cat run.log_tmp >> $RUNLOG ; fi
+     		    echo "---------->> $OUTPUTFILE is up to date " >> $RUNLOG 2>&1
+     		    echo " no previous $RUNLOG: assume NORMAL END" >> $RUNLOG 2>&1
 		fi
 	    fi
-	    ENDVAL=`cat run.log | grep -v 'ABNORMAL END' | grep -c 'NORMAL END'`
+	    ENDVAL=`cat $RUNLOG | grep -v 'ABNORMAL END' | grep -c 'NORMAL END'`
 	fi
-	rm -f run.log_00
+	rm -f run.log_tmp
 	#if test "x$RETVAL" = x0 ; then
 	if [ $RETVAL -eq 0 -a $ENDVAL -gt 0 ] ; then
-	    tail run.log
+	    tail $RUNLOG
 	    echo successful 1>&2
 	    # === Reduce the size of the testing emails!
 	    #cp $OUTPUTFILE $CDIR"/"$OUTPUTFILE
             if test -s STDERR.0000 ; then cp STDERR.0000 $CDIR"/STDERR.0000" ; fi
 	    return 0
 	else
-	    tail run.log
+	    tail $RUNLOG
 	    echo failed '(run:' $RETVAL ' end:' $ENDVAL ')' 1>&2
-	    cp run.log $CDIR"/run.log"
+	    cp $RUNLOG $CDIR"/"$RUNLOG
             if test -s STDERR.0000 ; then cp STDERR.0000 $CDIR"/STDERR.0000" ; fi
 	    return 1
 	fi
@@ -719,6 +722,35 @@
 }
 
 
+check_eedata()
+{
+    # check_eedata size.h eedata
+    if [ $# -eq 2 ] ; then
+      if test -f $1 -a -f $2 ; then
+      sx=`grep "^     & *nSx *=" $1 | sed "s/^     & *nSx *=//" | sed 's/, *$//'`
+      nx=`grep "^ *nTx *=" $2 | tail -1 | sed 's/^ *nTx *= *//' | sed "s/, *$//"`
+      if test "x$nx" = x ; then
+	rx=10
+      else
+	rx=`expr $sx % $nx`
+      fi
+      sy=`grep "^     & *nSy *=" $1 | sed "s/^     & *nSy *=//" | sed 's/, *$//'`
+      ny=`grep "^ *nTy *=" $2 | tail -1 | sed 's/^ *nTy *= *//' | sed "s/, *$//"`
+      if test "x$ny" = x ; then
+	ry=20
+      else
+	ry=`expr $sy % $ny`
+      fi
+      echo `expr $rx + $ry`
+    else
+      echo '-1'
+      fi
+    else
+      echo '-2'
+    fi
+
+}
+
 ###############################################################################
 
 
@@ -762,8 +794,9 @@
 DELDIR=
 
 ADM=
+OADM=
 
-# Additional monitor types
+# list of pTracers to check for monitor output
 PTRACERS_NUM="1 2 3 4 5"
 
 MATCH_CRIT=13
@@ -863,6 +896,8 @@
 
 	-adm | -ad) ADM=t ;;
 
+	-oad) OADM=t; NODEPEND=t ;;
+
 	-ieee) IEEE=true ;;
 	-noieee) IEEE= ;;
         -gsl) GSL=t ;;
@@ -905,6 +940,11 @@
     inputdir=input_ad
     ref_outp="output_adm.txt"
     EXECUTABLE="mitgcmuv_ad"
+elif test "x$OADM" = xt ; then
+    code_dir=code_oad
+    inputdir=input_oad
+    ref_outp="output_oadm.txt"
+    EXECUTABLE="mitgcmuv_ad"
 else
     code_dir=code
     inputdir=input
@@ -957,6 +997,7 @@
     OPTFILE=$MITGCM_OF
 fi
 
+RUNLOG="run.log"
 OUTPUTFILE=$ref_outp
 if test "x$COMMAND" = x ; then
     COMMAND="./$EXECUTABLE > $OUTPUTFILE"
@@ -971,7 +1012,7 @@
 #  (use default or load experiment-specific list from file "tr_checklist")
 # content : 1rst = main variable used to decide if it pass or FAIL
 #         others = number of matching digits to be printed in summary.txt
-if test "x$ADM" = x ; then
+if test "x$ADM" = x -a "x$OADM" = x; then
     DEF_CHECK_LIST='PS PS T+ S+ U+ V+ pt1+ pt2+ pt3+ pt4+ pt5+'
     EMPTY_RESULTS='.. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. ..'
     LEN_CHECK_LIST=`echo $DEF_CHECK_LIST | sed 's/ [a-zA-Z0-9]*+/&mn &mx &av &sd/g' | awk '{print NF-1}'`
@@ -995,13 +1036,20 @@
 
 #  Create a uniquely named directory to store results
 CMDLINE=$0
-for xx in "$@" ; do CMDLINE="$CMDLINE '$xx'" ; done
+for xx in "$@" ; do nw=`echo $xx | wc -w`
+    if test $nw = '1' ; then CMDLINE="$CMDLINE $xx"
+                        else CMDLINE="$CMDLINE '$xx'" ; fi
+done
+#for xx in "$@" ; do CMDLINE="$CMDLINE '$xx'" ; done
 MACH=`hostname`
 UNAMEA=`uname -a`
 DATE=`date +%Y%m%d`
 BASE="tr_"$MACH"_"$DATE"_"
 if test "x$OUTDIR" != x ; then
     BASE="tr_"$OUTDIR"_"$DATE"_"
+else
+    short_name=`hostname | sed 's/\..*$//'`
+    BASE="tr_"$short_name"_"$DATE"_"
 fi
 DNUM=0
 DRESULTS="$BASE$DNUM"
@@ -1048,7 +1096,7 @@
 fi
 echo
 echo >> $SUMMARY
-if test "x$ADM" = x ; then
+if test "x$ADM" = x -a "x$OADM" = x; then
     if [ $MATCH_CRIT -lt 10 ] ; 
     then line_0="default  "$MATCH_CRIT ; 
     else line_0="default "$MATCH_CRIT ; fi
@@ -1059,14 +1107,11 @@
     line_3="N n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d"
     line_4="2 d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  ."
     for ii in $PTRACERS_NUM ; do
-	#  tst=`eval 'echo $HAVE_PTR0'$ii`
-	#  if test "x$tst" = xt ; then
 	line_0="$line_0  --PTR 0"$ii"--"
 	line_1="$line_1        m  s"
 	line_2="$line_2  m  m  e  ."
 	line_3="$line_3  i  a  a  d"
 	line_4="$line_4  n  x  n  ."
-	#  fi
     done
     echo "$line_0" | tee -a $SUMMARY
     echo "$line_1" | tee -a $SUMMARY
@@ -1106,7 +1151,7 @@
     #  Cleanup only!
     if test "x$CLEANUP" = xt ; then
 	if test -r $BUILD_DIR/Makefile ; then
-            echo '  ------  clean dir:' $dir/build
+            echo '  ------  clean dir:' $dir/$builddir
 	    ( cd $BUILD_DIR ; make CLEAN )
 	fi
 	if test -d $dir/$rundir/CVS ; then
@@ -1127,11 +1172,8 @@
 	echo "can't read \"$fout\" -- skipping $dir"
 	continue
     fi
-    if test "x$ADM" = x ; then
-	check_for_add_mon_output  $fout
-    fi
 
-    # Check for additional types of monitor output
+    # Check for specific files for particular type of run
 
     if test ! -r $CODE_DIR"/SIZE.h_mpi" -a "x$MPI" = "xt" ; then
 	echo "can't find \"$CODE_DIR/SIZE.h_mpi\" -- skipping $dir"
@@ -1141,6 +1183,13 @@
 	echo "can't find \"$dir/input/eedata.mth\" -- skipping $dir"
 	continue
     fi
+    if test "x$MPI" = "xt" -a "x$MULTI_THREAD" = "xt" ; then
+      retv=`check_eedata $CODE_DIR"/SIZE.h_mpi" $dir"/input/eedata.mth"`
+      if test $retv != 0 ; then
+	echo "input/eedata.mth tiling misfit -- skipping $dir"
+	continue
+      fi
+    fi
 
     #  Check whether there are "extra runs" for this testdir
     extra_runs=
@@ -1151,14 +1200,23 @@
 	refExOut=`echo $ref_outp | sed "s/\./.${name}./"`
 	outf="$dir/results/$refExOut"
         if test -f $outf -a -r $outf ; then
-	    if test "x$MULTI_THREAD" = "xt" ; then 
-		if test -r $dir"/"$exd"/eedata.mth" ; then
+	  if test "x$MULTI_THREAD" = "xt" ; then 
+	    if test -r $dir"/"$exd"/eedata.mth" ; then
+	      if test "x$MPI" = "xt" ; then
+		retv=`check_eedata $CODE_DIR"/SIZE.h_mpi" $dir"/"$exd"/eedata.mth"`
+		if test $retv = 0 ; then
 		    extra_runs="$extra_runs $name"
-		#else echo $dir"/"$exd"/eedata.mth: not found"
+		else
+		    echo $exd"/eedata.mth tiling misfit -- skipping $dir"
 		fi
-	    else
+	      else
 		extra_runs="$extra_runs $name"
+	      fi
+	    #else echo $dir"/"$exd"/eedata.mth: not found"
 	    fi
+	  else
+	    extra_runs="$extra_runs $name"
+	  fi
 	fi
     done