anchor.diff:                      enddo
                       call TMscore(Lch,x1,y1,z1,nn1,Lch,x2,y2,z2,nn2,TM,
        &                           Rcomm,Lcomm)
  -                    write(*,*)'TMSCORE ',i,j,TM,Rcomm,Lcomm
       !                armsd=dsqrt(rms/Lch) !RMSD12
                       if (TM .eq. 0.0) then
                           rmsd = RMSD_max
  
positive.code: cjmht
        deallocate(x,y,z)
        deallocate(amat)
        deallocate(x1,y1,z1)
        deallocate(x2,y2,z2)
        deallocate(nn1,nn2)
        
        stop
        end
  
negative.code: cjmht dimension amat(nst,nst)    !RMSD matrics
        real, allocatable :: x(:,:), y(:,:), z(:,:)
        real, allocatable :: amat(:,:)
        dimension mark(nst)       !for removing used structures
        dimension n_str_near(nst) !numbef of neighboring structures
        dimension itra(nst),istr(nst),E(nst)
  