/[MITgcm]/MITgcm/verification/testreport
ViewVC logotype

Diff of /MITgcm/verification/testreport

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.32 by edhill, Sun Dec 7 07:06:13 2003 UTC revision 1.80 by jmc, Thu Jun 8 21:39:35 2006 UTC
# Line 13  usage() Line 13  usage()
13      echo "  (-help|-h)               print usage"      echo "  (-help|-h)               print usage"
14      echo "  (-mpi)                   use MPI input files"      echo "  (-mpi)                   use MPI input files"
15      echo "  (-ieee|-noieee)          if possible, use IEEE compiler flags"      echo "  (-ieee|-noieee)          if possible, use IEEE compiler flags"
16      echo "                             (DEF=\"noieee\")"      echo "                             (DEF=\"-ieee\")"
17      echo "  (-optfile=|-of=)STRING   list of optfiles to use"      echo "  (-optfile=|-of=)STRING   list of optfiles to use"
18      echo "  (-a|-addr) STRING        list of email recipients"      echo "  (-a|-addr) STRING        list of email recipients"
19      echo "                             (DEF=\"edhill@mitgcm.org\")"      echo "                             (DEF=\"edhill@mitgcm.org\")"
# Line 27  usage() Line 27  usage()
27      echo "                             (DEF=\"make output.txt\")"      echo "                             (DEF=\"make output.txt\")"
28      echo "  (-m|-make) STRING        command to use for \"make\""      echo "  (-m|-make) STRING        command to use for \"make\""
29      echo "                             (DEF=\"make\")"      echo "                             (DEF=\"make\")"
30        echo "  (-odir) STRING           used to build output directory name"
31        echo "                             (DEF=\"hostname\")"
32        echo "  (-ptracers|-ptr) STRING  specify which ptracers to test"
33        echo "                             (DEF=\"1 2 3 4 5\")"
34        echo "  (-j) JOBS                use \"make -j JOBS\" for parallel builds"
35      echo "  (-clean)                 *ONLY* run \"make CLEAN\""      echo "  (-clean)                 *ONLY* run \"make CLEAN\""
36      echo "  (-quick|-q)              same as \"-nogenmake -noclean -nodepend\""      echo "  (-quick|-q)              same as \"-nogenmake -noclean -nodepend\""
37      echo "  (-nogenmake|-ng)         skip the genmake stage"      echo "  (-nogenmake|-ng)         skip the genmake stage"
38      echo "  (-noclean|-nc)           skip the \"make clean\" stage"      echo "  (-noclean|-nc)           skip the \"make clean\" stage"
39      echo "  (-nodepend|-nd)          skip the \"make depend\" stage"      echo "  (-nodepend|-nd)          skip the \"make depend\" stage"
40        echo "  (-deldir|-dd)            on success, delete the output directory"
41      echo      echo
42      echo "and where STRING follows a whitespace-delimited format"      echo "and where STRING can be a whitespace-delimited list"
43      echo "such as:"      echo "such as:"
44        echo
45      echo "  -t 'exp0 exp2 exp3' "      echo "  -t 'exp0 exp2 exp3' "
46      echo "  -addr='abc@123.com testing@home.org'"      echo "  -addr='abc@123.com testing@home.org'"
47      echo      echo
48        echo "provided that the expression is properly quoted within the current"
49        echo "shell (note the use of single quotes to protect white space)."
50        echo
51      exit 1      exit 1
52  }  }
53    
54  #  build the mpack utility  #  build the mpack utility
55  build_mpack()  build_mpack()
56  {  {
57      echo -n "building the mpack utility...  "      printf "building the mpack utility...  "
58      if test ! -x "$MPACKDIR/mpack" ; then      if test ! -x "$MPACKDIR/mpack" ; then
59          if test ! -d $MPACKDIR ; then          if test ! -d $MPACKDIR ; then
60              echo              echo
# Line 54  build_mpack() Line 64  build_mpack()
64              echo              echo
65              HAVE_MPACK=f              HAVE_MPACK=f
66          fi          fi
67          echo -n "building mpack...  "          printf "building mpack...  "
68          if test "x$CC" = x ; then          if test "x$CC" = x ; then
69              export CC=cc              export CC=cc
70          fi          fi
# Line 77  build_mpack() Line 87  build_mpack()
87    
88  testoutput_for_prop()  testoutput_for_prop()
89  {  {
90      # testoutput_for_prop dir s1 label subdir      # testoutput_for_prop dir s1 label subdir extension
91      #      #
92      #  compares files in $dir/$subdir/output.txt and $dir/results/output.txt      #  compares files in $dir/$subdir/output.txt and $dir/results/output.txt
93      #  using search strings s1 and text label      #  using search strings s1 and text label
# Line 87  testoutput_for_prop() Line 97  testoutput_for_prop()
97      fi      fi
98      if [ -r $1/$4/output.txt ]; then      if [ -r $1/$4/output.txt ]; then
99          grep "$2" $1/$4/output.txt | sed 's/.*=//' | cat -n > tmp1.txt          grep "$2" $1/$4/output.txt | sed 's/.*=//' | cat -n > tmp1.txt
100          lncnt=`wc -l tmp1.txt | awk '{print $1}' `          lncntA=`wc -l tmp1.txt | awk '{print $1}' `
101          if [ $lncnt -lt 3 ]; then          if [ $lncntA -lt 3 ]; then
102              if [ $verbose -gt 0 ]; then              if [ $verbose -gt 0 ]; then
103                  echo Not enough lines of output when searching for "$2" 1>&2                  echo Not enough lines of output when searching for "$2" 1>&2
104              fi              fi
# Line 99  testoutput_for_prop() Line 109  testoutput_for_prop()
109          return 99          return 99
110      fi      fi
111      if [ $debug -gt 0 ]; then      if [ $debug -gt 0 ]; then
112          echo testoutput_for_prop: grep "$2" $1/results/output.txt 1>&2          echo testoutput_for_prop: grep "$2" $1/results/output.txt$5 1>&2
113      fi      fi
114      grep "$2" $1/results/output.txt | sed 's/.*=//' | cat -n > tmp2.txt      grep "$2" $1/results/output.txt$5 | sed 's/.*=//' | cat -n > tmp2.txt
115      lncnt=`wc -l tmp2.txt | awk '{print $1}' `      lncntB=`wc -l tmp2.txt | awk '{print $1}' `
116      if [ $lncnt -lt 3 ]; then      if [ $lncntB -lt 3 ]; then
117          if [ $verbose -gt 0 ]; then          if [ $verbose -gt 0 ]; then
118              echo Not enough lines of output when searching for "$2" 1>&2              echo Not enough lines of output when searching for "$2" 1>&2
119          fi          fi
120          return 99          return 99
121      fi      fi
122        if [ $lncntA -ne $lncntB ]; then
123            if [ $verbose -gt 0 ]; then
124                echo Not same Nb of lines when searching for "$2" ":" $lncntA $lncntB 1>&2
125            fi
126            return 99
127        fi
128        has_nan=`cat tmp1.txt | grep -i nan | wc -l`
129        if [ $has_nan -gt 0  ] ; then
130            echo testoutput_for_prop: output.txt contains $has_nan NaN values  1>&2
131            return 99
132        fi
133        has_inf=`cat tmp1.txt | grep -i inf | wc -l`
134        if [ $has_inf -gt 0  ] ; then
135            echo testoutput_for_prop: output.txt contains $has_inf Inf values  1>&2
136            return 99
137        fi
138      if [ $debug -gt 0 ]; then      if [ $debug -gt 0 ]; then
139          echo testoutput_for_prop: join tmp1.txt tmp2.txt 1>&2          echo testoutput_for_prop: join tmp1.txt tmp2.txt 1>&2
140      fi      fi
# Line 155  dashnum() Line 181  dashnum()
181    
182  testoutput_ad()  testoutput_ad()
183  {  {
184      grep $3 $1/results_ad/output.txt_adm | awk '{print NR " " $5}' > t05.txt      grep $3 $1/results_ad/output_adm.txt | awk '{print NR " " $5}' > t05.txt
185      grep $3 $1/$2/output.txt_adm | awk '{print NR " " $5}' > t15.txt      grep $3 $1/$2/output_adm.txt | awk '{print NR " " $5}' > t15.txt
186      grep $3 $1/results_ad/output.txt_adm | awk '{print NR " " $6}' > t06.txt      grep $3 $1/results_ad/output_adm.txt | awk '{print NR " " $6}' > t06.txt
187      grep $3 $1/$2/output.txt_adm | awk '{print NR " " $6}' > t16.txt      grep $3 $1/$2/output_adm.txt | awk '{print NR " " $6}' > t16.txt
188      join t05.txt t15.txt > t5.txt      join t05.txt t15.txt > t5.txt
189      join t06.txt t16.txt > t6.txt      join t06.txt t16.txt > t6.txt
190      echo "-1" >> t5.txt      echo "-1" >> t5.txt
# Line 169  testoutput_ad() Line 195  testoutput_ad()
195      rm -f t[01][56].txt t[56].txt      rm -f t[01][56].txt t[56].txt
196  }  }
197    
198    check_for_add_mon_output()
199    {
200        # Check for additional types of monitor output
201        if test "x$1" = x ; then
202            return
203        fi
204    
205        for ii in $PTRACERS_NUM ; do
206            eval "HAVE_PTR0"$ii"=f"
207        done
208    
209        ptr_add="trcstat_ptracerXX_min trcstat_ptracerXX_max"
210        ptr_add="$ptr_add trcstat_ptracerXX_mean trcstat_ptracerXX_sd"
211        for ii in $PTRACERS_NUM ; do
212            for jj in $ptr_add ; do
213                name=`eval "echo $jj | sed -e 's|XX|0"$ii"|g'"`
214                tst=`grep $name $1 | wc -l | awk '{print $1}'`
215                if test ! "x$tst" = x0 ; then
216                    eval "HAVE_PTR0"$ii"=t"
217                fi
218            done
219            #  eval 'echo "HAVE_PTR0'$ii' = $HAVE_PTR0'$ii'"'
220        done
221    }
222    
223  testoutput()  testoutput()
224  {  {
225      # testoutput directory subdir      # testoutput directory subdir extension
226      #      #
227      #  test output in "directory"      #  test output in "directory"
228      if test "x$ADM" = x ; then      if test "x$ADM" = x ; then
229          if [ $debug -gt 0 ]; then          if [ $debug -gt 0 ]; then
230              echo testoutput: testoutput_for_prop $1 cg2d_init_res 1>&2              echo testoutput: testoutput_for_prop $1 cg2d_init_res 1>&2
231          fi          fi
232          testoutput_for_prop $1 "cg2d_init_res" "cg2d init. residual" $2; cg2dres=$?          testoutput_for_prop $1 "cg2d_init_res" "cg2d init. residual" $2 $3; cg2dres=$?
233          if [ $debug -gt 0 ]; then          if [ $debug -gt 0 ]; then
234              echo testoutput: cg2dres=$cg2dres 1>&2              echo testoutput: cg2dres=$cg2dres 1>&2
235          fi          fi
236          testoutput_for_prop $1 "dynstat_theta_min" "theta minimum" $2; tmin=$?          testoutput_for_prop $1 "dynstat_theta_min" "theta minimum" $2 $3; tmin=$?
237          testoutput_for_prop $1 "dynstat_theta_max" "theta maximum" $2; tmax=$?          testoutput_for_prop $1 "dynstat_theta_max" "theta maximum" $2 $3; tmax=$?
238          testoutput_for_prop $1 "dynstat_theta_mean" "theta mean" $2; tmean=$?          testoutput_for_prop $1 "dynstat_theta_mean" "theta mean" $2 $3; tmean=$?
239          testoutput_for_prop $1 "dynstat_theta_sd" "theta s.d." $2; tsd=$?          testoutput_for_prop $1 "dynstat_theta_sd" "theta s.d." $2 $3; tsd=$?
240          testoutput_for_prop $1 "dynstat_salt_min" "salt minimum" $2; smin=$?          testoutput_for_prop $1 "dynstat_salt_min" "salt minimum" $2 $3; smin=$?
241          testoutput_for_prop $1 "dynstat_salt_max" "salt maximum" $2; smax=$?          testoutput_for_prop $1 "dynstat_salt_max" "salt maximum" $2 $3; smax=$?
242          testoutput_for_prop $1 "dynstat_salt_mean" "salt mean" $2; smean=$?          testoutput_for_prop $1 "dynstat_salt_mean" "salt mean" $2 $3; smean=$?
243          testoutput_for_prop $1 "dynstat_salt_sd" "salt s.d." $2; ssd=$?          testoutput_for_prop $1 "dynstat_salt_sd" "salt s.d." $2 $3; ssd=$?
244          testoutput_for_prop $1 "dynstat_uvel_min" "U minimum" $2; umin=$?          testoutput_for_prop $1 "dynstat_uvel_min" "U minimum" $2 $3; umin=$?
245          testoutput_for_prop $1 "dynstat_uvel_max" "U maximum" $2; umax=$?          testoutput_for_prop $1 "dynstat_uvel_max" "U maximum" $2 $3; umax=$?
246          testoutput_for_prop $1 "dynstat_uvel_mean" "U mean" $2; umean=$?          testoutput_for_prop $1 "dynstat_uvel_mean" "U mean" $2 $3; umean=$?
247          testoutput_for_prop $1 "dynstat_uvel_sd" "U s.d." $2; usd=$?          testoutput_for_prop $1 "dynstat_uvel_sd" "U s.d." $2 $3; usd=$?
248          testoutput_for_prop $1 "dynstat_vvel_min" "V minimum" $2; vmin=$?          testoutput_for_prop $1 "dynstat_vvel_min" "V minimum" $2 $3; vmin=$?
249          testoutput_for_prop $1 "dynstat_vvel_max" "V maximum" $2; vmax=$?          testoutput_for_prop $1 "dynstat_vvel_max" "V maximum" $2 $3; vmax=$?
250          testoutput_for_prop $1 "dynstat_vvel_mean" "V mean" $2; vmean=$?          testoutput_for_prop $1 "dynstat_vvel_mean" "V mean" $2 $3; vmean=$?
251          testoutput_for_prop $1 "dynstat_vvel_sd" "V s.d." $2; vsd=$?          testoutput_for_prop $1 "dynstat_vvel_sd" "V s.d." $2 $3; vsd=$?
252          dashnum $cg2dres $tmin $tmax $tmean $tsd $smin $smax $smean $ssd \  
253              $umin $umax $umean $usd $vmin $vmax $vmean $vsd          #  This is for PTRACERS
254            for ii in $PTRACERS_NUM ; do
255                eval `echo "p0"$ii"_min=99"`
256                eval `echo "p0"$ii"_max=99"`
257                eval `echo "p0"$ii"_mean=99"`
258                eval `echo "p0"$ii"_sd=99"`
259                tst=`eval 'echo "$HAVE_PTR0'$ii'"'`
260                #echo 'tst = '$tst
261                if test "x$tst" = xt ; then
262                    a="trcstat_ptracer0"
263                    testoutput_for_prop $1 "$a"$ii"_min"  "p0"$ii"_min"  $2 $3
264                    RETVAL=$? ; eval `echo "p0"$ii"_min="$RETVAL`
265                    testoutput_for_prop $1 "$a"$ii"_max"  "p0"$ii"_max"  $2 $3
266                    RETVAL=$? ; eval `echo "p0"$ii"_max="$RETVAL`
267                    testoutput_for_prop $1 "$a"$ii"_mean" "p0"$ii"_mean" $2 $3
268                    RETVAL=$? ; eval `echo "p0"$ii"_mean="$RETVAL`
269                    testoutput_for_prop $1 "$a"$ii"_sd"   "p0"$ii"_sd"   $2 $3
270                    RETVAL=$? ; eval `echo "p0"$ii"_sd="$RETVAL`
271                fi
272            done
273    
274            allargs="$cg2dres $tmin $tmax $tmean $tsd $smin $smax $smean $ssd"
275            allargs="$allargs $umin $umax $umean $usd $vmin $vmax $vmean $vsd"
276            allargs="$allargs $p01_min $p01_max $p01_mean $p01_sd"
277            allargs="$allargs $p02_min $p02_max $p02_mean $p02_sd"
278            allargs="$allargs $p03_min $p03_max $p03_mean $p03_sd"
279            allargs="$allargs $p04_min $p04_max $p04_mean $p04_sd"
280            allargs="$allargs $p05_min $p05_max $p05_mean $p05_sd"
281    
282            eval "dashnum $allargs"
283    
284      else      else
285          testoutput_ad $1 $2 "precision_grdchk_result"          testoutput_ad $1 $2 "precision_grdchk_result"
286      fi      fi
# Line 211  genmakemodel() Line 292  genmakemodel()
292      if test "x$NOGENMAKE" = xt ; then      if test "x$NOGENMAKE" = xt ; then
293          echo "genmake skipped!"          echo "genmake skipped!"
294      else      else
295          GENMAKE2="$BASH ../../../tools/genmake2"          if test "x$BASH" = x ; then
296                GENMAKE2="../../../tools/genmake2"
297            else
298                GENMAKE2="$BASH ../../../tools/genmake2 -bash $BASH"
299            fi
300          (          (
301              cd $1;              cd $1;
302              command="$GENMAKE2  -ds -m $MAKE"              command="$GENMAKE2  -ds -m $MAKE"
# Line 219  genmakemodel() Line 304  genmakemodel()
304                  command="$command --mods=../code"                  command="$command --mods=../code"
305              else              else
306                  command="$command --mods=../code_ad"                  command="$command --mods=../code_ad"
                 command="$command -adof=../../../tools/adjoint_options/adjoint_staf"  
307              fi              fi
308              if test "x$OPTFILE" != xNONE ; then              if test "x$OPTFILE" != xNONE ; then
309                  command="$command --optfile=$OPTFILE"                  command="$command --optfile=$OPTFILE"
# Line 227  genmakemodel() Line 311  genmakemodel()
311              if test "x$IEEE" != x ; then              if test "x$IEEE" != x ; then
312                  command="$command -ieee"                  command="$command -ieee"
313              fi              fi
314                if test "x$MPI" = xt ; then
315                    command="$command -mpi"
316                fi
317              printf 'genmake ... ' 1>&2              printf 'genmake ... ' 1>&2
318              $command > make.log 2>&1              $command > make.log 2>&1
319              RETVAL=$?              RETVAL=$?
320              cp Makefile $CDIR              #  Reduce the size of the testing emails!
321                head -100 Makefile > $CDIR/Makefile_head
322              if test "x$RETVAL" != x0 ; then              if test "x$RETVAL" != x0 ; then
323                  tail make.log                  tail make.log
324                  echo "genmakemodel: genmake failed" 1>&2                  echo "genmakemodel: genmake failed" 1>&2
# Line 271  makeclean() Line 359  makeclean()
359      fi      fi
360  }  }
361    
362    postclean()
363    {
364        # postclean directory
365        if test "x$POSTCLEAN" = xt ; then
366            (
367                cd $1;
368                if test -r Makefile ; then
369                    $MAKE CLEAN >> /dev/null 2>&1
370                fi
371                exit 0
372            )
373        fi
374    }
375    
376  makedependmodel()  makedependmodel()
377  {  {
378      # makedependmodel directory      # makedependmodel directory
# Line 302  makemodel() Line 404  makemodel()
404          if test -r Makefile ; then          if test -r Makefile ; then
405              printf 'make ... ' 1>&2              printf 'make ... ' 1>&2
406              if test "x$ADM" = x ; then              if test "x$ADM" = x ; then
407                  $MAKE >> make.log 2>&1                  if test "x$JOBS" = x ; then
408                        $MAKE >> make.log 2>&1
409                    else
410                        $MAKE -j $JOBS >> make.log 2>&1
411                    fi
412              else              else
413                  $MAKE adall >> make.log 2>&1                  $MAKE adall >> make.log 2>&1
414              fi              fi
# Line 329  symlink_mpifiles() Line 435  symlink_mpifiles()
435      code_dir=$2      code_dir=$2
436      BUILD_DIR=$dir/$3      BUILD_DIR=$dir/$3
437      CODE_DIR=$dir/$code_dir      CODE_DIR=$dir/$code_dir
438        
439      # These are files that should replace their counter-part when using -mpi      # These are files that should replace their counter-part when using -mpi
440      MPI_FILES=`(cd $CODE_DIR; find . -name "*_mpi")`      MPI_FILES=`(cd $CODE_DIR; find . -name "*_mpi")`
441    
# Line 344  symlink_mpifiles() Line 450  symlink_mpifiles()
450              RETVAL=$?              RETVAL=$?
451              if test "x$RETVAL" != x0 ; then              if test "x$RETVAL" != x0 ; then
452                  if ! test -f $BUILD_DIR/$i ; then                  if ! test -f $BUILD_DIR/$i ; then
453                  #echo Linking $name to $i                      #echo Linking $name to $i
454                      (cd $BUILD_DIR; ln -sf ../$code_dir/$i $name)                      (cd $BUILD_DIR; ln -sf ../$code_dir/$i $name)
455                  fi                  fi
456              fi              fi
457          done          done
458      else      else
459      # NO: We undo any _mpi symbolically linked files          # NO: We undo any _mpi symbolically linked files
460          for ii in $MPI_FILES ; do          for ii in $MPI_FILES ; do
461              i=`echo $ii | sed 's:^\./::'`              i=`echo $ii | sed 's:^\./::'`
462              name=`echo $i | sed 's:_mpi::' `              name=`echo $i | sed 's:_mpi::' `
463              if test -L $BUILD_DIR/$name ; then              if test -L $BUILD_DIR/$name ; then
464                  linktarg=`(cd $BUILD_DIR; readlink $name)`                  cmp $BUILD_DIR/$name "../$code_dir/$name"_mpi > /dev/null 2>&1
465                  if test $linktarg = "../$code_dir/$name"_mpi ; then                  RETVAL=$?
466                  #echo Un-linking $name from $linktarg                  if test "x$RETVAL" = x0 ; then
467                        #echo Un-linking $name from $linktarg
468                      rm -f $BUILD_DIR/$name                      rm -f $BUILD_DIR/$name
469                  fi                  fi
470              fi              fi
# Line 377  linkdata() Line 484  linkdata()
484              if test "x$ADM" = x ; then              if test "x$ADM" = x ; then
485                  files=`( cd ../input ; ls -1 | grep -v CVS )`                  files=`( cd ../input ; ls -1 | grep -v CVS )`
486                  for i in $files ; do                  for i in $files ; do
487                      if test ! -d "../input/"$i ; then                      if test ! -d "../input/"$i -a ! -f $i ; then
488                          ln -sf "../input/"$i $i                          ln -sf "../input/"$i $i
489                      fi                      fi
490                  done                  done
491              else              else
                 files=`( cd ../input ; ls -1 *.bin | grep -v CVS )`  
                 for i in $files ; do  
                     if test ! -d "../input/"$i ; then  
                         ln -sf "../input/"$i $i  
                     fi  
                 done  
492                  files=`( cd ../input_ad ; ls -1 | grep -v CVS )`                  files=`( cd ../input_ad ; ls -1 | grep -v CVS )`
493                  for i in $files ; do                  for i in $files ; do
494                      if test ! -d "../input_ad/"$i ; then                      if test ! -d "../input_ad/"$i ; then
495                          ln -sf "../input_ad/"$i $i                          ln -sf "../input_ad/"$i $i
496                      fi                      fi
497                  done                  done
498                    files=`( cd ../input ; ls -1 | grep -v CVS )`
499                    for i in $files ; do
500                        if test ! -d "../input/"$i ; then
501                         if test ! -r $i  ; then
502                            ln -sf "../input/"$i $i
503                         fi
504                        fi
505                    done
506              fi              fi
507          )          )
508      fi      fi
# Line 407  runmodel() Line 516  runmodel()
516      #  (where "$COMMAND" is relative to "directory")      #  (where "$COMMAND" is relative to "directory")
517      (      (
518          cd $1          cd $1
519          printf 'runmodel ... ' 1>&2          printf 'runmodel in %s ...' $1 1>&2
520          # make output.txt          # make output.txt
521          $COMMAND >> run.log 2>&1          echo
522          RETVAL=$?          rm -f run.log
523            if test $OUTPUTFILE -ot $EXECUTABLE ; then
524               ( eval $COMMAND ) > run.log 2>&1
525               RETVAL=$?
526            else
527               echo " $OUTPUTFILE is up to date " > run.log 2>&1
528               RETVAL=0
529            fi
530            # echo "COMMAND='$COMMAND'"
531            # echo "pwd='"`pwd`"'"
532          if test "x$RETVAL" = x0 ; then          if test "x$RETVAL" = x0 ; then
533                tail run.log
534              echo successful 1>&2              echo successful 1>&2
535              if test "x$ADM" = x ; then              # === Reduce the size of the testing emails!
536                  cp output.txt $CDIR"/output.txt"              # if test "x$ADM" = x ; then
537              else              #   cp output.txt $CDIR"/output.txt"
538                  cp output.txt_adm $CDIR"/output.txt_adm"              # else
539              fi              #   cp output_adm.txt $CDIR"/output_adm.txt"
540                # fi
541                if test -s STDERR.0000 ; then cp STDERR.0000 $CDIR"/STDERR.0000" ; fi
542              return 0              return 0
543          else          else
544              tail run.log              tail run.log
545              echo failed 1>&2              echo failed 1>&2
546              cp run.log $CDIR"/run.log"              cp run.log $CDIR"/run.log"
547                if test -s STDERR.0000 ; then cp STDERR.0000 $CDIR"/STDERR.0000" ; fi
548              return 1              return 1
549          fi          fi
550      )      )
# Line 432  createcodelet() Line 554  createcodelet()
554  {  {
555      # create codelet for comparing model output      # create codelet for comparing model output
556    
557      echo -n "creating the comparison code...  "      printf "creating the comparison code...  "
558      cat > tmp_cmpnum.c <<EOF      cat > tmp_cmpnum.c <<EOF
559  #include <stdio.h>  #include <stdio.h>
560  #include <math.h>  #include <math.h>
561  int main( int argc, char** argv )  {  int main( int argc, char** argv )  {
562    int linnum,best;    int linnum,best,lncnt;
563    double a,b,diff;    double a,b,abave,relerr;
564    best = -16;    best = -22;
565    while( 1 )  {    lncnt = 0;
566      while( 1 & (lncnt+=1) < 999 )  {
567      scanf("%d", &linnum);      scanf("%d", &linnum);
568      if (linnum == -1)  break;      if (linnum == -1)  break;
569      scanf("%lf", &a);  scanf("%lf", &b);      scanf("%lf", &a);  scanf("%lf", &b);
570      diff = 0.5*(fabs(a)+fabs(b));      abave = 0.5*(fabs(a)+fabs(b));
571      if (diff > 1.e-12) {      if (abave > 0.0) {
572        diff=fabs(a-b)/diff;        relerr=fabs(a-b)/abave;
573        if (diff > 0.0) {        if (relerr > 0.0) { linnum = (int)rint(log10(relerr)); }
574          linnum = (int)log10(diff);        else { linnum = -16 ; }
575          best = (best > linnum) ? best : linnum;        best = (best > linnum) ? best : linnum;
       }  
       else {  
         if (best == -16 && diff != 0)  best = -22;  
       }  
576      }      }
577    }    }
578      if (lncnt == 999) best=-29;
579    printf("%d\n", -best);    printf("%d\n", -best);
580    return 0;    return 0;
581  }  }
582  EOF  EOF
583      cc -o tmp_cmpnum tmp_cmpnum.c -lm      $CC -o tmp_cmpnum tmp_cmpnum.c -lm
584    
585      if [ -x ./tmp_cmpnum ]; then      if [ -x ./tmp_cmpnum ]; then
586          echo "OK"          echo "OK"
587          return 0          return 0
588      else      else
589          echo          echo
590          echo "ERROR: failed to compile comparison code"          echo "ERROR: failed to compile comparison code -- please specify"
591            echo "  a C compiler using the CC environment variable."
592          exit 1          exit 1
593      fi      fi
594  }  }
# Line 536  clean=0 Line 657  clean=0
657  expts=''  expts=''
658  # ieee=1  # ieee=1
659    
660  IEEE=  IEEE=true
661  if test "x$MITGCM_IEEE" != x ; then  if test "x$MITGCM_IEEE" != x ; then
662      IEEE=$MITGCM_IEEE      IEEE=$MITGCM_IEEE
663  fi  fi
# Line 547  QUICK=f Line 668  QUICK=f
668  NOGENMAKE=f  NOGENMAKE=f
669  NOCLEAN=f  NOCLEAN=f
670  NODEPEND=f  NODEPEND=f
671    POSTCLEAN=f
672    
673  BASH=  BASH=
674  OPTFILE=NONE  OPTFILE=NONE
# Line 556  MPACKDIR="../tools/mpack-1.6" Line 678  MPACKDIR="../tools/mpack-1.6"
678  HAVE_MPACK=  HAVE_MPACK=
679  MPACK="$MPACKDIR/mpack"  MPACK="$MPACKDIR/mpack"
680  COMMAND=  COMMAND=
681  MAKE=make  if test "x$MAKE" = x ; then
682        MAKE=make
683    fi
684    if test "x$CC" = x ; then
685        CC=cc
686    fi
687    JOBS=
688  MPI=f  MPI=f
689    OUTDIR=
690    DELDIR=
691    
692  ADM=  ADM=
693    
694  echo -n "parsing options...  "  # Additional monitor types
695    PTRACERS_NUM="1 2 3 4 5"
696    
697    printf "parsing options...  "
698    
699  ac_prev=  ac_prev=
700  for ac_option ; do  for ac_option ; do
# Line 610  for ac_option ; do Line 743  for ac_option ; do
743          -make=* | --make=*)          -make=* | --make=*)
744              MAKE=$ac_optarg ;;              MAKE=$ac_optarg ;;
745    
746            -odir | --odir)
747                ac_prev=OUTDIR ;;
748            -odir=* | --odir=*)
749                OUTDIR=$ac_optarg ;;
750    
751            -ptracers | --ptracers | -ptr | --ptr)
752                ac_prev=PTRACERS_NUM ;;
753            -ptracers=* | --ptracers=* | -ptr=* | --ptr=*)
754                PTRACERS_NUM=$ac_optarg ;;
755    
756            -j) ac_prev=JOBS ;;
757            -j=*) JOBS=$ac_optarg ;;
758    
759          -clean | --clean)          -clean | --clean)
760              CLEANUP=t ;;              CLEANUP=t ; DELDIR=t ;;
761    
762          -quick | --quick | -q | --q)          -quick | --quick | -q | --q)
763              QUICK=t ;;              QUICK=t ;;
# Line 622  for ac_option ; do Line 768  for ac_option ; do
768          -nodepend | --nodepend | -nd | --nd)          -nodepend | --nodepend | -nd | --nd)
769              NODEPEND=t ;;              NODEPEND=t ;;
770    
771            -postclean | --postclean | -pc | --pc)
772                POSTCLEAN=t ;;
773    
774          -mpi) MPI=t ;;          -mpi) MPI=t ;;
775    
776          -adm | -ad) ADM=t ;;          -adm | -ad) ADM=t ;;
# Line 633  for ac_option ; do Line 782  for ac_option ; do
782          -debug) debug=1 ;;          -debug) debug=1 ;;
783          -quiet) verbose=0 ;;          -quiet) verbose=0 ;;
784    
785            -deldir | -dd) DELDIR=t ;;
786    
787          -*)          -*)
788              echo "Error: unrecognized option: "$ac_option              echo "Error: unrecognized option: "$ac_option
789              usage              usage
# Line 661  if test "x$OPTFILE" = xNONE -a "x$MITGCM Line 812  if test "x$OPTFILE" = xNONE -a "x$MITGCM
812      OPTFILE=$MITGCM_OF      OPTFILE=$MITGCM_OF
813  fi  fi
814    
815    if test "x$ADM" = xt ; then
816        EXECUTABLE="mitgcmuv_ad"
817        OUTPUTFILE="output_adm.txt"
818    else
819        EXECUTABLE="mitgcmuv"
820        OUTPUTFILE="output.txt"
821    fi
822    
823  if test "x$ADM" = xt -a "x$COMMAND" = x ; then  if test "x$ADM" = xt -a "x$COMMAND" = x ; then
824      COMMAND="./mitgcmuv_ad > output.txt_adm 2>&1"      COMMAND="./$EXECUTABLE > $OUTPUTFILE"
825  fi  fi
826    
827  if test "x$COMMAND" = x ; then  if test "x$COMMAND" = x ; then
828      COMMAND="make output.txt"      COMMAND="./$EXECUTABLE > $OUTPUTFILE"
829  fi  fi
830    
831  echo "OK"  #echo "OK"
832    echo "OK (COMMAND= $COMMAND )"
833    
834  #  create the FORTRAN comparison code  #  create the FORTRAN comparison code
835  createcodelet  createcodelet
# Line 686  MACH=`hostname` Line 846  MACH=`hostname`
846  UNAMEA=`uname -a`  UNAMEA=`uname -a`
847  DATE=`date +%Y%m%d`  DATE=`date +%Y%m%d`
848  BASE="tr_"$MACH"_"$DATE"_"  BASE="tr_"$MACH"_"$DATE"_"
849    if test "x$OUTDIR" != x ; then
850        BASE="tr_"$OUTDIR"_"$DATE"_"
851    fi
852  DNUM=0  DNUM=0
853  DRESULTS="$BASE$DNUM"  DRESULTS="$BASE$DNUM"
854  while test -e $DRESULTS ; do  while test -e $DRESULTS ; do
# Line 699  if test "x$RETVAL" != x0 ; then Line 862  if test "x$RETVAL" != x0 ; then
862      exit 1      exit 1
863  fi  fi
864  SUMMARY="$DRESULTS/summary.txt"  SUMMARY="$DRESULTS/summary.txt"
865  echo -n "Start time:  " >> $SUMMARY  printf "Start time:  " >> $SUMMARY
866  start_date=`date`  start_date=`date`
867  echo $start_date > $SUMMARY  echo $start_date > $SUMMARY
868    
# Line 731  fi Line 894  fi
894  echo  echo
895  echo >> $SUMMARY  echo >> $SUMMARY
896  if test "x$ADM" = x ; then  if test "x$ADM" = x ; then
897      cat << EOF | tee -a $SUMMARY      line_0="            ----T-----  ----S-----  ----U-----  ----V-----"
898                  T           S           U           V      line_1="G D M    c        m  s        m  s        m  s        m  s"
899  G D M    c        m  s        m  s        m  s        m  s      line_2="E p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  ."
900  E p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .      line_3="N n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d"
901  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  ."
902  2 d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .      for ii in $PTRACERS_NUM ; do
903            #  tst=`eval 'echo $HAVE_PTR0'$ii`
904  EOF          #  if test "x$tst" = xt ; then
905            line_0="$line_0  --PTR 0"$ii"--"
906            line_1="$line_1        m  s"
907            line_2="$line_2  m  m  e  ."
908            line_3="$line_3  i  a  a  d"
909            line_4="$line_4  n  x  n  ."
910            #  fi
911        done
912        echo "$line_0" | tee -a $SUMMARY
913        echo "$line_1" | tee -a $SUMMARY
914        echo "$line_2" | tee -a $SUMMARY
915        echo "$line_3" | tee -a $SUMMARY
916        echo "$line_4" | tee -a $SUMMARY
917        echo " "       | tee -a $SUMMARY
918  else  else
919      echo "ADJOINT=true" >> $SUMMARY      echo "ADJOINT=true" >> $SUMMARY
920      echo >> $SUMMARY      echo >> $SUMMARY
# Line 762  for dir in $TESTDIRS ; do Line 938  for dir in $TESTDIRS ; do
938          if test -r $dir/input/Makefile ; then          if test -r $dir/input/Makefile ; then
939              ( cd $dir/input ; make CLEAN )              ( cd $dir/input ; make CLEAN )
940          fi          fi
941            (
942                cd $dir
943                rm -rf tr_run.*
944            )
945          continue          continue
946      fi      fi
947    
# Line 771  for dir in $TESTDIRS ; do Line 951  for dir in $TESTDIRS ; do
951      if test "x$ADM" = x ; then      if test "x$ADM" = x ; then
952          fout=$dir"/results/output.txt"          fout=$dir"/results/output.txt"
953      else      else
954          fout=$dir"/results_ad/output.txt_adm"          fout=$dir"/results_ad/output_adm.txt"
955      fi      fi
956      if test ! -r $fout ; then      if test ! -r $fout ; then
957          echo "can't read \"$fout\" -- skipping $dir"          echo "can't read \"$fout\" -- skipping $dir"
958          continue          continue
959      fi      fi
960        if test "x$ADM" = x ; then
961            check_for_add_mon_output  $fout
962        fi
963    
964        # Check for additional types of monitor output
965    
966      builddir="input"      builddir="input"
967      rundir="input"      rundir="input"
# Line 788  for dir in $TESTDIRS ; do Line 973  for dir in $TESTDIRS ; do
973          linkdata $use_seperate_build $dir/$rundir          linkdata $use_seperate_build $dir/$rundir
974      fi      fi
975            
976        #  Check whether there are "extra runs" for this testdir
977        extra_runs=
978        if test "x$ADM" = x -a "x$use_seperate_build" = x1 ; then
979            ex_run_dirs=`( cd $dir ; echo input.* )`
980            echo "ex_run_dirs='$ex_run_dirs'"
981            for exd in $ex_run_dirs ; do
982                name=`echo $exd | sed -e 's/input.//g'`
983                outf="$dir/results/output.txt.$name"
984                if test -f $outf -a -r $outf ; then
985                    extra_runs="$extra_runs $name"
986                fi
987            done
988        fi
989    
990      if test "x$ADM" = x ; then      if test "x$ADM" = x ; then
991          code_dir=code          code_dir=code
992          CODE_DIR=$dir/code          CODE_DIR=$dir/code
# Line 807  for dir in $TESTDIRS ; do Line 1006  for dir in $TESTDIRS ; do
1006      echo "Experiment:  $dir"      echo "Experiment:  $dir"
1007      echo      echo
1008      unset genmake makedepend make run      unset genmake makedepend make run
1009      results='-- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --'      results='-- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --'
1010    
1011      #  Create an output dir for each OPTFILE/tdir combination      #  Create an output dir for each OPTFILE/tdir combination
1012      rel_CDIR=$DRESULTS"/"$dir      rel_CDIR=$DRESULTS"/"$dir
# Line 830  for dir in $TESTDIRS ; do Line 1029  for dir in $TESTDIRS ; do
1029      echo      echo
1030      if test "x$ADM" = x ; then      if test "x$ADM" = x ; then
1031          fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results`          fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results`
1032            echo
1033            echo "$fres" >> $SUMMARY
1034            echo "fresults='$fres'" > $CDIR"/summary.txt"
1035            echo "MACH='$MACH'" >> $CDIR"/summary.txt"
1036            echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt"
1037            echo "DATE='$DATE'" >> $CDIR"/summary.txt"
1038            echo "tdir='$dir'" >> $CDIR"/summary.txt"
1039    
1040            for ex in $extra_runs ; do
1041                unset run
1042                results='-- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --'
1043                #  Create an output dir for each OPTFILE/tdir.ex combination
1044                rel_CDIR=$DRESULTS"/"$dir"."$ex
1045                mkdir $rel_CDIR
1046                CDIR=`pwd`"/$rel_CDIR"
1047                test ! -e "$dir/tr_run.$ex" && mkdir "$dir/tr_run.$ex"
1048                for ldir in input.$ex input ; do
1049                    (
1050                        cd "$dir/$ldir" > /dev/null 2>&1
1051                        ls -1 2>/dev/null \
1052                            | sed -e 's|^CVS$||g' | sed -e 's|^output.txt$||g'
1053                    ) > tr_exrun_links
1054                    (
1055                        cd "$dir/tr_run.$ex"
1056                        cat ../../tr_exrun_links | while read i ; do
1057                            if test ! "x$i" = x ; then
1058                                test ! -r $i  &&  ln -s "../"$ldir"/"$i $i
1059                            fi
1060                        done
1061                    )
1062                    test -e tr_exrun_links  &&  rm -f tr_exrun_links
1063                done
1064                ldir=build
1065                (
1066                    cd $dir/tr_run.$ex
1067                    test ! -e mitgcmuv  &&  ln -s "../"$ldir"/"mitgcmuv .
1068                )
1069                runmodel $dir/tr_run.$ex && run=Y \
1070                    && results=`testoutput $dir tr_run.$ex "."$ex`
1071                fres=`printf '%s %s %s %s' ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N}`
1072                fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results`
1073                fres="$fres.$ex"
1074                echo
1075                echo "$fres" >> $SUMMARY
1076                echo "fresults='$fres'" > $CDIR"/summary.txt"
1077                echo "MACH='$MACH'" >> $CDIR"/summary.txt"
1078                echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt"
1079                echo "DATE='$DATE'" >> $CDIR"/summary.txt"
1080                echo "tdir='$dir.$ex'" >> $CDIR"/summary.txt"
1081            done
1082      else      else
1083          fres=`printf '%s %s %s %s' ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N}`          fres=`printf '%s %s %s %s' ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N}`
1084          fres=$fres"$results   $dir"          fres=$fres"$results   $dir"
1085            echo
1086            echo "$fres" >> $SUMMARY
1087            echo "fresults='$fres'" > $CDIR"/summary.txt"
1088            echo "MACH='$MACH'" >> $CDIR"/summary.txt"
1089            echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt"
1090            echo "DATE='$DATE'" >> $CDIR"/summary.txt"
1091            echo "tdir='$dir'" >> $CDIR"/summary.txt"
1092            grep -A3 'Seconds in section "ALL' $dir/$rundir/$OUTPUTFILE \
1093                               >> $CDIR"/summary.txt"
1094      fi      fi
1095      echo  
1096      echo "$fres" >> $SUMMARY      postclean $dir/$builddir
     echo "fresults='$fres'" > $CDIR"/summary.txt"  
     echo "MACH='$MACH'" >> $CDIR"/summary.txt"  
     echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt"  
     echo "DATE='$DATE'" >> $CDIR"/summary.txt"  
     echo "tdir='$dir'" >> $CDIR"/summary.txt"  
1097            
1098      echo "-------------------------------------------------------------------------------"      echo "-------------------------------------------------------------------------------"
1099            
1100  done  done
1101    
1102  echo -n "Start time:  " >> $SUMMARY  printf "Start time:  " >> $SUMMARY
1103  echo $start_date >> $SUMMARY  echo $start_date >> $SUMMARY
1104  echo -n "End time:    " >> $SUMMARY  printf "End time:    " >> $SUMMARY
1105  date >> $SUMMARY  date >> $SUMMARY
1106    
1107  #  If addresses were supplied and mpack built successfully, then try  #  If addresses were supplied and mpack built successfully, then try
# Line 859  else Line 1112  else
1112      if test "x$HAVE_MPACK" = xt ; then      if test "x$HAVE_MPACK" = xt ; then
1113          tar -cf $DRESULTS".tar" $DRESULTS > /dev/null 2>&1 \          tar -cf $DRESULTS".tar" $DRESULTS > /dev/null 2>&1 \
1114              && gzip $DRESULTS".tar" \              && gzip $DRESULTS".tar" \
1115              && $MPACK -s MITgcm-test -m 1500000 $DRESULTS".tar.gz" $ADDRESSES              && $MPACK -s MITgcm-test -m 3555000 $DRESULTS".tar.gz" $ADDRESSES
1116          RETVAL=$?          RETVAL=$?
1117          if test "x$RETVAL" != x0 ; then          if test "x$RETVAL" != x0 ; then
1118              echo              echo
# Line 882  fi Line 1135  fi
1135  rm -f tmp_cmpnum.c tmp_cmpnum  rm -f tmp_cmpnum.c tmp_cmpnum
1136    
1137  if test "x$CLEANUP" != xt ; then  if test "x$CLEANUP" != xt ; then
1138      cat $SUMMARY      cat $SUMMARY | sed 's/ -- -- -- --//g'
1139      if test -e tr_out.txt ; then      if test -e tr_out.txt ; then
1140          mv tr_out.txt tr_out.txt.old          mv tr_out.txt tr_out.txt.old
1141      fi      fi
1142      cat $SUMMARY > tr_out.txt      cat $SUMMARY | sed 's/ -- -- -- --//g' > tr_out.txt
1143    fi
1144    
1145    if test "x$DELDIR" = xt ; then
1146        rm -rf $DRESULTS
1147  fi  fi
1148    

Legend:
Removed from v.1.32  
changed lines
  Added in v.1.80

  ViewVC Help
Powered by ViewVC 1.1.22