anchor.diff:        close(4)
   c^^^^^^^^^^^^^^^^^^^^^^^^ n_str_all done ^^^^^^^^^^^^^^^^^^^^^^^^
   
  +cjmht we can now allocate the arrays
  +      allocate(x(Lch,n_str_all))
  +      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
    11   read(3,1239,end=5)tx,ii,tx,tx,ii,tx,a1,a2,a3
  
positive.code:       enddo
  c^^^^^^^^^^^^^^^^^^^^output done ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  
        stop
        end
        
  
negative.code: ccc   
        dimension xt(ndim),yt(ndim),zt(ndim) !temporal coordinate
        dimension x_n(ndim),y_n(ndim),z_n(ndim) !native structure
  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)
  