apparently cray compilers use the device copies of integer variables for loop bounds

Author: Daniel Vickers

src/simulation/m_ib_patches.fpp

@@ -391,7 +391,6 @@ contains
         real(wp)                           :: corner_distance       !< Equation of state parameters
         real(wp), dimension(1:3)           :: xy_local              !< x and y coordinates in local IB frame
         real(wp), dimension(1:2)           :: length, center        !< x and y coordinates in local IB frame
-        real(wp), dimension(1:3,1:3)       :: inverse_rotation
 
         ! Transferring the rectangle's centroid and length information
 

src/simulation/m_ibm.fpp

@@ -76,6 +76,11 @@ contains
 
         call nvtxStartRange("SETUP-IBM-MODULE")
 
+        ! GPU routines require updated cell centers
+        $:GPU_UPDATE(device='[num_ibs, num_gbl_ibs, x_cc, y_cc, dx, dy, x_domain, y_domain, ib_bc_x%beg, ib_bc_y%beg]')
+        if (p /= 0) then
+            $:GPU_UPDATE(device='[z_cc, dz, z_domain, ib_bc_z%beg]')
+        end if
         $:GPU_UPDATE(device='[patch_ib(1:num_ibs)]')
 
         ! do all set up for moving immersed boundaries
@@ -98,11 +103,6 @@ contains
         end if
 #endif
 
-        ! GPU routines require updated cell centers
-        $:GPU_UPDATE(device='[num_ibs, num_gbl_ibs, x_cc, y_cc, dx, dy, x_domain, y_domain, ib_bc_x%beg, ib_bc_y%beg]')
-        if (p /= 0) then
-            $:GPU_UPDATE(device='[z_cc, dz, z_domain, ib_bc_z%beg]')
-        end if
         call s_update_ib_lookup()
 
         ! recompute the new ib_patch locations