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

Annotation of /MITgcm_contrib/jmc_script/tst1+1

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


Revision 1.1 - (hide 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 jmc 1.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