anchor.diff:              stage_uh_vh_elev(2:3) = 0.0_dp
           else
   
  -            ! Much experimentation was conducted here. This problem is sensitive to
  -            ! imperfections in our semi-transmissive boundary conditions which allow a 
  -            ! stage forcing. (Not too surprising, because the boundary is quite close to the
  -            ! coast, in 'not very deep' water).
   
               select case(boundary_type)
  -
  -            case('boundary_stage_transmissive_normal_momentum')
  -                ! Do nothing, because we do not need uh/vh.
  -                stage_uh_vh_elev(2:3) = 0.0_dp 
  -
  -            case('boundary_stage_transmissive_momentum')
  -                ! Do nothing, because we do not need uh/vh.
  -                stage_uh_vh_elev(2:3) = 0.0_dp 
   
               case('boundary_stage_radiation_momentum') 
                   ! These will never be used for this boundary
  -                stage_uh_vh_elev(2:3) = 0.0_dp
  -
  -            case('flather_with_vh_equal_zero') 
  -                ! This absorbs, but also distorts the stage at the boundary too much when
  -                ! wave frequencies are lower
                   stage_uh_vh_elev(2:3) = 0.0_dp
   
               case('flather_with_vh_from_continuity')
  
positive.code:     use global_mod, only: ip, dp, minimum_allowed_depth, default_nonlinear_timestepping_method
      use domain_mod, only: domain_type
      use multidomain_mod, only: multidomain_type, setup_multidomain, test_multidomain_mod
      use boundary_mod, only: boundary_stage_transmissive_normal_momentum, flather_boundary, &
          boundary_stage_transmissive_momentum, boundary_stage_radiation_momentum
      use local_routines
      use timer_mod
      use logging_mod, only: log_output_unit
  
negative.code:             call elevation_data%get_xy(x, y, domain%U(:,j,ELV), domain%nx(1), &
                  bilinear=1_ip)
  
              flatten_bathymetry_near_boundary = .true.
  
          end do
  
  