32#include "implicit_f.inc"
36 INTEGER,
INTENT(IN) :: IN(4)
38 my_real,
INTENT(INOUT) :: nod_normal(3,*)
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
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
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
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
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
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
70 nod_normal(1,in(3))=nod_normal(1,in(3))+nx
71 nod_normal(2,in(3))=nod_normal(2,in(3))+ny
72 nod_normal(3,in(3))=nod_normal(3,in(3))+nz