I am trying to work out the indexing for a 3D complex array in Cuda fortran. I\'m setting up the kernels via:
TYPE(dim3) :: block = dim3(Nx/blockx,Ny/blocky,Nz/bl