anchor.diff:  
       subroutine setup_multidomain_domains(domains, verbose, use_wetdry_limiting_nesting,&
           periodic_xs, periodic_ys, extra_halo_buffer, extra_cells_in_halo, &
  -        all_dx_md, all_timestepping_methods_md, md_label, md_p2p)
  +        all_dx_md, all_timestepping_methods_md, md_label, md_p2p, md_log_output_unit)
           !! Main routine for setting up the multidomain
           !!
           !! Determines which domains we need to have a two-way-nesting-comms relation
  
positive.code:             end do
  
          end do
          flush(log_output_unit)
  
          ! Free memory we don't need. Note some of these variables
          ! will grow in size rapidly with num_images, so good to be careful.
  
negative.code:                 ys=yc_tmp, &
                  periodic_xs=periodic_xs,&
                  periodic_ys=periodic_ys,&
                  max_parent_dx_ratio=nint(domains(j)%max_parent_dx_ratio),&
                  md_log_output_unit=md_log_output_unit)
  
              ! Keep track of the maximum number of 'boxes' in the box metadata
              tmp2 = shape(domains(j)%nesting%recv_metadata)
  