Isolating bug 3

Author: Diego Vaca

src/simulation/m_bubbles_EL.fpp

@@ -825,7 +825,7 @@ contains
         real(wp) :: charvol, charpres, charvol2, charpres2
         integer, dimension(3) :: cellaux
         integer :: i, j, k
-        integer :: smearGrid, smearGridz
+        integer :: mapCells_pinf, smearGrid, smearGridz
         logical :: celloutside
 
         scoord = mtn_s(bub_id, 1:3, 2)
@@ -912,9 +912,14 @@ contains
 
         else if (lag_params%cluster_type >= 2) then
             ! Bubble dynamic closure from Maeda and Colonius (2018)
+            if (lag_params%smooth_type == 1) then
+                mapCells_pinf = mapCells
+            else
+                stop "lag_params%cluster_type: 2 requires lag_params%smooth_type: 1."
+            end if
 
             ! Include the cell that contains the bubble (mapCells+1+mapCells)
-            smearGrid = mapCells - (-mapCells) + 1
+            smearGrid = mapCells_pinf - (-mapCells_pinf) + 1
             smearGridz = smearGrid
             if (p == 0) smearGridz = 1
 
@@ -988,9 +993,6 @@ contains
             f_pinfl = charpres2/charvol2
             vol = charvol
             dc = (3._wp*abs(vol)/(4._wp*pi))**(1._wp/3._wp)
-        else
-
-            print *, "Check cluterflag. Exiting."
 
         end if