diff --git a/Source/SourceTerms/ERF_buoyancy_utils.H b/Source/SourceTerms/ERF_buoyancy_utils.H index 061240ae5..e8db1015e 100644 --- a/Source/SourceTerms/ERF_buoyancy_utils.H +++ b/Source/SourceTerms/ERF_buoyancy_utils.H @@ -26,7 +26,7 @@ buoyancy_dry_anelastic (int& i, amrex::Real theta_d_0_wface = amrex::Real(0.5) * (theta_d_0_lo + theta_d_0_hi); - return (-grav_gpu * (theta_d_wface - theta_d_0_wface) / theta_d_0_wface); + return (-grav_gpu * (theta_d_wface - theta_d_0_wface) / theta_d_0_wface * 0.5 * (r0_arr(i,j,k) + r0_arr(i,j,k-1))); } AMREX_GPU_DEVICE @@ -87,7 +87,7 @@ buoyancy_moist_anelastic (int& i, amrex::Real theta_v_0_wface = amrex::Real(0.5) * (theta_v_0_lo + theta_v_0_hi); - return (-grav_gpu * (theta_v_wface - theta_v_0_wface) / theta_v_0_wface); + return (-grav_gpu * (theta_v_wface - theta_v_0_wface) / theta_v_0_wface * 0.5 * (r0_arr(i,j,k) + r0_arr(i,j,k-1))); } AMREX_GPU_DEVICE