anchor.diff:  cjmht dimension amat(nst,nst)    !RMSD matrics
         real, allocatable :: x(:,:), y(:,:), z(:,:)
         real, allocatable :: amat(:,:)
  -cjmht
         dimension mark(nst)       !for removing used structures
         dimension n_str_near(nst) !numbef of neighboring structures
         dimension itra(nst),istr(nst),E(nst)
  
positive.code:       double precision score,score_max
        dimension xa(nmax),ya(nmax),za(nmax)
  !$OMP THREADPRIVATE(/stru/,/nres/,/para/,/align/,/nscore/,/scores/)
  cjmht dimension x1(nmax),y1(nmax),z1(nmax),n1(nmax)
  cjmht dimension x2(nmax),y2(nmax),z2(nmax),n2(nmax)
        real, allocatable, intent(inout) :: x1(:),y1(:),z1(:)
        integer, allocatable, intent(in) :: n1(:)
        real, allocatable, intent(inout) :: x2(:),y2(:),z2(:)
        integer, allocatable, intent(in) :: n2(:)
  
  ccc   RMSD:
        double precision r_1(3,nmax),r_2(3,nmax),r_3(3,nmax),w(nmax)
  
negative.code:       allocate(y(Lch,n_str_all))
        allocate(z(Lch,n_str_all))
        allocate(amat(n_str_all,n_str_all))
        
  cccccccccccc read native structure cccccccccccccccccccc
        i=0
  