test including brunt updates during solver iterations

Author: evbauer

star/private/hydro_vars.f90

@@ -606,7 +606,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

star/private/solver_support.f90

@@ -829,7 +829,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.