17 |
esac |
esac |
18 |
flag=$1 |
flag=$1 |
19 |
|
|
20 |
|
#- input : |
|
#- input : |
|
21 |
# 1) data.tst = parameter file "data" for 2.Nit iter run |
# 1) data.tst = parameter file "data" for 2.Nit iter run |
22 |
# 2) in dir temp : pickup(_cd).0..0{num0}.001.001.data |
# 2) in dir temp : pickup(_cd).0..0{num0}.001.001.data |
23 |
|
|
24 |
num0=`sed -n 's/nIter0=//p' data.tst | sed 's/,//g' | sed 's/ //g'` |
num0=`sed -n 's/nIter0=//p' data.tst | sed 's/,//g' | sed 's/ //g'` |
25 |
#num0="72000" |
#num0="72000" |
26 |
|
|
|
gcmExc=mitgcmuv |
|
27 |
# list of pickup(s) that are needed for a restart : |
# list of pickup(s) that are needed for a restart : |
28 |
listP='pickup' |
listP='pickup' |
29 |
#listP='pickup pickup_cd' |
#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 |
|
|
37 |
#- dir where to put the results : |
#- dir where to put the results : |
38 |
dir1=res_2it |
dir1=res_2it |
39 |
dir2=res_1iA |
dir2=res_1iA |
40 |
dir3=res_1iB |
dir3=res_1iB |
41 |
|
|
42 |
|
#make sure that local dir is in the path: |
43 |
|
export PATH=${PATH}:. |
44 |
|
|
45 |
#--------------------------- |
#--------------------------- |
46 |
Dbl=`expr $Nit + $Nit` |
Dbl=`expr $Nit + $Nit` |
47 |
num1=`expr $num0 + $Nit` |
num1=`expr $num0 + $Nit` |
50 |
num1c=`add0upto10c $num1` |
num1c=`add0upto10c $num1` |
51 |
num2c=`add0upto10c $num2` |
num2c=`add0upto10c $num2` |
52 |
echo $num0 $num1 $num2 |
echo $num0 $num1 $num2 |
53 |
echo $num0c $num1c $num2c |
echo $num0c $num1c $num2c |
54 |
echo 'gcmExc='$gcmExc |
echo 'gcmExc='$gcmExc |
55 |
# exit |
# exit |
56 |
#------------------------ |
#------------------------ |
63 |
cp -p data.tst data |
cp -p data.tst data |
64 |
ln -s temp/pickup*.${num0c}* . |
ln -s temp/pickup*.${num0c}* . |
65 |
echo '==> START RUN 2 x' $Nit 'it' |
echo '==> START RUN 2 x' $Nit 'it' |
66 |
./$gcmExc > std_outp.2it 2>&1 |
if [ $mpi -ge 1 ] |
67 |
if test -f pickup.ckptA.001.001.data |
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 |
then out=`egrep -c 'STOP ABNORMAL END' std_outp.2it` |
then out=`egrep -c 'STOP ABNORMAL END' std_outp.2it` |
72 |
else out=9 ; fi |
else out=9 ; fi |
73 |
if test $out != 0 ; then echo 'RUN 2it STOP before END => exit' ; exit ; fi |
if test $out != 0 ; then echo 'RUN 2it STOP before END => exit' ; exit ; fi |
80 |
mv $dir1/pickup*.${num0c}* . |
mv $dir1/pickup*.${num0c}* . |
81 |
#-- |
#-- |
82 |
sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$Nit/g" data.tst > data |
sed "s/^ nTimeSteps=$Dbl/ nTimeSteps=$Nit/g" data.tst > data |
83 |
echo '==> START RUN 1iA' |
echo '==> START RUN 1iA' |
84 |
./$gcmExc > std_outp.1iA 2>&1 |
if [ $mpi -ge 1 ] |
85 |
if test -f pickup.ckptA.001.001.data |
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 |
then out=`egrep -c 'STOP ABNORMAL END' std_outp.1iA` |
then out=`egrep -c 'STOP ABNORMAL END' std_outp.1iA` |
90 |
else out=9 ; fi |
else out=9 ; fi |
91 |
if test $out != 0 ; then echo 'RUN 1iA STOP before END => exit' ; exit ; fi |
if test $out != 0 ; then echo 'RUN 1iA STOP before END => exit' ; exit ; fi |
97 |
mv *.0??.0??.* $dir2 |
mv *.0??.0??.* $dir2 |
98 |
rm $dir2/pickup*.${num0c}* |
rm $dir2/pickup*.${num0c}* |
99 |
#-- |
#-- |
100 |
ln -s $dir2/pickup*.$num1c.* . |
ln -s $dir2/pickup*.$num1c.* . |
101 |
mv data data.tmp |
mv data data.tmp |
102 |
sed "s/^ nIter0=$num0/ nIter0=$num1/g" data.tmp > data |
sed "s/^ nIter0=$num0/ nIter0=$num1/g" data.tmp > data |
103 |
echo '==> START RUN 1iB' |
echo '==> START RUN 1iB' |
104 |
./$gcmExc > std_outp.1iB 2>&1 |
if [ $mpi -ge 1 ] |
105 |
if test -f pickup.ckptA.001.001.data |
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 |
then out=`egrep -c 'STOP ABNORMAL END' std_outp.1iB` |
then out=`egrep -c 'STOP ABNORMAL END' std_outp.1iB` |
110 |
else out=9 ; fi |
else out=9 ; fi |
111 |
if test $out != 0 ; then echo 'RUN 1iB STOP before END => exit' ; exit ; fi |
if test $out != 0 ; then echo 'RUN 1iB STOP before END => exit' ; exit ; fi |
131 |
| sed 's/.* cg2d_init_res =//' |
| sed 's/.* cg2d_init_res =//' |
132 |
echo ' run 2it:' |
echo ' run 2it:' |
133 |
grep "cg2d_init_res" std_outp.2it \ |
grep "cg2d_init_res" std_outp.2it \ |
134 |
| sed 's/.* cg2d_init_res =//' |
| sed 's/.* cg2d_init_res =//' |
135 |
exit |
exit |
136 |
#------------ |
#------------ |
137 |
fi |
fi |