1! { dg-do compile }
2! { dg-options "-O3" }
3      implicit real*8 (a-h,o-z)
4      common clop6(3),dps(6),aml6(6,6)
5      dimension y1(3)
6      dimension dclo(3)
7      dimension dx(3),dy(3)
8      save
9      do 80 ii=1,itco
10        y1(3)=dps(1)
11        do 40 l=1,3
12          dy(l)=clop6(l)-y1(l)
13   40   continue
14        dczp=abs(dy(3))
15        if(dcx.le.c1m10.and.dcz.le.c1m10.and.dcxp.le.c1m10.and.dczp
16     +  .le.c1m10.and.dcy.le.c1m10.and.dcyp.le.c1m10) goto 90
17   80 continue
18      write(6) itco
19      ii=itco
20   90 continue
21      if(ii.ne.itco) then
22        do 65 k=1,3
23          do 55 j=1,3
24            jj=2*j
25            kk=2*k
26            dclo(k)=aml6(kk-1,jj-1)*dx(j)+dclo(k)
27            dclo(k)=aml6(kk-1,jj)*dy(j)+dclo(k)
28   55     continue
29   65   continue
30      endif
31      end
32
33