31
32
33
34#include "implicit_f.inc"
35
36
37
38 INTEGER NSN,FLAG,TAGKINE(*)
39 INTEGER NSV(*)
41 . a(3,*), frigap(*)
42
43
44
45 INTEGER II, I
47 . xx,yy,zz,
48 . s11,s21,s31,s12,s22,s32,s13,s23,s33
49
50
51
52
53 s11=frigap(5)
54 s21=frigap(6)
55 s31=frigap(7)
56 s12=frigap(8)
57 s22=frigap(9)
58 s32=frigap(10)
59 s13=frigap(12)
60 s23=frigap(13)
61 s33=frigap(14)
62 IF(flag==1)THEN
63 DO ii=1,nsn
64 IF(tagkine(ii)>0)THEN
65 i=nsv(ii)
66 xx=a(1,i)
67 yy=a(2,i)
68 zz=a(3,i)
69 a(1,i)=s11*xx+s12*yy+s13*zz
70 a(2,i)=s21*xx+s22*yy+s23*zz
71 a(3,i)=s31*xx+s32*yy+s33*zz
72 ENDIF
73 ENDDO
74 ELSE
75
76 DO ii=1,nsn
77 IF(tagkine(ii)>0)THEN
78 i=nsv(ii)
79 xx=a(1,i)
80 yy=a(2,i)
81 zz=a(3,i)
82 a(1,i)=s11*xx+s21*yy+s31*zz
83 a(2,i)=s12*xx+s22*yy+s32*zz
84 a(3,i)=s13*xx+s23*yy+s33*zz
85 ENDIF
86 ENDDO
87 ENDIF
88
89 RETURN