anchor.diff:  
           ! Stage
           domain%U(:,:,STG) = initial_stage_1 
  +        ! Depth integrated velocity
  +        domain%U(:,:, UH:VH) = 0.0_dp
   
           ! Elevation
           domain%U(:,:,ELV) = 0._dp
  
positive.code:             end do
          end do
  
          domain%manning_squared = 0.0_dp
  
      end subroutine
  
negative.code:     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
  