/[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.93 by jmc, Wed Sep 5 16:09:31 2007 UTC revision 1.94 by jmc, Thu Sep 6 22:31:45 2007 UTC
# Line 236  testoutput_run() Line 236  testoutput_run()
236      #      #
237      #  test output from 1 run in "directory"      #  test output from 1 run in "directory"
238      if test "x$ADM" = x ; then      if test "x$ADM" = x ; then
239          if [ $debug -gt 0 ]; then          # default list of output variables to be checked:
240              echo testoutput_run: testoutput_var $1 cg2d_init_res 1>&2          #  1rst : main variable used to decide if it pass or FAIL
241          fi          #  others : number of matching digits to be printed in summary.txt
242          testoutput_var $1 "cg2d_init_res" "cg2d init. residual" $2 $3; cg2dres=$?          listChk=$DEF_CHECK_LIST
243          if [ $debug -gt 0 ]; then          nbDef=`echo $listChk | awk '{print NF-1+3*gsub("+","E")}'`
244              echo testoutput_run: cg2dres=$cg2dres 1>&2          #  load experiment-specific list from file "tr_checklist" (if it exist)
245          fi          if test -r $1/$2/tr_checklist ; then listChk=`cat $1/$2/tr_checklist` ; fi
246          testoutput_var $1 "dynstat_theta_min" "theta minimum" $2 $3; tmin=$?          sVar=`echo $listChk | awk '{print $1}'`
247          testoutput_var $1 "dynstat_theta_max" "theta maximum" $2 $3; tmax=$?          # remove 1rst var and expand the list: + => min max mean s.d
248          testoutput_var $1 "dynstat_theta_mean" "theta mean" $2 $3; tmean=$?          listVar=`echo $listChk | awk '{ for(i=2;i<=NF;i++){if (sub("+","",$i)) \
249          testoutput_var $1 "dynstat_theta_sd" "theta s.d." $2 $3; tsd=$?              {printf " %s %s %s %s",$i"mn",$i"mx",$i"av",$i"sd"} else {printf " %s",$i} }}'`
250          testoutput_var $1 "dynstat_salt_min" "salt minimum" $2 $3; smin=$?          if [ $debug -gt 0 ]; then echo "testoutput_run: listVar(I)='$listVar'" 1>&2 ; fi
251          testoutput_var $1 "dynstat_salt_max" "salt maximum" $2 $3; smax=$?          for ii in 1 2 3 4 5 6 7 8 9 ; do
252          testoutput_var $1 "dynstat_salt_mean" "salt mean" $2 $3; smean=$?            tst=`eval 'echo "$HAVE_PTR0'$ii'"'`
253          testoutput_var $1 "dynstat_salt_sd" "salt s.d." $2 $3; ssd=$?           #echo "-- ptr test=" $tst "number of var=" `echo $listVar | awk '{print NF}'` 1>&2
254          testoutput_var $1 "dynstat_uvel_min" "U minimum" $2 $3; umin=$?            if test "x$tst" != xt ; then listVar=`echo "$listVar" | sed "s/ pt$ii..//g"` ; fi
255          testoutput_var $1 "dynstat_uvel_max" "U maximum" $2 $3; umax=$?          done
256          testoutput_var $1 "dynstat_uvel_mean" "U mean" $2 $3; umean=$?          tst=`echo $sVar $listVar | awk '{ for(i=2;i<=NF;i++){t+=($i==$1)}; print t }'`
257          testoutput_var $1 "dynstat_uvel_sd" "U s.d." $2 $3; usd=$?          if test $tst != 1 ; then
258          testoutput_var $1 "dynstat_vvel_min" "V minimum" $2 $3; vmin=$?            if test $tst = 0 ; then echo "==> WARNING: selected var >$sVar< not found" 1>&2
259          testoutput_var $1 "dynstat_vvel_max" "V maximum" $2 $3; vmax=$?                   else echo "==> WARNING: found selected var >$sVar< $tst times" 1>&2 ; fi
260          testoutput_var $1 "dynstat_vvel_mean" "V mean" $2 $3; vmean=$?            echo "==> WARNING: in checked list:" $listVar 1>&2
261          testoutput_var $1 "dynstat_vvel_sd" "V s.d." $2 $3; vsd=$?          #- put it back once:
262              listVar=" $sVar "`echo "$listVar " | sed "s/ $sVar / /g"`
263          #  This is for PTRACERS          fi
264          for ii in $PTRACERS_NUM ; do          if [ $debug -gt 0 ]; then echo "testoutput_run: listVar(M)='$listVar'" 1>&2 ; fi
265              eval `echo "p0"$ii"_min=99"`          echo "listVar='$listVar'" > $CDIR"/summary.txt"
266              eval `echo "p0"$ii"_max=99"`          allargs=""
267              eval `echo "p0"$ii"_mean=99"`          for xx in $listVar
268              eval `echo "p0"$ii"_sd=99"`          do
269              tst=`eval 'echo "$HAVE_PTR0'$ii'"'`            case $xx in
270              #echo 'tst = '$tst             'PS')  if [ $debug -gt 0 ]
271              if test "x$tst" = xt ; then                    then echo testoutput_run: testoutput_var $1 cg2d_init_res 1>&2 ; fi
272                  a="trcstat_ptracer0"                    testoutput_var $1 "cg2d_init_res" "cg2d init. residual" $2 $3; yy=$?
273                  testoutput_var $1 "$a"$ii"_min"  "p0"$ii"_min"  $2 $3                    if [ $debug -gt 0 ] ; then echo testoutput_run: cg2dres=$yy 1>&2 ; fi ;;
274                  RETVAL=$? ; eval `echo "p0"$ii"_min="$RETVAL`             'Tmn') testoutput_var $1 "dynstat_theta_min"  "theta minimum"  $2 $3; yy=$? ;;
275                  testoutput_var $1 "$a"$ii"_max"  "p0"$ii"_max"  $2 $3             'Tmx') testoutput_var $1 "dynstat_theta_max"  "theta maximum"  $2 $3; yy=$? ;;
276                  RETVAL=$? ; eval `echo "p0"$ii"_max="$RETVAL`             'Tav') testoutput_var $1 "dynstat_theta_mean" "theta mean"     $2 $3; yy=$? ;;
277                  testoutput_var $1 "$a"$ii"_mean" "p0"$ii"_mean" $2 $3             'Tsd') testoutput_var $1 "dynstat_theta_sd"   "theta s.d."     $2 $3; yy=$? ;;
278                  RETVAL=$? ; eval `echo "p0"$ii"_mean="$RETVAL`             'Smn') testoutput_var $1 "dynstat_salt_min"  "salt minimum"    $2 $3; yy=$? ;;
279                  testoutput_var $1 "$a"$ii"_sd"   "p0"$ii"_sd"   $2 $3             'Smx') testoutput_var $1 "dynstat_salt_max"  "salt maximum"    $2 $3; yy=$? ;;
280                  RETVAL=$? ; eval `echo "p0"$ii"_sd="$RETVAL`             'Sav') testoutput_var $1 "dynstat_salt_mean" "salt mean"       $2 $3; yy=$? ;;
281              fi             'Ssd') testoutput_var $1 "dynstat_salt_sd"   "salt s.d."       $2 $3; yy=$? ;;
282          done             'Umn') testoutput_var $1 "dynstat_uvel_min"  "U minimum"       $2 $3; yy=$? ;;
283               'Umx') testoutput_var $1 "dynstat_uvel_max"  "U maximum"       $2 $3; yy=$? ;;
284          allargs="$cg2dres $tmin $tmax $tmean $tsd $smin $smax $smean $ssd"             'Uav') testoutput_var $1 "dynstat_uvel_mean" "U mean"          $2 $3; yy=$? ;;
285          allargs="$allargs $umin $umax $umean $usd $vmin $vmax $vmean $vsd"             'Usd') testoutput_var $1 "dynstat_uvel_sd"   "U s.d."          $2 $3; yy=$? ;;
286          allargs="$allargs $p01_min $p01_max $p01_mean $p01_sd"             'Vmn') testoutput_var $1 "dynstat_vvel_min"  "V minimum"       $2 $3; yy=$? ;;
287          allargs="$allargs $p02_min $p02_max $p02_mean $p02_sd"             'Vmx') testoutput_var $1 "dynstat_vvel_max"  "V maximum"       $2 $3; yy=$? ;;
288          allargs="$allargs $p03_min $p03_max $p03_mean $p03_sd"             'Vav') testoutput_var $1 "dynstat_vvel_mean" "V mean"          $2 $3; yy=$? ;;
289          allargs="$allargs $p04_min $p04_max $p04_mean $p04_sd"             'Vsd') testoutput_var $1 "dynstat_vvel_sd"   "V s.d."          $2 $3; yy=$? ;;
290          allargs="$allargs $p05_min $p05_max $p05_mean $p05_sd"          'pt1mn'|'pt2mn'|'pt3mn'|'pt4mn'|'pt5mn') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'`
291               testoutput_var $1 "trcstat_ptracer0"$ii"_min"  "p0"$ii"_min"   $2 $3; yy=$? ;;
292          eval "dashnum $allargs"          'pt1mx'|'pt2mx'|'pt3mx'|'pt4mx'|'pt5mx') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'`
293               testoutput_var $1 "trcstat_ptracer0"$ii"_max"  "p0"$ii"_max"   $2 $3; yy=$? ;;
294            'pt1av'|'pt2av'|'pt3av'|'pt4av'|'pt5av') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'`
295               testoutput_var $1 "trcstat_ptracer0"$ii"_mean" "p0"$ii"_mean" $2 $3; yy=$? ;;
296            'pt1sd'|'pt2sd'|'pt3sd'|'pt4sd'|'pt5sd') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'`
297               testoutput_var $1 "trcstat_ptracer0"$ii"_sd"   "p0"$ii"_sd"    $2 $3; yy=$? ;;
298                  *) yy=99; echo "WARNING: asking for var=$xx : not recognized !" 1>&2 ;;
299              esac
300              if test $xx = $sVar
301              then allargs="$allargs > $yy <"
302              else allargs="$allargs $yy"
303              fi
304            done
305    
306            nbVar=`echo $listVar | awk '{print NF}'`
307            if [ $nbVar -lt $nbDef ] ; then
308            #-- fill line (up to standard length) with dot:
309              adNul=`expr $nbDef - $nbVar | awk '{for(i=1;i<=$1;i++){print "."}}'`
310              echo $allargs $adNul
311            else
312              echo $allargs
313            fi
314    
315      else      else
316          testoutput_ad $1 $2 "precision_grdchk_result"          testoutput_ad $1 $2 "precision_grdchk_result"
# Line 660  formatresults() Line 681  formatresults()
681      nm=$1      nm=$1
682      printf '%s %s %s %s' $2 $3 $4 $5      printf '%s %s %s %s' $2 $3 $4 $5
683      shift; shift; shift; shift; shift;      shift; shift; shift; shift; shift;
684      printf '%3s' $@      listPrt=$@
685        listRes=`echo $listPrt | sed 's/>//' | sed 's/<//'`
686        xx=`echo $listPrt | sed 's/.*>//' | sed 's/<.*//' | awk '{print $1}'`
687        printf '%3s' $listPrt | sed 's/ 99/ --/g' | sed 's/  > />/' | sed 's/  < /</'
688    #   printf '%3s' $listPrt | sed 's/ 99/ --/g' | sed 's/  > / /' | sed 's/  < / /'
689            
690      if [ $1 = '--' ]; then      if [ $xx = '..' ]; then
691            printf ' N/O '
692        elif [ $xx = '--' ]; then
693            printf ' N/O '
694        elif [ $xx = 99 ]; then
695          printf ' N/O '          printf ' N/O '
696      else      else
697          if [ $1 -ge $MATCH_CRIT ]; then          if [ $xx -ge $MATCH_CRIT ]; then
698              printf ' pass'              printf ' pass'
699          else          else
700              printf ' FAIL'              printf ' FAIL'
# Line 932  fi Line 961  fi
961  #echo "OK"  #echo "OK"
962  echo "OK (COMMAND= $COMMAND )"  echo "OK (COMMAND= $COMMAND )"
963    
964    # set the Default List of output variables to be checked:
965    #  (use default or load experiment-specific list from file "tr_checklist")
966    # content : 1rst = main variable used to decide if it pass or FAIL
967    #         others = number of matching digits to be printed in summary.txt
968    if test "x$ADM" != xt ; then
969        DEF_CHECK_LIST='PS PS T+ S+ U+ V+ pt1+ pt2+ pt3+ pt4+ pt5+'
970        EMPTY_RESULTS='.. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. ..'
971        nbDef=`echo $DEF_CHECK_LIST | awk '{print NF-1+3*gsub("+","E")}'`
972        ii=`echo $EMPTY_RESULTS | awk '{print NF}'`
973        EMPTY_RESULTS=$EMPTY_RESULTS`expr $nbDef - $ii | awk 'BEGIN{FS=":"}{for(i=1;i<=$1;i++){printf "  ."}}'`
974    fi
975    
976  #  create the FORTRAN comparison code  #  create the FORTRAN comparison code
977  createcodelet  createcodelet
978    
# Line 996  echo Line 1037  echo
1037  echo >> $SUMMARY  echo >> $SUMMARY
1038  if test "x$ADM" = x ; then  if test "x$ADM" = x ; then
1039      if [ $MATCH_CRIT -lt 10 ] ;      if [ $MATCH_CRIT -lt 10 ] ;
1040      then line_0="      "$MATCH_CRIT ;      then line_0="default  "$MATCH_CRIT ;
1041      else line_0="     "$MATCH_CRIT ; fi      else line_0="default "$MATCH_CRIT ; fi
1042      line_0="   $line_0  ----T-----  ----S-----  ----U-----  ----V-----"         line_0="$line_0  ----T-----  ----S-----  ----U-----  ----V-----"
1043  #   line_0="            ----T-----  ----S-----  ----U-----  ----V-----"  #   line_0="            ----T-----  ----S-----  ----U-----  ----V-----"
1044      line_1="G D M    c        m  s        m  s        m  s        m  s"      line_1="G D M    c        m  s        m  s        m  s        m  s"
1045      line_2="E p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  ."      line_2="E p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  ."
# Line 1128  for dir in $TESTDIRS ; do Line 1169  for dir in $TESTDIRS ; do
1169      fi      fi
1170      echo      echo
1171      unset genmake makedepend make run      unset genmake makedepend make run
1172      results='-- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --'      results=$EMPTY_RESULTS
1173    
1174      #  Create an output dir for each OPTFILE/tdir combination      #  Create an output dir for each OPTFILE/tdir combination
1175      rel_CDIR=$DRESULTS"/"$dir      rel_CDIR=$DRESULTS"/"$dir
# Line 1156  for dir in $TESTDIRS ; do Line 1197  for dir in $TESTDIRS ; do
1197          fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results`          fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results`
1198          echo          echo
1199          echo "$fres" >> $SUMMARY          echo "$fres" >> $SUMMARY
1200          echo "fresults='$fres'" > $CDIR"/summary.txt"          touch $CDIR"/summary.txt"
1201            echo "fresults='$fres'" >> $CDIR"/summary.txt"
1202          echo "MACH='$MACH'" >> $CDIR"/summary.txt"          echo "MACH='$MACH'" >> $CDIR"/summary.txt"
1203          echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt"          echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt"
1204          echo "DATE='$DATE'" >> $CDIR"/summary.txt"          echo "DATE='$DATE'" >> $CDIR"/summary.txt"
# Line 1164  for dir in $TESTDIRS ; do Line 1206  for dir in $TESTDIRS ; do
1206    
1207          for ex in $extra_runs ; do          for ex in $extra_runs ; do
1208              unset run              unset run
1209              results='-- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --'              results=$EMPTY_RESULTS
1210              #  Create an output dir for each OPTFILE/tdir.ex combination              #  Create an output dir for each OPTFILE/tdir.ex combination
1211              rel_CDIR=$DRESULTS"/"$dir"."$ex              rel_CDIR=$DRESULTS"/"$dir"."$ex
1212              mkdir $rel_CDIR              mkdir $rel_CDIR
# Line 1179  for dir in $TESTDIRS ; do Line 1221  for dir in $TESTDIRS ; do
1221              fres="$fres.$ex"              fres="$fres.$ex"
1222              echo              echo
1223              echo "$fres" >> $SUMMARY              echo "$fres" >> $SUMMARY
1224              echo "fresults='$fres'" > $CDIR"/summary.txt"              touch $CDIR"/summary.txt"
1225                echo "fresults='$fres'" >> $CDIR"/summary.txt"
1226              echo "MACH='$MACH'" >> $CDIR"/summary.txt"              echo "MACH='$MACH'" >> $CDIR"/summary.txt"
1227              echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt"              echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt"
1228              echo "DATE='$DATE'" >> $CDIR"/summary.txt"              echo "DATE='$DATE'" >> $CDIR"/summary.txt"
# Line 1193  for dir in $TESTDIRS ; do Line 1236  for dir in $TESTDIRS ; do
1236          fres=$fres"$results   $dir"          fres=$fres"$results   $dir"
1237          echo          echo
1238          echo "$fres" >> $SUMMARY          echo "$fres" >> $SUMMARY
1239          echo "fresults='$fres'" > $CDIR"/summary.txt"          touch $CDIR"/summary.txt"
1240            echo "fresults='$fres'" >> $CDIR"/summary.txt"
1241          echo "MACH='$MACH'" >> $CDIR"/summary.txt"          echo "MACH='$MACH'" >> $CDIR"/summary.txt"
1242          echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt"          echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt"
1243          echo "DATE='$DATE'" >> $CDIR"/summary.txt"          echo "DATE='$DATE'" >> $CDIR"/summary.txt"
# Line 1247  fi Line 1291  fi
1291  rm -f tr_cmpnum.c tr_cmpnum  rm -f tr_cmpnum.c tr_cmpnum
1292    
1293  if test "x$CLEANUP" != xt ; then  if test "x$CLEANUP" != xt ; then
1294      cat $SUMMARY | sed 's/ -- -- -- --//g'      cat $SUMMARY | sed 's/  \.  \.  \.  \.  \.  \.  \.  \.  \.  \.  \.  \.//'
1295    #   cat $SUMMARY | sed 's/  \.//g'
1296      if test -e tr_out.txt ; then      if test -e tr_out.txt ; then
1297          mv tr_out.txt tr_out.txt.old          mv tr_out.txt tr_out.txt.old
1298      fi      fi
1299      cat $SUMMARY | sed 's/ -- -- -- --//g' > tr_out.txt      cat $SUMMARY | sed 's/  \.//g' > tr_out.txt
1300  fi  fi
1301    
1302  if test "x$DELDIR" = xt ; then  if test "x$DELDIR" = xt ; then

Legend:
Removed from v.1.93  
changed lines
  Added in v.1.94

  ViewVC Help
Powered by ViewVC 1.1.22