33
34
35#include "implicit_f.inc"
36#include "comlock.inc"
37
38
39
40#include "mvsiz_p.inc"
41
42
43
44 INTEGER, INTENT(IN) :: NEL
45
47 . ms(*), rho(*),volu(*),msnf(*),off(*)
48 INTEGER NC1(*),NC2(*),NC3(*),NC4(*)
49
50
51
52
53
54
55 INTEGER I
56
58 . dmass(mvsiz)
59
60 DO i=1,nel
61 dmass(i)=fourth*rho(i)*volu(i)*off(i)
62 ENDDO
63
64#include "lockon.inc"
65 DO i=1,nel
66 ms(nc1(i))=ms(nc1(i)) + dmass(i)
67 msnf(nc1(i))=msnf(nc1(i)) + dmass(i)
68 ENDDO
69#include "lockoff.inc"
70
71#include "lockon.inc"
72 DO i=1,nel
73 msnf(nc2(i))=msnf(nc2(i)) + dmass(i)
74 ms(nc2(i))=ms(nc2(i)) + dmass(i)
75 ENDDO
76#include "lockoff.inc"
77
78#include "lockon.inc"
79 DO i=1,nel
80 ms(nc3(i))=ms(nc3(i)) + dmass(i)
81 msnf(nc3(i))=msnf(nc3(i)) + dmass(i)
82 ENDDO
83#include "lockoff.inc"
84
85#include "lockon.inc"
86 DO i=1,nel
87 ms(nc4(i))=ms(nc4(i)) + dmass(i)
88 msnf(nc4(i))=msnf(nc4(i)) + dmass(i)
89 ENDDO
90#include "lockoff.inc"
91
92
93 RETURN