anchor.diff:      ! grid size (number of x/y cells)
       integer(ip), parameter :: global_nx(2) = [945_ip, 615_ip] 
   
  -    ! Number of domains in model; loop variable
  -    integer(ip):: nd, j
  +    integer(ip):: j
   
       ! Set the model name
       md%output_basedir = './OUTPUTS/'
  
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: module local_routines 
      !! 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
  