anchor.diff:        enddo
   c^^^^^^^^^^^^^^^^^^^^output done ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
   
  -cjmht
  -      deallocate(x,y,z)
  -      deallocate(amat)
  -      
         stop
         end
         
  
positive.code:          y1(j)=t(2)+u(2,1)*xa(j)+u(2,2)*ya(j)+u(2,3)*za(j)
           z1(j)=t(3)+u(3,1)*xa(j)+u(3,2)*ya(j)+u(3,3)*za(j)
        enddo
  cjmht
  cjmht Attempt to get rmsd of TM-score component
  ccc   RMSD (d<5.0)-------->
        LL=0
        do i=1,n_cut
           m=i_ali(i)             ![1,nseqA]
           r_1(1,i)=xa(iA(m))
           r_1(2,i)=ya(iA(m))
           r_1(3,i)=za(iA(m))
           r_2(1,i)=xb(iB(m))
           r_2(2,i)=yb(iB(m))
           r_2(3,i)=zb(iB(m))
           LL=LL+1
        enddo
        call u3b(w,r_1,r_2,LL,0,rms,u,t,ier)
        Rcomm=dsqrt(rms/LL)
        Lcomm=LL
  cjmht END
        TM=score_max
  
  c^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  
negative.code: *******************************************************************************
  *     SPICKER_2.0, released on December 29, 2010
  *     
  