32#include "implicit_f.inc"
36 INTEGER,
INTENT(IN) :: IN(4)
38 my_real,
INTENT(INOUT) :: nod_normal(3,*)
44 my_real x01,y01,z01, x02,y02,z02, x03,y03,z03, x04,y04,z04
45 my_real xn1,yn1,zn1, xn2,yn2,zn2, xn3,yn3,zn3, xn4,yn4,zn4
49 IF(in(4) /= in(3))
THEN
50 nod_normal(1,in(1))=nod_normal(1,in(1))+nx
51 nod_normal(2,in(1))=nod_normal(2,in(1))+ny
52 nod_normal(3,in(1))=nod_normal(3,in(1))+nz
54 nod_normal(1,in(2))=nod_normal(1,in(2))+nx
55 nod_normal(2,in(2))=nod_normal(2,in(2))+ny
56 nod_normal(3,in(2))=nod_normal(3,in(2))+nz
58 nod_normal(1,in(3))=nod_normal(1,in(3))+nx
59 nod_normal(2,in(3))=nod_normal(2,in(3))+ny
60 nod_normal(3,in(3))=nod_normal(3,in(3))+nz
62 nod_normal(1,in(4))=nod_normal(1,in(4))+nx
63 nod_normal(2,in(4))=nod_normal(2,in(4))+ny
64 nod_normal(3,in(4))=nod_normal(3,in(4))+nz
66 nod_normal(1,in(1))=nod_normal(1,in(1))+nx
67 nod_normal(2,in(1))=nod_normal(2,in(1))+ny
68 nod_normal(3,in(1))=nod_normal(3,in(1))+nz
70 nod_normal(1,in(2))=nod_normal(1,in(2))+nx
71 nod_normal(2,in(2))=nod_normal(2,in(2))+ny
72 nod_normal(3,in(2))=nod_normal(3,in(2))+nz
74 nod_normal(1,in(3))=nod_normal(1,in(3))+nx
75 nod_normal(2,in(3))=nod_normal(2,in(3))+ny
76 nod_normal(3,in(3))=nod_normal(3,in(3))+nz