first draft
copied fix_momentum, work on the Kokkos View types for the arrays. the next step is parallel_for and parallel_reduce transformations.