anchor.diff:    return
   end subroutine SOR_solver
   
  -! -----------------------------------------------------------------------------
  +! ---------------------------------------------------------------------------
   !> Update velocities using the barotropic tendency due to the pressure
   !> gradient.
   
  -subroutine update_velocities_for_barotropic_tendency(array, etanew, g_vec, xstep, ystep, dspace, dt, nx, ny, layers)
  +subroutine update_velocities_for_barotropic_tendency(array, etanew, g_vec, &
  +    xstep, ystep, dspace, dt, nx, ny, layers)
     implicit none
   
     double precision, intent(inout) :: array(0:nx+1, 0:ny+1, layers)
  
positive.code:   return
  end subroutine evaluate_b_iso
  
  ! -----------------------------------------------------------------------------
  
  subroutine evaluate_b_RedGrav(b, h, u, v, nx, ny, layers, gr)
    implicit none
  
negative.code:       uav = uav/real(avwrite)
        vav = vav/real(avwrite)
  
        call write_output(hav, uav, vav, etaav, wind_x, wind_y, &
            nx, ny, layers, n, RedGrav, .false., 'av')
  
        ! Check if there are NaNs in the data
        call break_if_NaN(h, nx, ny, layers, n)
  