anchor.diff:              if(myid .eq. 0) write(*,*) "chi2_new = ", chi2_new
   
               ! Test if the move should be accepted or rejected based on del_chi
  -            if(del_chi <0.0)then
  +            if(del_chi < 0.0) then
                   ! Accept the move
                   call fem_accept_move(mpi_comm_world)
                   e1 = e2 ! eam
  
positive.code:     scale_fac = scale_fac_initial
      res = 0.61/Q
      nk = size(k)
      beta=1./((8.6171e-05)*temperature)
      square_pixel = .TRUE. ! HRMC uses square pixels, not round.
  
      call fem_initialize(m, res, k, nk, ntheta, nphi, npsi, scatfact_e, istat,  square_pixel)
  
negative.code:                 write(53,*) nphi, npsi, ntheta
                  write(53,*) Q
              close(53)
  
              ! Write final vk
              open(unit=54,file=trim(vkf_fn),form='formatted',status='unknown')
              do i=1, nk
                  write(54,*)k(i),vk(i)
              enddo
              close(54)
  
              ! Write final model
              open(unit=55,file=trim(final_model_fn),form='formatted',status='unknown')
              write(55,*) m%natoms
  