anchor.diff:  ccc   
         dimension xt(ndim),yt(ndim),zt(ndim) !temporal coordinate
         dimension x_n(ndim),y_n(ndim),z_n(ndim) !native structure
  -      dimension x(ndim,nst),y(ndim,nst),z(ndim,nst) !used structures
  -      dimension amat(nst,nst)    !RMSD matrics
  +cjmht dimension x(ndim,nst),y(ndim,nst),z(ndim,nst) !used structures
  +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)
  
positive.code: *******************************************************************************
  *     SPICKER_2.0, released on December 29, 2010
  *     
  
negative.code: c       -2: no result obtained because of negative weights w
  c           or all weights equal to zero.
  cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
        subroutine rmsd_score_all(ndim,w,n_str,Lch,x,y,z,amat, rmsd_delta,
       &                          rmsd_a)
  c     REM: I, J, K, L, M, or N are int
        implicit none
  c     Intent In
  cjmht integer ndim,nst
        integer ndim
        integer n_str,Lch
        double precision w(ndim)
  cjmht real x(ndim,nst),y(ndim,nst),z(ndim,nst)
        real, allocatable, intent(in) :: x(:,:), y(:,:), z(:,:)
  c     Intent Out
  cjmht real amat(nst,nst)
        real, allocatable, intent(inout) :: amat(:,:)
        real rmsd_delta
        double precision rmsd_a
  c     Local variables
  