Annotation of /MITgcm_contrib/jmc_script/tst1+1

#!/bin/sh

# Nit = Nb of iter of 2nd & 3rd run (1rst one is 2 x Nit iter long)
Dbl=`sed -n 's/nTimeSteps=//p' data.tst | sed 's/,//g'`
Nit=`expr $Dbl / 2`
#Nit=2

case $# in
  1) ;;
   *) echo "Usage: `basename $0` flag"
      echo " Check restart: compare 1 run of 2 x $Nit it long"
      echo "     with 2 consecutive runs of $Nit it long each"
      echo " where: flag = 1 -> do the 3 run"
      echo "        flag = 2 -> compare std_outp"
      echo "        flag = 3 -> compare pickup files"
      exit ;;
esac
flag=$1

#- input :
# 1) data.tst = parameter file "data" for 2.Nit iter run
# 2) in dir temp : pickup(_cd).0..0{num0}.001.001.data

num0=`sed -n 's/nIter0=//p' data.tst | sed 's/,//g' | sed 's/ //g'`
#num0="72000"

# list of pickup(s) that are needed for a restart :
listP='pickup'
#listP='pickup pickup_cd pickup_nh'

mpi=0
if [ $mpi -ge 1 ]
then gcmExc=run_ifc_mpi
else gcmExc=mitgcmuv
fi

#- dir where to put the results :
dir1=res_2it
dir2=res_1iA
dir3=res_1iB

#make sure that local dir is in the path:
export PATH=${PATH}:.

#---------------------------
Dbl=`expr $Nit + $Nit`
num1=`expr $num0 + $Nit`
num2=`expr $num0 + $Dbl`
num0c=`add0upto10c $num0`
num1c=`add0upto10c $num1`
num2c=`add0upto10c $num2`
echo $num0 $num1 $num2
echo $num0c $num1c $num2c
echo 'gcmExc='$gcmExc
# exit
#------------------------

if test $flag = 1 ; then
#------------
rm -r $dir1 $dir2 $dir3
mkdir $dir1 $dir2 $dir3

cp -p data.tst data
ln -s temp/pickup*.${num0c}* .
echo '==> START RUN 2 x' $Nit 'it'
if [ $mpi -ge 1 ]
then ./$gcmExc std_outp.2it > TTT.out.2it 2>&1
else ./$gcmExc > std_outp.2it 2>&1
fi
if test -f pickup.ckptA.001.001.data
then out=`egrep -c 'STOP ABNORMAL END' std_outp.2it`
else out=9 ; fi
if test $out != 0 ; then echo 'RUN 2it STOP before END => exit' ; exit ; fi
echo '==>  END  RUN 2 x' $Nit 'it'
for xx in $listP
do
 rnp -s $xx.ckptA $xx.$num2c
done
mv *.0??.0??.* $dir1
mv $dir1/pickup*.${num0c}* .
#--
sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$Nit/g" data.tst > data
echo '==> START RUN 1iA'
if [ $mpi -ge 1 ]
then ./$gcmExc std_outp.1iA > TTT.out.1iA 2>&1
else ./$gcmExc > std_outp.1iA 2>&1
fi
if test -f pickup.ckptA.001.001.data
then out=`egrep -c 'STOP ABNORMAL END' std_outp.1iA`
else out=9 ; fi
if test $out != 0 ; then echo 'RUN 1iA STOP before END => exit' ; exit ; fi
echo '==>  END  RUN 1iA'
for xx in $listP
do
 rnp -s $xx.ckptA $xx.$num1c
done
mv *.0??.0??.* $dir2
rm $dir2/pickup*.${num0c}*
#--
ln -s $dir2/pickup*.$num1c.* .
mv data data.tmp
sed "s/^ nIter0=$num0/ nIter0=$num1/g" data.tmp > data
echo '==> START RUN 1iB'
if [ $mpi -ge 1 ]
then ./$gcmExc std_outp.1iB > TTT.out.1iB 2>&1
else ./$gcmExc > std_outp.1iB 2>&1
fi
if test -f pickup.ckptA.001.001.data
then out=`egrep -c 'STOP ABNORMAL END' std_outp.1iB`
else out=9 ; fi
if test $out != 0 ; then echo 'RUN 1iB STOP before END => exit' ; exit ; fi
echo '==>  END  RUN 1iB'
for xx in $listP
do
 rnp -s $xx.ckptA $xx.$num2c
done
mv *.0??.0??.* $dir3
rm $dir3/pickup*.$num1c.*
exit
#------------
fi

if test $flag = 2 ; then
#------------
echo ' ' ; echo '-- compare cg2d_init_res :'
echo ' run 1iA:'
 grep "cg2d_init_res" std_outp.1iA \
      | sed 's/.* cg2d_init_res =//'
echo ' run 1iB:'
 grep "cg2d_init_res" std_outp.1iB \
      | sed 's/.* cg2d_init_res =//'
echo ' run 2it:'
 grep "cg2d_init_res" std_outp.2it \
      | sed 's/.* cg2d_init_res =//'
exit
#------------
fi

if test $flag = 3 ; then
#------------
difex=`ls ~jmc/bin/dif_data`
for xx in $listP
do
 $difex $dir1 $dir3 $xx.$num2c
done
exit
#------------
fi
1	jmc	1.1	#!/bin/sh
2
3	jmc	1.2	# Nit = Nb of iter of 2nd & 3rd run (1rst one is 2 x Nit iter long)
4	jmc	1.4	Dbl=`sed -n 's/nTimeSteps=//p' data.tst \| sed 's/,//g'`
5			Nit=`expr $Dbl / 2`
6			#Nit=2
7	jmc	1.2
8			case $# in
9			1) ;;
10			*) echo "Usage: `basename $0` flag"
11			echo " Check restart: compare 1 run of 2 x $Nit it long"
12			echo " with 2 consecutive runs of $Nit it long each"
13			echo " where: flag = 1 -> do the 3 run"
14			echo " flag = 2 -> compare std_outp"
15			echo " flag = 3 -> compare pickup files"
16			exit ;;
17			esac
18			flag=$1
19
20	jmc	1.5	#- input :
21	jmc	1.2	# 1) data.tst = parameter file "data" for 2.Nit iter run
22	jmc	1.1	# 2) in dir temp : pickup(_cd).0..0{num0}.001.001.data
23
24	jmc	1.4	num0=`sed -n 's/nIter0=//p' data.tst \| sed 's/,//g' \| sed 's/ //g'`
25			#num0="72000"
26	jmc	1.1
27	jmc	1.2	# list of pickup(s) that are needed for a restart :
28			listP='pickup'
29	jmc	1.5	#listP='pickup pickup_cd pickup_nh'
30
31			mpi=0
32			if [ $mpi -ge 1 ]
33			then gcmExc=run_ifc_mpi
34			else gcmExc=mitgcmuv
35			fi
36	jmc	1.1
37			#- dir where to put the results :
38			dir1=res_2it
39			dir2=res_1iA
40			dir3=res_1iB
41	jmc	1.5
42			#make sure that local dir is in the path:
43			export PATH=${PATH}:.
44
45	jmc	1.1	#---------------------------
46	jmc	1.2	Dbl=`expr $Nit + $Nit`
47			num1=`expr $num0 + $Nit`
48			num2=`expr $num0 + $Dbl`
49	jmc	1.1	num0c=`add0upto10c $num0`
50			num1c=`add0upto10c $num1`
51			num2c=`add0upto10c $num2`
52			echo $num0 $num1 $num2
53	jmc	1.5	echo $num0c $num1c $num2c
54	jmc	1.2	echo 'gcmExc='$gcmExc
55	jmc	1.1	# exit
56			#------------------------
57	jmc	1.2
58			if test $flag = 1 ; then
59			#------------
60			rm -r $dir1 $dir2 $dir3
61			mkdir $dir1 $dir2 $dir3
62
63	jmc	1.1	cp -p data.tst data
64			ln -s temp/pickup.${num0c} .
65	jmc	1.2	echo '==> START RUN 2 x' $Nit 'it'
66	jmc	1.5	if [ $mpi -ge 1 ]
67			then ./$gcmExc std_outp.2it > TTT.out.2it 2>&1
68			else ./$gcmExc > std_outp.2it 2>&1
69			fi
70			if test -f pickup.ckptA.001.001.data
71	jmc	1.2	then out=`egrep -c 'STOP ABNORMAL END' std_outp.2it`
72			else out=9 ; fi
73	jmc	1.1	if test $out != 0 ; then echo 'RUN 2it STOP before END => exit' ; exit ; fi
74	jmc	1.2	echo '==> END RUN 2 x' $Nit 'it'
75	jmc	1.1	for xx in $listP
76			do
77	jmc	1.4	rnp -s $xx.ckptA $xx.$num2c
78	jmc	1.1	done
79	jmc	1.3	mv .0??.0??. $dir1
80	jmc	1.1	mv $dir1/pickup.${num0c} .
81			#--
82	jmc	1.2	sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$Nit/g" data.tst > data
83	jmc	1.5	echo '==> START RUN 1iA'
84			if [ $mpi -ge 1 ]
85			then ./$gcmExc std_outp.1iA > TTT.out.1iA 2>&1
86			else ./$gcmExc > std_outp.1iA 2>&1
87			fi
88			if test -f pickup.ckptA.001.001.data
89	jmc	1.2	then out=`egrep -c 'STOP ABNORMAL END' std_outp.1iA`
90			else out=9 ; fi
91	jmc	1.1	if test $out != 0 ; then echo 'RUN 1iA STOP before END => exit' ; exit ; fi
92			echo '==> END RUN 1iA'
93			for xx in $listP
94			do
95	jmc	1.4	rnp -s $xx.ckptA $xx.$num1c
96	jmc	1.1	done
97	jmc	1.3	mv .0??.0??. $dir2
98	jmc	1.1	rm $dir2/pickup.${num0c}
99			#--
100	jmc	1.5	ln -s $dir2/pickup.$num1c. .
101	jmc	1.1	mv data data.tmp
102			sed "s/^ nIter0=$num0/ nIter0=$num1/g" data.tmp > data
103	jmc	1.5	echo '==> START RUN 1iB'
104			if [ $mpi -ge 1 ]
105			then ./$gcmExc std_outp.1iB > TTT.out.1iB 2>&1
106			else ./$gcmExc > std_outp.1iB 2>&1
107			fi
108			if test -f pickup.ckptA.001.001.data
109	jmc	1.2	then out=`egrep -c 'STOP ABNORMAL END' std_outp.1iB`
110			else out=9 ; fi
111	jmc	1.1	if test $out != 0 ; then echo 'RUN 1iB STOP before END => exit' ; exit ; fi
112			echo '==> END RUN 1iB'
113			for xx in $listP
114			do
115	jmc	1.4	rnp -s $xx.ckptA $xx.$num2c
116	jmc	1.1	done
117	jmc	1.3	mv .0??.0??. $dir3
118	jmc	1.1	rm $dir3/pickup.$num1c.
119	jmc	1.2	exit
120			#------------
121			fi
122
123			if test $flag = 2 ; then
124			#------------
125	jmc	1.1	echo ' ' ; echo '-- compare cg2d_init_res :'
126			echo ' run 1iA:'
127			grep "cg2d_init_res" std_outp.1iA \
128			\| sed 's/.* cg2d_init_res =//'
129			echo ' run 1iB:'
130			grep "cg2d_init_res" std_outp.1iB \
131			\| sed 's/.* cg2d_init_res =//'
132			echo ' run 2it:'
133			grep "cg2d_init_res" std_outp.2it \
134	jmc	1.5	\| sed 's/.* cg2d_init_res =//'
135	jmc	1.1	exit
136	jmc	1.2	#------------
137			fi
138
139			if test $flag = 3 ; then
140			#------------
141	jmc	1.4	difex=`ls ~jmc/bin/dif_data`
142	jmc	1.2	for xx in $listP
143			do
144			$difex $dir1 $dir3 $xx.$num2c
145			done
146			exit
147			#------------
148			fi