anchor.diff:              ! for each box, with the variables:
               !    domain_index, domain_image, box_left_i, box_right_i, box_bottom_j, box_top_j
               !
  -            ! FIXME: This would have to be changed to allow 'multi-time-level-communication',
  -            !        to deal with variable halow_width
               call get_domain_xs_ys(domains(j), xc_tmp, yc_tmp)
               call convert_priority_domain_info_to_boxes(&
                   priority_domain_index = domains(j)%nesting%priority_domain_index, &
  
positive.code: 
      end subroutine
  
      ! Routine to put all domains outputs in a single folder
      subroutine setup_multidomain_output_folder(md)
  
          class(multidomain_type), intent(inout) :: md
          integer:: i, date_time_values(8) ! For date_and_time
          character(len=charlen) :: datetime_char, mkdir_command
  
          call date_and_time(values=date_time_values)
  #ifdef COARRAY
  
negative.code: 
  #if defined(COARRAY_USE_MPI_FOR_INTENSIVE_COMMS)
  #ifdef REALFLOAT
      integer, parameter :: mympi_dp = MPI_REAL
  #else
      integer, parameter :: mympi_dp = MPI_DOUBLE_PRECISION
  #endif
  #endif
  
  