From 22401acfd92763530f8c26877f9276ddaa8acb68 Mon Sep 17 00:00:00 2001
From: evbauer <18405113+evbauer@users.noreply.github.com>
Date: Fri, 16 Jun 2023 13:55:43 -0400
Subject: [PATCH] test including brunt updates during solver iterations

---
 star/private/hydro_vars.f90     | 2 +-
 star/private/solver_support.f90 | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/star/private/hydro_vars.f90 b/star/private/hydro_vars.f90
index 00206213d..d9c9cf157 100644
--- a/star/private/hydro_vars.f90
+++ b/star/private/hydro_vars.f90
@@ -593,7 +593,7 @@ subroutine set_hydro_vars( &
             
          end if
 
-         if (.not. skip_brunt) then ! skip_brunt during solver iterations
+         if (.not. skip_brunt) then
             if (dbg) write(*,*) 'call do_brunt_N2'
             call do_brunt_N2(s, nzlo, nzhi, ierr)
             if (failed('do_brunt_N2')) return
diff --git a/star/private/solver_support.f90 b/star/private/solver_support.f90
index e76ff8d22..69d3d1b6f 100644
--- a/star/private/solver_support.f90
+++ b/star/private/solver_support.f90
@@ -841,7 +841,7 @@ subroutine set_vars_for_solver(s, nvar, nzlo, nzhi, dt, ierr)
             skip_grads = .true., &
             skip_rotation = .true., &
             skip_m_grav_and_grav = .true., &
-            skip_brunt = .true., &
+            skip_brunt = .false., & ! .true., & ! EXPERIMENTAL: tesing not skipping brunt updates during solver iterations
             skip_mixing_info = .true., &
             skip_set_cz_bdy_mass = .true., &
             skip_other_cgrav = .true.