From d38cf680944d4f41c99d948a8a528ad5f54e5d27 Mon Sep 17 00:00:00 2001 From: AMLattanzi Date: Wed, 20 Dec 2023 12:21:26 -0800 Subject: [PATCH] Fix qv_for_p in slow rhs pre. --- Source/TimeIntegration/ERF_slow_rhs_pre.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/Source/TimeIntegration/ERF_slow_rhs_pre.cpp b/Source/TimeIntegration/ERF_slow_rhs_pre.cpp index 496623354..2eaa94cb0 100644 --- a/Source/TimeIntegration/ERF_slow_rhs_pre.cpp +++ b/Source/TimeIntegration/ERF_slow_rhs_pre.cpp @@ -549,7 +549,7 @@ void erf_slow_rhs_pre (int level, int finest_level, //if (cell_data(i,j,k,RhoTheta_comp) < 0.) printf("BAD THETA AT %d %d %d %e %e \n", // i,j,k,cell_data(i,j,k,RhoTheta_comp),cell_data(i,j,k+1,RhoTheta_comp)); AMREX_ASSERT(cell_data(i,j,k,RhoTheta_comp) > 0.); - Real qv_for_p = (use_moisture) ? cell_data(i,j,k,RhoQ1_comp) : 0.0; + Real qv_for_p = (use_moisture) ? cell_data(i,j,k,RhoQ1_comp)/cell_data(i,j,k,Rho_comp) : 0.0; pp_arr(i,j,k) = getPgivenRTh(cell_data(i,j,k,RhoTheta_comp),qv_for_p) - p0_arr(i,j,k); }); } // end profile