/[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.38 by edhill, Fri Mar 12 19:55:47 2004 UTC revision 1.188 by jmc, Fri Jul 19 18:49:33 2013 UTC
# Line 11  usage() Line 11  usage()
11      echo      echo
12      echo "where possible OPTIONS are:"      echo "where possible OPTIONS are:"
13      echo "  (-help|-h)               print usage"      echo "  (-help|-h)               print usage"
14      echo "  (-mpi)                   use MPI input files"      echo "  (-mth)                   run multi-threaded (using eedata.mth)"
15      echo "  (-ieee|-noieee)          if possible, use IEEE compiler flags"      echo "  (-mpi)                   use MPI to compile and run on 2 processors"
16      echo "                             (DEF=\"noieee\")"      echo "  (-MPI)  NUMBER           use MPI to compile and run on max NUMBER procs"
17      echo "  (-optfile=|-of=)STRING   list of optfiles to use"      echo "  (-mfile|-mf) STRING      MPI: file with list of possible machines to run on"
18      echo "  (-a|-addr) STRING        list of email recipients"      echo "  (-fast)                  use optfile default for compiler flags (no '-ieee')"
19      echo "                             (DEF=\"edhill@mitgcm.org\")"      echo "                            DEF=off => use IEEE numerics option (if available)"
20      echo "  (-t|-tdir) STRING        list of test dirs to use"      echo "  (-devel)                 use optfile developement flags (if available)"
21      echo "                             (DEF=\"\" which builds all)"      echo "  (-gsl)                   compile with \"-gsl\" flag"
22      echo "  (-b|-bash) STRING        preferred location of a \"bash\" or"      echo "  (-use_r4|-ur4)           if allowed, use real*4 type for '_RS' variable"
23        echo "  (-optfile|-of) STRING    list of optfiles to use"
24        echo "  (-addr|-a) STRING        list of email recipients"
25        echo "                             (DEF=\"\" no email is sent)"
26        echo "  (-mpackdir|-mpd) DIR     location of the mpack utility"
27        echo "                             (DEF=\"../tools/mpack-1.6\")"
28        echo "  (-tdir|-t) STRING        list of group and/or exp. dirs to test"
29        echo "                             (recognized groups: basic, tutorials)"
30        echo "                             (DEF=\"\" which test all)"
31        echo "                             (if list= 'start_from THIS_EXP' then"
32        echo "                              test THIS_EXP + all the following)"
33        echo "  (-skipdir|-skd) STRING   list of exp. dirs to skip"
34        echo "                             (DEF=\"\" which test all)"
35        echo "  (-bash|-b) STRING        preferred location of a \"bash\" or"
36      echo "                             Bourne-compatible \"sh\" shell"      echo "                             Bourne-compatible \"sh\" shell"
37      echo "                             (DEF=\"\" for \"bash\")"      echo "                             (DEF=\"\" for \"bash\")"
38        echo "  (-tlm)                   perform a Tangent Linear run"
39      echo "  (-adm|-ad)               perform an adjoint run"      echo "  (-adm|-ad)               perform an adjoint run"
40      echo "  (-command) STRING        command to run"      echo "  (-oad)                   perform an OpenAD adjoint run"
41      echo "                             (DEF=\"make output.txt\")"      echo "  (-command|-c) STRING     command to run"
42      echo "  (-m|-make) STRING        command to use for \"make\""      echo "  (-makedepend|-md) STRING command to use for \"makedepend\""
43        echo "  (-make|-m) STRING        command to use for \"make\""
44      echo "                             (DEF=\"make\")"      echo "                             (DEF=\"make\")"
45      echo "  (-clean)                 *ONLY* run \"make CLEAN\""      echo "  (-odir) STRING           used to build output directory name"
46        echo "                             (DEF=\"hostname\")"
47    #   echo "  (-ptracers|-ptr) STRING  specify which ptracers to test"
48    #   echo "                             (DEF=\"1 2 3 4 5\")"
49        echo "  (-match) NUMBER          Matching Criteria (number of digits)"
50        echo "                             (DEF=\"$MATCH_CRIT\")"
51        echo "  (-j) JOBS                use \"make -j JOBS\" for parallel builds"
52        echo "  (-ef) STRING             used as genmake2 \"-extra_flag\" argument"
53        echo "  (-clean)                 *ONLY* run \"make CLEAN\" & clean run-dir"
54        echo "  (-norun|-nr)             skip the \"runmodel\" stage (stop after make)"
55        echo "  (-obj)                   only produces objects (=norun & no executable)"
56        echo "  (-runonly|-ro)           *ONLY* run stage (=\"-quick\" without make)"
57      echo "  (-quick|-q)              same as \"-nogenmake -noclean -nodepend\""      echo "  (-quick|-q)              same as \"-nogenmake -noclean -nodepend\""
58      echo "  (-nogenmake|-ng)         skip the genmake stage"      echo "  (-nogenmake|-ng)         skip the genmake stage"
59      echo "  (-noclean|-nc)           skip the \"make clean\" stage"      echo "  (-noclean|-nc)           skip the \"make clean\" stage"
60      echo "  (-nodepend|-nd)          skip the \"make depend\" stage"      echo "  (-nodepend|-nd)          skip the \"make depend\" stage"
61        echo "  (-postclean|-pc)         after each exp. test, clean build-dir & run-dir"
62        echo "  (-deloutp|-do)           delete output files after successful run"
63        echo "  (-deldir|-dd)            on success, delete the output directory"
64        echo "  (-ts)                    provide timing information per timestep"
65        echo "  (-papis)                 provide MFlop/s per timestep using PAPI"
66        echo "  (-pcls)                  provide MFlop/s per timestep using PCL"
67      echo      echo
68      echo "and where STRING follows a whitespace-delimited format"      echo "and where STRING can be a whitespace-delimited list"
69      echo "such as:"      echo "such as:"
70        echo
71      echo "  -t 'exp0 exp2 exp3' "      echo "  -t 'exp0 exp2 exp3' "
72      echo "  -addr='abc@123.com testing@home.org'"      echo "  -addr='abc@123.com testing@home.org'"
73      echo      echo
74        echo "provided that the expression is properly quoted within the current"
75        echo "shell (note the use of single quotes to protect white space)."
76        echo
77      exit 1      exit 1
78  }  }
79    
# Line 45  usage() Line 81  usage()
81  build_mpack()  build_mpack()
82  {  {
83      printf "building the mpack utility...  "      printf "building the mpack utility...  "
84      if test ! -x "$MPACKDIR/mpack" ; then      MPACK="$MPACKDIR/mpack"
85        if test ! -x $MPACK ; then
86          if test ! -d $MPACKDIR ; then          if test ! -d $MPACKDIR ; then
87              echo              echo
88              echo "Error: can't find \"$MPACKDIR\""              echo "Error: can't find \"$MPACKDIR\""
# Line 54  build_mpack() Line 91  build_mpack()
91              echo              echo
92              HAVE_MPACK=f              HAVE_MPACK=f
93          fi          fi
         printf "building mpack...  "  
94          if test "x$CC" = x ; then          if test "x$CC" = x ; then
95              export CC=cc              export CC=cc
96          fi          fi
97            printf "building mpack (using CC=$CC)...  "
98          ( cd $MPACKDIR && ./configure && $MAKE ) > tr_build_mpack.out 2>&1          ( cd $MPACKDIR && ./configure && $MAKE ) > tr_build_mpack.out 2>&1
99          RETVAL=$?          RETVAL=$?
100          if test "x$RETVAL" != x0 ; then          if test "x$RETVAL" != x0 ; then
# Line 68  build_mpack() Line 105  build_mpack()
105          else          else
106              rm -f tr_build_mpack.out              rm -f tr_build_mpack.out
107              HAVE_MPACK=t              HAVE_MPACK=t
108                echo "done"
109          fi          fi
110      else      else
111          HAVE_MPACK=t          HAVE_MPACK=t
112            echo "already exist"
113      fi      fi
     echo "OK"  
114  }  }
115    
116  testoutput_for_prop()  testoutput_var()
117  {  {
118      # testoutput_for_prop dir s1 label subdir extension      # testoutput_var dir s1 label subdir reference_output
119      #      #
120      #  compares files in $dir/$subdir/output.txt and $dir/results/output.txt      #  compares 1 variable output selected from file $dir/$subdir/$OUTPUTFILE
121      #  using search strings s1 and text label      #     with same output from reference file $dir/$reference_output
122        #  using search strings s1 and text label
123    
124      if [ $debug -gt 0 ]; then      if [ $debug -gt 0 ]; then
125          echo testoutput_for_prop: grep "$2" $1/$4/output.txt 1>&2          echo testoutput_var: grep "$2" $1/$4/$OUTPUTFILE 1>&2
126      fi      fi
127      if [ -r $1/$4/output.txt ]; then      if [ -r $1/$4/$OUTPUTFILE ]; then
128          grep "$2" $1/$4/output.txt | sed 's/.*=//' | cat -n > tmp1.txt          grep "$2" $1/$4/$OUTPUTFILE | sed 's/.*=//' | cat -n > tmp1.txt
129          lncnt=`wc -l tmp1.txt | awk '{print $1}' `          lncntA=`wc -l tmp1.txt | awk '{print $1}' `
130          if [ $lncnt -lt 3 ]; then          if [ $lncntA -lt 2 ]; then
131              if [ $verbose -gt 0 ]; then              if [ $verbose -gt 0 ]; then
132                  echo Not enough lines of output when searching for "$2" 1>&2                  echo Not enough lines of output when searching for "$2" 1>&2
133              fi              fi
134              return 99              rm -f tmp1.txt ; return 99
135          fi          fi
136      else      else
137          echo testoutput_for_prop: output.txt from model run was not readable 1>&2          echo testoutput_var: $OUTPUTFILE from model run was not readable 1>&2
138          return 99          return 99
139      fi      fi
140      if [ $debug -gt 0 ]; then      if [ $debug -gt 0 ]; then
141          echo testoutput_for_prop: grep "$2" $1/results/output.txt$5 1>&2          echo testoutput_var: grep "$2" $1/$5 1>&2
142      fi      fi
143      grep "$2" $1/results/output.txt$5 | sed 's/.*=//' | cat -n > tmp2.txt      grep "$2" $1/$5 | sed 's/.*=//' | cat -n > tmp2.txt
144      lncnt=`wc -l tmp2.txt | awk '{print $1}' `      lncntB=`wc -l tmp2.txt | awk '{print $1}' `
145      if [ $lncnt -lt 3 ]; then      if [ $lncntB -lt 2 ]; then
146          if [ $verbose -gt 0 ]; then          if [ $verbose -gt 0 ]; then
147              echo Not enough lines of output when searching for "$2" 1>&2              echo Not enough lines of output when searching for "$2" 1>&2
148          fi          fi
149          return 99          rm -f tmp1.txt tmp2.txt ; return 99
150        fi
151        if [ $lncntA -ne $lncntB ]; then
152            if [ $verbose -gt 0 ]; then
153                echo Not same Nb of lines when searching for "$2" ":" $lncntA $lncntB 1>&2
154            fi
155            rm -f tmp1.txt tmp2.txt ; return 99
156        fi
157        has_nan=`cat tmp1.txt | grep -i nan | wc -l`
158        if [ $has_nan -gt 0  ] ; then
159            echo testoutput_var: $OUTPUTFILE contains $has_nan NaN values  1>&2
160            rm -f tmp1.txt tmp2.txt ; return 99
161        fi
162        has_inf=`cat tmp1.txt | grep -i inf | wc -l`
163        if [ $has_inf -gt 0  ] ; then
164            echo testoutput_var: $OUTPUTFILE contains $has_inf Inf values  1>&2
165            rm -f tmp1.txt tmp2.txt ; return 99
166      fi      fi
167      if [ $debug -gt 0 ]; then      if [ $debug -gt 0 ]; then
168          echo testoutput_for_prop: join tmp1.txt tmp2.txt 1>&2          echo testoutput_var: join tmp1.txt tmp2.txt 1>&2
169      fi      fi
170      join tmp1.txt tmp2.txt | awk '{print $1 " " $2 " " $3}' > tmp3.txt      join tmp1.txt tmp2.txt | awk '{print $1 " " $2 " " $3}' > tmp3.txt
171      if [ $debug -gt 0 ]; then      if [ $debug -gt 0 ]; then
172          echo testoutput_for_prop: compare_lines 1>&2          echo testoutput_var: compare_lines 1>&2
173      fi      fi
174      if [ $verbose -gt 1 ]; then      if [ $verbose -gt 1 ]; then
175          cat tmp3.txt 1>&2          cat tmp3.txt 1>&2
# Line 122  testoutput_for_prop() Line 177  testoutput_for_prop()
177      echo "-1" >> tmp3.txt      echo "-1" >> tmp3.txt
178      # On the SGI O3K (*not* the O2K), "cat -n" inserts a ":" after the line number      # On the SGI O3K (*not* the O2K), "cat -n" inserts a ":" after the line number
179      cat tmp3.txt | sed -e 's|:||g' > tmp4.txt      cat tmp3.txt | sed -e 's|:||g' > tmp4.txt
180      digits_of_similarity=`./tmp_cmpnum < tmp4.txt`      digits_of_similarity=`./tr_cmpnum < tmp4.txt`
181      if [ $digits_of_similarity -eq 99 ]; then      if [ $digits_of_similarity -eq 99 ]; then
182          if [ $verbose -gt 0 ]; then          if [ $verbose -gt 0 ]; then
183              echo testoutput_for_prop: No comparison was available for \"$2\" 1>&2              echo testoutput_var: No comparison was available for \"$3\" 1>&2
184          fi          fi
185          digits_of_similarity=99          digits_of_similarity=99
186      else      else
187          if [ $verbose -gt 0 ]; then          if [ $verbose -gt 0 ]; then
188              echo There were $digits_of_similarity decimal places of similarity for \"$2\" 1>&2              echo There were $digits_of_similarity decimal places of similarity for \"$3\" 1>&2
189          fi          fi
190      fi      fi
191      rm -f tmp1.txt tmp2.txt tmp3.txt tmp4.txt      rm -f tmp1.txt tmp2.txt tmp3.txt tmp4.txt
192        
193      return $digits_of_similarity      return $digits_of_similarity
194  }  }
195    
196  dashnum()  testoutput_run()
197  {  {
198      # dashnum n1 n2 n3 ...      # testoutput_run directory subdir reference_output
199      #      #
200      #  print numbers using %3i format or "--" if number = 99      #  test output from 1 run in "directory"
201    # --> same processing for adjoint & forward test
202      for num in $@ ; do          # default list of output variables to be checked:
203          if [ $num = 99 ]; then          #  1rst : main variable used to decide if it pass or FAIL
204              printf ' --'          #  others : number of matching digits to be printed in summary.txt
205          else          listChk=$DEF_CHECK_LIST
206              printf '%3i' $num          #  load experiment-specific list from file "tr_checklist" (if it exist)
207            if test -r $1/$2/tr_checklist ; then listChk=`cat $1/$2/tr_checklist` ; fi
208            sVar=`echo $listChk | awk '{print $1}'`
209            # remove 1rst var and expand the list: + => min max mean s.d
210            listVar=`echo $listChk | sed 's/ [a-zA-Z0-9]*+/&mn &mx &av &sd/g' \
211                                   | sed 's/+//g' | sed "s/^$sVar//"`
212            if [ $debug -gt 0 ]; then echo "testoutput_run: listVar(I)='$listVar'" 1>&2 ; fi
213            # set ref.outp file-name and, if compressed, uncompress it into subdir
214            refoutp=results/$3
215            if test ! -r $1/results/$3 ; then
216              if test -r $1/results/${3}.gz ; then
217                tmpNam1=$3
218                tmpNam2=`echo $tmpNam1 | sed 's/\..*\././'`
219                if test "x$tmpNam1" = "x$tmpNam2" ; then
220                    tmpNam1=`echo $tmpNam2 | sed 's/\./.std./'`
221                fi
222                refoutp=$2/$tmpNam1
223                if [ $debug -gt 0 ]; then
224                    echo "testoutput_run: gunzip 'results/$3.gz' into '$refoutp'" 1>&2
225                fi
226                gzip -cd $1/results/${3}.gz > $1/$refoutp
227              else
228                echo "missing reference output file '$3' (.gz) in '$1/results/'"
229              fi
230          fi          fi
231      done          # check for ptracer output in reference_output file :
232  }          if test $KIND = 0 ; then
233              ptr_mon="trcstat_ptracerXX_min trcstat_ptracerXX_max"
234  testoutput_ad()            ptr_mon="$ptr_mon trcstat_ptracerXX_mean trcstat_ptracerXX_sd"
235  {            for ii in $PTRACERS_NUM ; do
236      grep $3 $1/results_ad/output.txt_adm | awk '{print NR " " $5}' > t05.txt              ptrfound=0
237      grep $3 $1/$2/output.txt_adm | awk '{print NR " " $5}' > t15.txt              for jj in $ptr_mon ; do
238      grep $3 $1/results_ad/output.txt_adm | awk '{print NR " " $6}' > t06.txt                  name=`eval "echo $jj | sed -e 's|XX|0"$ii"|g'"`
239      grep $3 $1/$2/output.txt_adm | awk '{print NR " " $6}' > t16.txt                  tst=`grep $name $1/$refoutp | wc -l | awk '{print $1}'`
240      join t05.txt t15.txt > t5.txt                  if test ! "x$tst" = x0 ; then ptrfound=1 ; fi
241      join t06.txt t16.txt > t6.txt              done
242      echo "-1" >> t5.txt              if test $ptrfound = '1' ; then
243      echo "-1" >> t6.txt                  eval "HAVE_PTR0"$ii"=t"
244      digits_5=`./tmp_cmpnum < t5.txt`              else
245      digits_6=`./tmp_cmpnum < t6.txt`                  eval "HAVE_PTR0"$ii"=f"
246      dashnum $digits_5 $digits_6                # remove this ptr from the list of output variable to check
247      rm -f t[01][56].txt t[56].txt                # echo "-- ptr test=" $tst "number of var=" `echo $listVar | awk '{print NF}'` 1>&2
248  }                  listVar=`echo "$listVar" | sed "s/ pt$ii..//g"`
249                fi
250              # eval 'echo "HAVE_PTR0'$ii' = $HAVE_PTR0'$ii'"' 1>&2
251              done
252            fi
253            tst=`echo $sVar $listVar | awk '{ for(i=2;i<=NF;i++){if($i==$1)t+=1}; print t }'`
254            if test $tst != 1 ; then
255              if test $tst = 0 ; then echo "==> WARNING: selected var >$sVar< not found" 1>&2
256                     else echo "==> WARNING: found selected var >$sVar< $tst times" 1>&2 ; fi
257              echo "==> WARNING: in checked list:" $listVar 1>&2
258            #- put it back once:
259              listVar=" $sVar "`echo "$listVar " | sed "s/ $sVar / /g"`
260            fi
261            if [ $debug -gt 0 ]; then echo "testoutput_run: listVar(M)='$listVar'" 1>&2 ; fi
262            echo "listVar='$listVar'" >> $locDIR"/summary.txt"
263            #---
264            allargs=""
265            for xx in $listVar
266            do
267              case $xx in
268               'PS')  if [ $debug -gt 0 ]
269                      then echo testoutput_run: testoutput_var $1 cg2d_init_res 1>&2 ; fi
270                      testoutput_var $1 "cg2d_init_res" "Press. Solver (cg2d)" $2 $refoutp ; yy=$?
271                      if [ $debug -gt 0 ] ; then echo testoutput_run: cg2dres=$yy 1>&2 ; fi ;;
272              'Cost') testoutput_var $1 "ADM  ref_cost_function" "ADM CostFct" $2 $refoutp ; yy=$? ;;
273             'AdGrd') testoutput_var $1 "ADM  adjoint_gradient"  "ADM Ad Grad" $2 $refoutp ; yy=$? ;;
274             'FDGrd') testoutput_var $1 "ADM  finite-diff_grad"  "ADM FD Grad" $2 $refoutp ; yy=$? ;;
275            'tlmCst') testoutput_var $1 "TLM  ref_cost_function" "TLM CostFct" $2 $refoutp ; yy=$? ;;
276            'tlmGrd') testoutput_var $1 "TLM  tangent-lin_grad"  "TLM TL Grad" $2 $refoutp ; yy=$? ;;
277            'fwdGrd') testoutput_var $1 "TLM  finite-diff_grad"  "TLM FD Grad" $2 $refoutp ; yy=$? ;;
278               'Tmn') testoutput_var $1 "dynstat_theta_min"  "Theta minimum"  $2 $refoutp ; yy=$? ;;
279               'Tmx') testoutput_var $1 "dynstat_theta_max"  "Theta maximum"  $2 $refoutp ; yy=$? ;;
280               'Tav') testoutput_var $1 "dynstat_theta_mean" "Theta mean"     $2 $refoutp ; yy=$? ;;
281               'Tsd') testoutput_var $1 "dynstat_theta_sd"   "Theta Std.Dev"  $2 $refoutp ; yy=$? ;;
282               'Smn') testoutput_var $1 "dynstat_salt_min"  "Salt minimum"    $2 $refoutp ; yy=$? ;;
283               'Smx') testoutput_var $1 "dynstat_salt_max"  "Salt maximum"    $2 $refoutp ; yy=$? ;;
284               'Sav') testoutput_var $1 "dynstat_salt_mean" "Salt mean"       $2 $refoutp ; yy=$? ;;
285               'Ssd') testoutput_var $1 "dynstat_salt_sd"   "Salt Std.Dev"    $2 $refoutp ; yy=$? ;;
286               'Umn') testoutput_var $1 "dynstat_uvel_min"  "U minimum"       $2 $refoutp ; yy=$? ;;
287               'Umx') testoutput_var $1 "dynstat_uvel_max"  "U maximum"       $2 $refoutp ; yy=$? ;;
288               'Uav') testoutput_var $1 "dynstat_uvel_mean" "U mean"          $2 $refoutp ; yy=$? ;;
289               'Usd') testoutput_var $1 "dynstat_uvel_sd"   "U Std.Dev"       $2 $refoutp ; yy=$? ;;
290               'Vmn') testoutput_var $1 "dynstat_vvel_min"  "V minimum"       $2 $refoutp ; yy=$? ;;
291               'Vmx') testoutput_var $1 "dynstat_vvel_max"  "V maximum"       $2 $refoutp ; yy=$? ;;
292               'Vav') testoutput_var $1 "dynstat_vvel_mean" "V mean"          $2 $refoutp ; yy=$? ;;
293               'Vsd') testoutput_var $1 "dynstat_vvel_sd"   "V Std.Dev"       $2 $refoutp ; yy=$? ;;
294            'pt1mn'|'pt2mn'|'pt3mn'|'pt4mn'|'pt5mn') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'`
295               testoutput_var $1 "trcstat_ptracer0"$ii"_min"  "p0"$ii"_min"   $2 $refoutp ; yy=$? ;;
296            'pt1mx'|'pt2mx'|'pt3mx'|'pt4mx'|'pt5mx') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'`
297               testoutput_var $1 "trcstat_ptracer0"$ii"_max"  "p0"$ii"_max"   $2 $refoutp ; yy=$? ;;
298            'pt1av'|'pt2av'|'pt3av'|'pt4av'|'pt5av') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'`
299               testoutput_var $1 "trcstat_ptracer0"$ii"_mean" "p0"$ii"_mean" $2 $refoutp ; yy=$? ;;
300            'pt1sd'|'pt2sd'|'pt3sd'|'pt4sd'|'pt5sd') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'`
301               testoutput_var $1 "trcstat_ptracer0"$ii"_sd"   "p0"$ii"_StDv"  $2 $refoutp ; yy=$? ;;
302             'Qntmn') testoutput_var $1 "forcing_qnet_min" "Qnet minimum"  $2 $refoutp ; yy=$? ;;
303             'Qntmx') testoutput_var $1 "forcing_qnet_max" "Qnet maximum"  $2 $refoutp ; yy=$? ;;
304             'Qntav') testoutput_var $1 "forcing_qnet_mean" "Qnet mean"       $2 $refoutp ; yy=$? ;;
305             'Qntsd') testoutput_var $1 "forcing_qnet_sd"  "Qnet Std.Dev"  $2 $refoutp ; yy=$? ;;
306             'aSImn') testoutput_var $1 "seaice_area_min"   "SIce Area min"   $2 $refoutp ; yy=$? ;;
307             'aSImx') testoutput_var $1 "seaice_area_max"   "SIce Area max"   $2 $refoutp ; yy=$? ;;
308             'aSIav') testoutput_var $1 "seaice_area_mean"  "SIce Area mean"  $2 $refoutp ; yy=$? ;;
309             'aSIsd') testoutput_var $1 "seaice_area_sd"    "SIce Area StDv"  $2 $refoutp ; yy=$? ;;
310             'hSImn') testoutput_var $1 "seaice_heff_min"   "SIce Heff min"   $2 $refoutp ; yy=$? ;;
311             'hSImx') testoutput_var $1 "seaice_heff_max"   "SIce Heff max"   $2 $refoutp ; yy=$? ;;
312             'hSIav') testoutput_var $1 "seaice_heff_mean"  "SIce Heff mean"  $2 $refoutp ; yy=$? ;;
313             'hSIsd') testoutput_var $1 "seaice_heff_sd"    "SIce Heff StDv"  $2 $refoutp ; yy=$? ;;
314             'uSImn') testoutput_var $1 "seaice_uice_min"   "SIce Uice min"   $2 $refoutp ; yy=$? ;;
315             'uSImx') testoutput_var $1 "seaice_uice_max"   "SIce Uice max"   $2 $refoutp ; yy=$? ;;
316             'uSIav') testoutput_var $1 "seaice_uice_mean"  "SIce Uice mean"  $2 $refoutp ; yy=$? ;;
317             'uSIsd') testoutput_var $1 "seaice_uice_sd"    "SIce Uice StDv"  $2 $refoutp ; yy=$? ;;
318             'vSImn') testoutput_var $1 "seaice_vice_min"   "SIce Vice min"   $2 $refoutp ; yy=$? ;;
319             'vSImx') testoutput_var $1 "seaice_vice_max"   "SIce Vice max"   $2 $refoutp ; yy=$? ;;
320             'vSIav') testoutput_var $1 "seaice_vice_mean"  "SIce Vice mean"  $2 $refoutp ; yy=$? ;;
321             'vSIsd') testoutput_var $1 "seaice_vice_sd"    "SIce Vice StDv"  $2 $refoutp ; yy=$? ;;
322            'AthSiG') testoutput_var $1 "thSI_Ice_Area_G" "thSIc Area Global" $2 $refoutp ; yy=$? ;;
323            'AthSiS') testoutput_var $1 "thSI_Ice_Area_S" "thSIc Area South"  $2 $refoutp ; yy=$? ;;
324            'AthSiN') testoutput_var $1 "thSI_Ice_Area_N" "thSIc Area North"  $2 $refoutp ; yy=$? ;;
325            'HthSiG') testoutput_var $1 "thSI_IceH_ave_G" "thSIc H Glob-ave"  $2 $refoutp ; yy=$? ;;
326            'HthSiS') testoutput_var $1 "thSI_IceH_ave_S" "thSIc H South-av"  $2 $refoutp ; yy=$? ;;
327            'HthSiN') testoutput_var $1 "thSI_IceH_ave_N" "thSIc H North-av"  $2 $refoutp ; yy=$? ;;
328            'HthMxS') testoutput_var $1 "thSI_IceH_max_S" "thSIc H South-max" $2 $refoutp ; yy=$? ;;
329            'HthMxN') testoutput_var $1 "thSI_IceH_max_N" "thSIc H North-max" $2 $refoutp ; yy=$? ;;
330            'StrmIc') testoutput_var $1 "STREAMICE_FP_ERR" "StreamIce Solver" $2 $refoutp ; yy=$? ;;
331                  *) yy=99; echo "WARNING: asking for var=$xx : not recognized !" 1>&2 ;;
332              esac
333              if test $xx = $sVar
334              then allargs="$allargs > $yy <"
335              else allargs="$allargs $yy"
336              fi
337            done
338    
339  testoutput()          nbVar=`echo $listVar | awk '{print NF}'`
340  {          if [ $nbVar -lt $LEN_CHECK_LIST ] ; then
341      # testoutput directory subdir extension          #-- fill line (up to standard length) with dot:
342      #            adNul=`expr $LEN_CHECK_LIST - $nbVar | awk '{for(i=1;i<=$1;i++){print "."}}'`
343      #  test output in "directory"            echo $allargs $adNul
344      if test "x$ADM" = x ; then          else
345          if [ $debug -gt 0 ]; then            echo $allargs
346              echo testoutput: testoutput_for_prop $1 cg2d_init_res 1>&2          fi
347          fi  # <-- same processing for adjoint & forward test
         testoutput_for_prop $1 "cg2d_init_res" "cg2d init. residual" $2 $3; cg2dres=$?  
         if [ $debug -gt 0 ]; then  
             echo testoutput: cg2dres=$cg2dres 1>&2  
         fi  
         testoutput_for_prop $1 "dynstat_theta_min" "theta minimum" $2 $3; tmin=$?  
         testoutput_for_prop $1 "dynstat_theta_max" "theta maximum" $2 $3; tmax=$?  
         testoutput_for_prop $1 "dynstat_theta_mean" "theta mean" $2 $3; tmean=$?  
         testoutput_for_prop $1 "dynstat_theta_sd" "theta s.d." $2 $3; tsd=$?  
         testoutput_for_prop $1 "dynstat_salt_min" "salt minimum" $2 $3; smin=$?  
         testoutput_for_prop $1 "dynstat_salt_max" "salt maximum" $2 $3; smax=$?  
         testoutput_for_prop $1 "dynstat_salt_mean" "salt mean" $2 $3; smean=$?  
         testoutput_for_prop $1 "dynstat_salt_sd" "salt s.d." $2 $3; ssd=$?  
         testoutput_for_prop $1 "dynstat_uvel_min" "U minimum" $2 $3; umin=$?  
         testoutput_for_prop $1 "dynstat_uvel_max" "U maximum" $2 $3; umax=$?  
         testoutput_for_prop $1 "dynstat_uvel_mean" "U mean" $2 $3; umean=$?  
         testoutput_for_prop $1 "dynstat_uvel_sd" "U s.d." $2 $3; usd=$?  
         testoutput_for_prop $1 "dynstat_vvel_min" "V minimum" $2 $3; vmin=$?  
         testoutput_for_prop $1 "dynstat_vvel_max" "V maximum" $2 $3; vmax=$?  
         testoutput_for_prop $1 "dynstat_vvel_mean" "V mean" $2 $3; vmean=$?  
         testoutput_for_prop $1 "dynstat_vvel_sd" "V s.d." $2 $3; vsd=$?  
         dashnum $cg2dres $tmin $tmax $tmean $tsd $smin $smax $smean $ssd \  
             $umin $umax $umean $usd $vmin $vmax $vmean $vsd  
     else  
         testoutput_ad $1 $2 "precision_grdchk_result"  
     fi  
348  }  }
349    
350  genmakemodel()  genmakemodel()
# Line 219  genmakemodel() Line 361  genmakemodel()
361          (          (
362              cd $1;              cd $1;
363              command="$GENMAKE2  -ds -m $MAKE"              command="$GENMAKE2  -ds -m $MAKE"
364              if test "x$ADM" = x ; then              if test "x$MKDEPEND" != x ; then
365                  command="$command --mods=../code"                  command="$command -makedepend=$MKDEPEND"
366              else              fi
367                if test $KIND = 1 -o $KIND = 2 ; then
368                  command="$command --mods=../code_ad"                  command="$command --mods=../code_ad"
369                  command="$command -adof=../../../tools/adjoint_options/adjoint_staf"              elif test $KIND = 4 ; then
370                    command="$command -adof ../../../tools/adjoint_options/adjoint_oad -mods '../code_oad ../../OpenAD/code_oad_all'"
371                else
372                    command="$command -mods=../code"
373              fi              fi
374              if test "x$OPTFILE" != xNONE ; then              if test "x$OPTFILE" != xNONE ; then
375                  command="$command --optfile=$OPTFILE"                  command="$command -optfile=$OPTFILE"
376              fi              fi
377              if test "x$IEEE" != x ; then              if test $OptLev = 1 ; then
378                  command="$command -ieee"                  command="$command -ieee"
379              fi              fi
380              printf 'genmake ... ' 1>&2              if test $OptLev = 0 ; then
381              $command > make.log 2>&1                  command="$command -devel"
382                fi
383                if test "x$GSL" = xt ; then
384                    command="$command -gsl"
385                fi
386                if test "x$MPI" != x0 ; then
387                    command="$command -mpi"
388                fi
389                if test "x$MULTI_THREAD" = xt ; then
390                #- run multi-threaded using OpenMP:
391                    command="$command -omp"
392                fi
393                if test "x$USE_R4" = xt ; then
394                    command="$command -use_r4"
395                fi
396                if test "x$EXTRFLG" != x ; then
397                    command="$command -extra_flag $EXTRFLG"
398                fi
399                if test "x$TS" = xt ; then
400                    command="$command -ts"
401                fi
402                if test "x$PAPIS" = xt ; then
403                    command="$command -papis"
404                else
405                if test "x$PCLS" = xt ; then
406                    command="$command -pcls"
407                fi
408                fi
409                printf 'genmake ... '
410                eval $command > genmake.tr_log 2>&1
411              RETVAL=$?              RETVAL=$?
412              cp Makefile $CDIR              #  Reduce the size of the testing emails!
413                head -100 Makefile > $CDIR/Makefile_head
414              if test "x$RETVAL" != x0 ; then              if test "x$RETVAL" != x0 ; then
415                  tail make.log                  tail genmake.tr_log
416                  echo "genmakemodel: genmake failed" 1>&2                  echo "genmakemodel: genmake failed"
417                  cp genmake_* make.log $CDIR                  cp genmake.log genmake_* genmake.tr_log $CDIR
418                  return 1                  return 1
419              else              else
420                  echo "successful" 1>&2                  echo "successful"
421              fi              fi
422          )          )
423      fi      fi
# Line 250  genmakemodel() Line 426  genmakemodel()
426  makeclean()  makeclean()
427  {  {
428      # makeclean directory      # makeclean directory
429        if test "x$NODEPEND" = xf ; then rm -f $1/make.tr_log ; fi
430      if test "x$NOCLEAN" = xt ; then      if test "x$NOCLEAN" = xt ; then
431          echo "make CLEAN skipped!"          echo "make Clean skipped!"
432      else      else
433          (          (
434              cd $1;              cd $1;
435              if test -e output.txt ; then              #if test -e $OUTPUTFILE ; then rm -f $OUTPUTFILE ; fi
                 rm -f output.txt  
             fi  
             printf 'make CLEAN ... ' 2>&1  
436              if test -r Makefile ; then              if test -r Makefile ; then
437                  $MAKE CLEAN >> make.log 2>&1                  printf 'clean build-dir: make Clean ... '
438                    $MAKE Clean >> make.tr_log 2>&1
439                  RETVAL=$?                  RETVAL=$?
440                  if test "x$RETVAL" != x0 ; then                  if test "x$RETVAL" != x0 ; then
441                      tail make.log                      tail make.tr_log
442                      echo "makeclean: \"make CLEAN\" failed" 1>&2                      echo "makeclean: \"make Clean\" failed"
443                      cp make.log $CDIR"/make.log"                      cp make.tr_log genmake.log genmake.tr_log $CDIR
444                      return 1                      return 1
445                  fi                  fi
446                    echo successful
447                else
448                    echo ''
449              fi              fi
450              echo successful 1>&2              exit 0
451            )
452        fi
453    }
454    
455    run_clean()
456    {
457        # run_clean directory
458        if test "x$NOCLEAN" = xt ; then
459            echo "run_clean skipped!"
460        else
461            (
462                cd $1;
463                printf 'clean run-dir ... '
464                # part of what is done after "make clean" when doing "make CLEAN"
465                find . -name "*.meta" -exec rm {} \;
466                find . -name "*.data" -exec rm {} \;
467                find . -name "fort.*" -exec rm {} \;
468                find . -type l -exec rm {} \;
469                #- should remove executable only if sym-link (alredy done above)
470                rm -f $RUNLOG *.txt STD* w2_tile_topology.????.log *diagnostics.log datetime
471                rm -rf mnc_test_*
472                rm -f *_MIT_CE_000.opt0000 costfunction*0000
473                rm -f oad_cp.000.[0-9][0-9][0-9][0-9][0-9]
474                echo successful
475              exit 0              exit 0
476          )          )
477      fi      fi
# Line 283  makedependmodel() Line 485  makedependmodel()
485      else      else
486          (          (
487              cd $1;              cd $1;
488              printf 'make depend ... ' 1>&2              printf 'make depend ... '
489              $MAKE depend >> make.log 2>&1              $MAKE depend >> make.tr_log 2>&1
490              RETVAL=$?              RETVAL=$?
491              if test "x$RETVAL" != x0 ; then              if test "x$RETVAL" != x0 ; then
492                  tail make.log                  tail make.tr_log
493                  echo "makedependmodel: make depend failed" 1>&2                  echo "makedependmodel: make depend failed"
494                  cp make.log $CDIR"/make.log"                  cp make.tr_log genmake.log genmake.tr_log $CDIR
495                  return 1                  return 1
496              else              else
497                  echo successful 1>&2                  echo successful
498              fi              fi
499          )          )
500      fi      fi
# Line 302  makemodel() Line 504  makemodel()
504  {  {
505      # makemodel directory      # makemodel directory
506      (      (
507        mk_fail=0
508        if test "x$NOMAKE" = xt ; then
509            cd $1;
510            if test -x $EXECUTABLE ; then
511                echo "make skipped!"
512            else
513                echo "no executable!"
514                mk_fail=3
515            fi
516        else
517          cd $1;          cd $1;
518          if test -r Makefile ; then          if test -r Makefile ; then
519              printf 'make ... ' 1>&2              printf 'make ... '
520              if test "x$ADM" = x ; then              if test "x$REPLMAKE" = x ; then
521                  $MAKE >> make.log 2>&1                  $MAKE $TARG >> make.tr_log 2>&1
522              else              else
523                  $MAKE adall >> make.log 2>&1                  $REPLMAKE $TARG >> make.tr_log 2>&1
524              fi              fi
525              RETVAL=$?              RETVAL=$?
526              if test "x$RETVAL" != x0 ; then              if test "x$RETVAL" != x0 ; then
527                  tail make.log                  tail make.tr_log
528                  echo failed 1>&2                  echo failed
529                  cp make.log $CDIR"/make.log"                  cp genmake.log genmake.tr_log $CDIR
530                  return 1                  tail -$NBLINES_MKLOG make.tr_log > $CDIR"/make.tr_log_tail"
531                    rm -f $EXECUTABLE
532                    mk_fail=1
533              else              else
534                  echo successful 1>&2                  echo successful
535              fi              fi
536            else
537                echo "no Makefile !"
538                mk_fail=2
539            fi
540        fi
541        if test $KIND = 1 -a -f taf_ftl.log ; then
542            head -1 taf_ftl.log >> $CDIR"/summary.txt"
543        fi
544        if test $KIND = 2 -a -f taf_ad.log ; then
545            head -1 taf_ad.log >> $CDIR"/summary.txt"
546            nerr=`grep -c 'TAF *.* ERROR ' taf_ad.log`
547            nwar=`grep -c 'TAF RECOMPUTATION *.* WARNING ' taf_ad.log`
548            if test -f taf_output ; then
549                n2er=`grep -c 'TAF *.* ERROR ' taf_output`
550                nerr=`expr $nerr + $n2er`
551          fi          fi
552            echo " TAF reports $nerr Errors and $nwar Recomputation Warnings" \
553                                    >> $CDIR"/summary.txt"
554        fi
555        if test $mk_fail != 0 ; then return $mk_fail ; fi
556      )      )
557  }  }
558    
559    mk_mpi_size()
560    {
561        # mk_mpi_size input_file output_file proc_Nb threads_Nb_X threads_Nb_Y
562        #
563        #  make new SIZE.h (=output_file) from SIZE.h_mpi (=input_file)
564        #     for an MPI build with no more than proc_Nb processors ;
565        #  ensure that enough tiles per proc (nSx,nSy) remain for the given
566        #     number of threads (nTx,nTy) ;
567        #  return the effective number of processors.
568    
569        inp=$1
570        out=$2
571        np=$3
572        tx=$4
573        ty=$5
574        tmp=TTT.$$
575    
576        # dirX : select with direction to favor in MPI process repartition
577        #   dirX=1 : prefer to put more proc in X direction
578        #   dirX=0 : prefer to put more proc in Y direction
579        dirX=0
580    
581        px=`grep "^     & *nPx *=" $inp | sed "s/^     & *nPx *= *//" | sed 's/, *$//'`
582        py=`grep "^     & *nPy *=" $inp | sed "s/^     & *nPy *= *//" | sed 's/, *$//'`
583        sx=`grep "^     & *nSx *=" $inp | sed "s/^     & *nSx *= *//" | sed 's/, *$//'`
584        sy=`grep "^     & *nSy *=" $inp | sed "s/^     & *nSy *= *//" | sed 's/, *$//'`
585    
586        #- for each direction, assume # of threads is a multiple of total number of tiles
587        nx=$px
588        if [ `expr $sx % $tx` -ne 0 -a `expr $sx \* $px % $tx` -eq 0 ] ; then
589            nx=`expr $sx \* $px / $tx`
590            if [ $verbose -gt 1 ]; then
591                echo " change px from $px to $nx to accommodate $tx threads"
592            fi
593        fi
594        ny=$py
595        if [ `expr $sy % $ty` -ne 0 -a `expr $sy \* $py % $ty` -eq 0 ] ; then
596            ny=`expr $sy \* $py / $ty`
597            if [ $verbose -gt 1 ]; then
598                echo " change py from $py to $ny to accommodate $ty threads"
599            fi
600        fi
601        #- find the largest divisor of input_file proc Nb, but not larger than $np
602        pp=0
603        i=1
604        while [ $i -le $nx ] ; do
605          if [ `expr $nx % $i` -eq 0 ] ; then
606            j=1
607            while [ $j -le $ny ] ; do
608              if [ `expr $ny % $j` -eq 0 ] ; then
609                ij=`expr $i \* $j`
610                if [ $ij -gt $pp ] ; then
611                    flag=1
612                elif [ $ij -eq $pp ] ; then
613                    flag=$dirX
614                else
615                    flag=0
616                fi
617                if test $flag = 1 ; then
618                  if [ $ij -le $np ] ; then
619                    ix=$i ; jy=$j ; pp=$ij
620                    #echo "  ix,jy= $ix,$jy"
621                  fi
622                fi
623              fi
624              j=`expr $j + 1`
625            done
626          fi
627          i=`expr $i + 1`
628        done
629    
630        #- create new SIZE.h type file:
631        sx=`expr $sx \* $px / $ix`
632        sy=`expr $sy \* $py / $jy`
633        if [ $verbose -gt 1 ]; then
634            echo " px,py,np= $px,$py,$np : New MPI size: px,py= $ix,$jy : sx,sy= $sx,$sy"
635        fi
636        sed "/^     \& *nPx *=/s/[0-9]*,/$ix,/" $inp > $tmp
637        sed "/^     \& *nPy *=/s/[0-9]*,/$jy,/" $tmp > $out
638        sed "/^     \& *nSx *=/s/[0-9]*,/$sx,/" $out > $tmp
639        sed "/^     \& *nSy *=/s/[0-9]*,/$sy,/" $tmp > $out
640        rm -f $tmp
641        return $pp
642    }
643    
644  symlink_mpifiles()  symlink_mpifiles()
645  {  {
646      # Put special links so that MPI specific files are used      # Put special links so that MPI specific files are used
# Line 331  symlink_mpifiles() Line 649  symlink_mpifiles()
649    
650      dir=$1      dir=$1
651      code_dir=$2      code_dir=$2
652      BUILD_DIR=$dir/$3      build_dir=$dir/$3
     CODE_DIR=$dir/$code_dir  
653    
654      # These are files that should replace their counter-part when using -mpi      # These are files that should replace their counter-part when using -mpi
655      MPI_FILES=`(cd $CODE_DIR; find . -name "*_mpi")`      MPI_FILES=`(cd $dir/$code_dir; find . -name "*_mpi" -print)`
656    
657      #  Is this an MPI run?      for ii in $MPI_FILES ; do
658      if test "x$MPI" = xt ; then          i=`echo $ii | sed 's:^\./::'`
659          # YES: We symbolically link these files to the build          name=`echo $i | sed 's:_mpi::'`
660          # dir so long as there is no real file in place          file="../$code_dir/$i"
661          for ii in $MPI_FILES ; do          if test $name = 'SIZE.h' ; then file="SIZE.h.mpi" ; fi
662              i=`echo $ii | sed 's:^\./::'`  
663              name=`echo $i | sed 's:_mpi::' `          #  Is this an MPI run?
664              cmp $CODE_DIR/$i $BUILD_DIR/$name > /dev/null 2>&1          if test "x$MPI" = x0 ; then
665              RETVAL=$?              # NO: We undo any _mpi symbolically linked files
666              if test "x$RETVAL" != x0 ; then              if test -L $build_dir/$name ; then
667                  if ! test -f $BUILD_DIR/$i ; then                  ( cd $build_dir ; cmp $name $file > /dev/null 2>&1 )
668                  #echo Linking $name to $i                  RETVAL=$?
669                      (cd $BUILD_DIR; ln -sf ../$code_dir/$i $name)                  if test "x$RETVAL" = x0 ; then
670                        if [ $verbose -gt 1 ]; then
671                            echo "  Un-linking $name from ../$code_dir" ; fi
672                        rm -f $build_dir/$name
673                  fi                  fi
674              fi              fi
675          done          else
676      else              # YES: We symbolically link these files to the build
677      # NO: We undo any _mpi symbolically linked files              # dir so long as there is no real file in place
678          for ii in $MPI_FILES ; do              ( cd $build_dir ; cmp $name $file > /dev/null 2>&1 )
679              i=`echo $ii | sed 's:^\./::'`              RETVAL=$?
680              name=`echo $i | sed 's:_mpi::' `              if [ $verbose -gt 1 ]; then echo "  cmp $name $file returns: $RETVAL" ; fi
681              if test -L $BUILD_DIR/$name ; then              if test "x$RETVAL" != x0 ; then
682                  linktarg=`(cd $BUILD_DIR; readlink $name)`                  if test -h $build_dir/$name ; then rm -f $build_dir/$name ; fi
683                  if test $linktarg = "../$code_dir/$name"_mpi ; then                  if test ! -r $build_dir/$name ; then
684                  #echo Un-linking $name from $linktarg                      if [ $verbose -gt 1 ]; then echo "  Linking $name to $file" ; fi
685                      rm -f $BUILD_DIR/$name                      (cd $build_dir; ln -sf $file $name)
686                  fi                  fi
687              fi              fi
688          done          fi
689      fi      done
       
690  }  }
691    
692  linkdata()  linkdata()
693  {  {
694      # linkdata flag      # linkdata run_dir input_dir_1 input_dir_2 ...
695      #      #
696      # symbolically link data files to run directory      # symbolically link data files to run directory
697      if test "x$1" = x1 ; then      if test -d $1 ; then
698          (          (
699              cd $2              cd $1 ; shift
700              if test "x$ADM" = x ; then              echo 'linkdata from dirs:' $*
701                  files=`( cd ../input ; ls -1 | grep -v CVS )`              inpMPI=`(cd ../$1 ; find . -name "*.mpi" -print | sed 's:^\./::')`
702                  for i in $files ; do              for xx in $inpMPI ; do
703                      if test ! -d "../input/"$i ; then                if test -r "../"$1"/"$xx ; then
704                          ln -sf "../input/"$i $i                  # found 1 _mpi sfx file in 1rst input dir and it is readable
705                      fi                  yy=`echo $xx | sed 's:\.mpi$::'`
706                  done                  if test "x$MPI" = "x0" ; then
707              else                      # not mpi test: remove symbolic link
708                  files=`( cd ../input ; ls -1 *.bin | grep -v CVS )`                      if test -h $yy ; then rm -f $yy ; fi
709                  for i in $files ; do                  else
710                      if test ! -d "../input/"$i ; then                      # mpi test: remove symbolic link & link .mpi sfx file
711                          ln -sf "../input/"$i $i                      if test -h $yy ; then rm -f $yy ; fi
712                        if test ! -r $yy ; then
713                            ln -sf "../"$1"/"$xx $yy ;
714                            printf " $xx" 1>&2
715                      fi                      fi
716                  done                  fi
717                  files=`( cd ../input_ad ; ls -1 | grep -v CVS )`                fi
718                  for i in $files ; do              done
719                      if test ! -d "../input_ad/"$i ; then              if test -r "../"$1"/eedata.mth" ; then
720                          ln -sf "../input_ad/"$i $i              # found eedata.mth in 1rst input dir and it is readable
721                    if test "x$MULTI_THREAD" = "xt" ; then
722                    # multi-threaded test: remove symbolic link & link eedata.mth
723                        if test -h eedata ; then rm -f eedata ; fi
724                        if test ! -r eedata ; then
725                            ln -sf "../"$1"/eedata.mth" eedata ;
726                            printf ' eedata.mth' 1>&2
727                      fi                      fi
728                  done                  else
729                    # not multi-threaded test: remove eedata symbolic link
730                        if test -h eedata ; then rm -f eedata ; fi
731                    fi
732              fi              fi
733                prevDir='NONE'
734                for ldir in $* ; do
735                    if test -d "../"$ldir -a $ldir != $prevDir ; then
736                        printf " ldir=${ldir}:" 1>&2
737                        files=`( cd "../"$ldir ; ls -1 | grep -v CVS )`
738                        for i in $files ; do
739                            if test ! -d "../"$ldir/$i ; then
740                                if test ! -r $i  ; then
741                                    printf ' '$i 1>&2
742                                    ln -sf "../"$ldir"/"$i $i
743                                fi
744                            fi
745                        done
746                        printf ' ;' 1>&2
747                        if test -x "../"$ldir"/"prepare_run ; then
748                            "../"$ldir"/"prepare_run 1>&2
749                        else
750                            echo '' 1>&2
751                        fi
752                    fi
753                    prevDir=$ldir
754                done
755          )          )
756      fi      fi
757  }  }
# Line 407  runmodel() Line 760  runmodel()
760  {  {
761      # runmodel directory      # runmodel directory
762      #      #
763      #  runs "$COMMAND in "directory"      #  runs $COMMAND in "directory"
764      #  (where "$COMMAND" is relative to "directory")      #  (where "$COMMAND" is relative to "directory")
765      (      (
766          cd $1          cd $1
767          printf 'runmodel ... ' 1>&2          printf 'runmodel in %s ... ' $1
768          # make output.txt          if test "x$MPI" != x0 ; then
769          echo              #- adjust the MPI run command with the right number of Procs
770          echo "COMMAND='$COMMAND'"              #echo '' ; echo "  COMMAND='$COMMAND'"
771          echo "pwd='"`pwd`"'"              COMMAND=`echo $COMMAND | sed "s/ TR_NPROC / $LOC_NPROC /"`
772          ( eval $COMMAND ) >> run.log 2>&1              if test "x$MPI_MFILE" != x ; then
773          RETVAL=$?                COMMAND=`echo $COMMAND | sed "s/ TR_MFILE / ..\/..\/$LOC_MFILE /"`
774          if test "x$RETVAL" = x0 ; then              fi
775              echo successful 1>&2              #echo "  COMMAND='$COMMAND'"
776              if test "x$ADM" = x ; then          fi
777                  cp output.txt $CDIR"/output.txt"          if test -L $EXECUTABLE ; then
778              else            if test -x "../"$builddir"/"$EXECUTABLE ; then
779                  cp output.txt_adm $CDIR"/output.txt_adm"              cmp $EXECUTABLE "../"$builddir"/"$EXECUTABLE > /dev/null 2>&1
780                outD=$? ; if test "x$outD" != x0 ; then rm -f $EXECUTABLE ; fi
781              else rm -f $EXECUTABLE
782              fi
783            fi
784            if test ! -x $EXECUTABLE -a -x "../"$builddir"/"$EXECUTABLE ; then
785                echo " link" $EXECUTABLE "from dir ../"$builddir > run.log_tmp
786                ln -sf "../"$builddir"/"$EXECUTABLE .
787            fi
788            if test ! -x $EXECUTABLE ; then
789                rm -f $RUNLOG ; touch $RUNLOG
790                if test -f run.log_tmp ; then cat run.log_tmp >> $RUNLOG ; fi
791                echo " no executable:" $EXECUTABLE >> $RUNLOG
792                RETVAL=8
793                ENDVAL=-1
794            else
795              if test ! -f $OUTPUTFILE -o $OUTPUTFILE -ot $EXECUTABLE ; then
796                # output do not exist or is older than executable:
797                rm -f $OUTPUTFILE $RUNLOG ; touch $RUNLOG
798                if test -f run.log_tmp ; then cat run.log_tmp >> $RUNLOG ; fi
799            #- Divided Adjoint Run:
800            #  get the number of additional runs (add_DIVA_runs) from file "run_ADM_DIVA"
801                if test $KIND = 2 -a -f run_ADM_DIVA ; then
802                  adm_diva_nb=`sed -n '/^ *add_DIVA_runs\>.*=/p' run_ADM_DIVA | sed 's/ //g'`
803                  echo " Divided Adjoint Run: $adm_diva_nb" >> $RUNLOG
804                  eval "let $adm_diva_nb"
805                  if [ $add_DIVA_runs -ge 1 ] ; then
806                    rm -f costf* divided.ctrl snapshot*
807                    echo -n "(add_DIVA_runs=$add_DIVA_runs) ... "
808                    for ii in `seq 1 $add_DIVA_runs` ; do
809                      ( eval $COMMAND ) >> $RUNLOG 2>&1
810                      echo " additional DIVA run # $ii : done" >> $RUNLOG
811                      mv -f $OUTPUTFILE ${OUTPUTFILE}.diva_${ii}
812                    done
813                  fi
814                fi
815            #- special DIVA processing ends here
816                ( eval $COMMAND ) >> $RUNLOG 2>&1
817                RETVAL=$?
818                ENDVAL=`tail $OUTPUTFILE | grep -c 'PROGRAM MAIN: Execution ended Normally'`
819                if [ $POSTCLEAN -eq 1 -a $ENDVAL -gt 0 ] ; then
820                    find . -name "*.meta" -exec rm {} \;
821                    find . -name "*.data" -exec rm {} \;
822                    rm -rf mnc_test_*
823              fi              fi
824              else
825                RETVAL=0
826                ENDVAL=`tail $OUTPUTFILE | grep -c 'PROGRAM MAIN: Execution ended Normally'`
827                touch $RUNLOG
828                if test -f run.log_tmp ; then cat run.log_tmp >> $RUNLOG ; fi
829                echo "---------->> $OUTPUTFILE is up to date " >> $RUNLOG 2>&1
830              fi
831            fi
832            rm -f run.log_tmp
833            #- in all cases where OutputFile exists, report SIZE and time
834            if test -f $OUTPUTFILE ; then
835              grep '(PID\.TID 0000\.0001)      n.. =' $OUTPUTFILE \
836                    | sed 's/(PID.TID 0000.0001)   //' >> $CDIR"/summary.txt"
837              grep -A3 'Seconds in section "ALL' $OUTPUTFILE \
838                    | sed 's/(PID.TID 0000.0001)   //' >> $CDIR"/summary.txt"
839            fi
840            if test -s STDERR.0000 ; then cp STDERR.0000 $CDIR"/STDERR.0000" ; fi
841            if [ $RETVAL -eq 0 -a $ENDVAL -gt 0 ] ; then
842                echo successful
843                printf '=> output from running in %s :\n' $1 1>&2
844                tail $RUNLOG | sed 's/^.*/> &/g' 1>&2
845                return 0
846            elif [ $RETVAL -ne 0 -a $ENDVAL -gt 0 ] ; then
847                #-- for some weird cases (run is finihed but with error code)
848                echo 'finished with error (run:' $RETVAL ' end:' $ENDVAL ')'
849                printf '=> output from running in %s :\n' $1 1>&2
850                tail $RUNLOG | sed 's/^.*/> &/g' 1>&2
851              return 0              return 0
852          else          else
853              tail run.log              echo 'failed (run:' $RETVAL ' end:' $ENDVAL ')'
854              echo failed 1>&2              printf '=> output from running in %s :\n' $1 1>&2
855              cp run.log $CDIR"/run.log"              tail $RUNLOG | sed 's/^.*/> &/g' 1>&2
856                cp $RUNLOG $CDIR"/"$RUNLOG
857              return 1              return 1
858          fi          fi
859      )      )
# Line 439  createcodelet() Line 863  createcodelet()
863  {  {
864      # create codelet for comparing model output      # create codelet for comparing model output
865    
866      printf "creating the comparison code...  "      printf "creating the comparison code (using CC=$CC)...  "
867      cat > tmp_cmpnum.c <<EOF      cat > tr_cmpnum.c <<EOF
868  #include <stdio.h>  #include <stdio.h>
869  #include <math.h>  #include <math.h>
870  int main( int argc, char** argv )  {  int main( int argc, char** argv )  {
871    int linnum,best,lncnt;    int linnum,cmplin,best,lncnt;
872    double a,b,diff;    double a,b,abave,relerr;
873    best = -16;    best = -22;
874    lncnt = 0;    lncnt = 0;
875    while( 1 & (lncnt+=1) < 999 )  {    while( 1 & ( (lncnt+=1) < 999 ) )  {
876      scanf("%d", &linnum);      scanf("%d", &linnum);
877      if (linnum == -1)  break;      if (linnum == -1)  break;
878      scanf("%lf", &a);  scanf("%lf", &b);      scanf("%lf", &a);  scanf("%lf", &b);
879      diff = 0.5*(fabs(a)+fabs(b));      abave = 0.5*(fabs(a)+fabs(b));
880      if (diff > 1.e-12) {      if ( abave == abave ) {
881        diff=fabs(a-b)/diff;        if (abave > 0.0) {
882        if (diff > 0.0) {          relerr=fabs(a-b)/abave;
883          linnum = (int)log10(diff);          if (relerr > 0.0) { cmplin = (int)rint(log10(relerr)); }
884          best = (best > linnum) ? best : linnum;          else { cmplin = -16 ; }
885            best = (best > cmplin) ? best : cmplin; }
886          else { cmplin = -22 ; }
887       /* printf("%d ; %lf ; %lf\n",cmplin,a,b); */
888        }        }
889        else {     else {
890          if (best == -16 && diff != 0)  best = -22;     /* printf("%lf ; %lf ; %lf\n",abave,a,b); */
891        }        break; }
     }  
892    }    }
893    if (lncnt == 999) best=-29;    if (lncnt == 999) best=-29;
894      if (linnum != -1) best=-99;
895    printf("%d\n", -best);    printf("%d\n", -best);
896    return 0;    return 0;
897  }  }
898  EOF  EOF
899      cc -o tmp_cmpnum tmp_cmpnum.c -lm      $CC -o tr_cmpnum tr_cmpnum.c -lm
900    
901      if [ -x ./tmp_cmpnum ]; then      if [ -x ./tr_cmpnum ]; then
902          echo "OK"          echo "OK"
903          return 0          return 0
904      else      else
905          echo          echo
906          echo "ERROR: failed to compile comparison code"          echo "ERROR: failed to compile comparison code -- please specify"
907            echo "  a C compiler using the CC environment variable."
908          exit 1          exit 1
909      fi      fi
910  }  }
# Line 488  formatresults() Line 916  formatresults()
916      nm=$1      nm=$1
917      printf '%s %s %s %s' $2 $3 $4 $5      printf '%s %s %s %s' $2 $3 $4 $5
918      shift; shift; shift; shift; shift;      shift; shift; shift; shift; shift;
919      printf '%3s' $@      listPrt=$@
920            listRes=`echo $listPrt | sed 's/>//' | sed 's/<//'`
921      if [ $1 = '--' ]; then      xx=`echo $listPrt | sed 's/.*>//' | sed 's/<.*//' | awk '{print $1}'`
922        printf '%3s' $listPrt
923    #   line below does not work on hp-ux_ia64 : do those substitutions later on
924    #   printf '%3s' $listPrt | sed 's/ 99/ --/g' | sed 's/  > />/' | sed 's/  < /</'
925    
926        if [ $xx = '..' ]; then
927            printf ' N/O '
928        elif [ $xx = '--' ]; then
929            printf ' N/O '
930        elif [ $xx = 99 ]; then
931          printf ' N/O '          printf ' N/O '
932      else      else
933          if [ $1 -gt 12 ]; then          if [ $xx -ge $MATCH_CRIT ]; then
934              printf ' pass'              printf ' pass'
935          else          else
936              printf ' FAIL'              printf ' FAIL'
937          fi          fi
938      fi      fi
939      printf '  %s' $nm      printf '  %s' $nm
940        if test $KIND = 2 ; then
941            #-- apppend taf repport summary:
942            tafrep=`grep -c '^ TAF reports ' $CDIR/summary.txt`
943            if test $tafrep = 1 ; then
944                grep '^ TAF reports ' $CDIR/summary.txt | awk '{printf "  (e=%i, w=%i)",$3,$6}'
945            fi
946        fi
947      printf '\n'      printf '\n'
       
 }  
   
 show_help()  
 {  
     cat - << EOF  
 $0 [-help] [-quick] [-verbose] dir1 [dir2] [...]  
   
  -help|-h      Show this help message  
  -quiet     Reduce the amount of output  
  -verbose   Produce copious amounts of output  
  -debug     Produce even more output which will mean nothing to most  
  -force     Do "make CLEAN" before compiling. This forces a complete rebuild.  
  -clean     Do "make CLEAN" after compiling and testing.  
  -cleanup   Aggresively removes all model output, executables and object files  
             and then exits. Use with care.  
948    
 Normal usage:  
  $0 *       Configure, compile, run and analyze in all experiment directories  
 EOF  
949  }  }
950    
951  scandirs()  scandirs()
952  {  {
953      if [ $# -eq 0 ]; then      if [ $# -eq 1 ]; then
954          for arg in * ; do          for arg in * ; do
955              test -d $arg/input && echo $arg             #test -f $arg/$1 && echo $arg
956          done              test -f $arg/$1 -o -f $arg/$1.gz && echo $arg
957            done
958      else      else
959          echo $*          echo $*
960      fi      fi
961  }  }
962    
963    
964    check_eedata()
965    {
966        # check_eedata eedata size.h
967        if [ $# -eq 2 ] ; then
968         if test -f $1 -a -f $2 ; then
969          nx=`grep "^ *nTx *=" $1 | tail -1 | sed 's/^ *nTx *= *//' | sed "s/, *$//"`
970          sx=`grep "^     & *nSx *=" $2 | sed "s/^     & *nSx *=//" | sed 's/, *$//'`
971          if test "x$nx" = x ; then
972            rx=10
973          else
974            rx=`expr $sx % $nx`
975          fi
976          ny=`grep "^ *nTy *=" $1 | tail -1 | sed 's/^ *nTy *= *//' | sed "s/, *$//"`
977          sy=`grep "^     & *nSy *=" $2 | sed "s/^     & *nSy *=//" | sed 's/, *$//'`
978          if test "x$ny" = x ; then
979            ry=20
980          else
981            ry=`expr $sy % $ny`
982          fi
983          echo `expr $rx + $ry`
984         else
985          echo '-1'
986         fi
987        elif [ $# -eq 1 ] ; then
988         if test -f $1 ; then
989          nx=`grep "^ *nTx *=" $1 | tail -1 | sed 's/^ *nTx *= *//' | sed "s/, *$//"`
990          if test "x$nx" = x ; then nx=1 ; fi
991          ny=`grep "^ *nTy *=" $1 | tail -1 | sed 's/^ *nTy *= *//' | sed "s/, *$//"`
992          if test "x$ny" = x ; then ny=1 ; fi
993         #echo $nx $ny
994          echo $nx
995         else
996          echo '-1'
997         fi
998        else
999          echo '-2'
1000        fi
1001    
1002    }
1003    
1004  ###############################################################################  ###############################################################################
1005    
1006    
1007  #  Default properties  #  Default properties
1008  debug=0  debug=0
1009  verbose=1  verbose=1
1010  clean=0  NBLINES_MKLOG=16000
 expts=''  
 # ieee=1  
   
 IEEE=  
 if test "x$MITGCM_IEEE" != x ; then  
     IEEE=$MITGCM_IEEE  
 fi  
1011    
1012    OptLev=1
1013    GSL=f
1014    
1015  CLEANUP=f  CLEANUP=f
1016    NORUN=f
1017  QUICK=f  QUICK=f
1018    NOMAKE=f
1019  NOGENMAKE=f  NOGENMAKE=f
1020  NOCLEAN=f  NOCLEAN=f
1021  NODEPEND=f  NODEPEND=f
1022    POSTCLEAN=0
1023    
1024  BASH=  BASH=
1025  OPTFILE=NONE  OPTFILE=NONE
1026  ADDRESSES=  ADDRESSES=
1027  TESTDIRS=  TESTDIRS=
1028    SKIPDIRS=
1029  MPACKDIR="../tools/mpack-1.6"  MPACKDIR="../tools/mpack-1.6"
1030  HAVE_MPACK=  HAVE_MPACK=
1031  MPACK="$MPACKDIR/mpack"  MPACK=
1032  COMMAND=  COMMAND=
1033  MAKE=make  MKDEPEND=
1034  MPI=f  if test "x$MAKE" = x ; then
1035        MAKE=make
1036    fi
1037    REPLMAKE=
1038    if test "x$CC" = x ; then
1039        CC=cc
1040    fi
1041    JOBS=
1042    TARG=
1043    MPI=0
1044    MPI_MFILE=
1045    MULTI_THREAD=f
1046    OUTDIR=
1047    DELDIR=
1048    USE_R4=
1049    EXTRFLG=
1050    
1051    #- type of testing (KIND):
1052    #   KIND=0 : forward (= default) ;  KIND=1 : Tangent Linear with TAF ;
1053    #   KIND=2 : Adjoint with TAF ;     KIND=4 : Adjoint with OpenAD
1054    KIND=0
1055    
1056  ADM=  # list of pTracers to check for monitor output
1057    PTRACERS_NUM="1 2 3 4 5"
1058    
1059    MATCH_CRIT=10
1060    
1061  printf "parsing options...  "  printf "parsing options...  "
1062    
# Line 583  for ac_option ; do Line 1071  for ac_option ; do
1071      fi      fi
1072    
1073      ac_optarg=`expr "x$ac_option" : 'x[^=]*=\(.*\)'`      ac_optarg=`expr "x$ac_option" : 'x[^=]*=\(.*\)'`
1074        
1075      case $ac_option in      case $ac_option in
1076                    
1077          -help | --help | -h | --h)          -help | --help | -h | --h) usage ;;
             usage ;;  
           
         -optfile | --optfile | -of | --of)  
             ac_prev=OPTFILE ;;  
         -optfile=* | --optfile=* | -of=* | --of=*)  
             OPTFILE=$ac_optarg ;;  
           
         -addr | --addr | -a | --a)  
             ac_prev=ADDRESSES ;;  
         -addr=* | --addr=*)  
             ADDRESSES=$ac_optarg ;;  
   
         -tdir | --tdir | -t | --t)  
             ac_prev=TESTDIRS ;;  
         -tdir=* | --tdir=*)  
             TESTDIRS=$ac_optarg ;;  
   
         -bash | --bash | -b | --b)  
             ac_prev=BASH ;;  
         -bash=* | --bash=*)  
             BASH=$ac_optarg ;;  
   
         -command | --command | -c | --c)  
             ac_prev=COMMAND ;;  
         -command=* | --command=*)  
             COMMAND=$ac_optarg ;;  
   
         -make | --make | -m | --m)  
             ac_prev=MAKE ;;  
         -make=* | --make=*)  
             MAKE=$ac_optarg ;;  
   
         -clean | --clean)  
             CLEANUP=t ;;  
   
         -quick | --quick | -q | --q)  
             QUICK=t ;;  
         -nogenmake | --nogenmake | -ng | --ng)  
             NOGENMAKE=t ;;  
         -noclean | --noclean | -nc | --nc)  
             NOCLEAN=t ;;  
         -nodepend | --nodepend | -nd | --nd)  
             NODEPEND=t ;;  
   
         -mpi) MPI=t ;;  
1078    
1079          -adm | -ad) ADM=t ;;          -optfile | --optfile | -of | --of) ac_prev=OPTFILE ;;
1080            -optfile=* | --optfile=* | -of=* | --of=*) OPTFILE=$ac_optarg ;;
1081    
1082          -ieee) IEEE=true ;;          -addr | --addr | -a | --a) ac_prev=ADDRESSES ;;
1083          -noieee) IEEE= ;;          -addr=* | --addr=* | -a=* | --a=*) ADDRESSES=$ac_optarg ;;
1084            -mpackdir | --mpackdir | -mpd | --mpd) ac_prev=MPACKDIR ;;
1085            -mpackdir=* | --mpackdir=* | -mpd=* | --mpd=*) MPACKDIR=$ac_optarg ;;
1086    
1087            -tdir | --tdir | -t | --t) ac_prev=TESTDIRS ;;
1088            -tdir=* | --tdir=* | -t=* | --t=*) TESTDIRS=$ac_optarg ;;
1089            -skipdir | --skipdir | -skd | --skd) ac_prev=SKIPDIRS ;;
1090            -skipdir=* | --skipdir=* | -skd=* | --skd=*) SKIPDIRS=$ac_optarg ;;
1091    
1092            -bash | --bash | -b | --b) ac_prev=BASH ;;
1093            -bash=* | --bash=* | -b=* | --b=*) BASH=$ac_optarg ;;
1094    
1095            -command | --command | -c | --c) ac_prev=COMMAND ;;
1096            -command=* | --command=* | -c=* | --c=*) COMMAND=$ac_optarg ;;
1097    
1098            -makedepend | --makedepend | -md | --md) ac_prev=MKDEPEND ;;
1099            -makedepend=* | --makedepend=* | -md=* | --md=*) MKDEPEND=$ac_optarg ;;
1100    
1101            -make | --make | -m | --m) ac_prev=MAKE ;;
1102            -make=* | --make=* | -m=* | --m=*) MAKE=$ac_optarg ;;
1103    
1104            -repl_mk | --repl_mk ) ac_prev=REPLMAKE ;;
1105            -repl_mk=* | --repl_mk=*) REPLMAKE=$ac_optarg ;;
1106    
1107            -odir | --odir) ac_prev=OUTDIR ;;
1108            -odir=* | --odir=*) OUTDIR=$ac_optarg ;;
1109    
1110            -ptracers | --ptracers | -ptr | --ptr) ac_prev=PTRACERS_NUM ;;
1111            -ptracers=* | --ptracers=* | -ptr=* | --ptr=*) PTRACERS_NUM=$ac_optarg ;;
1112    
1113            -match | --match ) ac_prev=MATCH_CRIT ;;
1114            -match=* | --match=* ) MATCH_CRIT=$ac_optarg ;;
1115    
1116            -j | --j) ac_prev=JOBS ;;
1117            -j=* | --j=*) JOBS=$ac_optarg ;;
1118    
1119            -ef | --ef) ac_prev=EXTRFLG ;;
1120            -ef=* | --ef=*) EXTRFLG=$ac_optarg ;;
1121    
1122            -clean | --clean) CLEANUP=t ; DELDIR=t ;;
1123    
1124            -norun | --norun | -nr | --nr) NORUN=t ;;
1125            -obj | --obj ) TARG='obj' ; NORUN=t ;;
1126            -runonly | --runonly | -ro | --ro) QUICK=t ; NOMAKE=t ;;
1127            -quick | --quick | -q | --q) QUICK=t ;;
1128            -nogenmake | --nogenmake | -ng | --ng) NOGENMAKE=t ;;
1129            -noclean | --noclean | -nc | --nc) NOCLEAN=t ;;
1130            -nodepend | --nodepend | -nd | --nd) NODEPEND=t ;;
1131    
1132            -postclean | --postclean | -pc | --pc) POSTCLEAN=2 ;;
1133            -deloutp | --deloutp | -do | --do) POSTCLEAN=1 ;;
1134    
1135            -mpi | --mpi) MPI=2 ;;
1136            -MPI | --MPI) ac_prev=MPI ;;
1137            -MPI=* | --MPI=*) MPI=$ac_optarg ;;
1138    
1139            -mfile | --mfile | -mf | --mf) ac_prev=MPI_MFILE ;;
1140            -mfile=* | --mfile=* | -mf=* | --mf=*) MPI_MFILE=$ac_optarg ;;
1141    
1142            -mth) MULTI_THREAD=t ;;
1143    
1144            -tlm) if test $KIND = 0 ; then KIND=1 ; else
1145                    echo "Error: '-tlm', '-adm' and '-oad' are exclusive" ; usage
1146                  fi ;;
1147            -adm | -ad) if test $KIND = 0 ; then KIND=2 ; else
1148                    echo "Error: '-tlm', '-adm' and '-oad' are exclusive" ; usage
1149                  fi ;;
1150            -oad) if test $KIND = 0 ; then KIND=4 ; NODEPEND=t ; else
1151                    echo "Error: '-tlm', '-adm' and '-oad' are exclusive" ; usage
1152                  fi ;;
1153    
1154            -ieee)   echo "Warning: ignore option '-ieee' (already the default)"
1155                     printf " ... " ;;
1156            -noieee) echo "Warning: will use option '-fast' instead of '-noieee' (obsolete)"
1157                     printf " ... " ; OptLev=`expr $OptLev \* 2` ;;
1158            -fast)  OptLev=`expr $OptLev \* 2` ;;
1159            -devel) OptLev=0 ;;
1160            -gsl) GSL=t ;;
1161    
1162          -verbose) verbose=2 ;;          -verbose) verbose=2 ;;
1163          -debug) debug=1 ;;          -debug) debug=1 ;;
1164          -quiet) verbose=0 ;;          -quiet) verbose=0 ;;
1165    
1166          -*)          -deldir | -dd) DELDIR=t ;;
1167              echo "Error: unrecognized option: "$ac_option  
1168              usage          -use_r4|-ur4) USE_R4=t ;;
1169              ;;  
1170                    -ts) TS=t;;
1171          *)          -papis) PAPIS=t;;
1172              echo "Error: unrecognized argument: "$ac_option          -pcls) PCL=t;;
1173              usage  
1174              ;;          -*) echo "Error: unrecognized option: "$ac_option
1175                        usage ;;
1176            *)  echo "Error: unrecognized argument: "$ac_option
1177                usage ;;
1178    
1179       esac       esac
1180        
1181  done  done
1182    
1183  if test "x$QUICK" = xt ; then  if test "x$QUICK" = xt ; then
# Line 662  if test "x$QUICK" = xt ; then Line 1186  if test "x$QUICK" = xt ; then
1186      NODEPEND=t      NODEPEND=t
1187  fi  fi
1188    
1189    #- check length of MPI machine file:
1190    if test "x$MPI" != x0 -a "x$MPI_MFILE" != x ; then
1191        if test -r $MPI_MFILE ; then
1192            nl=`wc -l $MPI_MFILE | awk '{print $1}'`
1193            if [ $nl -lt $MPI ] ; then
1194              echo "Error: need at least $MPI nodes (currently only $nl) in MPI_MFILE=$MPI_FILE"
1195              usage
1196            fi
1197            if [ $verbose -gt 1 ]; then
1198                echo " MPI_MFILE=$MPI_MFILE : $nl procs for MPI=$MPI run"
1199            fi
1200        else
1201              echo "Error: cannot access MPI_MFILE=$MPI_FILE"
1202              usage
1203        fi
1204    fi
1205    
1206    #- setting for forward or ADM testing
1207    if test $KIND = 1 ; then
1208        TARG=ftlall
1209        code_dir=code_ad
1210        inputdir=input_ad
1211        ref_outp="output_tlm.txt"
1212        EXECUTABLE="mitgcmuv_ftl"
1213    elif test $KIND = 2 ; then
1214        if test "x$TARG" = xobj ; then TARG=adobj ; else TARG=adall ; fi
1215        code_dir=code_ad
1216        inputdir=input_ad
1217        ref_outp="output_adm.txt"
1218        EXECUTABLE="mitgcmuv_ad"
1219    elif test $KIND = 4 ; then
1220        TARG=adAll
1221        code_dir=code_oad
1222        inputdir=input_oad
1223        ref_outp="output_oadm.txt"
1224        EXECUTABLE="mitgcmuv_ad"
1225    else
1226        if test "x$JOBS" != x ; then TARG="-j $JOBS $TARG" ; fi
1227        code_dir=code
1228        inputdir=input
1229        ref_outp="output.txt"
1230        EXECUTABLE="mitgcmuv"
1231    fi
1232    
1233    xx=`echo $TESTDIRS | awk '{print $1}'`
1234  if test "x$TESTDIRS" = x ; then  if test "x$TESTDIRS" = x ; then
1235      TESTDIRS=`scandirs`      LIST=`scandirs results/$ref_outp`
1236    elif test $xx = 'start_from' ; then
1237        xx=`echo $TESTDIRS | awk '{print $2}'`
1238        LIST=`scandirs results/$ref_outp | sed -n "/$xx/,$ p"`
1239    else
1240        #- expand group of experiments:
1241        LIST=" "
1242        for xx in $TESTDIRS
1243        do
1244          case $xx in
1245            'basic') LIST=${LIST}" aim.5l_cs hs94.128x64x5 ideal_2D_oce"
1246                     LIST=${LIST}" lab_sea tutorial_baroclinic_gyre"
1247                     LIST=${LIST}" tutorial_global_oce_latlon tutorial_plume_on_slope"
1248                    ;;
1249            'tutorials')
1250                     LIST=${LIST}" "`ls | grep 'tutorial_'` ;;
1251            *)       LIST=${LIST}" "$xx ;;
1252          esac
1253        done
1254  fi  fi
1255    #echo 'LIST='${LIST}'<'
1256    #- skip dirs, remove duplicate and non-directory:
1257    TESTDIRS=" "
1258    count=0
1259    for xx in $LIST
1260    do
1261        yy=`echo $SKIPDIRS | grep -c $xx`
1262        if test $yy = 0 ; then
1263            if test -d $xx ; then
1264                yy=`echo $TESTDIRS | grep -c $xx`
1265                if test $yy = 0 ; then TESTDIRS=${TESTDIRS}" "$xx ; fi
1266            else count=1 ;
1267                echo ""; echo -n " -- skip \"$xx\" (not a directory !)"
1268            fi
1269        else
1270            if test $count = 1 ; then echo -n ", \"$xx\""
1271            else count=1 ; echo "" ;  echo -n " skip: \"$xx\""
1272            fi
1273        fi
1274    done
1275    if test $count = 1 ; then echo "" ; echo -n " ... " ; fi
1276    #echo 'TESTDIRS='${TESTDIRS}'<'
1277    
1278  if test "x$OPTFILE" = xNONE -a "x$MITGCM_OF" != x ; then  if test "x$OPTFILE" = xNONE -a "x$MITGCM_OF" != x ; then
1279      OPTFILE=$MITGCM_OF      OPTFILE=$MITGCM_OF
1280  fi  fi
1281    
1282  if test "x$ADM" = xt -a "x$COMMAND" = x ; then  LOC_MFILE='tr_mpi_mfile'
1283      COMMAND="./mitgcmuv_ad > output.txt_adm 2>&1"  RUNLOG="run.tr_log"
1284  fi  OUTPUTFILE=$ref_outp
   
1285  if test "x$COMMAND" = x ; then  if test "x$COMMAND" = x ; then
1286      COMMAND="make output.txt"      COMMAND="./$EXECUTABLE > $OUTPUTFILE"
1287    fi
1288    if test "x$MPI" != x0 ; then
1289        OUTPUTFILE="STDOUT.0000"
1290  fi  fi
1291    
1292  echo "OK"  echo "OK (COMMAND= $COMMAND )"
1293    
1294    # set the Default List of output variables to be checked:
1295    #  (use default or load experiment-specific list from file "tr_checklist")
1296    # content : 1rst = main variable used to decide if it pass or FAIL
1297    #         others = number of matching digits to be printed in summary.txt
1298    if test $KIND = 1 ; then
1299        DEF_CHECK_LIST='tlmGrd tlmCst tlmGrd fwdGrd'
1300        EMPTY_RESULTS='.. .. ..'
1301        LEN_CHECK_LIST=`echo $DEF_CHECK_LIST | sed 's/ [a-zA-Z0-9]*+/&mn &mx &av &sd/g' | awk '{print NF-1}'`
1302    elif test $KIND = 2 -o $KIND = 4 ; then
1303        DEF_CHECK_LIST='AdGrd Cost AdGrd FDGrd'
1304        EMPTY_RESULTS='.. .. ..'
1305        LEN_CHECK_LIST=`echo $DEF_CHECK_LIST | sed 's/ [a-zA-Z0-9]*+/&mn &mx &av &sd/g' | awk '{print NF-1}'`
1306    else
1307        DEF_CHECK_LIST='PS PS T+ S+ U+ V+ pt1+ pt2+ pt3+ pt4+ pt5+'
1308        EMPTY_RESULTS='.. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. ..'
1309        LEN_CHECK_LIST=`echo $DEF_CHECK_LIST | sed 's/ [a-zA-Z0-9]*+/&mn &mx &av &sd/g' | awk '{print NF-1}'`
1310        ii=`echo $EMPTY_RESULTS | awk '{print NF}'`
1311        EMPTY_RESULTS=$EMPTY_RESULTS`expr $LEN_CHECK_LIST - $ii | awk 'BEGIN{FS=":"}{for(i=1;i<=$1;i++){printf "  ."}}'`
1312    fi
1313    
1314  #  create the FORTRAN comparison code  #  create the FORTRAN comparison code
1315  createcodelet  if test "x$CLEANUP" = xt -o -x tr_cmpnum ; then
1316        echo "skipping comparison code build"
1317    else
1318        createcodelet
1319    fi
1320    
1321  #  build the mpack utility  #  build the mpack utility (if ADDRESSES = NONE, do it to test the build)
1322  if test "x$ADDRESSES" = xNONE -o "x$ADDRESSES" = x ; then  if test "x$ADDRESSES" = x ; then
1323      echo "skipping mpack build"      echo "skipping mpack build"
1324  else  else
1325      build_mpack      build_mpack
1326  fi  fi
1327    
1328  #  Create a uniquely named directory to store results  #  Create a uniquely named directory to store results
1329    CMDLINE=$0
1330    for xx in "$@" ; do nw=`echo $xx | wc -w`
1331        if test $nw = '1' ; then CMDLINE="$CMDLINE $xx" ; else
1332          nb=`echo $xx | grep -c '='`
1333          if test $nb = 0 ; then CMDLINE="$CMDLINE '$xx'"
1334          else yy=`echo "$xx'" | sed "s/=/='/"` ;
1335                                 CMDLINE="$CMDLINE $yy" ; fi
1336        fi
1337    done
1338    #for xx in "$@" ; do CMDLINE="$CMDLINE '$xx'" ; done
1339  MACH=`hostname`  MACH=`hostname`
1340  UNAMEA=`uname -a`  UNAMEA=`uname -a`
1341  DATE=`date +%Y%m%d`  DATE=`date +%Y%m%d`
1342  BASE="tr_"$MACH"_"$DATE"_"  BASE="tr_"$MACH"_"$DATE"_"
1343    if test "x$OUTDIR" != x ; then
1344        BASE="tr_"$OUTDIR"_"$DATE"_"
1345    else
1346        short_name=`hostname | sed 's/\..*$//'`
1347        BASE="tr_"$short_name"_"$DATE"_"
1348    fi
1349  DNUM=0  DNUM=0
1350  DRESULTS="$BASE$DNUM"  DRESULTS="$BASE$DNUM"
1351  while test -e $DRESULTS ; do  while test -e $DRESULTS ; do
# Line 708  if test "x$RETVAL" != x0 ; then Line 1359  if test "x$RETVAL" != x0 ; then
1359      exit 1      exit 1
1360  fi  fi
1361  SUMMARY="$DRESULTS/summary.txt"  SUMMARY="$DRESULTS/summary.txt"
 printf "Start time:  " >> $SUMMARY  
1362  start_date=`date`  start_date=`date`
1363  echo $start_date > $SUMMARY  echo $start_date > $SUMMARY
1364    echo 'run:' $CMDLINE >> $SUMMARY
1365    echo 'on :' $UNAMEA  >> $SUMMARY
1366    
1367  of_path=  of_path=
1368  if test "x$OPTFILE" != xNONE ; then  if test "x$OPTFILE" != xNONE ; then
# Line 728  if test "x$OPTFILE" != xNONE ; then Line 1380  if test "x$OPTFILE" != xNONE ; then
1380          echo >> $SUMMARY          echo >> $SUMMARY
1381          echo "  OPTFILE=$OPTFILE" >> $SUMMARY          echo "  OPTFILE=$OPTFILE" >> $SUMMARY
1382      else      else
1383          echo | tee $SUMMARY          echo | tee -a $SUMMARY
1384          echo "ERROR: can't read OPTFILE=\"$OPTFILE\"" | tee $SUMMARY          echo "ERROR: can't read OPTFILE=\"$OPTFILE\"" | tee -a $SUMMARY
1385          exit 1          exit 1
1386      fi      fi
1387  else  else
1388      echo >> $SUMMARY      echo >> $SUMMARY
1389      echo "No \"OPTFILE\" was explicitly specified by testreport," >> $SUMMARY      echo "No \"OPTFILE\" was specified ; genmake2 found and uses:" >> $SUMMARY
1390      echo "   so the genmake default will be used." >> $SUMMARY      #-note: to be filled later after 1rst run
1391  fi  fi
1392  echo  echo
1393  echo >> $SUMMARY  echo >> $SUMMARY
1394  if test "x$ADM" = x ; then  if test $KIND = 0 ; then
1395      cat << EOF | tee -a $SUMMARY      line_0=`printf '%s %2i' 'default' $MATCH_CRIT`
1396                  T           S           U           V         line_0="$line_0  ----T-----  ----S-----  ----U-----  ----V-----"
1397  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"
1398  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  ."
1399  N n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d      line_3="n n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d"
1400  2 d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .      line_4="2 d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  ."
1401        for ii in $PTRACERS_NUM ; do
1402  EOF          line_0="$line_0  --PTR 0"$ii"--"
1403            line_1="$line_1        m  s"
1404            line_2="$line_2  m  m  e  ."
1405            line_3="$line_3  i  a  a  d"
1406            line_4="$line_4  n  x  n  ."
1407        done
1408  else  else
1409      echo "ADJOINT=true" >> $SUMMARY      line_0=`printf '%s %2i' 'default   ' $MATCH_CRIT`
1410      if test $KIND = 1 ; then
1411       #echo "TANGLIN=true" >> $SUMMARY
1412        echo "TangLin generated by TAF" >> $SUMMARY
1413      elif test $KIND = 3 ; then
1414        echo "TangLin generated by OpenAD" >> $SUMMARY
1415      elif test $KIND = 2 ; then
1416       #echo "ADJOINT=true" >> $SUMMARY
1417        echo "Adjoint generated by TAF" >> $SUMMARY
1418      else
1419        echo "Adjoint generated by OpenAD" >> $SUMMARY
1420      fi
1421      if test $KIND = 1 -o $KIND = 3 ; then
1422        line_1="G D M    C  T  F"
1423        line_2="e p a R  o  L  D"
1424      else
1425        line_1="G D M    C  A  F"
1426        line_2="e p a R  o  d  D"
1427      fi
1428        line_3="n n k u  s  G  G"
1429        line_4="2 d e n  t  r  r"
1430      echo >> $SUMMARY      echo >> $SUMMARY
     cat << EOF | tee -a $SUMMARY  
 G D M    C  G  
 E p a R  o  r  
 N n k u  s  a  
 2 d e n  t  d  
   
 EOF  
1431  fi  fi
1432    if test "x$CLEANUP" != xt ; then
1433        echo "$line_0" | tee -a $SUMMARY
1434        echo "$line_1" | tee -a $SUMMARY
1435        echo "$line_2" | tee -a $SUMMARY
1436        echo "$line_3" | tee -a $SUMMARY
1437        echo "$line_4" | tee -a $SUMMARY
1438        echo ""        | tee -a $SUMMARY
1439    fi
1440    echo "-------------------------------------------------------------------------------"
1441    
1442  #  ...and each test directory...  #  ...and each test directory...
1443  for dir in $TESTDIRS ; do  for dir in $TESTDIRS ; do
1444        
1445        # set builddir & rundir:
1446        builddir="build"
1447        if test ! -d $dir/$builddir ; then mkdir $dir/$builddir ; fi
1448        rundir="run"
1449        pfxdir="tr_$rundir"
1450        if test ! -d $dir/$rundir ; then
1451            rundir=$builddir
1452        fi
1453        CODE_DIR=$dir/$code_dir
1454        BUILD_DIR=$dir/$builddir
1455    
1456      #  Cleanup only!      #  Cleanup only!
1457      if test "x$CLEANUP" = xt ; then      if test "x$CLEANUP" = xt ; then
1458          if test -r $dir/build/Makefile ; then          echo -n '  --- dir:' $BUILD_DIR ': '
1459              ( cd $dir/build ; make CLEAN )          makeclean $BUILD_DIR
1460            (   cd $BUILD_DIR
1461                rm -f $EXECUTABLE *.bak
1462                rm -f genmake_state genmake_*optfile genmake.log
1463                rm -f SIZE.h.mpi genmake.tr_log make.tr_log
1464            )
1465            if test -d $dir/$rundir/CVS ; then
1466                echo -n '  --- dir:' $dir/$rundir ': '
1467                run_clean $dir/$rundir
1468          fi          fi
1469          if test -r $dir/input/Makefile ; then          trdir=`( cd $dir ; find . -type d -name "$pfxdir.*" -print | sed 's/^.\///')`
1470              ( cd $dir/input ; make CLEAN )          ttd=`echo $trdir | wc -w`
1471            if test $ttd != 0 ; then
1472                echo '  --- rm dir:' $trdir
1473                ( cd $dir ; rm -rf $trdir )
1474          fi          fi
1475          continue          continue
1476      fi      fi
1477    
1478      #  Verify that the testdir exists and contains previous      #  Verify that the testdir exists and contains previous
1479      #  results in the correct location--or skip this directory!      #  results in the correct location--or skip this directory!
1480      fout=      fout=$dir"/results/"$ref_outp
1481      if test "x$ADM" = x ; then     #if test ! -r $fout ; then
1482          fout=$dir"/results/output.txt"      if test ! -r $fout -a ! -r ${fout}.gz ; then
     else  
         fout=$dir"/results_ad/output.txt_adm"  
     fi  
     if test ! -r $fout ; then  
1483          echo "can't read \"$fout\" -- skipping $dir"          echo "can't read \"$fout\" -- skipping $dir"
1484          continue          continue
1485      fi      fi
1486    
1487      builddir="input"      # Check for specific files for particular type of run
1488      rundir="input"  
1489      use_seperate_build=0      if test ! -r $CODE_DIR"/SIZE.h_mpi"  -a "x$MPI" != "x0" ; then
1490      if test -d $dir/build -a -r $dir/build ; then          echo "can't find \"$CODE_DIR/SIZE.h_mpi\" -- skipping $dir"
1491          builddir="build"          continue
         rundir="build"  
         use_seperate_build=1  
         linkdata $use_seperate_build $dir/$rundir  
1492      fi      fi
1493            if test ! -r $dir"/input/eedata.mth" -a "x$MULTI_THREAD" = "xt" ; then
1494      #  Check whether there are "extra runs" for this testdir          echo "can't find \"$dir/input/eedata.mth\" -- skipping $dir"
1495      extra_runs=          continue
     if test "x$ADM" = x -a "x$use_seperate_build" = x1 ; then  
         ex_run_dirs=`( cd $dir ; echo input.* )`  
         echo "ex_run_dirs='$ex_run_dirs'"  
         for exd in $ex_run_dirs ; do  
             name=`echo $exd | sed -e 's/input.//g'`  
             outf="$dir/results/output.txt.$name"  
             if test -f $outf -a -r $outf ; then  
                 extra_runs="$extra_runs $name"  
             fi  
         done  
1496      fi      fi
1497    
1498      if test "x$ADM" = x ; then      if test "x$MPI" != "x0" ; then
1499          code_dir=code          ntx=1 ; nty=1
1500          CODE_DIR=$dir/code          if test "x$MULTI_THREAD" = "xt" ; then
1501      else            ff=$dir"/input/eedata.mth"
1502          code_dir=code_ad            ntx=`grep "^ *nTx *=" $ff | tail -1 | sed 's/^ *nTx *= *//' | sed "s/, *$//"`
1503          CODE_DIR=$dir/code_ad            nty=`grep "^ *nTy *=" $ff | tail -1 | sed 's/^ *nTy *= *//' | sed "s/, *$//"`
1504              if test "x$ntx" = x ; then ntx=1 ; fi
1505              if test "x$nty" = x ; then nty=1 ; fi
1506            fi
1507            #- create new SIZE.h with no more than '$MPI' Procs
1508            mk_mpi_size $CODE_DIR"/SIZE.h_mpi" $BUILD_DIR"/tr_size.mpi" $MPI $ntx $nty
1509            LOC_NPROC=$?
1510            (   cd $BUILD_DIR
1511                if test -r SIZE.h.mpi ; then
1512                    cmp tr_size.mpi SIZE.h.mpi > /dev/null 2>&1 ; RETVAL=$?
1513                else RETVAL=1
1514                fi
1515                if test "x$RETVAL" = x0 ; then
1516                    rm -f tr_size.mpi
1517                else
1518                    rm -f SIZE.h.mpi ; mv tr_size.mpi SIZE.h.mpi
1519                fi
1520            )
1521            if test "x$MPI_MFILE" != x ; then
1522                #- create new MPI machine-file with the right number of Procs
1523                rm -f $LOC_MFILE
1524                cat $MPI_MFILE | sort | uniq | head -$LOC_NPROC > $LOC_MFILE
1525                nl=`wc -l $LOC_MFILE | awk '{print $1}'`
1526                if [ $nl -lt $LOC_NPROC ] ; then
1527                    rm -f $LOC_MFILE
1528                    cat $MPI_MFILE | head -$LOC_NPROC > $LOC_MFILE
1529                    #sed -n "1,$LOC_NPROC p" $MPI_MFILE > $LOC_MFILE
1530                fi
1531                if [ $verbose -gt 1 ]; then
1532                    nl=`wc -l $LOC_MFILE | awk '{print $1}'`
1533                    echo " new LOC_MFILE=$LOC_MFILE : $nl procs for LOC_NPROC=$LOC_NPROC"
1534                fi
1535            fi
1536            if test "x$MULTI_THREAD" = "xt" ; then
1537                retv=`check_eedata $dir"/input/eedata.mth" $BUILD_DIR"/SIZE.h.mpi"`
1538                if test $retv != 0 ; then
1539                    echo "input/eedata.mth tiling misfit -- skipping $dir"
1540                    continue
1541                fi
1542            fi
1543      fi      fi
     BUILD_DIR=$dir/$builddir  
1544    
1545      if test ! -r $CODE_DIR"/SIZE.h_mpi" -a "x$MPI" = "xt" ; then      #  Check whether there are "extra runs" for this testdir
1546          echo "can't find \"$CODE_DIR/SIZE.h_mpi\" -- skipping $dir"      extra_runs=
1547          continue      if test "x$NORUN" = xf ; then
1548            ex_run_dirs=`( cd $dir ; ls -d $inputdir.* 2> /dev/null )`
1549      fi      fi
1550        #echo "ex_run_dirs='$ex_run_dirs'"
1551        for exd in $ex_run_dirs ; do
1552            name=`echo $exd | sed -e "s/$inputdir\.//"`
1553            refExOut=`echo $ref_outp | sed "s/\./.${name}./"`
1554            outf="$dir/results/$refExOut"
1555            if test -r $outf -o -r ${outf}.gz ; then
1556              if test "x$MULTI_THREAD" = "xt" ; then
1557                if test -r $dir"/"$exd"/eedata.mth" ; then
1558                  if test "x$MPI" = "x0" ; then
1559                    extra_runs="$extra_runs $name"
1560                  else
1561                    retv=`check_eedata $dir"/"$exd"/eedata.mth" $BUILD_DIR"/SIZE.h.mpi"`
1562                    if test $retv = 0 ; then
1563                        extra_runs="$extra_runs $name"
1564                    else
1565                        echo $exd"/eedata.mth tiling misfit -- skipping $dir"
1566                    fi
1567                  fi
1568                #else echo $dir"/"$exd"/eedata.mth: not found"
1569                fi
1570              else
1571                extra_runs="$extra_runs $name"
1572              fi
1573            fi
1574        done
1575    
     echo "-------------------------------------------------------------------------------"  
1576      echo      echo
1577      echo "Experiment:  $dir"      if test "x$extra_runs" = "x" ; then
1578            echo "Experiment:  $dir"
1579        else
1580            echo "Experiment:  $dir ; extra_runs=$extra_runs"
1581        fi
1582      echo      echo
1583      unset genmake makedepend make run      unset genmake makedepend make run
1584      results='-- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --'      results=$EMPTY_RESULTS
1585    
1586      #  Create an output dir for each OPTFILE/tdir combination      #  Create an output dir & summary.txt file for each tested experiment (tdir)
1587      rel_CDIR=$DRESULTS"/"$dir      locDIR=$DRESULTS"/"$dir
1588      mkdir $rel_CDIR      mkdir $locDIR
1589      CDIR=`pwd`"/$rel_CDIR"      #- report to this experiment local summary file ---
1590            echo "DATE='$DATE' ; tdir='$dir'" > $locDIR"/summary.txt"
1591      if test "x$CLEANUP" = xt ; then      echo "MACH='$MACH'" >> $locDIR"/summary.txt"
1592          makeclean $dir/$builddir      echo "UNAMEA='$UNAMEA'" >> $locDIR"/summary.txt"
1593        CDIR=`pwd`"/$locDIR"
1594    
1595        if test "x$NORUN" = xt ; then
1596                run=N
1597            genmakemodel $dir/$builddir && genmake=Y \
1598                && makeclean $dir/$builddir \
1599                && symlink_mpifiles $dir $code_dir $builddir \
1600                && makedependmodel $dir/$builddir && makedepend=Y \
1601                && makemodel $dir/$builddir && make=Y
1602      else      else
1603          genmakemodel $dir/$builddir && genmake=Y \          genmakemodel $dir/$builddir && genmake=Y \
1604              && makeclean $dir/$builddir \              && makeclean $dir/$builddir \
1605              && symlink_mpifiles $dir $code_dir $builddir \              && symlink_mpifiles $dir $code_dir $builddir \
1606              && makedependmodel $dir/$builddir && makedepend=Y \              && makedependmodel $dir/$builddir && makedepend=Y \
1607              && makemodel $dir/$builddir && make=Y \              && makemodel $dir/$builddir && make=Y \
1608              && linkdata $use_seperate_build $dir/$rundir \              && run_clean $dir/$rundir \
1609                && linkdata $dir/$rundir $inputdir \
1610              && runmodel $dir/$rundir && run=Y \              && runmodel $dir/$rundir && run=Y \
1611              && results=`testoutput $dir $rundir`              && results=`testoutput_run $dir $rundir $ref_outp`
1612      fi      fi
1613            #echo "results='$results'"
1614      echo  
     if test "x$ADM" = x ; then  
1615          fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results`          fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results`
1616          echo          echo 1>&2
1617          echo "$fres" >> $SUMMARY          echo "$fres" | sed 's/ 99/ --/g' | sed 's/  > />/' | sed 's/  < /</' >> $SUMMARY
1618          echo "fresults='$fres'" > $CDIR"/summary.txt"          echo "fresults='$fres'" | sed 's/ 99/ --/g' >> $locDIR"/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"  
1619    
         OLD_COMMAND=$COMMAND  
         COMMAND="./mitgcmuv > output.txt"  
1620          for ex in $extra_runs ; do          for ex in $extra_runs ; do
1621              test -e "$dir/tr_run" && rm -rf "$dir/tr_run"              unset run
1622              mkdir "$dir/tr_run"              results=$EMPTY_RESULTS
1623              links=`( cd "$dir/input" > /dev/null 2>&1 ; ls -1 | grep -v CVS )`              #  reference output file
1624              (              refExOut=`echo $ref_outp | sed "s/\./.${ex}./g"`
1625                  cd "$dir/tr_run"              #  Create an output dir & summary.txt file for each extra run (tdir.ex)
1626                  for i in $links; do              locDIR=$DRESULTS"/"$dir"."$ex
1627                      ln -s ../input/$i $i              mkdir $locDIR
1628                  done              #- report to this experiment local summary file ---
1629              )              echo "DATE='$DATE' ; tdir='$dir.$ex'" > $locDIR"/summary.txt"
1630              links=`( cd "$dir/input.$ex" > /dev/null 2>&1 ; ls -1 | grep -v CVS )`              #echo "MACH='$MACH'" >> $locDIR"/summary.txt"
1631              (              #echo "UNAMEA='$UNAMEA'" >> $locDIR"/summary.txt"
1632                  cd "$dir/tr_run"              CDIR=`pwd`"/$locDIR"
1633                  for i in $links; do              test ! -e "$dir/$pfxdir.$ex" && mkdir "$dir/$pfxdir.$ex"
1634                      test -e $i  &&  rm -f $i              run_clean $dir/$pfxdir.$ex
1635                      ln -s ../input.$ex/$i $i              linkdata $dir/$pfxdir.$ex $inputdir.$ex $inputdir
1636                  done              runmodel $dir/$pfxdir.$ex && run=Y \
1637                  ln -s ../$builddir/mitgcmuv mitgcmuv              && results=`testoutput_run $dir $pfxdir.$ex $refExOut`
1638              )              fres=`formatresults $dir.$ex ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results`
1639              runmodel $dir/tr_run && run=Y \              echo 1>&2
1640                  && results=`testoutput $dir tr_run "."$ex`              echo "$fres" | sed 's/ 99/ --/g' | sed 's/  > />/' | sed 's/  < /</' >> $SUMMARY
1641              fres=`printf '%s %s %s %s' ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N}`              echo "fresults='$fres'" | sed 's/ 99/ --/g' >> $locDIR"/summary.txt"
1642              fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results`              if test "x$POSTCLEAN" = x2 ; then
1643              fres="$fres.$ex"                  run_clean $dir/$pfxdir.$ex
1644              echo              fi
             echo "$fres" >> $SUMMARY  
             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"  
1645          done          done
1646          COMMAND=$OLD_COMMAND  
1647      else      if test ! -f $DRESULTS"/"genmake_state ; then
1648          fres=`printf '%s %s %s %s' ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N}`          if test -f $dir/$builddir/Makefile ; then
1649          fres=$fres"$results   $dir"              mkOpt=`grep '^# OPTFILE=' $dir/$builddir/Makefile 2>/dev/null | head -1 | sed 's/^# //'`
1650          echo              echo "from '$dir/$builddir/Makefile', extract:" > $DRESULTS/genmake_state
1651          echo "$fres" >> $SUMMARY              sed -n '/^# executed by:/,+1 p' $dir/$builddir/Makefile >> $DRESULTS/genmake_state
1652          echo "fresults='$fres'" > $CDIR"/summary.txt"              echo " $mkOpt" >> $DRESULTS/genmake_state
1653          echo "MACH='$MACH'" >> $CDIR"/summary.txt"              if test "x$OPTFILE" = xNONE ; then
1654          echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt"                  eval $mkOpt
1655          echo "DATE='$DATE'" >> $CDIR"/summary.txt"                  sed "/^No \"OPTFILE\" was specified ; genmake2/a\  OPTFILE=${OPTFILE}"\
1656          echo "tdir='$dir'" >> $CDIR"/summary.txt"                          $SUMMARY > tmp.tr_log
1657                    RETVAL=$?
1658                    if test "x$RETVAL" = x0 ; then rm -f $SUMMARY
1659                            cp tmp.tr_log $SUMMARY
1660                    else rm -f tmp.tr_log
1661                    fi
1662                fi
1663                gmkLog=$dir/$builddir/genmake.log
1664                if test -r $gmkLog ; then
1665                    grep '^Get compiler version using:' $gmkLog > /dev/null 2>&1
1666                    RETVAL=$?
1667                    if test "x$RETVAL" = x0 ; then
1668                      echo -n "from '$gmkLog', " >> $DRESULTS/genmake_state
1669                      echo "extract compiler version:"   >> $DRESULTS/genmake_state
1670                      sed -n '/Get compiler version/,/<-- compiler version/p' \
1671                         $gmkLog | grep -v '^... compiler version ' > tmp.tr_log
1672                      sed -n '1p' tmp.tr_log >> $DRESULTS/genmake_state
1673                      sed -n '2,/^$/p' tmp.tr_log | sed '/^$/d' | sed 's/^./ &/' \
1674                                             >> $DRESULTS/genmake_state
1675                      rm -f tmp.tr_log
1676                    fi
1677                    echo -n "from '$gmkLog', " >> $DRESULTS/genmake_state
1678                    echo "get NETCDF & LAPACK settings:" >> $DRESULTS/genmake_state
1679                    sed -n '/set HAVE_NETCDF=/p' $gmkLog >> $DRESULTS/genmake_state
1680                    sed -n '/set HAVE_LAPACK=/p' $gmkLog >> $DRESULTS/genmake_state
1681                fi
1682            fi
1683      fi      fi
1684            #postclean $dir/$builddir
1685        if test "x$POSTCLEAN" = x2 ; then
1686            makeclean $dir/$builddir \
1687                && run_clean $dir/$rundir
1688        fi
1689        if test "x$MPI" != x0 -a "x$MPI_MFILE" != x ; then rm -f $LOC_MFILE ; fi
1690    
1691      echo "-------------------------------------------------------------------------------"      echo "-------------------------------------------------------------------------------"
1692        
1693  done  done
1694    
1695  printf "Start time:  " >> $SUMMARY  printf "Start time:  " >> $SUMMARY
1696  echo $start_date >> $SUMMARY  echo "$start_date" >> $SUMMARY
1697  printf "End time:    " >> $SUMMARY  printf "End time:    " >> $SUMMARY
1698  date >> $SUMMARY  date >> $SUMMARY
1699    
# Line 925  else Line 1705  else
1705      if test "x$HAVE_MPACK" = xt ; then      if test "x$HAVE_MPACK" = xt ; then
1706          tar -cf $DRESULTS".tar" $DRESULTS > /dev/null 2>&1 \          tar -cf $DRESULTS".tar" $DRESULTS > /dev/null 2>&1 \
1707              && gzip $DRESULTS".tar" \              && gzip $DRESULTS".tar" \
1708              && $MPACK -s MITgcm-test -m 1500000 $DRESULTS".tar.gz" $ADDRESSES              && $MPACK -s MITgcm-test -m 3555000 $DRESULTS".tar.gz" $ADDRESSES
1709          RETVAL=$?          RETVAL=$?
1710          if test "x$RETVAL" != x0 ; then          if test "x$RETVAL" != x0 ; then
1711              echo              echo
# Line 938  else Line 1718  else
1718              echo "An email containing results was sent to the following addresses:"              echo "An email containing results was sent to the following addresses:"
1719              echo "  \"$ADDRESSES\""              echo "  \"$ADDRESSES\""
1720              echo              echo
1721                test -f $DRESULTS".tar"  &&  rm -f $DRESULTS".tar"
1722                test -f $DRESULTS".tar.gz"  &&  rm -f $DRESULTS".tar.gz"
1723          fi          fi
         test -f $DRESULTS".tar"  &&  rm -f $DRESULTS".tar"  
         test -f $DRESULTS".tar.gz"  &&  rm -f $DRESULTS".tar.gz"  
1724      fi      fi
1725  fi  fi
1726    
1727  # rm -f tmp_cmpnum.f a.out  if test "x$QUICK" = xf -a "x$NORUN" = xf ; then
1728  rm -f tmp_cmpnum.c tmp_cmpnum      rm -f tr_cmpnum.c tr_cmpnum
1729    fi
1730    
1731  if test "x$CLEANUP" != xt ; then  if test "x$CLEANUP" != xt ; then
1732      cat $SUMMARY      cat $SUMMARY | sed 's/ \.  \.  \.  \.  \.  \.  \.  \.  \.  \.  \.  \. //'
1733      if test -e tr_out.txt ; then      if test -e tr_out.txt ; then
1734          mv tr_out.txt tr_out.txt.old          mv tr_out.txt tr_out.txt.old
1735      fi      fi
1736      cat $SUMMARY > tr_out.txt      cat $SUMMARY | sed '/^[YN] [YN] [YN] [YN]/ s/ \. //g' > tr_out.txt
1737  fi  fi
1738    
1739    if test "x$DELDIR" = xt ; then
1740        rm -rf $DRESULTS
1741    fi
1742    echo "======== End of testreport execution ========"

Legend:
Removed from v.1.38  
changed lines
  Added in v.1.188

  ViewVC Help
Powered by ViewVC 1.1.22