consolidate calculations for vel_liq

Exawind · Jan 16, 2025 · fbfd78d · fbfd78d
1 parent d64f260
commit fbfd78d
Showing 1 changed file with 2 additions and 5 deletions.
diff --git a/amr-wind/ocean_waves/relaxation_zones/relaxation_zones_ops.cpp b/amr-wind/ocean_waves/relaxation_zones/relaxation_zones_ops.cpp
@@ -183,11 +183,8 @@ void apply_relaxation_zones(CFDSim& sim, const RelaxZonesBaseData& wdata)
                             } else {
                                 out_vel = 0.0;
                             }
-                            amrex::Real vel_liq = vel(i, j, k, n);
-                            const amrex::Real dvel_liq =
-                                utils::combine_linear(Gamma, out_vel, vel_liq) -
-                                vel_liq;
-                            vel_liq += dvel_liq;
+                            const amrex::Real vel_liq = utils::combine_linear(
+                                Gamma, out_vel, vel(i, j, k, n));
                             amrex::Real integrated_vel_liq =
                                 volfrac(i, j, k) * vel_liq;
                             integrated_vel_liq += amrex::max(0.0, dvf) *