anchor.diff:          ! HRMC loop begins. The loop never stops.
           do while (i .le. step_end)
   
  -            i=i+1
  +            i = i+1
               if(myid .eq. 0) write(*,*) "Starting step", i
   
  -            call random_move(m,w,xx_cur,yy_cur,zz_cur,xx_new,yy_new,zz_new, max_move)
  +            call random_move(m, w, xx_cur, yy_cur, zz_cur, xx_new, yy_new, zz_new, max_move)
               ! check_curoffs returns false if the new atom placement is too close to
               ! another atom. Returns true if the move is okay. (hard shere cutoff)
               do while( .not. check_cutoffs(m,cutoff_r,w) )
  
positive.code:     real, allocatable, dimension(:,:) :: cutoff_r 
      real, pointer, dimension(:,:) :: scatfact_e
      real :: xx_cur, yy_cur, zz_cur, xx_new, yy_new, zz_new
      real :: scale_fac, scale_fac_initial, beta
      double precision :: chi2_old, chi2_new, del_chi, chi2_gr, chi2_vk, chi2_no_energy, chi2_initial
      real :: R
      integer :: i, j, step_start, step_end
  
negative.code:     character (len=256) :: param_filename
      character (len=256) :: eam_filename
      character (len=256) :: jobID, c, step_str
      character (len=256) :: vki_fn, vkf_fn, output_model_fn, final_model_fn, chi_squared_file, acceptance_rate_fn, femfile
      character (len=256) :: paramfile_restart
      real :: temperature
      real :: max_move
  