diff --git a/amr-wind/incflo_advance.cpp b/amr-wind/incflo_advance.cpp index d87bd7d841..3890fcd063 100644 --- a/amr-wind/incflo_advance.cpp +++ b/amr-wind/incflo_advance.cpp @@ -57,9 +57,9 @@ void incflo::advance(int inonlin) ApplyPredictor(false, inonlin); if (!m_use_godunov) { - if (inonlin > 0) { - amrex::Abort("Advection iterations are not supported for MOL"); - } + // if (inonlin > 0) { + // amrex::Abort("Advection iterations are not supported for MOL"); + // } ApplyCorrector(); } @@ -354,15 +354,18 @@ void incflo::ApplyPredictor(bool incremental_projection, int inonlin) } // With scalars computed, compute advection of momentum - if (inonlin == 0) { - icns().compute_advection_term(amr_wind::FieldState::Old); - } else { - icns().compute_advection_term(amr_wind::FieldState::NPH); - } - - // ************************************************************************************* - // Define (or if use_godunov, re-define) the forcing terms and viscous terms - // independently for the right hand side, without 1/rho term + // if (inonlin == 0) { + // icns().compute_advection_term(amr_wind::FieldState::Old); + //} else { + // icns().compute_advection_term(amr_wind::FieldState::NPH); + //} + const auto fstate = + (inonlin == 0) ? amr_wind::FieldState::Old : amr_wind::FieldState::NPH; + icns().compute_advection_term(fstate); + + // ************************************************************************************* + // Define (or if use_godunov, re-define) the forcing terms and viscous + // terms independently for the right hand side, without 1/rho term // ************************************************************************************* icns().compute_source_term(amr_wind::FieldState::NPH); icns().compute_diffusion_term(amr_wind::FieldState::Old); diff --git a/amr-wind/utilities/console_io.cpp b/amr-wind/utilities/console_io.cpp index e33653016a..0f37287794 100644 --- a/amr-wind/utilities/console_io.cpp +++ b/amr-wind/utilities/console_io.cpp @@ -224,8 +224,6 @@ void print_nonlinear_residual( const auto& velocity_new = sim.pde_manager().icns().fields().field; const auto& oset_mask = sim.repo().get_int_field("mask_cell"); - const int ncomp = AMREX_SPACEDIM; - for (int lev = 0; lev < nlevels; ++lev) { amrex::iMultiFab level_mask;