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:       dimension order(nst)
        dimension ires(ndim)
        double precision rmsd_a,rmsd2_a
        
  cjmht dimension x1(ndim),y1(ndim),z1(ndim),nn1(ndim)
  cjmht dimension x2(ndim),y2(ndim),z2(ndim),nn2(ndim)
        real, allocatable :: x1(:),y1(:),z1(:)
        real, allocatable :: x2(:),y2(:),z2(:)
        integer, allocatable :: nn1(:), nn2(:)
  
  cjmht
        integer score_type
  
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
  