anchor.diff:           end do
            end if
          end do
  -       do j=1,nb
  -         write(55,*) gr(j)
  -       end do
  +c       do j=1,nb
  +c         write(55,*) gr(j)
  +c       end do
          i1 = count(oc1)
          i2 = count(oc2)
   c
  -       write(57,*) i1,i2
          fourpi = 16.0*atan(1.0)
          const = fourpi*real(i1*i2)/(vol*3.0d0)
          do k=1,nb
  
positive.code:              ang = SUM(rd*rd1)/(d1*d)
               ocons = ang.ge.ang0
               ocons = ocons.and.(ang.le.ang1)
               write(45,*) i2,i1,ic,ang,ocons
             end if
             if(d.gt.0.01) k = int(d/dr)+1
             if(k.le.nb.and.ocons ) then
  
negative.code:       if(string.eq." ") return
        l  = len_trim(string)
        do i=1,l
         if(string(i:i).ne." ") exit
        end do
        is_integer = verify(string(i:l), anumber).eq.0
  