/[MITgcm]/MITgcm/verification/testscript
ViewVC logotype

Contents of /MITgcm/verification/testscript

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


Revision 1.29 - (show annotations) (download)
Wed Jun 9 14:03:35 2004 UTC (19 years, 9 months ago) by adcroft
Branch: MAIN
CVS Tags: checkpoint57t_post, checkpoint58l_post, checkpoint53f_post, checkpoint54a_pre, checkpoint55c_post, checkpoint57m_post, checkpoint58e_post, checkpoint57v_post, checkpoint57f_post, checkpoint57s_post, checkpoint57j_post, checkpoint58b_post, checkpoint58m_post, checkpoint57b_post, checkpoint53d_post, checkpoint57f_pre, checkpoint57k_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint57a_post, checkpoint58r_post, checkpoint55j_post, checkpoint56b_post, checkpoint59c, checkpoint57y_post, checkpoint58g_post, checkpoint57x_post, checkpoint54a_post, checkpoint55h_post, checkpoint58n_post, checkpoint58x_post, checkpoint57g_pre, checkpoint54b_post, checkpoint58h_post, checkpoint57e_post, checkpoint58w_post, checkpoint54d_post, checkpoint56c_post, checkpoint54e_post, checkpoint58j_post, checkpoint55b_post, checkpoint57h_post, checkpoint57y_pre, checkpoint55, checkpoint55a_post, checkpoint57c_pre, checkpoint57o_post, checkpoint55g_post, checkpoint57r_post, checkpoint59e, checkpoint59d, checkpoint59f, checkpoint59a, checkpoint55f_post, checkpoint57i_post, checkpoint59b, checkpoint57h_pre, checkpoint59, checkpoint58, checkpoint57a_pre, checkpoint54, checkpoint57, checkpoint56, checkpoint57h_done, checkpoint58f_post, checkpoint53g_post, checkpoint57n_post, checkpoint58d_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint58a_post, checkpoint58i_post, checkpoint57q_post, checkpoint58o_post, checkpoint57z_post, checkpoint54f_post, eckpoint57e_pre, checkpoint58c_post, checkpoint58k_post, checkpoint57c_post, checkpoint58u_post, checkpoint58y_post, checkpoint58v_post, checkpoint58s_post, checkpoint55e_post, checkpoint54c_post, checkpoint58p_post, checkpoint58t_post, checkpoint58q_post, checkpoint55i_post, checkpoint57l_post, checkpoint56a_post, checkpoint55d_post
Changes since 1.28: +1 -0 lines
Added vertical diffusivity profile (T/S) due to Bryan and Lewis, 1979.
New parameters:
 diffKrBL79surf - surface diffusivity
 diffKrBL79deep - deep diffusivity
 diffKrBL79Ho   - turning depth for arctan function
 diffKrBL79scl  - depth scale for arctan function
This diffusivity is added to all other diffusivities. The defaults are
set so as to give zero diffusivity.

1 #!/bin/sh
2
3 compare_lines()
4 {
5 # use codelet to compare lines
6 if [ $verbose -gt 1 ]; then
7 cat tmp3.txt 1>&2
8 fi
9 return `./a.out < tmp3.txt`
10 }
11
12 testoutput_for_prop()
13 {
14 # testoutput_for_prop dir s1 label subdir
15 #
16 # compares files in $dir/$subdir/output.txt and $dir/results/output.txt
17 # using search strings s1 and text label
18
19 if [ $debug -gt 0 ]; then
20 echo testoutput_for_prop: grep "$2" $1/$4/output.txt 1>&2
21 fi
22 if [ -r $1/$4/output.txt ]; then
23 grep "$2" $1/$4/output.txt | sed 's/.*=//' | nl > tmp1.txt
24 lncnt=`wc -l tmp1.txt | awk '{print $1}' `
25 if [ $lncnt -lt 3 ]; then
26 if [ $verbose -gt 0 ]; then
27 echo Not enough lines of output when searching for "$2" 1>&2
28 fi
29 return 99
30 fi
31 else
32 echo testoutput_for_prop: output.txt from model run was not readable 1>&2
33 return 99
34 fi
35 if [ $debug -gt 0 ]; then
36 echo testoutput_for_prop: grep "$2" $1/results/output.txt 1>&2
37 fi
38 grep "$2" $1/results/output.txt | sed 's/.*=//' | nl > tmp2.txt
39 lncnt=`wc -l tmp2.txt | awk '{print $1}' `
40 if [ $lncnt -lt 3 ]; then
41 if [ $verbose -gt 0 ]; then
42 echo Not enough lines of output when searching for "$2" 1>&2
43 fi
44 return 99
45 fi
46 if [ $debug -gt 0 ]; then
47 echo testoutput_for_prop: join tmp1.txt tmp2.txt 1>&2
48 fi
49 join tmp1.txt tmp2.txt | awk '{print $1 " " $2 " " $3}' > tmp3.txt
50 if [ $debug -gt 0 ]; then
51 echo testoutput_for_prop: compare_lines 1>&2
52 fi
53 compare_lines
54 digits_of_similarity=$?
55 if [ $digits_of_similarity -eq 99 ]; then
56 if [ $verbose -gt 0 ]; then
57 echo testoutput_for_prop: No comparison was available for \"$2\" 1>&2
58 fi
59 digits_of_similarity=99
60 else
61 if [ $verbose -gt 0 ]; then
62 echo There were $digits_of_similarity decimal places of similarity for \"$2\" 1>&2
63 fi
64 fi
65 rm tmp1.txt tmp2.txt tmp3.txt
66
67 return $digits_of_similarity
68 }
69
70 dashnum()
71 {
72 # dashnum n1 n2 n3 ...
73 #
74 # print numbers using %3i format or "--" if number = 99
75 for num in $@
76 do
77 if [ $num = 99 ]; then
78 printf ' --'
79 else
80 printf '%3i' $num
81 fi
82 done
83 }
84
85 testoutput()
86 {
87 # testoutput diretory subdir
88 #
89 # test output in "directory"
90
91 if [ $debug -gt 0 ]; then
92 echo testoutput: testoutput_for_prop $1 cg2d_init_res 1>&2
93 fi
94 testoutput_for_prop $1 "cg2d_init_res" "cg2d init. residual" $2; cg2dres=$?
95 if [ $debug -gt 0 ]; then
96 echo testoutput: cg2dres=$cg2dres 1>&2
97 fi
98
99 if [ $longtest -gt 0 ]; then
100 testoutput_for_prop $1 "dynstat_theta_min" "theta minimum" $2; tmin=$?
101 testoutput_for_prop $1 "dynstat_theta_max" "theta maximum" $2; tmax=$?
102 testoutput_for_prop $1 "dynstat_theta_mean" "theta mean" $2; tmean=$?
103 testoutput_for_prop $1 "dynstat_theta_sd" "theta s.d." $2; tsd=$?
104 testoutput_for_prop $1 "dynstat_salt_min" "salt minimum" $2; smin=$?
105 testoutput_for_prop $1 "dynstat_salt_max" "salt maximum" $2; smax=$?
106 testoutput_for_prop $1 "dynstat_salt_mean" "salt mean" $2; smean=$?
107 testoutput_for_prop $1 "dynstat_salt_sd" "salt s.d." $2; ssd=$?
108 testoutput_for_prop $1 "dynstat_uvel_min" "U minimum" $2; umin=$?
109 testoutput_for_prop $1 "dynstat_uvel_max" "U maximum" $2; umax=$?
110 testoutput_for_prop $1 "dynstat_uvel_mean" "U mean" $2; umean=$?
111 testoutput_for_prop $1 "dynstat_uvel_sd" "U s.d." $2; usd=$?
112 testoutput_for_prop $1 "dynstat_vvel_min" "V minimum" $2; vmin=$?
113 testoutput_for_prop $1 "dynstat_vvel_max" "V maximum" $2; vmax=$?
114 testoutput_for_prop $1 "dynstat_vvel_mean" "V mean" $2; vmean=$?
115 testoutput_for_prop $1 "dynstat_vvel_sd" "V s.d." $2; vsd=$?
116 else
117 testoutput_for_prop $1 "dynstat_theta_min" "theta minimum" $2; tmin=$?
118 testoutput_for_prop $1 "dynstat_theta_max" "theta maximum" $2; tmax=$?
119 testoutput_for_prop $1 "dynstat_salt_min" "salt minimum" $2; smin=$?
120 testoutput_for_prop $1 "dynstat_salt_max" "salt maximum" $2; smax=$?
121 testoutput_for_prop $1 "dynstat_uvel_min" "U minimum" $2; umin=$?
122 testoutput_for_prop $1 "dynstat_uvel_max" "U maximum" $2; umax=$?
123 testoutput_for_prop $1 "dynstat_vvel_min" "V minimum" $2; vmin=$?
124 testoutput_for_prop $1 "dynstat_vvel_max" "V maximum" $2; vmax=$?
125 fi
126
127 dashnum $cg2dres $tmin $tmax $tmean $tsd $smin $smax $smean $ssd \
128 $umin $umax $umean $usd $vmin $vmax $vmean $vsd
129 #printf '%3i' $cg2dres $tmin $tmax $tmean $tsd $smin $smax $smean $ssd \
130 # $umin $umax $umean $usd $vmin $vmax $vmean $vsd
131 }
132
133 genmakemodel()
134 {
135 # genmakemodel directory
136 ( cd $1;
137 if [ $quick -eq 0 -o ! -r Makefile ]; then
138 printf 'genmake ... ' 1>&2
139 find . -type l -exec rm {} \;
140 if [ $ieee -eq 0 ]; then
141 ../../../tools/genmake2 -mods=../code > make.log 2>&1
142 else
143 ../../../tools/genmake2 -ieee -mods=../code > make.log 2>&1
144 fi
145 if [ $? -ne 0 ]; then
146 tail make.log
147 echo genmakemodel: $genmake failed 1>&2
148 return 1
149 else
150 echo succesful 1>&2
151 fi
152 fi
153 )
154 }
155
156 makecleancompile()
157 {
158 # makecleancompile directory
159 ( cd $1;
160 if [ $force -gt 0 ]; then
161 rm -f output.txt
162 printf 'make clean ... ' 2>&1
163 make CLEAN >> make.log 2>&1
164 if [ $? -ne 0 ]; then
165 tail make.log
166 echo makecleancompile: make clean failed 1>&2
167 return 1
168 else
169 echo succesful 1>&2
170 fi
171 fi
172 )
173 }
174
175 makecleanupafter()
176 {
177 # makeupafter directory
178 ( cd $1;
179 if [ $clean -gt 0 ]; then
180 rm -f output.txt
181 printf 'make clean ... ' 2>&1
182 make CLEAN >> make.log 2>&1
183 if [ $? -ne 0 ]; then
184 tail make.log
185 echo makeupafter: make clean failed 1>&2
186 return 1
187 else
188 echo succesful 1>&2
189 fi
190 fi
191 )
192 }
193
194 makedependmodel()
195 {
196 # makedependmodel directory
197 ( cd $1;
198 if [ $quick -eq 0 -o ! -r Makefile ]; then
199 printf 'make depend ... ' 1>&2
200 make cleanlinks >> make.log 2>&1
201 make depend >> make.log 2>&1
202 if [ $? -ne 0 ]; then
203 tail make.log
204 echo makemodel: make depend failed 1>&2
205 return 1
206 else
207 echo succesful 1>&2
208 fi
209 fi
210 )
211 }
212
213 makemodel()
214 {
215 # makemodel directory
216 ( cd $1;
217 if [ -r Makefile ]; then
218 printf 'make ... ' 1>&2
219 make >> make.log 2>&1
220 if [ $? -ne 0 ]; then
221 tail make.log
222 echo failed 1>&2
223 return 1
224 else
225 echo succesful 1>&2
226 fi
227 fi
228 )
229 }
230
231 linkdata()
232 {
233 # linkdata flag
234 #
235 # symbolically link data files to run directory
236 if test "x$1" = x1 ; then
237 (
238 cd $2
239 if test "x$ADM" = x ; then
240 files=`( cd ../input ; ls -1 | grep -v CVS )`
241 for i in $files ; do
242 if test ! -d "../input/"$i ; then
243 ln -sf "../input/"$i $i
244 fi
245 done
246 else
247 files=`( cd ../input ; ls -1 *.bin | grep -v CVS )`
248 for i in $files ; do
249 if test ! -d "../input/"$i ; then
250 ln -sf "../input/"$i $i
251 fi
252 done
253 files=`( cd ../input_ad ; ls -1 | grep -v CVS )`
254 for i in $files ; do
255 if test ! -d "../input_ad/"$i ; then
256 ln -sf "../input_ad/"$i $i
257 fi
258 done
259 fi
260 )
261 fi
262 }
263
264 runmodel()
265 {
266 # runmodel directory exe
267 #
268 # runs the model "exe" in "directory" (exe is relative to directory)
269
270 ( cd $1
271 if [ -x $2 ]; then
272 if [ $quick -eq 0 ]; then
273 rm -f output.txt
274 fi
275 if [ $debug -gt 0 ]; then
276 echo -n runmodel: pwd= 1>&2
277 pwd 1>&2
278 fi
279 printf 'runmodel: ' 1>&2
280 rm -f *.nc
281 make output.txt && return 0
282 return 1
283 fi
284 )
285 }
286
287 copynewresults()
288 {
289 # copynewresults directory
290 #
291 # runs the model "exe" in "directory" (exe is relative to directory)
292
293 ( cd $1
294 if [ $newresults -eq 1 ]; then
295 if [ -r output.txt ]; then
296 cp -f output.txt ../results/
297 return 0
298 fi
299 fi
300 )
301 }
302
303 createcodelet()
304 {
305 # create codelet for comparing model output
306 cat > tmp_cmpnum.f <<EOFA
307 program cmpnum
308 implicit none
309 real*8 a,b,diff
310 integer linnum,best
311 best=-16
312 99 read(*,*,end=70,err=60) linnum,a,b
313 diff=0.5*(abs(a)+abs(b))
314 c print *,a,b,diff,abs(a-b)/diff
315 if (diff.gt.1.e-12) then
316 diff=abs(a-b)/diff
317 if (diff.gt.0.) then
318 c print *,int(log10(diff)),diff
319 linnum=int(log10(diff))
320 best=max(best,linnum)
321 endif
322 else
323 if (best.eq.-16.and.diff.ne.0.) best=-22
324 endif
325 goto 99
326 60 stop 'cmpnum: An error occured reading a,b'
327 70 print *,-best
328 end
329 EOFA
330 f77 tmp_cmpnum.f
331 if [ -x ./a.out ]; then
332 return 0
333 else
334 echo createcodelet: failed to compile codelet 1>&2
335 exit 1
336 fi
337 }
338
339 formatresults()
340 {
341 # formatresults expt genmake depend make run results*
342
343 nm=$1
344 printf '%s %s %s %s' $2 $3 $4 $5
345 shift; shift; shift; shift; shift;
346 printf '%3s' $@
347
348 if [ $1 = '--' ]; then
349 printf ' N/O '
350 else
351 if [ $1 -gt 12 ]; then
352 printf ' pass'
353 else
354 printf ' FAIL'
355 fi
356 fi
357 printf ' %s' $nm
358 printf '\n'
359
360 }
361
362 show_help()
363 {
364 cat - << EOF
365 $0 [-help] [-quick] [-verbose] dir1 [dir2] [...]
366 -help Show this help message
367 -quick Skip "genmake" and "make depend" if the Makefile exists
368 -newresults Copy the output.txt into results/
369 -quiet Reduce the amount of output
370 -verbose Produce copious amounts of output
371 -debug Produce even more output which will mean nothing to most
372 -force Do "make CLEAN" before compiling. This forces a complete rebuild.
373 -clean Do "make CLEAN" after compiling and testing.
374 -shorttest Don't compare numeric output for mean and s.d. of variables.
375 -noieee By default, $0 uses the -ieee option for genmake. This turns it off.
376 -cleanup Aggresively removes all model output, executables and object files
377 and then exits. Use with care.
378
379 Normal usage:
380 $0 * Configure, compile, run and analyze in all experiment directories
381 EOF
382 }
383
384 scandirs()
385 {
386 if [ $# -eq 0 ]; then
387 for arg in *
388 do
389 test -d $arg/input && echo $arg
390 done
391 else
392 echo $*
393 fi
394 }
395
396 clean_up()
397 {
398 # Find all executables, object files, CPP'd source and model output
399 # and DELETE it.
400 for opt in '-name "mitgcmuv*"' '-name "*.o"' '-name "*.f"' '-name "*.c"' '-name "fort.*"' '-name "make.log"' '-path *results -prune -o -name "*.meta"' '-path *results -prune -o -name "*.data"' '-type l'
401 do
402 echo Cleaning: find . $opt -exec rm {} \;
403 find . $opt -exec rm {} \;
404 done
405 }
406
407 ###############################################################################
408 ###############################################################################
409 ###############################################################################
410
411 # Main function
412
413 # Default properties
414 debug=0
415 verbose=1
416 quick=0
417 newresults=0
418 force=0
419 clean=0
420 ieee=1
421 longtest=1
422 expts=''
423 genmake=genmake
424
425 # Process arguments
426 for arg in $@
427 do
428 case $arg in
429 -cleanup) clean_up; exit 0;;
430 -quick) quick=1;;
431 -newresults) newresults=1;;
432 -verbose) verbose=2;;
433 -debug) debug=1;;
434 -force) force=1;;
435 -clean) clean=1;;
436 -noieee) ieee=0;;
437 -shorttest) longtest=0;;
438 -quiet) verbose=0;;
439 -help) show_help; exit 0;;
440 -*) echo Unrecognized option:$arg; exit 9;;
441 *) test -d $arg && expts=`echo $expts $arg`;;
442 esac
443 done
444
445 if [ $force -gt 0 -a $quick -gt 0 ]; then
446 echo You specified -quick and -force together which conflict.
447 echo Please specify either -quick or -force or neither but not both.
448 exit 1
449 fi
450
451 #if [ ${#expts} -eq 0 ]; then
452 # echo Scanning all directories
453 # for arg in *
454 # do
455 # test -d $arg && expts=`echo $expts $arg`
456 # done
457 #fi
458 expts=`scandirs $expts`
459
460 createcodelet
461
462 date > summary.txt
463 if [ $longtest -gt 0 ]; then
464 cat << EOF >> summary.txt
465 T S U V
466 C D M c m s m s m s m s
467 n p a R g m m e . m m e . m m e . m m e .
468 f n k u 2 i a a d i a a d i a a d i a a d
469 g d e n d n x n . n x n . n x n . n x n .
470
471 EOF
472 else
473 cat << EOF >> summary.txt
474 T S U V
475 C D M c
476 n p a R g m m m m m m m m
477 f n k u 2 i a i a i a i a
478 g d e n d n x n x n x n x
479
480 EOF
481 fi
482
483 # Now configue, make, run and test in each directory
484 for dir in $expts
485 do
486 if [ -r $dir/results/output.txt ]; then
487 echo -------------------------------------------------------------------------------
488 echo
489 echo Experiment: $dir
490 echo
491 unset genmakepass makedepend make run
492 if [ $longtest -gt 0 ]; then
493 results='-- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --'
494 else
495 results='-- -- -- -- -- -- -- -- --'
496 fi
497 if [ -r $dir/build ]; then
498 seperatebuilddir=1
499 builddir=build
500 rundir=build
501 (cd $dir/input; rm -f *.{o,f,c,F,h} work* output.txt Make* make.log; )
502 else
503 seperatebuilddir=0
504 builddir=input
505 rundir=input
506 fi
507 genmakemodel $dir/$builddir && genmakepass=Y \
508 && makecleancompile $dir/$builddir \
509 && makedependmodel $dir/$builddir && makedepend=Y \
510 && makemodel $dir/$builddir && make=Y \
511 && linkdata $seperatebuilddir $dir/$rundir \
512 && runmodel $dir/$builddir mitgcmuv && run=Y \
513 && results=`testoutput $dir $rundir` \
514 && copynewresults $dir/$rundir \
515 && makecleanupafter $dir/$builddir
516 echo
517 formatresults $dir ${genmakepass:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results
518 echo
519 formatresults $dir ${genmakepass:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results >> summary.txt
520 fi
521 done
522
523 rm tmp_cmpnum.f a.out
524
525 echo -------------------------------------------------------------------------------
526 echo
527 date >> summary.txt
528 cat summary.txt

  ViewVC Help
Powered by ViewVC 1.1.22