anchor.diff:      use logging_mod, only: log_output_unit
       implicit none
   
  +
       contains 
   
  -    subroutine set_initial_conditions(domain)            
  +    subroutine set_initial_conditions(domain, grid_smooth_iterations)            
           class(domain_type), intent(inout):: domain
  -        integer(ip):: i, j
  +        integer(ip), intent(in) :: grid_smooth_iterations
  +
  +        integer(ip):: j
           real(dp), allocatable:: x(:), y(:)
           type(multi_raster_type) :: elevation_rast
   
  
positive.code:     ! 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)%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
  
negative.code:     character(len=charlen) :: timestepping_method
  
      ! Approx timestep between outputs
      real(dp), parameter :: approximate_writeout_frequency = 1000
      real(dp), parameter :: final_time = 3600.0_dp
      real(dp), parameter :: mesh_refine = 1.0_dp
      real(dp) :: my_dt 
  
      ! Lower-left corner coordinate
  