1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
    #!/bin/bash
    rm -r ./ACC
    rm -r ./PSA
    rm -r ./other
    rm ACC.out
    lat=30.025
    lon=100.13
    stk=140.0
    dip=90.0
    depth=2
    L=95.0
    W=25.0
    dl=5.0
    dw=5.0
    hypoL=47.5 # position center 2.5,8
    hypoW=12.5
    num=0
    outfile=input
    inputfile=input.params
    for stress in $(seq 30 5 100)
    do
        for kappa in 0.002 0.003 0.004 0.005 0.006 0.007 0.008 
        do 
            num=`expr $num + 1`
            mkdir $num
            sed -i '7c \ 8.0 '$stress' 1 '$kappa' ' $inputfile
            sed -i '9c \ '$lat'  '$lon' ' $inputfile
            sed -i '11c \ '$stk'  '$dip' '$depth' ' $inputfile
            sed -i '29c \ '$L'  '$W' '$dl' '$dw' 70.0' $inputfile
            sed -i '36c \ '$hypoL' '$hypoW' 1' $inputfile
            sed -i '97c \  1' $inputfile  # 1 random 0 specify
            sed -i '93c \ '$RANDOM' 5' $inputfile
           ./exsim14 $inputfile
            cp $inputfile $num
            mv ACC PSA other ACC.out psa.out $num
            mv $num /media/aicesun/Document/litang/fault1/mw75/random/center/
        done
    done