anchor.diff:              end do
           end do
   
  +        ! Manning friction squared
           domain%manning_squared = 0.0_dp
   
       end subroutine
  
positive.code:     real(dp), parameter, dimension(2):: global_ll = [0._dp, 0._dp]
      ! grid size (number of x/y cells)
      integer(ip), parameter, dimension(2):: global_nx = [400, 400] * 2 + 1
  
      ! analytical solution
      real(dp), allocatable :: analytical_solution(:,:)
      character(len=charlen):: analytical_solution_file
  
      ! Misc
      integer :: j, nd
  
negative.code: 
          ! Stage
          domain%U(:,:,STG) = initial_stage_1 
          ! Depth integrated velocity
          domain%U(:,:, UH:VH) = 0.0_dp
  
          ! Elevation
          domain%U(:,:,ELV) = 0._dp
  