29
30
31
32#include "implicit_f.inc"
33
34
35
36 INTEGER,INTENT(IN) :: IN(4)
38 my_real,
INTENT(INOUT) :: nod_normal(3,*)
39
40
41
42
43
44
45 IF(in(4) /= in(3))THEN
46 nod_normal(1,in(1))=nod_normal(1,in(1))+nx
47 nod_normal(2,in(1))=nod_normal(2,in(1))+ny
48 nod_normal(3,in(1))=nod_normal(3,in(1))+nz
49
50 nod_normal(1,in(2))=nod_normal(1,in(2))+nx
51 nod_normal(2,in(2))=nod_normal(2,in(2))+ny
52 nod_normal(3,in(2))=nod_normal(3,in(2))+nz
53
54 nod_normal(1,in(3))=nod_normal(1,in(3))+nx
55 nod_normal(2,in(3))=nod_normal(2,in(3))+ny
56 nod_normal(3,in(3))=nod_normal(3,in(3))+nz
57
58 nod_normal(1,in(4))=nod_normal(1,in(4))+nx
59 nod_normal(2,in(4))=nod_normal(2,in(4))+ny
60 nod_normal(3,in(4))=nod_normal(3,in(4))+nz
61 ELSE
62 nod_normal(1,in(1))=nod_normal(1,in(1))+nx
63 nod_normal(2,in(1))=nod_normal(2,in(1))+ny
64 nod_normal(3,in(1))=nod_normal(3,in(1))+nz
65
66 nod_normal(1,in(2))=nod_normal(1,in(2))+nx
67 nod_normal(2,in(2))=nod_normal(2,in(2))+ny
68 nod_normal(3,in(2))=nod_normal(3,in(2))+nz
69
70 nod_normal(1,in(3))=nod_normal(1,in(3))+nx
71 nod_normal(2,in(3))=nod_normal(2,in
72 nod_normal(3,in(3))=nod_normal(3,in(3))+nz
73 ENDIF
74
75 RETURN