fix rho-weighting of anelastic buoyancy (erf-model#1924)

AMLattanzi · Nov 3, 2024 · 08ff588 · 08ff588
1 parent 84efdc3
commit 08ff588
Showing 1 changed file with 2 additions and 2 deletions.
diff --git 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