anchor.diff:      ! Setup basic metadata
       md%domains(1)%lw = global_lw
       md%domains(1)%lower_left =global_ll
  -    md%domains(1)%nx = global_nx
  +    md%domains(1)%nx = nint(global_nx * mesh_refine)
       md%domains(1)%dx = md%domains(1)%lw/md%domains(1)%nx
       md%domains(1)%dx_refinement_factor = 1.0_dp
       md%domains(1)%timestepping_refinement_factor = 1_ip
  
positive.code:     use logging_mod, only: log_output_unit
      implicit none
  
      contains 
  
      subroutine set_initial_conditions(domain)            
          class(domain_type), intent(inout):: domain
          integer(ip):: i, j
          real(dp), allocatable:: x(:), y(:)
          type(multi_raster_type) :: elevation_rast
  
  
negative.code:         domain%U(:,:,ELV) = min(domain%U(:,:,ELV), -100.0_dp)
  
          ! Ensure topography is smooth
          do j = 1, grid_smooth_iterations
              call domain%smooth_elevation
          end do
  
  