Fix missing END_GPU_PARALLEL_LOOP, swap v/g example properties, document beta_vars

- Add missing END_GPU_PARALLEL_LOOP in s_enforce_EL_bubbles_boundary_conditions
- Fix swapped vapor/gas properties in 2D_moving_lag_bubs and
  3D_moving_lag_particles examples (gam_v/gam_g, M_v/M_g, k_v/k_g,
  cp_v/cp_g, R_v/R_g were all inverted vs reference examples)
- Document beta_vars indices: 1=void fraction, 2=d(beta)/dt, 5=energy

Co-Authored-By: Claude Opus 4.6 (1M context) <noreply@anthropic.com>

Author: sbryngelson

examples/2D_moving_lag_bubs/case.py

@@ -132,16 +132,16 @@
     "bub_pp%pv": pv / p0,
     "bub_pp%vd": diffVapor / (L0 * c0),
     "bub_pp%mu_l": mu_host / (rho0 * L0 * c0),
-    "bub_pp%gam_v": gamma_g,
-    "bub_pp%gam_g": gamma_v,
-    "bub_pp%M_v": MW_g,
-    "bub_pp%M_g": MW_v,
-    "bub_pp%k_v": k_g * (T0 / (L0 * rho0 * c0 * c0 * c0)),
-    "bub_pp%k_g": k_v * (T0 / (L0 * rho0 * c0 * c0 * c0)),
-    "bub_pp%cp_v": cp_g * (T0 / (c0 * c0)),
-    "bub_pp%cp_g": cp_v * (T0 / (c0 * c0)),
-    "bub_pp%R_v": (R_uni / MW_g) * (T0 / (c0 * c0)),
-    "bub_pp%R_g": (R_uni / MW_v) * (T0 / (c0 * c0)),
+    "bub_pp%gam_v": gamma_v,
+    "bub_pp%gam_g": gamma_g,
+    "bub_pp%M_v": MW_v,
+    "bub_pp%M_g": MW_g,
+    "bub_pp%k_v": k_v * (T0 / (L0 * rho0 * c0 * c0 * c0)),
+    "bub_pp%k_g": k_g * (T0 / (L0 * rho0 * c0 * c0 * c0)),
+    "bub_pp%cp_v": cp_v * (T0 / (c0 * c0)),
+    "bub_pp%cp_g": cp_g * (T0 / (c0 * c0)),
+    "bub_pp%R_v": (R_uni / MW_v) * (T0 / (c0 * c0)),
+    "bub_pp%R_g": (R_uni / MW_g) * (T0 / (c0 * c0)),
     # Viscosity
     "viscous": "T",
     "fluid_pp(1)%Re(1)": 1.0 / (mu_host / (rho0 * c0 * L0)),

examples/3D_moving_lag_particles/case.py

@@ -127,16 +127,16 @@
     "bub_pp%pv": pv / p0,
     "bub_pp%vd": diffVapor / (L0 * c0),
     "bub_pp%mu_l": mu_host / (rho0 * L0 * c0),
-    "bub_pp%gam_v": gamma_g,
-    "bub_pp%gam_g": gamma_v,
-    "bub_pp%M_v": MW_g,
-    "bub_pp%M_g": MW_v,
-    "bub_pp%k_v": k_g * (T0 / (L0 * rho0 * c0 * c0 * c0)),
-    "bub_pp%k_g": k_v * (T0 / (L0 * rho0 * c0 * c0 * c0)),
-    "bub_pp%cp_v": cp_g * (T0 / (c0 * c0)),
-    "bub_pp%cp_g": cp_v * (T0 / (c0 * c0)),
-    "bub_pp%R_v": (R_uni / MW_g) * (T0 / (c0 * c0)),
-    "bub_pp%R_g": (R_uni / MW_v) * (T0 / (c0 * c0)),
+    "bub_pp%gam_v": gamma_v,
+    "bub_pp%gam_g": gamma_g,
+    "bub_pp%M_v": MW_v,
+    "bub_pp%M_g": MW_g,
+    "bub_pp%k_v": k_v * (T0 / (L0 * rho0 * c0 * c0 * c0)),
+    "bub_pp%k_g": k_g * (T0 / (L0 * rho0 * c0 * c0 * c0)),
+    "bub_pp%cp_v": cp_v * (T0 / (c0 * c0)),
+    "bub_pp%cp_g": cp_g * (T0 / (c0 * c0)),
+    "bub_pp%R_v": (R_uni / MW_v) * (T0 / (c0 * c0)),
+    "bub_pp%R_g": (R_uni / MW_g) * (T0 / (c0 * c0)),
     # Viscosity
     "viscous": "T",
     "fluid_pp(1)%Re(1)": 1.0 / (mu_host / (rho0 * c0 * L0)),

src/common/m_constants.fpp

@@ -61,7 +61,7 @@ module m_constants
 
     ! Lagrange bubbles constants
     integer, parameter :: mapCells = 3         !< Number of cells around the bubble where the smoothening function will have effect
-    integer, parameter :: beta_vars(1:3) = [1, 2, 5] !< Indices of variables to communicate for bubble coupling
+    integer, parameter :: beta_vars(1:3) = [1, 2, 5] !< q_beta indices to communicate: 1=void fraction, 2=d(beta)/dt, 5=energy source
     real(wp), parameter :: R_uni = 8314._wp    !< Universal gas constant - J/kmol/K
     integer, parameter :: lag_io_vars = 21     !< Number of variables per particle for MPI_IO
 

src/simulation/m_bubbles_EL.fpp

@@ -1589,6 +1589,7 @@ contains
             cell = fd_number - buff_size
             call s_locate_cell(mtn_pos(k, 1:3, 2), cell, mtn_s(k, 1:3, 2))
         end do
+        $:END_GPU_PARALLEL_LOOP()
 
         call nvtxEndRange ! LAG-BC