OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
assem_s20.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine assem_s20 (ixs, ixs20, nel, iddl, ndof, k_diag, k_lt, iadk, jdik, off, k11, k12, k13, k14, k15, k16, k17, k18, k19, k110, k111, k112, k113, k114, k115, k116, k117, k118, k119, k120, k22, k23, k24, k25, k26, k27, k28, k29, k210, k211, k212, k213, k214, k215, k216, k217, k218, k219, k220, k33, k34, k35, k36, k37, k38, k39, k310, k311, k312, k313, k314, k315, k316, k317, k318, k319, k320, k44, k45, k46, k47, k48, k49, k410, k411, k412, k413, k414, k415, k416, k417, k418, k419, k420, k55, k56, k57, k58, k59, k510, k511, k512, k513, k514, k515, k516, k517, k518, k519, k520, k66, k67, k68, k69, k610, k611, k612, k613, k614, k615, k616, k617, k618, k619, k620, k77, k78, k79, k710, k711, k712, k713, k714, k715, k716, k717, k718, k719, k720, k88, k89, k810, k811, k812, k813, k814, k815, k816, k817, k818, k819, k820, k99, k910, k911, k912, k913, k914, k915, k916, k917, k918, k919, k920, k1010, k1011, k1012, k1013, k1014, k1015, k1016, k1017, k1018, k1019, k1020, k1111, k1112, k1113, k1114, k1115, k1116, k1117, k1118, k1119, k1120, k1212, k1213, k1214, k1215, k1216, k1217, k1218, k1219, k1220, k1313, k1314, k1315, k1316, k1317, k1318, k1319, k1320, k1414, k1415, k1416, k1417, k1418, k1419, k1420, k1515, k1516, k1517, k1518, k1519, k1520, k1616, k1617, k1618, k1619, k1620, k1717, k1718, k1719, k1720, k1818, k1819, k1820, k1919, k1920, k2020)

Function/Subroutine Documentation

◆ assem_s20()

subroutine assem_s20 ( integer, dimension(nixs,*) ixs,
integer, dimension(12,*) ixs20,
integer nel,
integer, dimension(*) iddl,
integer, dimension(*) ndof,
k_diag,
k_lt,
integer, dimension(*) iadk,
integer, dimension(*) jdik,
off,
k11,
k12,
k13,
k14,
k15,
k16,
k17,
k18,
k19,
k110,
k111,
k112,
k113,
k114,
k115,
k116,
k117,
k118,
k119,
k120,
k22,
k23,
k24,
k25,
k26,
k27,
k28,
k29,
k210,
k211,
k212,
k213,
k214,
k215,
k216,
k217,
k218,
k219,
k220,
k33,
k34,
k35,
k36,
k37,
k38,
k39,
k310,
k311,
k312,
k313,
k314,
k315,
k316,
k317,
k318,
k319,
k320,
k44,
k45,
k46,
k47,
k48,
k49,
k410,
k411,
k412,
k413,
k414,
k415,
k416,
k417,
k418,
k419,
k420,
k55,
k56,
k57,
k58,
k59,
k510,
k511,
k512,
k513,
k514,
k515,
k516,
k517,
k518,
k519,
k520,
k66,
k67,
k68,
k69,
k610,
k611,
k612,
k613,
k614,
k615,
k616,
k617,
k618,
k619,
k620,
k77,
k78,
k79,
k710,
k711,
k712,
k713,
k714,
k715,
k716,
k717,
k718,
k719,
k720,
k88,
k89,
k810,
k811,
k812,
k813,
k814,
k815,
k816,
k817,
k818,
k819,
k820,
k99,
k910,
k911,
k912,
k913,
k914,
k915,
k916,
k917,
k918,
k919,
k920,
k1010,
k1011,
k1012,
k1013,
k1014,
k1015,
k1016,
k1017,
k1018,
k1019,
k1020,
k1111,
k1112,
k1113,
k1114,
k1115,
k1116,
k1117,
k1118,
k1119,
k1120,
k1212,
k1213,
k1214,
k1215,
k1216,
k1217,
k1218,
k1219,
k1220,
k1313,
k1314,
k1315,
k1316,
k1317,
k1318,
k1319,
k1320,
k1414,
k1415,
k1416,
k1417,
k1418,
k1419,
k1420,
k1515,
k1516,
k1517,
k1518,
k1519,
k1520,
k1616,
k1617,
k1618,
k1619,
k1620,
k1717,
k1718,
k1719,
k1720,
k1818,
k1819,
k1820,
k1919,
k1920,
k2020 )

Definition at line 33 of file assem_s20.F.

57 use element_mod , only : nixs
58C----6---------------------------------------------------------------7---------8
59C I m p l i c i t T y p e s
60C-----------------------------------------------
61#include "implicit_f.inc"
62C-----------------------------------------------
63C G l o b a l P a r a m e t e r s
64C-----------------------------------------------
65#include "mvsiz_p.inc"
66C-----------------------------------------------
67C C o m m o n B l o c k s
68C-----------------------------------------------
69C-----------------------------------------------------------------
70C D u m m y A r g u m e n t s
71C-----------------------------------------------
72 INTEGER IXS(NIXS,*),IXS20(12,*),NEL ,IDDL(*) ,NDOF(*) ,
73 . IADK(*) ,JDIK(*)
74C REAL
75 my_real
76 . k_diag(*) ,k_lt(*) ,off(*)
77 my_real
78 . k11(3,3,*) ,k12(3,3,*) ,k13(3,3,*) ,k14(3,3,*) ,k15(3,3,*) ,
79 . k16(3,3,*) ,k17(3,3,*) ,k18(3,3,*) ,k19(3,3,*) ,k110(3,3,*),
80 . k111(3,3,*),k112(3,3,*),k113(3,3,*),k114(3,3,*),k115(3,3,*),
81 . k116(3,3,*),k117(3,3,*),k118(3,3,*),k119(3,3,*),k120(3,3,*),
82 . k22(3,3,*) ,k23(3,3,*) ,k24(3,3,*) ,k25(3,3,*) ,k26(3,3,*) ,
83 . k27(3,3,*) ,k28(3,3,*) ,k29(3,3,*) ,k210(3,3,*),k211(3,3,*),
84 . k212(3,3,*),k213(3,3,*),k214(3,3,*),k215(3,3,*),k216(3,3,*),
85 . k217(3,3,*),k218(3,3,*),k219(3,3,*),k220(3,3,*),k33(3,3,*) ,
86 . k34(3,3,*) ,k35(3,3,*) ,k36(3,3,*) ,k37(3,3,*) ,k38(3,3,*) ,
87 . k39(3,3,*) ,k310(3,3,*),k311(3,3,*),k312(3,3,*),k313(3,3,*),
88 . k314(3,3,*),k315(3,3,*),k316(3,3,*),k317(3,3,*),k318(3,3,*),
89 . k319(3,3,*),k320(3,3,*),k44(3,3,*) ,k45(3,3,*) ,k46(3,3,*) ,
90 . k47(3,3,*) ,k48(3,3,*) ,k49(3,3,*) ,k410(3,3,*),k411(3,3,*),
91 . k412(3,3,*),k413(3,3,*),k414(3,3,*),k415(3,3,*),k416(3,3,*),
92 . k417(3,3,*),k418(3,3,*),k419(3,3,*),k420(3,3,*),k55(3,3,*) ,
93 . k56(3,3,*) ,k57(3,3,*) ,k58(3,3,*) ,k59(3,3,*) ,k510(3,3,*),
94 . k511(3,3,*),k512(3,3,*),k513(3,3,*),k514(3,3,*),k515(3,3,*),
95 . k516(3,3,*),k517(3,3,*),k518(3,3,*),k519(3,3,*),k520(3,3,*),
96 . k66(3,3,*) ,k67(3,3,*) ,k68(3,3,*) ,k69(3,3,*) ,k610(3,3,*),
97 . k611(3,3,*),k612(3,3,*),k613(3,3,*),k614(3,3,*),k615(3,3,*),
98 . k616(3,3,*),k617(3,3,*),k618(3,3,*),k619(3,3,*),k620(3,3,*),
99 . k77(3,3,*) ,k78(3,3,*) ,k79(3,3,*) ,k710(3,3,*),k711(3,3,*),
100 . k712(3,3,*),k713(3,3,*),k714(3,3,*),k715(3,3,*),k716(3,3,*),
101 . k717(3,3,*),k718(3,3,*),k719(3,3,*),k720(3,3,*),k88(3,3,*) ,
102 . k89(3,3,*) ,k810(3,3,*),k811(3,3,*),k812(3,3,*),k813(3,3,*),
103 . k814(3,3,*),k815(3,3,*),k816(3,3,*),k817(3,3,*),k818(3,3,*),
104 . k819(3,3,*),k820(3,3,*),k99(3,3,*) ,k910(3,3,*),k911(3,3,*),
105 . k912(3,3,*),k913(3,3,*),k914(3,3,*),k915(3,3,*),k916(3,3,*),
106 . k917(3,3,*),k918(3,3,*),k919(3,3,*),k920(3,3,*),k1010(3,3,*),
107 . k1011(3,3,*),k1012(3,3,*),k1013(3,3,*),k1014(3,3,*),k1015(3,3,*),
108 . k1016(3,3,*),k1017(3,3,*),k1018(3,3,*),k1019(3,3,*),k1020(3,3,*),
109 . k1111(3,3,*),k1112(3,3,*),k1113(3,3,*),k1114(3,3,*),k1115(3,3,*),
110 . k1116(3,3,*),k1117(3,3,*),k1118(3,3,*),k1119(3,3,*),k1120(3,3,*),
111 . k1212(3,3,*),k1213(3,3,*),k1214(3,3,*),k1215(3,3,*),k1216(3,3,*),
112 . k1217(3,3,*),k1218(3,3,*),k1219(3,3,*),k1220(3,3,*),k1313(3,3,*),
113 . k1314(3,3,*),k1315(3,3,*),k1316(3,3,*),k1317(3,3,*),k1318(3,3,*),
114 . k1319(3,3,*),k1320(3,3,*),k1414(3,3,*),k1415(3,3,*),k1416(3,3,*),
115 . k1417(3,3,*),k1418(3,3,*),k1419(3,3,*),k1420(3,3,*),k1515(3,3,*),
116 . k1516(3,3,*),k1517(3,3,*),k1518(3,3,*),k1519(3,3,*),k1520(3,3,*),
117 . k1616(3,3,*),k1617(3,3,*),k1618(3,3,*),k1619(3,3,*),k1620(3,3,*),
118 . k1717(3,3,*),k1718(3,3,*),k1719(3,3,*),k1720(3,3,*),k1818(3,3,*),
119 . k1819(3,3,*),k1820(3,3,*),k1919(3,3,*),k1920(3,3,*),k2020(3,3,*)
120C-----------------------------------------------
121C L o c a l V a r i a b l e s
122C-----------------------------------------------
123 INTEGER K,EP,NN(MVSIZ,20),J,I
124C----6---------------------------------------------------------------7---------8
125 DO k = 1,8
126 DO ep = 1,nel
127 nn(ep,k)=ixs(k+1,ep)
128 ENDDO
129 ENDDO
130
131 DO k = 9,20
132 DO ep = 1,nel
133 nn(ep,k)=ixs20(k-8,ep)
134 ENDDO
135 ENDDO
136
137C-------block diagonal k11,k22...k88-----
138 CALL assem_kii(nn(1,1),nel,iddl,iadk,k_diag,k_lt ,k11,3,off)
139 CALL assem_kii(nn(1,2),nel,iddl,iadk,k_diag,k_lt ,k22,3,off)
140 CALL assem_kii(nn(1,3),nel,iddl,iadk,k_diag,k_lt ,k33,3,off)
141 CALL assem_kii(nn(1,4),nel,iddl,iadk,k_diag,k_lt ,k44,3,off)
142 CALL assem_kii(nn(1,5),nel,iddl,iadk,k_diag,k_lt ,k55,3,off)
143 CALL assem_kii(nn(1,6),nel,iddl,iadk,k_diag,k_lt ,k66,3,off)
144 CALL assem_kii(nn(1,7),nel,iddl,iadk,k_diag,k_lt ,k77,3,off)
145 CALL assem_kii(nn(1,8),nel,iddl,iadk,k_diag,k_lt ,k88,3,off)
146 CALL assem_kii(nn(1,9),nel,iddl,iadk,k_diag,k_lt ,k99,3,off)
147 CALL assem_kii(nn(1,10),nel,iddl,iadk,k_diag,k_lt ,k1010,3,off)
148 CALL assem_kii(nn(1,11),nel,iddl,iadk,k_diag,k_lt ,k1111,3,off)
149 CALL assem_kii(nn(1,12),nel,iddl,iadk,k_diag,k_lt ,k1212,3,off)
150 CALL assem_kii(nn(1,13),nel,iddl,iadk,k_diag,k_lt ,k1313,3,off)
151 CALL assem_kii(nn(1,14),nel,iddl,iadk,k_diag,k_lt ,k1414,3,off)
152 CALL assem_kii(nn(1,15),nel,iddl,iadk,k_diag,k_lt ,k1515,3,off)
153 CALL assem_kii(nn(1,16),nel,iddl,iadk,k_diag,k_lt ,k1616,3,off)
154 CALL assem_kii(nn(1,17),nel,iddl,iadk,k_diag,k_lt ,k1717,3,off)
155 CALL assem_kii(nn(1,18),nel,iddl,iadk,k_diag,k_lt ,k1818,3,off)
156 CALL assem_kii(nn(1,19),nel,iddl,iadk,k_diag,k_lt ,k1919,3,off)
157 CALL assem_kii(nn(1,20),nel,iddl,iadk,k_diag,k_lt ,k2020,3,off)
158C-------non diag kij-----
159 CALL assem_kij(nn(1,1),nn(1,2),nel,iddl,iadk,jdik,
160 . k_diag,k_lt,k12,3,off)
161 CALL assem_kij(nn(1,1),nn(1,3),nel,iddl,iadk,jdik,
162 . k_diag,k_lt,k13,3,off)
163 CALL assem_kij(nn(1,1),nn(1,4),nel,iddl,iadk,jdik,
164 . k_diag,k_lt,k14,3,off)
165 CALL assem_kij(nn(1,1),nn(1,5),nel,iddl,iadk,jdik,
166 . k_diag,k_lt,k15,3,off)
167 CALL assem_kij(nn(1,1),nn(1,6),nel,iddl,iadk,jdik,
168 . k_diag,k_lt,k16,3,off)
169 CALL assem_kij(nn(1,1),nn(1,7),nel,iddl,iadk,jdik,
170 . k_diag,k_lt,k17,3,off)
171 CALL assem_kij(nn(1,1),nn(1,8),nel,iddl,iadk,jdik,
172 . k_diag,k_lt,k18,3,off)
173 CALL assem_kij(nn(1,1),nn(1,9),nel,iddl,iadk,jdik,
174 . k_diag,k_lt,k19,3,off)
175 CALL assem_kij(nn(1,1),nn(1,10),nel,iddl,iadk,jdik,
176 . k_diag,k_lt,k110,3,off)
177 CALL assem_kij(nn(1,1),nn(1,11),nel,iddl,iadk,jdik,
178 . k_diag,k_lt,k111,3,off)
179 CALL assem_kij(nn(1,1),nn(1,12),nel,iddl,iadk,jdik,
180 . k_diag,k_lt,k112,3,off)
181 CALL assem_kij(nn(1,1),nn(1,13),nel,iddl,iadk,jdik,
182 . k_diag,k_lt,k113,3,off)
183 CALL assem_kij(nn(1,1),nn(1,14),nel,iddl,iadk,jdik,
184 . k_diag,k_lt,k114,3,off)
185 CALL assem_kij(nn(1,1),nn(1,15),nel,iddl,iadk,jdik,
186 . k_diag,k_lt,k115,3,off)
187 CALL assem_kij(nn(1,1),nn(1,16),nel,iddl,iadk,jdik,
188 . k_diag,k_lt,k116,3,off)
189 CALL assem_kij(nn(1,1),nn(1,17),nel,iddl,iadk,jdik,
190 . k_diag,k_lt,k117,3,off)
191 CALL assem_kij(nn(1,1),nn(1,18),nel,iddl,iadk,jdik,
192 . k_diag,k_lt,k118,3,off)
193 CALL assem_kij(nn(1,1),nn(1,19),nel,iddl,iadk,jdik,
194 . k_diag,k_lt,k119,3,off)
195 CALL assem_kij(nn(1,1),nn(1,20),nel,iddl,iadk,jdik,
196 . k_diag,k_lt,k120,3,off)
197
198
199
200 CALL assem_kij(nn(1,2),nn(1,3),nel,iddl,iadk,jdik,
201 . k_diag,k_lt,k23,3,off)
202 CALL assem_kij(nn(1,2),nn(1,4),nel,iddl,iadk,jdik,
203 . k_diag,k_lt,k24,3,off)
204 CALL assem_kij(nn(1,2),nn(1,5),nel,iddl,iadk,jdik,
205 . k_diag,k_lt,k25,3,off)
206 CALL assem_kij(nn(1,2),nn(1,6),nel,iddl,iadk,jdik,
207 . k_diag,k_lt,k26,3,off)
208 CALL assem_kij(nn(1,2),nn(1,7),nel,iddl,iadk,jdik,
209 . k_diag,k_lt,k27,3,off)
210 CALL assem_kij(nn(1,2),nn(1,8),nel,iddl,iadk,jdik,
211 . k_diag,k_lt,k28,3,off)
212 CALL assem_kij(nn(1,2),nn(1,9),nel,iddl,iadk,jdik,
213 . k_diag,k_lt,k29,3,off)
214 CALL assem_kij(nn(1,2),nn(1,10),nel,iddl,iadk,jdik,
215 . k_diag,k_lt,k210,3,off)
216 CALL assem_kij(nn(1,2),nn(1,11),nel,iddl,iadk,jdik,
217 . k_diag,k_lt,k211,3,off)
218 CALL assem_kij(nn(1,2),nn(1,12),nel,iddl,iadk,jdik,
219 . k_diag,k_lt,k212,3,off)
220 CALL assem_kij(nn(1,2),nn(1,13),nel,iddl,iadk,jdik,
221 . k_diag,k_lt,k213,3,off)
222 CALL assem_kij(nn(1,2),nn(1,14),nel,iddl,iadk,jdik,
223 . k_diag,k_lt,k214,3,off)
224 CALL assem_kij(nn(1,2),nn(1,15),nel,iddl,iadk,jdik,
225 . k_diag,k_lt,k215,3,off)
226 CALL assem_kij(nn(1,2),nn(1,16),nel,iddl,iadk,jdik,
227 . k_diag,k_lt,k216,3,off)
228 CALL assem_kij(nn(1,2),nn(1,17),nel,iddl,iadk,jdik,
229 . k_diag,k_lt,k217,3,off)
230 CALL assem_kij(nn(1,2),nn(1,18),nel,iddl,iadk,jdik,
231 . k_diag,k_lt,k218,3,off)
232 CALL assem_kij(nn(1,2),nn(1,19),nel,iddl,iadk,jdik,
233 . k_diag,k_lt,k219,3,off)
234 CALL assem_kij(nn(1,2),nn(1,20),nel,iddl,iadk,jdik,
235 . k_diag,k_lt,k220,3,off)
236
237
238 CALL assem_kij(nn(1,3),nn(1,4),nel,iddl,iadk,jdik,
239 . k_diag,k_lt,k34,3,off)
240 CALL assem_kij(nn(1,3),nn(1,5),nel,iddl,iadk,jdik,
241 . k_diag,k_lt,k35,3,off)
242 CALL assem_kij(nn(1,3),nn(1,6),nel,iddl,iadk,jdik,
243 . k_diag,k_lt,k36,3,off)
244 CALL assem_kij(nn(1,3),nn(1,7),nel,iddl,iadk,jdik,
245 . k_diag,k_lt,k37,3,off)
246 CALL assem_kij(nn(1,3),nn(1,8),nel,iddl,iadk,jdik,
247 . k_diag,k_lt,k38,3,off)
248 CALL assem_kij(nn(1,3),nn(1,9),nel,iddl,iadk,jdik,
249 . k_diag,k_lt,k39,3,off)
250 CALL assem_kij(nn(1,3),nn(1,10),nel,iddl,iadk,jdik,
251 . k_diag,k_lt,k310,3,off)
252 CALL assem_kij(nn(1,3),nn(1,11),nel,iddl,iadk,jdik,
253 . k_diag,k_lt,k311,3,off)
254 CALL assem_kij(nn(1,3),nn(1,12),nel,iddl,iadk,jdik,
255 . k_diag,k_lt,k312,3,off)
256 CALL assem_kij(nn(1,3),nn(1,13),nel,iddl,iadk,jdik,
257 . k_diag,k_lt,k313,3,off)
258 CALL assem_kij(nn(1,3),nn(1,14),nel,iddl,iadk,jdik,
259 . k_diag,k_lt,k314,3,off)
260 CALL assem_kij(nn(1,3),nn(1,15),nel,iddl,iadk,jdik,
261 . k_diag,k_lt,k315,3,off)
262 CALL assem_kij(nn(1,3),nn(1,16),nel,iddl,iadk,jdik,
263 . k_diag,k_lt,k316,3,off)
264 CALL assem_kij(nn(1,3),nn(1,17),nel,iddl,iadk,jdik,
265 . k_diag,k_lt,k317,3,off)
266 CALL assem_kij(nn(1,3),nn(1,18),nel,iddl,iadk,jdik,
267 . k_diag,k_lt,k318,3,off)
268 CALL assem_kij(nn(1,3),nn(1,19),nel,iddl,iadk,jdik,
269 . k_diag,k_lt,k319,3,off)
270 CALL assem_kij(nn(1,3),nn(1,20),nel,iddl,iadk,jdik,
271 . k_diag,k_lt,k320,3,off)
272
273
274
275 CALL assem_kij(nn(1,4),nn(1,5),nel,iddl,iadk,jdik,
276 . k_diag,k_lt,k45,3,off)
277 CALL assem_kij(nn(1,4),nn(1,6),nel,iddl,iadk,jdik,
278 . k_diag,k_lt,k46,3,off)
279 CALL assem_kij(nn(1,4),nn(1,7),nel,iddl,iadk,jdik,
280 . k_diag,k_lt,k47,3,off)
281 CALL assem_kij(nn(1,4),nn(1,8),nel,iddl,iadk,jdik,
282 . k_diag,k_lt,k48,3,off)
283 CALL assem_kij(nn(1,4),nn(1,9),nel,iddl,iadk,jdik,
284 . k_diag,k_lt,k49,3,off)
285 CALL assem_kij(nn(1,4),nn(1,10),nel,iddl,iadk,jdik,
286 . k_diag,k_lt,k410,3,off)
287 CALL assem_kij(nn(1,4),nn(1,11),nel,iddl,iadk,jdik,
288 . k_diag,k_lt,k411,3,off)
289 CALL assem_kij(nn(1,4),nn(1,12),nel,iddl,iadk,jdik,
290 . k_diag,k_lt,k412,3,off)
291 CALL assem_kij(nn(1,4),nn(1,13),nel,iddl,iadk,jdik,
292 . k_diag,k_lt,k413,3,off)
293 CALL assem_kij(nn(1,4),nn(1,14),nel,iddl,iadk,jdik,
294 . k_diag,k_lt,k414,3,off)
295 CALL assem_kij(nn(1,4),nn(1,15),nel,iddl,iadk,jdik,
296 . k_diag,k_lt,k415,3,off)
297 CALL assem_kij(nn(1,4),nn(1,16),nel,iddl,iadk,jdik,
298 . k_diag,k_lt,k416,3,off)
299 CALL assem_kij(nn(1,4),nn(1,17),nel,iddl,iadk,jdik,
300 . k_diag,k_lt,k417,3,off)
301 CALL assem_kij(nn(1,4),nn(1,18),nel,iddl,iadk,jdik,
302 . k_diag,k_lt,k418,3,off)
303 CALL assem_kij(nn(1,4),nn(1,19),nel,iddl,iadk,jdik,
304 . k_diag,k_lt,k419,3,off)
305 CALL assem_kij(nn(1,4),nn(1,20),nel,iddl,iadk,jdik,
306 . k_diag,k_lt,k420,3,off)
307
308
309
310 CALL assem_kij(nn(1,5),nn(1,6),nel,iddl,iadk,jdik,
311 . k_diag,k_lt,k56,3,off)
312 CALL assem_kij(nn(1,5),nn(1,7),nel,iddl,iadk,jdik,
313 . k_diag,k_lt,k57,3,off)
314 CALL assem_kij(nn(1,5),nn(1,8),nel,iddl,iadk,jdik,
315 . k_diag,k_lt,k58,3,off)
316 CALL assem_kij(nn(1,5),nn(1,9),nel,iddl,iadk,jdik,
317 . k_diag,k_lt,k59,3,off)
318 CALL assem_kij(nn(1,5),nn(1,10),nel,iddl,iadk,jdik,
319 . k_diag,k_lt,k510,3,off)
320 CALL assem_kij(nn(1,5),nn(1,11),nel,iddl,iadk,jdik,
321 . k_diag,k_lt,k511,3,off)
322 CALL assem_kij(nn(1,5),nn(1,12),nel,iddl,iadk,jdik,
323 . k_diag,k_lt,k512,3,off)
324 CALL assem_kij(nn(1,5),nn(1,13),nel,iddl,iadk,jdik,
325 . k_diag,k_lt,k513,3,off)
326 CALL assem_kij(nn(1,5),nn(1,14),nel,iddl,iadk,jdik,
327 . k_diag,k_lt,k514,3,off)
328 CALL assem_kij(nn(1,5),nn(1,15),nel,iddl,iadk,jdik,
329 . k_diag,k_lt,k515,3,off)
330 CALL assem_kij(nn(1,5),nn(1,16),nel,iddl,iadk,jdik,
331 . k_diag,k_lt,k516,3,off)
332 CALL assem_kij(nn(1,5),nn(1,17),nel,iddl,iadk,jdik,
333 . k_diag,k_lt,k517,3,off)
334 CALL assem_kij(nn(1,5),nn(1,18),nel,iddl,iadk,jdik,
335 . k_diag,k_lt,k518,3,off)
336 CALL assem_kij(nn(1,5),nn(1,19),nel,iddl,iadk,jdik,
337 . k_diag,k_lt,k519,3,off)
338 CALL assem_kij(nn(1,5),nn(1,20),nel,iddl,iadk,jdik,
339 . k_diag,k_lt,k520,3,off)
340
341
342
343 CALL assem_kij(nn(1,6),nn(1,7),nel,iddl,iadk,jdik,
344 . k_diag,k_lt,k67,3,off)
345 CALL assem_kij(nn(1,6),nn(1,8),nel,iddl,iadk,jdik,
346 . k_diag,k_lt,k68,3,off)
347 CALL assem_kij(nn(1,6),nn(1,9),nel,iddl,iadk,jdik,
348 . k_diag,k_lt,k69,3,off)
349 CALL assem_kij(nn(1,6),nn(1,10),nel,iddl,iadk,jdik,
350 . k_diag,k_lt,k610,3,off)
351 CALL assem_kij(nn(1,6),nn(1,11),nel,iddl,iadk,jdik,
352 . k_diag,k_lt,k611,3,off)
353 CALL assem_kij(nn(1,6),nn(1,12),nel,iddl,iadk,jdik,
354 . k_diag,k_lt,k612,3,off)
355 CALL assem_kij(nn(1,6),nn(1,13),nel,iddl,iadk,jdik,
356 . k_diag,k_lt,k613,3,off)
357 CALL assem_kij(nn(1,6),nn(1,14),nel,iddl,iadk,jdik,
358 . k_diag,k_lt,k614,3,off)
359 CALL assem_kij(nn(1,6),nn(1,15),nel,iddl,iadk,jdik,
360 . k_diag,k_lt,k615,3,off)
361 CALL assem_kij(nn(1,6),nn(1,16),nel,iddl,iadk,jdik,
362 . k_diag,k_lt,k616,3,off)
363 CALL assem_kij(nn(1,6),nn(1,17),nel,iddl,iadk,jdik,
364 . k_diag,k_lt,k617,3,off)
365 CALL assem_kij(nn(1,6),nn(1,18),nel,iddl,iadk,jdik,
366 . k_diag,k_lt,k618,3,off)
367 CALL assem_kij(nn(1,6),nn(1,19),nel,iddl,iadk,jdik,
368 . k_diag,k_lt,k619,3,off)
369 CALL assem_kij(nn(1,6),nn(1,20),nel,iddl,iadk,jdik,
370 . k_diag,k_lt,k620,3,off)
371
372
373
374
375 CALL assem_kij(nn(1,7),nn(1,8),nel,iddl,iadk,jdik,
376 . k_diag, k_lt,k78,3,off)
377 CALL assem_kij(nn(1,7),nn(1,9),nel,iddl,iadk,jdik,
378 . k_diag,k_lt,k79,3,off)
379 CALL assem_kij(nn(1,7),nn(1,10),nel,iddl,iadk,jdik,
380 . k_diag,k_lt,k710,3,off)
381 CALL assem_kij(nn(1,7),nn(1,11),nel,iddl,iadk,jdik,
382 . k_diag,k_lt,k711,3,off)
383 CALL assem_kij(nn(1,7),nn(1,12),nel,iddl,iadk,jdik,
384 . k_diag,k_lt,k712,3,off)
385 CALL assem_kij(nn(1,7),nn(1,13),nel,iddl,iadk,jdik,
386 . k_diag,k_lt,k713,3,off)
387 CALL assem_kij(nn(1,7),nn(1,14),nel,iddl,iadk,jdik,
388 . k_diag,k_lt,k714,3,off)
389 CALL assem_kij(nn(1,7),nn(1,15),nel,iddl,iadk,jdik,
390 . k_diag,k_lt,k715,3,off)
391 CALL assem_kij(nn(1,7),nn(1,16),nel,iddl,iadk,jdik,
392 . k_diag,k_lt,k716,3,off)
393 CALL assem_kij(nn(1,7),nn(1,17),nel,iddl,iadk,jdik,
394 . k_diag,k_lt,k717,3,off)
395 CALL assem_kij(nn(1,7),nn(1,18),nel,iddl,iadk,jdik,
396 . k_diag,k_lt,k718,3,off)
397 CALL assem_kij(nn(1,7),nn(1,19),nel,iddl,iadk,jdik,
398 . k_diag,k_lt,k719,3,off)
399 CALL assem_kij(nn(1,7),nn(1,20),nel,iddl,iadk,jdik,
400 . k_diag,k_lt,k720,3,off)
401
402
403 CALL assem_kij(nn(1,8),nn(1,9),nel,iddl,iadk,jdik,
404 . k_diag,k_lt,k89,3,off)
405 CALL assem_kij(nn(1,8),nn(1,10),nel,iddl,iadk,jdik,
406 . k_diag,k_lt,k810,3,off)
407 CALL assem_kij(nn(1,8),nn(1,11),nel,iddl,iadk,jdik,
408 . k_diag,k_lt,k811,3,off)
409 CALL assem_kij(nn(1,8),nn(1,12),nel,iddl,iadk,jdik,
410 . k_diag,k_lt,k812,3,off)
411 CALL assem_kij(nn(1,8),nn(1,13),nel,iddl,iadk,jdik,
412 . k_diag,k_lt,k813,3,off)
413 CALL assem_kij(nn(1,8),nn(1,14),nel,iddl,iadk,jdik,
414 . k_diag,k_lt,k814,3,off)
415 CALL assem_kij(nn(1,8),nn(1,15),nel,iddl,iadk,jdik,
416 . k_diag,k_lt,k815,3,off)
417 CALL assem_kij(nn(1,8),nn(1,16),nel,iddl,iadk,jdik,
418 . k_diag,k_lt,k816,3,off)
419 CALL assem_kij(nn(1,8),nn(1,17),nel,iddl,iadk,jdik,
420 . k_diag,k_lt,k817,3,off)
421 CALL assem_kij(nn(1,8),nn(1,18),nel,iddl,iadk,jdik,
422 . k_diag,k_lt,k818,3,off)
423 CALL assem_kij(nn(1,8),nn(1,19),nel,iddl,iadk,jdik,
424 . k_diag,k_lt,k819,3,off)
425 CALL assem_kij(nn(1,8),nn(1,20),nel,iddl,iadk,jdik,
426 . k_diag,k_lt,k820,3,off)
427
428
429 CALL assem_kij(nn(1,9),nn(1,10),nel,iddl,iadk,jdik,
430 . k_diag,k_lt,k910,3,off)
431 CALL assem_kij(nn(1,9),nn(1,11),nel,iddl,iadk,jdik,
432 . k_diag,k_lt,k911,3,off)
433 CALL assem_kij(nn(1,9),nn(1,12),nel,iddl,iadk,jdik,
434 . k_diag,k_lt,k912,3,off)
435 CALL assem_kij(nn(1,9),nn(1,13),nel,iddl,iadk,jdik,
436 . k_diag,k_lt,k913,3,off)
437 CALL assem_kij(nn(1,9),nn(1,14),nel,iddl,iadk,jdik,
438 . k_diag,k_lt,k914,3,off)
439 CALL assem_kij(nn(1,9),nn(1,15),nel,iddl,iadk,jdik,
440 . k_diag,k_lt,k915,3,off)
441 CALL assem_kij(nn(1,9),nn(1,16),nel,iddl,iadk,jdik,
442 . k_diag,k_lt,k916,3,off)
443 CALL assem_kij(nn(1,9),nn(1,17),nel,iddl,iadk,jdik,
444 . k_diag,k_lt,k917,3,off)
445 CALL assem_kij(nn(1,9),nn(1,18),nel,iddl,iadk,jdik,
446 . k_diag,k_lt,k918,3,off)
447 CALL assem_kij(nn(1,9),nn(1,19),nel,iddl,iadk,jdik,
448 . k_diag,k_lt,k919,3,off)
449 CALL assem_kij(nn(1,9),nn(1,20),nel,iddl,iadk,jdik,
450 . k_diag,k_lt,k920,3,off)
451
452
453 CALL assem_kij(nn(1,10),nn(1,11),nel,iddl,iadk,jdik,
454 . k_diag,k_lt,k1011,3,off)
455 CALL assem_kij(nn(1,10),nn(1,12),nel,iddl,iadk,jdik,
456 . k_diag,k_lt,k1012,3,off)
457 CALL assem_kij(nn(1,10),nn(1,13),nel,iddl,iadk,jdik,
458 . k_diag,k_lt,k1013,3,off)
459 CALL assem_kij(nn(1,10),nn(1,14),nel,iddl,iadk,jdik,
460 . k_diag,k_lt,k1014,3,off)
461 CALL assem_kij(nn(1,10),nn(1,15),nel,iddl,iadk,jdik,
462 . k_diag,k_lt,k1015,3,off)
463 CALL assem_kij(nn(1,10),nn(1,16),nel,iddl,iadk,jdik,
464 . k_diag,k_lt,k1016,3,off)
465 CALL assem_kij(nn(1,10),nn(1,17),nel,iddl,iadk,jdik,
466 . k_diag,k_lt,k1017,3,off)
467 CALL assem_kij(nn(1,10),nn(1,18),nel,iddl,iadk,jdik,
468 . k_diag,k_lt,k1018,3,off)
469 CALL assem_kij(nn(1,10),nn(1,19),nel,iddl,iadk,jdik,
470 . k_diag,k_lt,k1019,3,off)
471 CALL assem_kij(nn(1,10),nn(1,20),nel,iddl,iadk,jdik,
472 . k_diag,k_lt,k1020,3,off)
473
474
475 CALL assem_kij(nn(1,11),nn(1,12),nel,iddl,iadk,jdik,
476 . k_diag,k_lt,k1112,3,off)
477 CALL assem_kij(nn(1,11),nn(1,13),nel,iddl,iadk,jdik,
478 . k_diag,k_lt,k1113,3,off)
479 CALL assem_kij(nn(1,11),nn(1,14),nel,iddl,iadk,jdik,
480 . k_diag,k_lt,k1114,3,off)
481 CALL assem_kij(nn(1,11),nn(1,15),nel,iddl,iadk,jdik,
482 . k_diag,k_lt,k1115,3,off)
483 CALL assem_kij(nn(1,11),nn(1,16),nel,iddl,iadk,jdik,
484 . k_diag,k_lt,k1116,3,off)
485 CALL assem_kij(nn(1,11),nn(1,17),nel,iddl,iadk,jdik,
486 . k_diag,k_lt,k1117,3,off)
487 CALL assem_kij(nn(1,11),nn(1,18),nel,iddl,iadk,jdik,
488 . k_diag,k_lt,k1118,3,off)
489 CALL assem_kij(nn(1,11),nn(1,19),nel,iddl,iadk,jdik,
490 . k_diag,k_lt,k1119,3,off)
491 CALL assem_kij(nn(1,11),nn(1,20),nel,iddl,iadk,jdik,
492 . k_diag,k_lt,k1120,3,off)
493
494
495
496 CALL assem_kij(nn(1,12),nn(1,13),nel,iddl,iadk,jdik,
497 . k_diag,k_lt,k1213,3,off)
498 CALL assem_kij(nn(1,12),nn(1,14),nel,iddl,iadk,jdik,
499 . k_diag,k_lt,k1214,3,off)
500 CALL assem_kij(nn(1,12),nn(1,15),nel,iddl,iadk,jdik,
501 . k_diag,k_lt,k1215,3,off)
502 CALL assem_kij(nn(1,12),nn(1,16),nel,iddl,iadk,jdik,
503 . k_diag,k_lt,k1216,3,off)
504 CALL assem_kij(nn(1,12),nn(1,17),nel,iddl,iadk,jdik,
505 . k_diag,k_lt,k1217,3,off)
506 CALL assem_kij(nn(1,12),nn(1,18),nel,iddl,iadk,jdik,
507 . k_diag,k_lt,k1218,3,off)
508 CALL assem_kij(nn(1,12),nn(1,19),nel,iddl,iadk,jdik,
509 . k_diag,k_lt,k1219,3,off)
510 CALL assem_kij(nn(1,12),nn(1,20),nel,iddl,iadk,jdik,
511 . k_diag,k_lt,k1220,3,off)
512
513
514 CALL assem_kij(nn(1,13),nn(1,14),nel,iddl,iadk,jdik,
515 . k_diag,k_lt,k1314,3,off)
516 CALL assem_kij(nn(1,13),nn(1,15),nel,iddl,iadk,jdik,
517 . k_diag,k_lt,k1315,3,off)
518 CALL assem_kij(nn(1,13),nn(1,16),nel,iddl,iadk,jdik,
519 . k_diag,k_lt,k1316,3,off)
520 CALL assem_kij(nn(1,13),nn(1,17),nel,iddl,iadk,jdik,
521 . k_diag,k_lt,k1317,3,off)
522 CALL assem_kij(nn(1,13),nn(1,18),nel,iddl,iadk,jdik,
523 . k_diag,k_lt,k1318,3,off)
524 CALL assem_kij(nn(1,13),nn(1,19),nel,iddl,iadk,jdik,
525 . k_diag,k_lt,k1319,3,off)
526 CALL assem_kij(nn(1,13),nn(1,20),nel,iddl,iadk,jdik,
527 . k_diag,k_lt,k1320,3,off)
528
529
530 CALL assem_kij(nn(1,14),nn(1,15),nel,iddl,iadk,jdik,
531 . k_diag,k_lt,k1415,3,off)
532 CALL assem_kij(nn(1,14),nn(1,16),nel,iddl,iadk,jdik,
533 . k_diag,k_lt,k1416,3,off)
534 CALL assem_kij(nn(1,14),nn(1,17),nel,iddl,iadk,jdik,
535 . k_diag,k_lt,k1417,3,off)
536 CALL assem_kij(nn(1,14),nn(1,18),nel,iddl,iadk,jdik,
537 . k_diag,k_lt,k1418,3,off)
538 CALL assem_kij(nn(1,14),nn(1,19),nel,iddl,iadk,jdik,
539 . k_diag,k_lt,k1419,3,off)
540 CALL assem_kij(nn(1,14),nn(1,20),nel,iddl,iadk,jdik,
541 . k_diag,k_lt,k1420,3,off)
542
543
544 CALL assem_kij(nn(1,15),nn(1,16),nel,iddl,iadk,jdik,
545 . k_diag,k_lt,k1516,3,off)
546 CALL assem_kij(nn(1,15),nn(1,17),nel,iddl,iadk,jdik,
547 . k_diag,k_lt,k1517,3,off)
548 CALL assem_kij(nn(1,15),nn(1,18),nel,iddl,iadk,jdik,
549 . k_diag,k_lt,k1518,3,off)
550 CALL assem_kij(nn(1,15),nn(1,19),nel,iddl,iadk,jdik,
551 . k_diag,k_lt,k1519,3,off)
552 CALL assem_kij(nn(1,15),nn(1,20),nel,iddl,iadk,jdik,
553 . k_diag,k_lt,k1520,3,off)
554
555
556 CALL assem_kij(nn(1,16),nn(1,17),nel,iddl,iadk,jdik,
557 . k_diag,k_lt,k1617,3,off)
558 CALL assem_kij(nn(1,16),nn(1,18),nel,iddl,iadk,jdik,
559 . k_diag,k_lt,k1618,3,off)
560 CALL assem_kij(nn(1,16),nn(1,19),nel,iddl,iadk,jdik,
561 . k_diag,k_lt,k1619,3,off)
562 CALL assem_kij(nn(1,16),nn(1,20),nel,iddl,iadk,jdik,
563 . k_diag,k_lt,k1620,3,off)
564
565
566 CALL assem_kij(nn(1,17),nn(1,18),nel,iddl,iadk,jdik,
567 . k_diag,k_lt,k1718,3,off)
568 CALL assem_kij(nn(1,17),nn(1,19),nel,iddl,iadk,jdik,
569 . k_diag,k_lt,k1719,3,off)
570 CALL assem_kij(nn(1,17),nn(1,20),nel,iddl,iadk,jdik,
571 . k_diag,k_lt,k1720,3,off)
572
573
574 CALL assem_kij(nn(1,18),nn(1,19),nel,iddl,iadk,jdik,
575 . k_diag,k_lt,k1819,3,off)
576 CALL assem_kij(nn(1,18),nn(1,20),nel,iddl,iadk,jdik,
577 . k_diag,k_lt,k1820,3,off)
578
579
580 CALL assem_kij(nn(1,19),nn(1,20),nel,iddl,iadk,jdik,
581 . k_diag,k_lt,k1920,3,off)
582
583C
584C----6---------------------------------------------------------------7---------8
585 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
assem_kii
subroutine assem_kii(ni, nel, iddl, iadk, k_diag, k_lt, kii, nd, off)
Definition imp_glob_k.F:964
assem_kij
subroutine assem_kij(ni, nj, nel, iddl, iadk, jdik, k_diag, k_lt, kij, nd, off)
Definition imp_glob_k.F:1043