almost working on GPUs

Author: Benjamin Wilfong

src/simulation/m_viscous.f90

@@ -107,12 +107,13 @@ subroutine s_compute_viscous_stress_tensor(q_prim_vf, grad_x_vf, grad_y_vf, grad
             do l = iz%beg, iz%end
                 do k = -1, 1
                     do j = ix%beg, ix%end
+                        print*, j, k, l
     !$acc loop seq
                         do i = 1, num_fluids
                             alpha_rho_visc(i) = q_prim_vf(i)%sf(j, k, l)
                             alpha_visc(i) = q_prim_vf(E_idx + i)%sf(j, k, l)
                         end do
-                        
+
                         if (bubbles) then
                             rho_visc = 0d0
                             gamma_visc = 0d0