anchor.diff:  
                   end if
   
  -            end subroutine 
  +            end subroutine
   
       end subroutine
   
       ! Make elevation constant in nesting send_regions that go to a single coarser cell, if the maximum elevation is above
       ! elevation_threshold
  -    ! 
  +    !
       ! This was done (in the past) to avoid wet-dry instabilities, caused by aggregating over wet-and-dry cells on a finer domain,
       ! which is then sent to a coarser domain.  Such an operation will break the hydrostatic balance, unless the elevation in the
       ! fine cells is constant
  
positive.code:         if(present(static_var)) then
              allocate(sv(size(static_var, kind=ip)))
              sv = static_var
          else 
              ! Default case -- store elevation once
              allocate(sv(1))
              sv = [ELV]
  
negative.code:     use cliffs_tolkova_mod, only: cliffs, setSSLim
      use extrapolation_limiting_mod, only: limited_gradient_dx_vectorized
  
  #ifdef SPHERICAL
      ! Compile with -DSPHERICAL to get the code to run in spherical coordinates
      use spherical_mod, only: area_lonlat_rectangle, deg2rad, earth_angular_freq
      use global_mod, only: radius_earth
  