/[MITgcm]/MITgcm_contrib/jmc_script/tst1+1
ViewVC logotype

Contents of /MITgcm_contrib/jmc_script/tst1+1

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


Revision 1.1 - (show annotations) (download)
Sun Oct 12 14:46:35 2003 UTC (21 years, 9 months ago) by jmc
Branch: MAIN
Some basic scripts are gather here:
a) the restart test : "1+1=2"
b) PBS batch script that resubmit itself (loop).

1 #!/bin/sh
2
3 #- input :
4 # 1) data.tst = parameter file "data" for 2 iter run
5 # 2) in dir temp : pickup(_cd).0..0{num0}.001.001.data
6
7 #- works only when start from iter that is a multiple of 10.
8 num0="3600"
9
10 exec=mitgcmuv
11 # list of pickup needed for a restart :
12 #listP='pickup pickup_cd'
13 listP='pickup pickup_ic'
14
15 #- dir where to put the results :
16 dir1=res_2it
17 dir2=res_1iA
18 dir3=res_1iB
19 #---------------------------
20 rm -r $dir1 $dir2 $dir3
21 mkdir $dir1 $dir2 $dir3
22 num1=`expr $num0 + 1`
23 num2=`expr $num0 + 2`
24 num0c=`add0upto10c $num0`
25 num1c=`add0upto10c $num1`
26 num2c=`add0upto10c $num2`
27 echo $num0 $num1 $num2
28 echo $num0c $num1c $num2c
29 echo 'exec='$exec
30 # exit
31 #------------------------
32 cp -p data.tst data
33 ln -s temp/pickup*.${num0c}* .
34 echo '==> START RUN 2it'
35 ./$exec > std_outp.2it 2>&1
36 out=`egrep -c 'STOP ABNORMAL END' std_outp.2it`
37 if test $out != 0 ; then echo 'RUN 2it STOP before END => exit' ; exit ; fi
38 echo '==> END RUN 2it'
39 for xx in $listP
40 do
41 rnp $xx.ckptA $xx.$num2c
42 done
43 mv *.00?.00?.* $dir1
44 mv $dir1/pickup*.${num0c}* .
45 #--
46 sed 's/^ nTimeSteps=2/ nTimeSteps=1/g' data.tst > data
47 echo '==> START RUN 1iA'
48 ./$exec > std_outp.1iA 2>&1
49 out=`egrep -c 'STOP ABNORMAL END' std_outp.1iA`
50 if test $out != 0 ; then echo 'RUN 1iA STOP before END => exit' ; exit ; fi
51 echo '==> END RUN 1iA'
52 for xx in $listP
53 do
54 rnp $xx.ckptA $xx.$num1c
55 done
56 mv *.00?.00?.* $dir2
57 rm $dir2/pickup*.${num0c}*
58 #--
59 ln -s $dir2/pickup*.$num1c.* .
60 mv data data.tmp
61 sed "s/^ nIter0=$num0/ nIter0=$num1/g" data.tmp > data
62 echo '==> START RUN 1iB'
63 ./$exec > std_outp.1iB 2>&1
64 out=`egrep -c 'STOP ABNORMAL END' std_outp.1iB`
65 if test $out != 0 ; then echo 'RUN 1iB STOP before END => exit' ; exit ; fi
66 echo '==> END RUN 1iB'
67 for xx in $listP
68 do
69 rnp $xx.ckptA $xx.$num2c
70 done
71 mv *.00?.00?.* $dir3
72 rm $dir3/pickup*.$num1c.*
73 #--
74 echo ' ' ; echo '-- compare cg2d_init_res :'
75 echo ' run 1iA:'
76 grep "cg2d_init_res" std_outp.1iA \
77 | sed 's/.* cg2d_init_res =//'
78 echo ' run 1iB:'
79 grep "cg2d_init_res" std_outp.1iB \
80 | sed 's/.* cg2d_init_res =//'
81 echo ' run 2it:'
82 grep "cg2d_init_res" std_outp.2it \
83 | sed 's/.* cg2d_init_res =//'
84 exit

  ViewVC Help
Powered by ViewVC 1.1.22