r2 on GPU

Author: conraddelgado

examples/2D_ibm_circle_periodic/case.py

@@ -35,7 +35,8 @@
 Nt = 5000
 t_save = 20
 
-Nx = 511  # to fully resolve requires ~ 40-60 cells across sphere diameter
+# to fully resolve requires 40-60 cells across sphere diameter
+Nx = 511
 Ny = Nx
 Nz = 0
 

examples/3D_ibm_sphere_periodic/case.py

@@ -35,7 +35,8 @@
 Nt = 2000
 t_save = Nt // 250
 
-Nx = 127  # to fully resolve requires ~ 40-60 cells across sphere diameter
+# to fully resolve requires 40-60 cells across sphere diameter
+Nx = 127
 Ny = Nx
 Nz = Ny
 

src/common/m_ib_patches.fpp

@@ -666,7 +666,7 @@ contains
             $:END_GPU_PARALLEL_LOOP()
         else
             $:GPU_PARALLEL_LOOP(private='[i,j,k,cart_y,cart_z]', copy='[ib_markers_sf]',&
-                      & copyin='[patch_id,center,radius]', collapse=3)
+                      & copyin='[patch_id,center,r2]', collapse=3)
             do k = 0, p
                 do j = 0, n
                     do i = 0, m