fixed fatal bug in how threads return if outside of the grid
ported Darcy routines to CUDA, compiles, needs testing