anchor.diff:          !
           do j = 1, nd_local
   
  -            if(verbose1) write(log_output_unit,*) ' '
  -            if(verbose1) write(log_output_unit,*) ' #################'
  -            if(verbose1) write(log_output_unit,*) ' Domain ID: ', domains(j)%myid
  +            if(verbose1) write(md_log_output_unit,*) ' '
  +            if(verbose1) write(md_log_output_unit,*) ' #################'
  +            if(verbose1) write(md_log_output_unit,*) ' Domain ID: ', domains(j)%myid
   
               domains(j)%boundary_exterior = (.NOT. domains(j)%is_nesting_boundary)
   
  
positive.code:                 gauge_ids = point_gauges(3,:))
          end do
  
          write(log_output_unit,*) '    Gauges are setup'
          flush(log_output_unit)
  
          deallocate(point_gauges, time_var_local, static_var_local)
  
  
negative.code:     !
      ! @param domains the array of domains
      !
      subroutine create_all_bbox_and_dx(domains, periodic_xs, periodic_ys, md_log_output_unit)
  
          type(domain_type), intent(inout) :: domains(:)
          real(dp), intent(in) :: periodic_xs(2), periodic_ys(2)
          integer(ip), intent(in) :: md_log_output_unit
  
          integer(ip):: nd, i, j, k, nd_local, ierr
  
  