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
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
|
o1 = mm**2
o2 = 5.0000000000000d1*o1
o3 = -8.0000000000000d1 + o2
o4 = 1/o3
o5 = 2.00000000000000d0*q(i,j,k)
o6 = cos(o5)
o7 = 5.0000000000000d0*o6
o8 = 3.00000000000000d0 + o7
o9 = dnnnnngtil
o10 = dnnngtil
o11 = -2.60000000000000d1*o10
o12 = dngtil
o13 = 2.50000000000000d1*o12
o14 = o11 + o13 + o9
o15 = -3.00000000000000d0*eta(i,j,k)
o16 = exp(o15)
o17 = sin(phi(i,j,k))
o18 = o17**2
o19 = mm*phi(i,j,k)
o20 = sin(o19)
o21 = sin(q(i,j,k))
o22 = o21**2
o23 = 2.50000000000000d1*o1
o24 = -4.0000000000000d1 + o23
o25 = 1/o24
o26 = cos(phi(i,j,k))
o27 = o26**2
o28 = -1.60000000000000d1 + o1
o29 = 3.00000000000000d0*o6
o30 = 2.00000000000000d0 + o29
o31 = 2.00000000000000d0*o28*o30*o9
o32 = -1.07000000000000d2*o1
o33 = 6.0900000000000d2*o6
o34 = -1.59000000000000d2*o1*o6
o35 = 4.0700000000000d2 + o32 + o33 + o34
o36 = o10*o35
o37 = 3.5000000000000d1*o1
o38 = 5.0000000000000d0*o1
o39 = -1.30000000000000d1 + o38
o40 = 9.0000000000000d0*o39*o6
o41 = -8.3000000000000d1 + o37 + o40
o42 = 5.0000000000000d0*o12*o41
o43 = o31 + o36 + o42
o44 = o22**2
o45 = cos(q(i,j,k))
o46 = o45**2
o47 = -3.00000000000000d0*mm*o14*o20*o22*o4*o8
o48 = -(mm*o20*o22*o25*o43)
o49 = o47 + o48
o50 = -8.0000000000000d0 + o38
o51 = 1/o50
o52 = cos(o19)
o53 = -9.0000000000000d0*o1
o54 = 5.0000000000000d0*o50*o6
o55 = 8.0000000000000d0 + o53 + o54
o56 = -7.5000000000000d1*o12*o55
o57 = 2.88000000000000d2*o9
o58 = -1.20000000000000d1*o1*o9
o59 = -9.4000000000000d1*o1
o60 = -4.1500000000000d2*o6
o61 = 1.75000000000000d2*o1*o6
o62 = 1.15000000000000d2 + o59 + o60 + o61
o63 = 4.0000000000000d0*o10*o62
o64 = dnnnnnnngtil
o65 = 3.2000000000000d1*o22*o64
o66 = -2.00000000000000d0*o1*o22*o64
o67 = -9.2800000000000d2*o22*o9
o68 = 1.03000000000000d2*o1*o22*o9
o69 = o57 + o58 + o63 + o65 + o66 + o67 + o68
o70 = 2.00000000000000d0*o69
o71 = o56 + o70
o72 = sin(o5)
o73 = dnnnnnngtil
o74 = 3.2000000000000d1*o73
o75 = -2.00000000000000d0*o1*o73
o76 = -3.2000000000000d1*o6*o73
o77 = 2.00000000000000d0*o1*o6*o73
o78 = dnnnngtil
o79 = -3.4000000000000d2*o78
o80 = 5.5000000000000d1*o1*o78
o81 = -8.1200000000000d2*o6*o78
o82 = 1.70000000000000d1*o1*o6*o78
o83 = 7.0000000000000d0*o1
o84 = -2.09000000000000d2*o6
o85 = 8.9000000000000d1*o1*o6
o86 = -3.10000000000000d1 + o83 + o84 + o85
o87 = dnngtil
o88 = -2.00000000000000d1*o86*o87
o89 = 7.0000000000000d0*o6
o90 = 1.00000000000000d0 + o89
o91 = gtil
o92 = 7.5000000000000d1*o50*o90*o91
o93 = o74 + o75 + o76 + o77 + o79 + o80 + o81 + o82 + o88 + o92
o94 = -2.40000000000000d2*byJ
o95 = 1.50000000000000d2*byJ*o1
o96 = -1.15200000000000d3*o46*o52*o78
o97 = 7.2000000000000d1*o1*o46*o52*o78
o98 = 3.9000000000000d1*o1
o99 = -2.68000000000000d2*o6
o100 = 1.00000000000000d2*o1*o6
o101 = 4.0000000000000d0*q(i,j,k)
o102 = cos(o101)
o103 = -5.8100000000000d2*o102
o104 = 2.45000000000000d2*o1*o102
o105 = -1.11000000000000d2 + o100 + o103 + o104 + o98 + o99
o106 = -5.0000000000000d0*o105*o52*o87
o107 = -1.92000000000000d2*o22*o46*o52*o73
o108 = 1.20000000000000d1*o1*o22*o46*o52*o73
o109 = 2.88000000000000d2*o22*o52*o78
o110 = -1.80000000000000d1*o1*o22*o52*o78
o111 = 5.5680000000000d3*o22*o46*o52*o78
o112 = -6.1800000000000d2*o1*o22*o46*o52*o78
o113 = 3.2000000000000d1*o44*o52*o73
o114 = -2.00000000000000d0*o1*o44*o52*o73
o115 = -9.2800000000000d2*o44*o52*o78
o116 = 1.03000000000000d2*o1*o44*o52*o78
o117 = o106 + o107 + o108 + o109 + o110 + o111 + o112 + o113 + o
& 114 + o115 + o116 + o94 + o95 + o96 + o97
o118 = 2.00000000000000d0*o117
o119 = 2.00000000000000d1*o6
o120 = 3.5000000000000d1*o102
o121 = 9.0000000000000d0 + o119 + o120
o122 = 1.87500000000000d1*o121*o50*o52*o91
o123 = o118 + o122
o124 = o21*o22
kxx(i,j,k) = mm*o16*o20*o25*o27*o43*o44 + o16*o27*o46*o49 - 1.00
& 000000000000d-1*o123*o16*o17*o22*o26*o51 - 1.00000000000000d-1*o
& 16*o17*o22*o26*o46*o51*o52*o71 + 3.00000000000000d0*mm*o14*o16*o
& 18*o20*o22*o4*o8 + 5.0000000000000d-2*mm*o16*o20*o21*o27*o45*o51
& *o72*o93
kxy(i,j,k) = mm*o16*o17*o20*o25*o26*o43*o44 + o16*o17*o26*o46*o4
& 9 - 5.0000000000000d-2*o123*o16*o18*o22*o51 + 5.0000000000000d-2
& *o123*o16*o22*o27*o51 - 5.0000000000000d-2*o16*o18*o22*o46*o51*o
& 52*o71 + 5.0000000000000d-2*o16*o22*o27*o46*o51*o52*o71 - 3.0000
& 0000000000d0*mm*o14*o16*o17*o20*o22*o26*o4*o8 + 5.0000000000000d
& -2*mm*o16*o17*o20*o21*o26*o45*o51*o72*o93
kxz(i,j,k) = mm*o124*o16*o20*o25*o26*o43*o45 - o16*o21*o26*o45*o
& 49 - 5.0000000000000d-2*o123*o16*o17*o21*o45*o51 + 5.00000000000
& 00d-2*o124*o16*o17*o45*o51*o52*o71 - 2.50000000000000d-2*mm*o16*
& o20*o22*o26*o51*o72*o93 + 2.50000000000000d-2*mm*o16*o20*o26*o46
& *o51*o72*o93
kyy(i,j,k) = mm*o16*o18*o20*o25*o43*o44 + o16*o18*o46*o49 + 1.00
& 000000000000d-1*o123*o16*o17*o22*o26*o51 + 1.00000000000000d-1*o
& 16*o17*o22*o26*o46*o51*o52*o71 + 3.00000000000000d0*mm*o14*o16*o
& 20*o22*o27*o4*o8 + 5.0000000000000d-2*mm*o16*o18*o20*o21*o45*o51
& *o72*o93
kyz(i,j,k) = mm*o124*o16*o17*o20*o25*o43*o45 - o16*o17*o21*o45*o
& 49 + 5.0000000000000d-2*o123*o16*o21*o26*o45*o51 - 5.00000000000
& 00d-2*o124*o16*o26*o45*o51*o52*o71 - 2.50000000000000d-2*mm*o16*
& o17*o20*o22*o51*o72*o93 + 2.50000000000000d-2*mm*o16*o17*o20*o46
& *o51*o72*o93
kzz(i,j,k) = mm*o16*o20*o22*o25*o43*o46 + o16*o22*o49 - 5.000000
& 0000000d-2*mm*o16*o20*o21*o45*o51*o72*o93
|
