Isolating bug 2

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 :: mapCells_pinf, smearGrid, smearGridz
+        integer :: smearGrid, smearGridz
         logical :: celloutside
 
         scoord = mtn_s(bub_id, 1:3, 2)
@@ -913,15 +913,8 @@ contains
         else if (lag_params%cluster_type >= 2) then
             ! Bubble dynamic closure from Maeda and Colonius (2018)
 
-            ! Range of cells included in Omega
-            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_pinf - (-mapCells_pinf) + 1
+            smearGrid = mapCells - (-mapCells) + 1
             smearGridz = smearGrid
             if (p == 0) smearGridz = 1