anchor.diff:      !------------------- Start HRMC. -----------------!
   
   #ifndef SERIAL
  -    call mpi_barrier(colored_comm, mpierr)
  +    call mpi_barrier(communicator, mpierr)
   #endif
   
           ! Calculate initial chi2
  
positive.code:     !------------------- Read inputs and initialize. -----------------!
  
      ! Set input filenames.
  #ifndef SERIAL
      call mpi_bcast(param_filename, 256, MPI_CHARACTER, 0, colored_comm, mpierr)
  #endif
      if(myid.eq.0) write(*,*) "Paramfile: ", trim(param_filename)
      
      ! Read input parameters
      call read_inputs(param_filename, femsim, model_filename, femfile, Q, ntheta, nphi, npsi, scale_fac_initial, eam_filename, step_start, step_end, temp_move_decrement, temperature, max_move, cutoff_r, seed, alpha, vk_exp, k, vk_exp_err)
  
negative.code:             endif
  
              if(energy_accepted) then
                  call fem_update(m, w, res, k, vk, scatfact_e, communicator, istat)
  
                  chi2_no_energy = chi_square(alpha, vk_exp, vk_exp_err, vk, scale_fac, nk)
                  chi2_new = chi2_no_energy + te2
                  del_chi = chi2_new - chi2_old
  #ifndef SERIAL
                  call mpi_bcast(del_chi, 1, mpi_double, 0, communicator, mpierr)
  #endif
  
                  ! Test if the move should be accepted or rejected based on del_chi
  