2 subroutine derv
(xx
,b
,bv
,det
,r
,s
,t
,ndopt
,cosxy
,thick
,edis
,
4 implicit real*8 (a
-h
,o
-z
)
7 common /shell1
/ disd
(9),ield
,ielp
,npt
,idw
,ndrot
8 common /shell4
/xji
(3,3),p
(3,32),h
(32)
10 dimension xx
(3,*),ndopt
(*),bv
(*),vni
(*),cosxy
(6,*),vnt
(*),
11 1 edis
(*),thick
(*),b
(*)
17 if (ndopt
(k
)) 127,127,130
20 b
(k2
+i
)=b
(k2
+i
) + (xji
(i
,1)*p
(1,k
) + xji
(i
,2)*p
(2,k
))*t
23 b
(k2
+i
+3)=b
(k2
+i
+3) - th*cosxy
(i
+3,kk
)
24 125 b
(k2
+i
+6)=b
(k2
+i
+6) + th*cosxy
(i
,kk
)