anchor.diff:      contains 
   
       subroutine set_initial_conditions_radial_dam(domain)            
  +        !
  +        ! This function sets the initial conditions in a domain
  +        !
           class(domain_type), target, intent(inout):: domain
           integer(ip):: i,j
           real(dp):: x, y, cx, cy, initial_stage_1, initial_stage_2, radius
  
positive.code: 
          ! Stage
          domain%U(:,:,STG) = initial_stage_1 
  
          ! Elevation
          domain%U(:,:,ELV) = 0._dp
  
negative.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
  
      ! Misc
      integer :: j, nd
  