anchor.diff:  module local_routines 
  -    !! Test the convergence of the code in a periodic domain
  +    !! Test flow algorithms in a large area spherical box with complex topography
       use global_mod, only: dp, ip, wall_elevation, pi, charlen
       use domain_mod, only: domain_type, STG, UH, VH, ELV
       use read_raster_mod, only: multi_raster_type
  
positive.code:         do j = 1, domain%nx(2)
              y = domain%y(j)
  
              call elevation_rast%get_xy(x, y, domain%U(:,j,ELV), domain%nx(1), bilinear=1_ip)
  
              domain%U(:,j,STG) = exp(sin(2*pi*x/10.0_dp))*cos(2*pi*y/10.0_dp)
          end do
  
          deallocate(x,y)
         
          ! Wet everywhere 
          domain%U(:,:,ELV) = min(domain%U(:,:,ELV), -100.0_dp)
  
          ! Ensure topography is smooth
  
negative.code: !@!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  
  program spherical_box
      !! Test flow algorithms in a large area spherical box with complex topography
  
      use global_mod, only: ip, dp, charlen
      use multidomain_mod, only: multidomain_type
      use timer_mod, only: timer_type
      use logging_mod, only: log_output_unit
  