1c { dg-do compile } 2 subroutine derv (xx,b,bv,det,r,s,t,ndopt,cosxy,thick,edis, 3 1 vni,vnt) 4 implicit real*8 (a-h,o-z) 5 save 6c 7 common /shell1/ disd(9),ield,ielp,npt,idw,ndrot 8 common /shell4/xji(3,3),p(3,32),h(32) 9c 10 dimension xx(3,*),ndopt(*),bv(*),vni(*),cosxy(6,*),vnt(*), 11 1 edis(*),thick(*),b(*) 12c 13 kk=0 14 k2=0 15 do 130 k=1,ield 16 k2=k2 + 3 17 if (ndopt(k)) 127,127,130 18 127 kk=kk + 1 19 do 125 i=1,3 20 b(k2+i)=b(k2+i) + (xji(i,1)*p(1,k) + xji(i,2)*p(2,k))*t 21 1 + xji(i,3)*h(k) 22 th=0.5*thick(kk) 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) 25 k2=k2 + 9 26 130 continue 27 return 28 end 29