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
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
|
o1 = deta**2
o2 = 1/o1
o3 = dq**2
o4 = 1/o3
o5 = 1/dq
o6 = tan(q)
o7 = 1/o6
o8 = dphi**2
o9 = 1/o8
o10 = eta**2
o11 = sigma**2
o12 = 1/o11
o13 = -(o10*o12)
o14 = -2.00000000000000d0*eta*eta0*o12
o15 = eta0**2
o16 = -(o12*o15)
o17 = o13 + o14 + o16
o18 = exp(o17)
o19 = sin(q)
o20 = o19**n
o21 = 2.00000000000000d0*amp*o18*o20
o22 = cos(phi)
o23 = o22**2
o24 = 2.00000000000000d0*amp*c*o18*o20*o23
o25 = 2.00000000000000d0*eta*eta0*o12
o26 = o13 + o16 + o25
o27 = exp(o26)
o28 = 2.00000000000000d0*amp*o20*o27
o29 = 2.00000000000000d0*amp*c*o20*o23*o27
o30 = o21 + o24 + o28 + o29
o31 = exp(o30)
o32 = o19**2
o33 = 1/o32
o34 = 1.00000000000000d0*o31*o33*o9
o35 = 1/dphi
o36 = o13 + o14 + o16 + o21 + o24 + o28 + o29
o37 = exp(o36)
o38 = sin(phi)
o39 = -2.00000000000000d0 + n
o40 = o19**o39
o41 = o13 + o16 + o21 + o24 + o25 + o28 + o29
o42 = exp(o41)
o43 = cos(q)
o44 = o43**2
o45 = n**2
o46 = o38**2
o47 = o11**2
o48 = 1/o47
o49 = amp**2
o50 = c**2
o51 = -2.00000000000000d0*o10*o12
o52 = -2.00000000000000d0*o12*o15
o53 = o21 + o24 + o28 + o29 + o51 + o52
o54 = exp(o53)
o55 = 2.00000000000000d0*n
o56 = -2.00000000000000d0 + o55
o57 = o19**o56
o58 = -4.0000000000000d0*eta*eta0*o12
o59 = o21 + o24 + o28 + o29 + o51 + o52 + o58
o60 = exp(o59)
o61 = 4.0000000000000d0*eta*eta0*o12
o62 = o21 + o24 + o28 + o29 + o51 + o52 + o61
o63 = exp(o62)
o64 = o6**2
o65 = 1/o64
o66 = 5.0000000000000d-1*eta
o67 = cosh(o66)
Ae(i,j,k) = o2
Aw(i,j,k) = o2
An(i,j,k) = o4 + 5.0000000000000d-1*o5*o7
As(i,j,k) = o4 - 5.0000000000000d-1*o5*o7
Ab(i,j,k) = o34 - 2.00000000000000d0*amp*c*o22*o35*o37*o38*o40 -
& 2.00000000000000d0*amp*c*o22*o35*o38*o40*o42
Aq(i,j,k) = o34 + 2.00000000000000d0*amp*c*o22*o35*o37*o38*o40 +
& 2.00000000000000d0*amp*c*o22*o35*o38*o40*o42
Ac(i,j,k) = -1.25000000000000d-1 - 2.00000000000000d0*o2 - 2.500
& 00000000000d-1*amp*n*o18*o20 - 5.0000000000000d-1*amp*o12*o18*o2
& 0 - 2.50000000000000d-1*amp*c*n*o18*o20*o23 - 5.0000000000000d-1
& *amp*c*o12*o18*o20*o23 - 2.50000000000000d-1*amp*n*o20*o27 - 5.0
& 000000000000d-1*amp*o12*o20*o27 - 2.50000000000000d-1*amp*c*n*o2
& 0*o23*o27 - 5.0000000000000d-1*amp*c*o12*o20*o23*o27 - 1.2500000
& 0000000d-1*o33 - 2.00000000000000d0*o4 - amp*c*o23*o37*o40 - amp
& *c*o23*o40*o42 - 2.50000000000000d-1*amp*n*o18*o40*o44 - 2.50000
& 000000000d-1*amp*c*n*o18*o23*o40*o44 - 2.50000000000000d-1*amp*n
& *o27*o40*o44 - 2.50000000000000d-1*amp*c*n*o23*o27*o40*o44 + 2.5
& 0000000000000d-1*amp*o18*o40*o44*o45 + 2.50000000000000d-1*amp*c
& *o18*o23*o40*o44*o45 + 2.50000000000000d-1*amp*o27*o40*o44*o45 +
& 2.50000000000000d-1*amp*c*o23*o27*o40*o44*o45 + amp*c*o37*o40*o
& 46 + amp*c*o40*o42*o46 + 2.00000000000000d0*amp*eta*eta0*o18*o20
& *o48 + amp*o10*o18*o20*o48 + amp*o15*o18*o20*o48 + 2.00000000000
& 000d0*amp*c*eta*eta0*o18*o20*o23*o48 + amp*c*o10*o18*o20*o23*o48
& + amp*c*o15*o18*o20*o23*o48 - 2.00000000000000d0*amp*eta*eta0*o
& 20*o27*o48 + amp*o10*o20*o27*o48 + amp*o15*o20*o27*o48 - 2.00000
& 000000000d0*amp*c*eta*eta0*o20*o23*o27*o48 + amp*c*o10*o20*o23*o
& 27*o48 + amp*c*o15*o20*o23*o27*o48 + 6.0000000000000d0*o23*o46*o
& 49*o50*o54*o57 + 3.00000000000000d0*o23*o46*o49*o50*o57*o60 + 3.
& 00000000000000d0*o23*o46*o49*o50*o57*o63 + 1.25000000000000d-1*o
& 65 - 2.00000000000000d0*o31*o33*o9
t1 = -2.50000000000000d-1*o67 + 5.0000000000000d-1*amp*n*o18*o20
& *o67 + amp*o12*o18*o20*o67 + 5.0000000000000d-1*amp*c*n*o18*o20*
& o23*o67 + amp*c*o12*o18*o20*o23*o67 + 5.0000000000000d-1*amp*n*o
& 20*o27*o67 + amp*o12*o20*o27*o67 + 5.0000000000000d-1*amp*c*n*o2
& 0*o23*o27*o67 + amp*c*o12*o20*o23*o27*o67 + 2.50000000000000d-1*
& o33*o67 + 2.00000000000000d0*amp*c*o23*o37*o40*o67 + 2.000000000
& 00000d0*amp*c*o23*o40*o42*o67 + 5.0000000000000d-1*amp*n*o18*o40
& *o44*o67 + 5.0000000000000d-1*amp*c*n*o18*o23*o40*o44*o67 + 5.00
& 00000000000d-1*amp*n*o27*o40*o44*o67 + 5.0000000000000d-1*amp*c*
& n*o23*o27*o40*o44*o67 - 5.0000000000000d-1*amp*o18*o40*o44*o45*o
& 67 - 5.0000000000000d-1*amp*c*o18*o23*o40*o44*o45*o67 - 5.000000
& 0000000d-1*amp*o27*o40*o44*o45*o67
t2 = -5.0000000000000d-1*amp*c*o23*o27*o40*o44*o45*o67 - 2.00000
& 000000000d0*amp*c*o37*o40*o46*o67 - 2.00000000000000d0*amp*c*o40
& *o42*o46*o67 - 4.0000000000000d0*amp*eta*eta0*o18*o20*o48*o67 -
& 2.00000000000000d0*amp*o10*o18*o20*o48*o67 - 2.00000000000000d0*
& amp*o15*o18*o20*o48*o67 - 4.0000000000000d0*amp*c*eta*eta0*o18*o
& 20*o23*o48*o67 - 2.00000000000000d0*amp*c*o10*o18*o20*o23*o48*o6
& 7 - 2.00000000000000d0*amp*c*o15*o18*o20*o23*o48*o67 + 4.0000000
& 000000d0*amp*eta*eta0*o20*o27*o48*o67 - 2.00000000000000d0*amp*o
& 10*o20*o27*o48*o67 - 2.00000000000000d0*amp*o15*o20*o27*o48*o67
& + 4.0000000000000d0*amp*c*eta*eta0*o20*o23*o27*o48*o67 - 2.00000
& 000000000d0*amp*c*o10*o20*o23*o27*o48*o67 - 2.00000000000000d0*a
& mp*c*o15*o20*o23*o27*o48*o67 - 1.20000000000000d1*o23*o46*o49*o5
& 0*o54*o57*o67 - 6.0000000000000d0*o23*o46*o49*o50*o57*o60*o67 -
& 6.0000000000000d0*o23*o46*o49*o50*o57*o63*o67 - 2.50000000000000
& d-1*o65*o67
Rhs(i,j,k) = t1 + t2
|
