fixing PR part 4

Author: Diego Vaca

.gitignore

@@ -56,6 +56,8 @@ examples/*/workloads/
 examples/*/run-*/
 examples/*/logs/
 examples/**/*.f90
+!examples/3D_lag_bubbles_shbubcollapse/input/lag_bubbles.dat
+!examples/3D_lag_bubbles_bubblescreen/input/lag_bubbles.dat
 workloads/
 
 benchmarks/*batch/*/

examples/3D_lagrange_bubblescreen/input/lag_bubbles.dat

@@ -152,7 +152,7 @@
     # Lagrangian Bubbles ===========================
      'bubbles_lagrange'                 : 'T',
      'bubble_model'                     : 2,    # Keller-Miksis model
-     'rkck_adap_dt'                     : 'F',  # Activate adaptive time stepper
+     'rkck_adap_dt'                     : 'T',  # Activate adaptive time stepper
      'rkck_tolerance'                   : 1.0e-05,
      'lag_params%nBubs_glb'             : 1,
      'lag_params%solver_approach'       : 2,    # Two-way coupled

examples/3D_lagrange_shbubcollapse/case.py

@@ -0,0 +1 @@
+   0.000000E+00   0.000000E+00    0.000000E+00    0.000000E+00    0.000000E+00    0.000000E+00    0.080000E-01    0.000000E+00

examples/3D_lagrange_shbubcollapse/input/lag_bubbles.dat

@@ -124,7 +124,7 @@ contains
         smearGridz = smearGrid
         if (p == 0) smearGridz = 1
 
-        !$acc parallel loop gang vector default(present) private(l, s_coord, cell, center) copyin(smearGrid, smearGridz)
+        !$acc parallel loop gang vector default(present) private(nodecoord, l, s_coord, cell, center) copyin(smearGrid, smearGridz)
         do l = 1, nBubs
             nodecoord(1:3) = 0
             center(1:3) = 0.0d0
@@ -138,7 +138,7 @@ contains
             strength_vol = volpart
             strength_vel = 4.0d0*pi*lbk_rad(l, 2)**2*lbk_vel(l, 2)
 
-            !$acc loop collapse(3) private(i, j, k, cellaux, nodecoord)
+            !$acc loop collapse(3) private(cellaux, nodecoord)
             do i = 1, smearGrid
                 do j = 1, smearGrid
                     do k = 1, smearGridz