anchor.diff:              stage_uh_vh_elev(2:3) = 0.0_dp
           else
   
  -            ! Much experimentation was conducted here. This problem seems to be sensitive to
  +            ! 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).
  
positive.code:     ! affects the tsunami in this case.
      !
      ! Possible values are:
      !   'boundary_stage_transmissive_normal_momentum' -- tsunami waves are over-amplified, probably 
      !       because of reflections from the boundary. tide is OK
      !   'boundary_stage_transmissive_momentum' -- Similar to above
      !   'flather_with_uh_equal_zero' -- tsunami waves under-amplified (i.e. too small). Tide is OK
      !   'flather_with_vh_from_continuity' -- about right? But depends on a scale for the boundary 
      !       VH term, which we can estimate from continuity, but needs tuning.
      !   
      character(len=charlen), parameter :: boundary_type = 'flather_with_vh_from_continuity'
      !character(len=charlen), parameter :: boundary_type = 'boundary_stage_transmissive_momentum'
  
      ! This will hold the information -- is seen by other parts of the module
  
negative.code:             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
  