anchor.diff:         integer, dimension(:,:), allocatable :: p2
          integer, dimension(:,:), allocatable :: p3
          integer, dimension(:,:), allocatable :: p4
  +       integer, dimension(:,:), allocatable :: p5
  +       integer, dimension(:,:), allocatable :: p6
  +       integer, dimension(:,:), allocatable :: p7
  +       integer, dimension(:,:), allocatable :: p8
  +       integer, dimension(:,:), allocatable :: p9
          character*16 , dimension(:), allocatable :: ta,tr
          integer i,il
          integer n2
          integer n3
          integer n4
  +       integer n5
  +       integer n6
  +       integer n7
  +       integer n8
  +       integer n9
    
  +       integer nmax
  +       nmax = 500
   c
          call smd_pdb_natoms(infile,ntot)
          allocate(c(3,ntot))
  
positive.code: 18     continue
         write(*,*) "input file is",trim(infile)
         call cluster(infile)
         stop
  c
         end
  
  
negative.code:       integer i2,i3,j3,ik
        integer fn_out
        logical ofound,omatch
        character*73 message
  
        message = "opening file"
        fn_out = 33
         open(unit=fn_out,file="fourimer.dat",
       $            form='formatted',status='unknown',err=911)
  