96 |
|
|
97 |
testoutput_var() |
testoutput_var() |
98 |
{ |
{ |
99 |
# testoutput_var dir s1 label subdir extension |
# testoutput_var dir s1 label subdir reference_output |
100 |
# |
# |
101 |
# compares 1 variable output selected from file $dir/$subdir/$OUTPUTFILE |
# compares 1 variable output selected from file $dir/$subdir/$OUTPUTFILE |
102 |
# with same output from reference file $dir/results/output.$extension |
# with same output from reference file $dir/results/$reference_output |
103 |
# using search strings s1 and text label |
# using search strings s1 and text label |
104 |
|
|
105 |
if [ $debug -gt 0 ]; then |
if [ $debug -gt 0 ]; then |
119 |
return 99 |
return 99 |
120 |
fi |
fi |
121 |
if [ $debug -gt 0 ]; then |
if [ $debug -gt 0 ]; then |
122 |
echo testoutput_var: grep "$2" $1/results/output.$5 1>&2 |
echo testoutput_var: grep "$2" $1/results/$5 1>&2 |
123 |
fi |
fi |
124 |
grep "$2" $1/results/output.$5 | sed 's/.*=//' | cat -n > tmp2.txt |
grep "$2" $1/results/$5 | sed 's/.*=//' | cat -n > tmp2.txt |
125 |
lncntB=`wc -l tmp2.txt | awk '{print $1}' ` |
lncntB=`wc -l tmp2.txt | awk '{print $1}' ` |
126 |
if [ $lncntB -lt 3 ]; then |
if [ $lncntB -lt 3 ]; then |
127 |
if [ $verbose -gt 0 ]; then |
if [ $verbose -gt 0 ]; then |
161 |
digits_of_similarity=`./tr_cmpnum < tmp4.txt` |
digits_of_similarity=`./tr_cmpnum < tmp4.txt` |
162 |
if [ $digits_of_similarity -eq 99 ]; then |
if [ $digits_of_similarity -eq 99 ]; then |
163 |
if [ $verbose -gt 0 ]; then |
if [ $verbose -gt 0 ]; then |
164 |
echo testoutput_var: No comparison was available for \"$2\" 1>&2 |
echo testoutput_var: No comparison was available for \"$3\" 1>&2 |
165 |
fi |
fi |
166 |
digits_of_similarity=99 |
digits_of_similarity=99 |
167 |
else |
else |
168 |
if [ $verbose -gt 0 ]; then |
if [ $verbose -gt 0 ]; then |
169 |
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 |
170 |
fi |
fi |
171 |
fi |
fi |
172 |
rm -f tmp1.txt tmp2.txt tmp3.txt tmp4.txt |
rm -f tmp1.txt tmp2.txt tmp3.txt tmp4.txt |
232 |
|
|
233 |
testoutput_run() |
testoutput_run() |
234 |
{ |
{ |
235 |
# testoutput_run directory subdir extension |
# testoutput_run directory subdir reference_output |
236 |
# |
# |
237 |
# test output from 1 run in "directory" |
# test output from 1 run in "directory" |
238 |
if test "x$ADM" = x ; then |
if test "x$ADM" = x ; then |
268 |
case $xx in |
case $xx in |
269 |
'PS') if [ $debug -gt 0 ] |
'PS') if [ $debug -gt 0 ] |
270 |
then echo testoutput_run: testoutput_var $1 cg2d_init_res 1>&2 ; fi |
then echo testoutput_run: testoutput_var $1 cg2d_init_res 1>&2 ; fi |
271 |
testoutput_var $1 "cg2d_init_res" "cg2d init. residual" $2 $3; yy=$? |
testoutput_var $1 "cg2d_init_res" "Press. Solver (cg2d)" $2 $3; yy=$? |
272 |
if [ $debug -gt 0 ] ; then echo testoutput_run: cg2dres=$yy 1>&2 ; fi ;; |
if [ $debug -gt 0 ] ; then echo testoutput_run: cg2dres=$yy 1>&2 ; fi ;; |
273 |
'Tmn') testoutput_var $1 "dynstat_theta_min" "theta minimum" $2 $3; yy=$? ;; |
'Cost') testoutput_var $1 "ADM precision_derivative_cost" "ADM Cost" $2 $3; yy=$? ;; |
274 |
'Tmx') testoutput_var $1 "dynstat_theta_max" "theta maximum" $2 $3; yy=$? ;; |
'Grad') testoutput_var $1 "ADM precision_derivative_grad" "ADM Grad" $2 $3; yy=$? ;; |
275 |
'Tav') testoutput_var $1 "dynstat_theta_mean" "theta mean" $2 $3; yy=$? ;; |
'Tmn') testoutput_var $1 "dynstat_theta_min" "Theta minimum" $2 $3; yy=$? ;; |
276 |
'Tsd') testoutput_var $1 "dynstat_theta_sd" "theta s.d." $2 $3; yy=$? ;; |
'Tmx') testoutput_var $1 "dynstat_theta_max" "Theta maximum" $2 $3; yy=$? ;; |
277 |
'Smn') testoutput_var $1 "dynstat_salt_min" "salt minimum" $2 $3; yy=$? ;; |
'Tav') testoutput_var $1 "dynstat_theta_mean" "Theta mean" $2 $3; yy=$? ;; |
278 |
'Smx') testoutput_var $1 "dynstat_salt_max" "salt maximum" $2 $3; yy=$? ;; |
'Tsd') testoutput_var $1 "dynstat_theta_sd" "Theta Std.Dev" $2 $3; yy=$? ;; |
279 |
'Sav') testoutput_var $1 "dynstat_salt_mean" "salt mean" $2 $3; yy=$? ;; |
'Smn') testoutput_var $1 "dynstat_salt_min" "Salt minimum" $2 $3; yy=$? ;; |
280 |
'Ssd') testoutput_var $1 "dynstat_salt_sd" "salt s.d." $2 $3; yy=$? ;; |
'Smx') testoutput_var $1 "dynstat_salt_max" "Salt maximum" $2 $3; yy=$? ;; |
281 |
|
'Sav') testoutput_var $1 "dynstat_salt_mean" "Salt mean" $2 $3; yy=$? ;; |
282 |
|
'Ssd') testoutput_var $1 "dynstat_salt_sd" "Salt Std.Dev" $2 $3; yy=$? ;; |
283 |
'Umn') testoutput_var $1 "dynstat_uvel_min" "U minimum" $2 $3; yy=$? ;; |
'Umn') testoutput_var $1 "dynstat_uvel_min" "U minimum" $2 $3; yy=$? ;; |
284 |
'Umx') testoutput_var $1 "dynstat_uvel_max" "U maximum" $2 $3; yy=$? ;; |
'Umx') testoutput_var $1 "dynstat_uvel_max" "U maximum" $2 $3; yy=$? ;; |
285 |
'Uav') testoutput_var $1 "dynstat_uvel_mean" "U mean" $2 $3; yy=$? ;; |
'Uav') testoutput_var $1 "dynstat_uvel_mean" "U mean" $2 $3; yy=$? ;; |
286 |
'Usd') testoutput_var $1 "dynstat_uvel_sd" "U s.d." $2 $3; yy=$? ;; |
'Usd') testoutput_var $1 "dynstat_uvel_sd" "U Std.Dev" $2 $3; yy=$? ;; |
287 |
'Vmn') testoutput_var $1 "dynstat_vvel_min" "V minimum" $2 $3; yy=$? ;; |
'Vmn') testoutput_var $1 "dynstat_vvel_min" "V minimum" $2 $3; yy=$? ;; |
288 |
'Vmx') testoutput_var $1 "dynstat_vvel_max" "V maximum" $2 $3; yy=$? ;; |
'Vmx') testoutput_var $1 "dynstat_vvel_max" "V maximum" $2 $3; yy=$? ;; |
289 |
'Vav') testoutput_var $1 "dynstat_vvel_mean" "V mean" $2 $3; yy=$? ;; |
'Vav') testoutput_var $1 "dynstat_vvel_mean" "V mean" $2 $3; yy=$? ;; |
290 |
'Vsd') testoutput_var $1 "dynstat_vvel_sd" "V s.d." $2 $3; yy=$? ;; |
'Vsd') testoutput_var $1 "dynstat_vvel_sd" "V Std.Dev" $2 $3; yy=$? ;; |
291 |
'pt1mn'|'pt2mn'|'pt3mn'|'pt4mn'|'pt5mn') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
'pt1mn'|'pt2mn'|'pt3mn'|'pt4mn'|'pt5mn') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
292 |
testoutput_var $1 "trcstat_ptracer0"$ii"_min" "p0"$ii"_min" $2 $3; yy=$? ;; |
testoutput_var $1 "trcstat_ptracer0"$ii"_min" "p0"$ii"_min" $2 $3; yy=$? ;; |
293 |
'pt1mx'|'pt2mx'|'pt3mx'|'pt4mx'|'pt5mx') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
'pt1mx'|'pt2mx'|'pt3mx'|'pt4mx'|'pt5mx') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
295 |
'pt1av'|'pt2av'|'pt3av'|'pt4av'|'pt5av') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
'pt1av'|'pt2av'|'pt3av'|'pt4av'|'pt5av') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
296 |
testoutput_var $1 "trcstat_ptracer0"$ii"_mean" "p0"$ii"_mean" $2 $3; yy=$? ;; |
testoutput_var $1 "trcstat_ptracer0"$ii"_mean" "p0"$ii"_mean" $2 $3; yy=$? ;; |
297 |
'pt1sd'|'pt2sd'|'pt3sd'|'pt4sd'|'pt5sd') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
'pt1sd'|'pt2sd'|'pt3sd'|'pt4sd'|'pt5sd') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
298 |
testoutput_var $1 "trcstat_ptracer0"$ii"_sd" "p0"$ii"_sd" $2 $3; yy=$? ;; |
testoutput_var $1 "trcstat_ptracer0"$ii"_sd" "p0"$ii"_StDv" $2 $3; yy=$? ;; |
299 |
*) yy=99; echo "WARNING: asking for var=$xx : not recognized !" 1>&2 ;; |
*) yy=99; echo "WARNING: asking for var=$xx : not recognized !" 1>&2 ;; |
300 |
esac |
esac |
301 |
if test $xx = $sVar |
if test $xx = $sVar |
712 |
{ |
{ |
713 |
if [ $# -eq 1 ]; then |
if [ $# -eq 1 ]; then |
714 |
for arg in * ; do |
for arg in * ; do |
715 |
test -d $arg/$1 && echo $arg |
test -f $arg/$1 && echo $arg |
716 |
done |
done |
717 |
else |
else |
718 |
echo $* |
echo $* |
897 |
NODEPEND=t |
NODEPEND=t |
898 |
fi |
fi |
899 |
|
|
900 |
|
#- setting for forward or ADM testing |
901 |
|
if test "x$ADM" = xt ; then |
902 |
|
code_dir=code_ad |
903 |
|
input_dirs='input_ad input' |
904 |
|
ref_outp="output_adm.txt" |
905 |
|
EXECUTABLE="mitgcmuv_ad" |
906 |
|
else |
907 |
|
code_dir=code |
908 |
|
input_dirs='input' |
909 |
|
ref_outp="output.txt" |
910 |
|
EXECUTABLE="mitgcmuv" |
911 |
|
fi |
912 |
|
|
913 |
if test "x$TESTDIRS" = x ; then |
if test "x$TESTDIRS" = x ; then |
914 |
if test "x$ADM" = xt ; then |
if test "x$ADM" = xt ; then |
915 |
LIST=`scandirs results_ad` |
LIST=`scandirs results_ad/$ref_outp` |
916 |
else |
else |
917 |
LIST=`scandirs results` |
LIST=`scandirs results/$ref_outp` |
918 |
fi |
fi |
919 |
else |
else |
920 |
#- expand group of experiments: |
#- expand group of experiments: |
959 |
OPTFILE=$MITGCM_OF |
OPTFILE=$MITGCM_OF |
960 |
fi |
fi |
961 |
|
|
962 |
if test "x$ADM" = xt ; then |
OUTPUTFILE=$ref_outp |
|
EXECUTABLE="mitgcmuv_ad" |
|
|
OUTPUTFILE="output_adm.txt" |
|
|
else |
|
|
EXECUTABLE="mitgcmuv" |
|
|
OUTPUTFILE="output.txt" |
|
|
fi |
|
|
|
|
963 |
if test "x$COMMAND" = x ; then |
if test "x$COMMAND" = x ; then |
964 |
COMMAND="./$EXECUTABLE > $OUTPUTFILE" |
COMMAND="./$EXECUTABLE > $OUTPUTFILE" |
965 |
fi |
fi |
974 |
# (use default or load experiment-specific list from file "tr_checklist") |
# (use default or load experiment-specific list from file "tr_checklist") |
975 |
# content : 1rst = main variable used to decide if it pass or FAIL |
# content : 1rst = main variable used to decide if it pass or FAIL |
976 |
# others = number of matching digits to be printed in summary.txt |
# others = number of matching digits to be printed in summary.txt |
977 |
if test "x$ADM" != xt ; then |
if test "x$ADM" = x ; then |
978 |
DEF_CHECK_LIST='PS PS T+ S+ U+ V+ pt1+ pt2+ pt3+ pt4+ pt5+' |
DEF_CHECK_LIST='PS PS T+ S+ U+ V+ pt1+ pt2+ pt3+ pt4+ pt5+' |
979 |
EMPTY_RESULTS='.. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. ..' |
EMPTY_RESULTS='.. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. ..' |
980 |
LEN_CHECK_LIST=`echo $DEF_CHECK_LIST | sed 's/ [a-zA-Z0-9]*+/&mn &mx &av &sd/g' | awk '{print NF-1}'` |
LEN_CHECK_LIST=`echo $DEF_CHECK_LIST | sed 's/ [a-zA-Z0-9]*+/&mn &mx &av &sd/g' | awk '{print NF-1}'` |
981 |
ii=`echo $EMPTY_RESULTS | awk '{print NF}'` |
ii=`echo $EMPTY_RESULTS | awk '{print NF}'` |
982 |
EMPTY_RESULTS=$EMPTY_RESULTS`expr $LEN_CHECK_LIST - $ii | awk 'BEGIN{FS=":"}{for(i=1;i<=$1;i++){printf " ."}}'` |
EMPTY_RESULTS=$EMPTY_RESULTS`expr $LEN_CHECK_LIST - $ii | awk 'BEGIN{FS=":"}{for(i=1;i<=$1;i++){printf " ."}}'` |
983 |
|
else |
984 |
|
DEF_CHECK_LIST='Cost Cost Grad' |
985 |
|
EMPTY_RESULTS='.. ..' |
986 |
|
LEN_CHECK_LIST=`echo $DEF_CHECK_LIST | sed 's/ [a-zA-Z0-9]*+/&mn &mx &av &sd/g' | awk '{print NF-1}'` |
987 |
fi |
fi |
988 |
|
|
989 |
# create the FORTRAN comparison code |
# create the FORTRAN comparison code |
1110 |
# results in the correct location--or skip this directory! |
# results in the correct location--or skip this directory! |
1111 |
fout= |
fout= |
1112 |
if test "x$ADM" = x ; then |
if test "x$ADM" = x ; then |
1113 |
fout=$dir"/results/output.txt" |
fout=$dir"/results/"$ref_outp |
1114 |
else |
else |
1115 |
fout=$dir"/results_ad/output_adm.txt" |
fout=$dir"/results_ad/"$ref_outp |
1116 |
fi |
fi |
1117 |
if test ! -r $fout ; then |
if test ! -r $fout ; then |
1118 |
echo "can't read \"$fout\" -- skipping $dir" |
echo "can't read \"$fout\" -- skipping $dir" |
1127 |
builddir="build" |
builddir="build" |
1128 |
if test ! -d $dir/$builddir ; then mkdir $dir/$builddir ; fi |
if test ! -d $dir/$builddir ; then mkdir $dir/$builddir ; fi |
1129 |
rundir="run" |
rundir="run" |
|
#rundir=$builddir |
|
1130 |
if test ! -d $dir/$rundir ; then |
if test ! -d $dir/$rundir ; then |
1131 |
rundir=$builddir |
rundir=$builddir |
1132 |
fi |
fi |
1133 |
|
CODE_DIR=$dir/$code_dir |
|
if test "x$ADM" = x ; then |
|
|
code_dir=code |
|
|
CODE_DIR=$dir/code |
|
|
input_dirs='input' |
|
|
else |
|
|
code_dir=code_ad |
|
|
CODE_DIR=$dir/code_ad |
|
|
input_dirs='input_ad input' |
|
|
fi |
|
1134 |
BUILD_DIR=$dir/$builddir |
BUILD_DIR=$dir/$builddir |
1135 |
|
|
1136 |
if test ! -r $CODE_DIR"/SIZE.h_mpi" -a "x$MPI" = "xt" ; then |
if test ! -r $CODE_DIR"/SIZE.h_mpi" -a "x$MPI" = "xt" ; then |
1192 |
&& makemodel $dir/$builddir && make=Y \ |
&& makemodel $dir/$builddir && make=Y \ |
1193 |
&& linkdata $dir/$rundir $input_dirs \ |
&& linkdata $dir/$rundir $input_dirs \ |
1194 |
&& runmodel $dir/$rundir && run=Y \ |
&& runmodel $dir/$rundir && run=Y \ |
1195 |
&& results=`testoutput_run $dir $rundir "txt"` |
&& results=`testoutput_run $dir $rundir $ref_outp` |
1196 |
fi |
fi |
1197 |
|
|
1198 |
echo |
echo |
1218 |
run_clean $dir/tr_run.$ex |
run_clean $dir/tr_run.$ex |
1219 |
linkdata $dir/tr_run.$ex input.$ex input |
linkdata $dir/tr_run.$ex input.$ex input |
1220 |
runmodel $dir/tr_run.$ex && run=Y \ |
runmodel $dir/tr_run.$ex && run=Y \ |
1221 |
&& results=`testoutput_run $dir tr_run.$ex ${ex}".txt"` |
&& results=`testoutput_run $dir tr_run.$ex "output.${ex}.txt"` |
1222 |
fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results` |
fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results` |
1223 |
fres="$fres.$ex" |
fres="$fres.$ex" |
1224 |
echo |
echo |
1244 |
echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt" |
echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt" |
1245 |
echo "DATE='$DATE'" >> $CDIR"/summary.txt" |
echo "DATE='$DATE'" >> $CDIR"/summary.txt" |
1246 |
echo "tdir='$dir'" >> $CDIR"/summary.txt" |
echo "tdir='$dir'" >> $CDIR"/summary.txt" |
1247 |
|
fi |
1248 |
|
if test "x$ADM" = xt ; then |
1249 |
grep -A3 'Seconds in section "ALL' $dir/$rundir/$OUTPUTFILE \ |
grep -A3 'Seconds in section "ALL' $dir/$rundir/$OUTPUTFILE \ |
1250 |
>> $CDIR"/summary.txt" |
>> $CDIR"/summary.txt" |
1251 |
fi |
fi |