anchor.diff:      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
  
positive.code: 
          ! Stage
          domain%U(:,:,STG) = initial_stage_1 
  
          ! Elevation
          domain%U(:,:,ELV) = 0._dp
  
negative.code:             end do
          end do
  
          ! Manning friction squared
          domain%manning_squared = 0.0_dp
  
      end subroutine
  