anchor.diff:      character(len=charlen) :: timestepping_method
   
       ! Approx timestep between outputs
  -    real(dp) :: approximate_writeout_frequency = 1000
  -    real(dp) :: final_time = 3600.0_dp 
  +    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
  
positive.code:         domain%U(:,:,ELV) = min(domain%U(:,:,ELV), -100.0_dp)
  
          ! Ensure topography is smooth
          do j = 1,10
              call domain%smooth_elevation
          end do
  
  
negative.code:     use logging_mod, only: log_output_unit
      implicit none
  
  
      contains 
  
      subroutine set_initial_conditions(domain, grid_smooth_iterations)            
          class(domain_type), intent(inout):: domain
          integer(ip), intent(in) :: grid_smooth_iterations
  
          integer(ip):: j
          real(dp), allocatable:: x(:), y(:)
          type(multi_raster_type) :: elevation_rast
  
  