anchor.diff:      md%domains(1)%timestepping_refinement_factor = 1_ip
       md%domains(1)%dx_refinement_factor = 1.0_dp
       md%domains(1)%timestepping_method = default_nonlinear_timestepping_method
  -    !md%domains(1)%compute_fluxes_inner_method = 'DE1'
  -    !md%domains(1)%theta = 4.0_dp !! No big impact for this problem.
  -    !md%domains(1)%timestepping_method = 'leapfrog_linear_plus_nonlinear_friction'
  -    !md%domains(1)%linear_solver_is_truely_linear = .false.
  -
  +    !md%domains(1)%use_eddy_viscosity = .true.
   
       print*, 1, ' lw: ', md%domains(1)%lw, ' ll: ', md%domains(1)%lower_left, ' dx: ', md%domains(1)%dx, &
           ' nx: ', md%domains(1)%nx
   
  -    ! A detailed domain [Cannot partially share a physical boundary with the outer domain]
  +    ! A detailed domain
       call md%domains(2)%match_geometry_to_parent(&
           parent_domain=md%domains(1), &
           lower_left=[28500.0_dp, 11500.0_dp], &
  
positive.code:     ! nd domains in this model
      nd = 2
      allocate(md%domains(nd))
      !md%load_balance_file = 'load_balance_partition.txt'
  
      !
      ! Setup basic metadata
  
negative.code:         dx_refinement_factor=nest_ratio, &
          timestepping_refinement_factor=nest_ratio,&
          rounding_method='nearest')
      md%domains(2)%timestepping_method = default_nonlinear_timestepping_method
      !md%domains(2)%use_eddy_viscosity = .true.
  
      print*, 2, ' lw: ', md%domains(2)%lw, ' ll: ', md%domains(2)%lower_left, ' dx: ', md%domains(2)%dx, &
          ' nx: ', md%domains(2)%nx
  