Fix qv_for_p in slow rhs pre.

erf-model · Dec 20, 2023 · d38cf68 · d38cf68
1 parent 6490d20
commit d38cf68
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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