34 . WEIGHT,STIFN,STIFR,MCP ,CONDN ,
35 . FTHE,INTBUF_TAB,ITHERM_FE,NODADT_THERM)
43#include "implicit_f.inc"
47 INTEGER IPARI(*), WEIGHT(*)
48 INTEGER ,
INTENT(IN) :: ITHERM_FE
49 INTEGER ,
INTENT(IN) :: NODADT_THERM
53 . ms(*),in(*),ar(*),vr(*),stifn(*),stifr(*),
54 . condn(*), mcp(*),fthe(*)
56 TYPE(intbuf_struct_) INTBUF_TAB
65 . k10, k11, k12, k13, k14, kfi, j10, j11, j12, j21, j22,
66 . jfi,nsn,nmn,ilev,nrts,nrtm
88 1 nmn ,a ,intbuf_tab%MSR ,ms ,weight ,
89 2 ar ,in ,stifn ,stifr ,ilev ,
90 3 mcp ,condn ,intbuf_tab%NMAS ,fthe,itherm_fe,nodadt_therm)
105 . A ,VR ,AR ,MS ,IN ,
106 . WEIGHT,STIFN,STIFR ,TAGNOD,INTBUF_TAB)
114#include "implicit_f.inc"
118 INTEGER IPARI(*), WEIGHT(*), TAGNOD(*)
122 . MS(*),IN(*),AR(*),(*),STIFN(*),STIFR(*)
124 TYPE(intbuf_struct_) INTBUF_TAB
138 1 nmn ,a ,intbuf_tab%MSR,ms ,weight ,
139 2 ar ,in ,stifn ,stifr ,ilev ,
140 3 intbuf_tab%NMAS,tagnod)
151 2 AR ,IN ,STIFN ,STIFR ,ILEV ,
152 3 MCP ,CONDN ,MMASS ,FTHE ,ITHERM_FE,NODADT_THERM)
156#include "implicit_f.inc"
160 INTEGER ,
INTENT(IN) :: ITHERM_FE
161 INTEGER ,
INTENT(IN) :: NODADT_THERM
166 . A(3,*) ,AR(3,*) ,MMASS(*) ,
167 . ms(*) ,in(*) ,stifn(*) ,stifr(*),
168 . mcp(*) ,condn(*),fthe(*)
172#include "com01_c.inc"
173#include "scr18_c.inc"
180 IF((ILEV == 1.OR.ilev == 3).AND.iroddl/=0)
THEN
181#include "vectorize.inc"
184 a(1,i)=a(1,i)*weight(i)
185 a(2,i)=a(2,i)*weight(i)
186 a(3,i)=a(3,i)*weight(i)
187 ar(1,i)=ar(1,i)*weight(i)
188 ar(2,i)=ar(2,i)*weight(i)
189 ar(3,i)=ar(3,i)*weight(i)
190 ms(i)=ms(i)*weight(i)
191 in(i)=in(i)*weight(i)
192 stifn(i)=stifn(i)*weight(i)
193 stifr(i)=stifr(i)*weight(i)
195 ELSEIF((ilev == 1.OR.ilev == 3).AND.iroddl == 0)
THEN
196#include "vectorize.inc"
199 a(1,i)=a(1,i)*weight(i)
200 a(2,i)=a(2,i)*weight(i)
201 a(3,i)=a(3,i)*weight(i)
202 ms(i)=ms(i)*weight(i)
203 stifn(i)=stifn(i)*weight(i)
205 ELSEIF((ilev == 0.OR.ilev == 2.OR.ilev == 4).AND.iroddl /= 0)
THEN
206#include "vectorize.inc"
209 a(1,i)=a(1,i)*weight(i)
210 a(2,i)=a(2,i)*weight(i)
211 a(3,i)=a(3,i)*weight(i)
213 ms(i)=ms(i)*weight(i)
214 stifn(i)=stifn(i)*weight(i)
215 ar(1,i)=ar(1,i)*weight(i)
216 ar(2,i)=ar(2,i)*weight(i)
217 ar(3,i)=ar(3,i)*weight(i)
218 in(i)=
max(em20,in(i))
219 in(i)=in(i)*weight(i)
220 stifr(i)=stifr(i)*weight(i)
222 ELSEIF((ilev == 0.OR.ilev == 2.OR.ilev == 4).AND.iroddl == 0)
THEN
223#include "vectorize.inc"
226 a(1,i)=a(1,i)*weight(i)
227 a(2,i)=a(2,i)*weight(i)
228 a(3,i)=a(3,i)*weight(i)
230 ms(i)=ms(i)*weight(i)
231 stifn(i)=stifn(i)*weight(i)
233 ELSEIF ((ilev == 10.OR.ilev == 11.OR.ilev == 12 .OR.
234 . ilev == 20.OR.ilev == 21.OR.ilev == 22).AND.
236#include "vectorize.inc"
239 a(1,i)=a(1,i)*weight(i)
240 a(2,i)=a(2,i)*weight(i)
241 a(3,i)=a(3,i)*weight(i)
242 ms(i)=ms(i)*weight(i)
243 in(i)=
max(em20,in(i))
244 in(i)=in(i)*weight(i)
245 ar(1,i)=ar(1,i)*weight(i)
246 ar(2,i)=ar(2,i)*weight(i)
247 ar(3,i)=ar(3,i)*weight(i)
248 stifn(i)=stifn(i)*weight(i)
249 stifr(i)=stifr(i)*weight(i)
251 ELSEIF ((ilev == 25 .or. ilev == 26) .AND. iroddl /= 0)
THEN
252#include "vectorize.inc"
255 a(1,i)=a(1,i)*weight(i)
256 a(2,i)=a(2,i)*weight(i)
257 a(3,i)=a(3,i)*weight(i)
258 ms(i)=ms(i)*weight(i)
259 in(i)=in(i)*weight(i)
260 ar(1,i)=ar(1,i)*weight(i)
261 ar(2,i)=ar(2,i)*weight(i)
262 ar(3,i)=ar(3,i)*weight(i)
263 stifn(i)=stifn(i)*weight(i)
266 ELSEIF ((ilev == 27) .AND. iroddl /= 0)
THEN
267#include "vectorize.inc"
270 a(1,i)=a(1,i)*weight(i)
271 a(2,i)=a(2,i)*weight(i)
272 a(3,i)=a(3,i)*weight(i)
274 ms(i)=ms(i)*weight(i)
275 in(i)=in(i)*weight(i)
276 ar(1,i)=ar(1,i)*weight(i)
277 ar(2,i)=ar(2,i)*weight(i)
278 ar(3,i)=ar(3,i)*weight(i)
279 stifn(i)=stifn(i)*weight(i)
280 stifr(i)=stifr(i)*weight(i)
282 ELSEIF ((ilev == 28) .AND. iroddl /= 0)
THEN
283#include "vectorize.inc"
286 a(1,i)=a(1,i)*weight(i)
287 a(2,i)=a(2,i)*weight(i)
288 a(3,i)=a(3,i)*weight(i)
290 ms(i)=ms(i)*weight(i)
291 in(i)=in(i)*weight(i)
292 ar(1,i)=ar(1,i)*weight(i)
293 ar(2,i)=ar(2,i)*weight(i)
294 ar(3,i)=ar(3,i)*weight(i)
295 stifn(i)=stifn(i)*weight(i)
296 stifr(i)=stifr(i)*weight(i)
298 ELSEIF ((ilev == 10 .OR. ilev == 11 .OR. ilev == 12 .OR.
299 . ilev == 20 .OR. ilev == 21 .OR. ilev == 22 .OR.
300 . ilev == 25 .or. ilev == 26 .or. ilev == 28 .OR.
301 . ilev == 27 ) .AND. iroddl == 0)
THEN
302#include "vectorize.inc"
305 a(1,i)=a(1,i)*weight(i)
306 a(2,i)=a(2,i)*weight(i)
307 a(3,i)=a(3,i)*weight(i)
308 ms(i)=ms(i)*weight(i)
309 stifn(i)=stifn(i)*weight(i)
311 ELSEIF(ilev == 30 .AND. iroddl /= 0)
THEN
312#include "vectorize.inc"
315 a(1,i)=a(1,i)*weight(i)
316 a(2,i)=a(2,i)*weight(i)
317 a(3,i)=a(3,i)*weight(i)
318 ar(1,i)=ar(1,i)*weight(i)
319 ar(2,i)=ar(2,i)*weight(i)
322 ms(i)=ms(i)*weight(i)
323 in(i)=
max(em20,in(i))
324 in(i)=in(i)*weight(i)
325 stifn(i)=stifn(i)*weight(i)
326 stifr(i)=stifr(i)*weight(i)
333 IF(itherm_fe > 0)
THEN
336 fthe(i) = fthe(i) * weight(i)
339 IF(nodadt_therm == 1)
THEN
342 condn(i) = condn(i) * weight(i)
355 2 AR ,IN ,STIFN,STIFR ,ILEV ,
360#include "implicit_f.inc"
365 . MSR(*), WEIGHT(*), TAGNOD(*)
368 . A(3,*),AR(3,*), MMASS(*),
369 . MS(*), IN(*),STIFN(*),(*)
373#include "com01_c.inc"
379 IF((ILEV == 1.OR.ilev == 3).AND.iroddl/=0)
THEN
380#include "vectorize.inc"
383 a(1,i)=a(1,i)*weight(i)
384 a(2,i)=a(2,i)*weight(i)
385 a(3,i)=a(3,i)*weight(i)
386 ar(1,i)=ar(1,i)*weight(i)
387 ar(2,i)=ar(2,i)*weight(i)
388 ar(3,i)=ar(3,i)*weight(i)
389 ms(i)=ms(i)*weight(i)
390 in(i)=in(i)*weight(i)
391 stifn(i)=stifn(i)*weight(i)
392 stifr(i)=stifr(i)*weight(i)
395 ELSEIF((ilev == 1.OR.ilev == 3).AND.iroddl == 0)
THEN
396#include "vectorize.inc"
399 a(1,i)=a(1,i)*weight(i)
400 a(2,i)=a(2,i)*weight(i)
401 a(3,i)=a(3,i)*weight(i)
402 ms(i)=ms(i)*weight(i)
403 stifn(i)=stifn(i)*weight(i)
406 ELSEIF((ilev == 0.OR.ilev == 2.OR.ilev == 4).AND.iroddl /= 0)
THEN
407#include "vectorize.inc"
410 a(1,i) =a(1,i) *weight(i)
411 a(2,i) =a(2,i) *weight(i)
412 a(3,i) =a(3,i) *weight(i)
413 ar(1,i)=ar(1,i)*weight(i)
414 ar(2,i)=ar(2,i)*weight(i)
415 ar(3,i)=ar(3,i)*weight(i)
418 in(i)=
max(em20,in(i))
419 in(i)=in(i)*weight(i)
420 stifn(i)=stifn(i)*weight(i)
421 stifr(i)=stifr(i)*weight(i)
424 ELSEIF((ilev == 0.OR.ilev == 2.OR.ilev == 4).AND.iroddl == 0)
THEN
425#include "vectorize.inc"
428 a(1,i)=a(1,i)*weight(i)
429 a(2,i)=a(2,i)*weight(i)
430 a(3,i)=a(3,i)*weight(i)
432 ms(i)=ms(i)*weight(i)
433 stifn(i)=stifn(i)*weight(i)
436 ELSEIF ((ilev == 10.OR.ilev == 11.OR.ilev == 12 .OR.
437 . ilev == 20.OR.ilev == 21.OR.ilev == 22) .AND.
439#include "vectorize.inc"
442 a(1,i)=a(1,i)*weight(i)
443 a(2,i)=a(2,i)*weight(i)
444 a(3,i)=a(3,i)*weight(i)
445 ms(i)=ms(i)*weight(i)
446 in(i)=
max(em20,in(i))
447 in(i)=in(i)*weight(i)
448 ar(1,i)=ar(1,i)*weight(i)
449 ar(2,i)=ar(2,i)*weight(i)
450 ar(3,i)=ar(3,i)*weight(i)
451 stifn(i)=stifn(i)*weight(i)
452 stifr(i)=stifr(i)*weight(i)
455 ELSEIF ((ilev == 25 .or. ilev == 26) .AND. iroddl /= 0)
THEN
456#include "vectorize.inc"
459 a(1,i)=a(1,i)*weight(i)
460 a(2,i)=a(2,i)*weight(i)
461 a(3,i)=a(3,i)*weight(i)
462 ms(i)=ms(i)*weight(i)
463 in(i)=in(i)*weight(i)
464 ar(1,i)=ar(1,i)*weight(i)
465 ar(2,i)=ar(2,i)*weight(i)
466 ar(3,i)=ar(3,i)*weight(i)
467 stifn(i)=stifn(i)*weight(i)
468 stifr(i)=stifr(i)*weight(i)
471 ELSEIF (ilev == 27 .AND. iroddl /= 0)
THEN
472#include "vectorize.inc"
475 a(1,i)=a(1,i)*weight(i)
476 a(2,i)=a(2,i)*weight(i)
477 a(3,i)=a(3,i)*weight(i)
479 ms(i)=ms(i)*weight(i)
480 in(i)=in(i)*weight(i)
481 ar(1,i)=ar(1,i)*weight(i)
482 ar(2,i)=ar(2,i)*weight(i)
483 ar(3,i)=ar(3,i)*weight(i)
484 stifn(i)=stifn(i)*weight(i)
485 stifr(i)=stifr(i)*weight(i)
488 ELSEIF (ilev == 28 .AND. iroddl /= 0)
THEN
489#include "vectorize.inc"
492 a(1,i)=a(1,i)*weight(i)
493 a(2,i)=a(2,i)*weight(i)
494 a(3,i)=a(3,i)*weight(i)
496 ms(i)=ms(i)*weight(i)
497 in(i)=in(i)*weight(i)
498 ar(1,i)=ar(1,i)*weight(i)
499 ar(2,i)=ar(2,i)*weight(i)
500 ar(3,i)=ar(3,i)*weight(i)
501 stifn(i)=stifn(i)*weight(i)
502 stifr(i)=stifr(i)*weight(i)
505 ELSEIF ((ilev == 10.OR.ilev == 11.OR.ilev == 12 .OR.
506 . ilev == 20.OR.ilev == 21.OR.ilev == 22 .OR.
507 . ilev == 25.or.ilev == 26.or.ilev == 28)
508 . .AND. iroddl == 0)
THEN
509#include "vectorize.inc"
512 a(1,i)=a(1,i)*weight(i)
513 a(2,i)=a(2,i)*weight(i)
514 a(3,i)=a(3,i)*weight(i)
516 stifn(i)=stifn(i)*weight(i)
519 ELSEIF(ilev == 30 .AND. iroddl /= 0)
THEN
520#include "vectorize.inc"
523 a(1,i)=a(1,i)*weight(i)
524 a(2,i)=a(2,i)*weight(i)
525 a(3,i)=a(3,i)*weight(i)
526 ar(1,i)=ar(1,i)*weight(i)
527 ar(2,i)=ar(2,i)*weight(i)
528 ar(3,i)=ar(3,i)*weight(i)
530 ms(i)=ms(i)*weight(i)
531 in(i)=
max(em20,in(i))
532 in(i)=in(i)*weight(i)
533 stifn(i)=stifn(i)*weight(i)
534 stifr(i)=stifr(i)*weight(i)