anchor.diff:          !
           do j = 1, nd_local
   
  -            if(verbose1) write(log_output_unit, *) 'Domain ', j
  +            if(verbose1) write(md_log_output_unit, *) 'Domain ', j
   
               ! Initial value to be refined in loop
               max_parent_dx_ratio = 1.0_dp
  
positive.code:         call co_sum(global_energy_total_on_rho)
  #endif
  
          write(log_output_unit, "(A)"         ) ''
          write(log_output_unit, "(A)"         ) '-----------'
          write(log_output_unit, "(A)"         ) 'Global stage range (over all domains and images): '
          write(log_output_unit, "(A, ES25.12E3)") '        ', global_max_stage
          write(log_output_unit, "(A, ES25.12E3)") '        ', global_min_stage
          write(log_output_unit, "(2A)"        ) 'Global speed range (over all domains and images): '
          write(log_output_unit, "(A, ES25.12E3)") '        ', global_max_speed
          write(log_output_unit, "(A, ES25.12E3)") '        ', global_min_speed
          write(log_output_unit, "(2A)"        ) &
              'Global energy-potential / rho (over all domains and images), zero when stage=domain%msl_linear: '
          write(log_output_unit, "(A, ES25.12E3)") '        ', global_energy_potential_on_rho
          write(log_output_unit, "(2A)"        ) 'Global energy-kinetic / rho (over all domains and images): '
          write(log_output_unit, "(A, ES25.12E3)") '        ', global_energy_kinetic_on_rho
          write(log_output_unit, "(2A)"        ) &
              'Global energy-total / rho (over all domains and images): '
          write(log_output_unit, "(A, ES25.12E3)") '        ', global_energy_total_on_rho
          write(log_output_unit, "(A)") '-----------'
          call md%report_mass_conservation_statistics()
  
          ! Unstable models may produce NaN energy -- this is useful to detect.
  
negative.code:     ! @param verbose logical
      !
      subroutine compute_multidomain_nesting_layer_width(domains, verbose,&
          periodic_xs, periodic_ys, extra_halo_buffer, extra_cells_in_halo,&
          md_log_output_unit)
          type(domain_type), intent(inout) :: domains(:)
          logical, optional, intent(in) :: verbose
          real(dp), intent(in) :: periodic_xs(2), periodic_ys(2)
          integer(ip), intent(in) :: extra_halo_buffer
          integer(ip), intent(in) :: extra_cells_in_halo
          integer(ip), intent(in) :: md_log_output_unit
  
          integer(ip):: i, j, ii, nd_local, nest_layer_width, boundary_flag, n1
  
  