/[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.50 by edhill, Wed Nov 24 17:06:07 2004 UTC revision 1.55 by edhill, Thu Dec 9 15:52:33 2004 UTC
# Line 519  createcodelet() Line 519  createcodelet()
519  #include <math.h>  #include <math.h>
520  int main( int argc, char** argv )  {  int main( int argc, char** argv )  {
521    int linnum,best,lncnt;    int linnum,best,lncnt;
522    double a,b,diff;    double a,b,abave,relerr;
523    best = -16;    best = -16;
524    lncnt = 0;    lncnt = 0;
525    while( 1 & (lncnt+=1) < 999 )  {    while( 1 & (lncnt+=1) < 999 )  {
526      scanf("%d", &linnum);      scanf("%d", &linnum);
527      if (linnum == -1)  break;      if (linnum == -1)  break;
528      scanf("%lf", &a);  scanf("%lf", &b);      scanf("%lf", &a);  scanf("%lf", &b);
529      diff = 0.5*(fabs(a)+fabs(b));      abave = 0.5*(fabs(a)+fabs(b));
530      if (diff > 1.e-12) {      if (abave > 0.0) {
531        diff=fabs(a-b)/diff;        relerr=fabs(a-b)/abave;
532        if (diff > 0.0) {        if (relerr > 0.0) {
533          linnum = (int)log10(diff);          linnum = (int)log10(relerr);
534          best = (best > linnum) ? best : linnum;          best = (best > linnum) ? best : linnum;
535        }        }
536        else {      }
537          if (best == -16 && diff != 0)  best = -22;      else {
538        }        if (best == -16)  best = -22;
539      }      }
540    }    }
541    if (lncnt == 999) best=-29;    if (lncnt == 999) best=-29;
# Line 881  for dir in $TESTDIRS ; do Line 881  for dir in $TESTDIRS ; do
881      fout=      fout=
882      if test "x$ADM" = x ; then      if test "x$ADM" = x ; then
883          fout=$dir"/results/output.txt"          fout=$dir"/results/output.txt"
         check_for_add_mon_output  $fout  
884      else      else
885          fout=$dir"/results_ad/output.txt_adm"          fout=$dir"/results_ad/output.txt_adm"
886      fi      fi
# Line 889  for dir in $TESTDIRS ; do Line 888  for dir in $TESTDIRS ; do
888          echo "can't read \"$fout\" -- skipping $dir"          echo "can't read \"$fout\" -- skipping $dir"
889          continue          continue
890      fi      fi
891        if test "x$ADM" = x ; then
892            check_for_add_mon_output  $fout
893        fi
894    
895      # Check for additional types of monitor output      # Check for additional types of monitor output
896    
# Line 967  for dir in $TESTDIRS ; do Line 969  for dir in $TESTDIRS ; do
969          echo "tdir='$dir'" >> $CDIR"/summary.txt"          echo "tdir='$dir'" >> $CDIR"/summary.txt"
970    
971          OLD_COMMAND=$COMMAND          OLD_COMMAND=$COMMAND
972          COMMAND="./mitgcmuv > output.txt"          # COMMAND="./mitgcmuv > output.txt"
973          for ex in $extra_runs ; do          for ex in $extra_runs ; do
974              test -e "$dir/tr_run.$ex" && rm -rf "$dir/tr_run.$ex"              test -e "$dir/tr_run.$ex" && rm -rf "$dir/tr_run.$ex"
975              mkdir "$dir/tr_run.$ex"              mkdir "$dir/tr_run.$ex"
Legend:
Removed from v.1.50  
changed lines
  Added in v.1.55
  ViewVC Help
Powered by ViewVC 1.1.22