36
37
38
39#include "implicit_f.inc"
40
41
42
43 INTEGER NEL
44
46 . bxy1(*),bxy2(*),bxy3(*),bxy4(*),
47 . bxy5(*),bxy6(*),bxy7(*),bxy8(*),
48 . byx1(*),byx2(*),byx3(*),byx4(*),
49 . byx5(*),byx6(*),byx7(*),byx8(*),
50 . bxz1(*),bxz2(*),bxz3(*),bxz4(*),
51 . bxz5(*),bxz6(*),bxz7(*),bxz8(*),
52 . bzx1(*),bzx2(*),bzx3(*),bzx4(*),
53 . bzx5(*),bzx6(*),bzx7(*),bzx8(*),
54 . byz1(*),byz2(*),byz3(*),byz4(*),
55 . byz5(*),byz6(*),byz7(*),byz8(*),
56 . bzy1(*),bzy2(*),bzy3(*),bzy4(*),
57 . bzy5(*),bzy6(*),bzy7(*),bzy8(*),
58 . pij(nel,*)
59
60
61
62 INTEGER I,J,K
63
64 DO i=1,nel
65
66 bxy1(i) = pij(i,25)
67 byx1(i) = pij(i,26)
68 bxz1(i) = pij(i,27)
69 bzx1(i) = pij(i,28)
70 byz1(i) = pij(i,29)
71 bzy1(i) = pij(i,30)
72
73 bxy2(i) = pij(i,31)
74 byx2(i) = pij(i,32)
75 bxz2(i) = pij(i,33)
76 bzx2(i) = pij(i,34)
77 byz2(i) = pij(i,35)
78 bzy2(i) = pij(i,36)
79
80 bxy3(i) = pij(i,37)
81 byx3(i) = pij(i,38)
82 bxz3(i) = pij(i,39)
83 bzx3(i) = pij(i,40)
84 byz3(i) = pij(i,41)
85 bzy3(i) = pij(i,42)
86
87 bxy4(i) = pij(i,43)
88 byx4(i) = pij(i,44)
89 bxz4(i) = pij(i,45)
90 bzx4(i) = pij(i,46)
91 byz4(i) = pij(i,47)
92 bzy4(i) = pij(i,48)
93
94 bxy5(i) = pij(i,49)
95 byx5(i) = pij(i,50)
96 bxz5(i) = pij(i,51)
97 bzx5(i) = pij(i,52)
98 byz5(i) = pij(i,53)
99 bzy5(i) = pij(i,54)
100
101 bxy6(i) = pij(i,55)
102 byx6(i) = pij(i,56)
103 bxz6(i) = pij(i,57)
104 bzx6(i) = pij(i,58)
105 byz6(i) = pij(i,59)
106 bzy6(i) = pij(i,60)
107
108 bxy7(i) = pij(i,61)
109 byx7(i) = pij(i,62)
110 bxz7(i) = pij(i,63)
111 bzx7(i) = pij(i,64)
112 byz7(i) = pij(i,65)
113 bzy7(i) = pij(i,66)
114
115 bxy8(i) = pij(i,67)
116 byx8(i) = pij(i,68)
117 bxz8(i) = pij(i,69)
118 bzx8(i) = pij(i,70)
119 byz8(i) = pij(i,71)
120 bzy8(i) = pij(i,72)
121 END DO
122
123 RETURN