fixup! Prepare for OpenACC porting in mpas_reconstruct_2d

Fix spacing and whitespace for review comments

Author: G. Dylan Dickerson

src/operators/mpas_vector_reconstruction.F

@@ -218,7 +218,7 @@ subroutine mpas_reconstruct_2d(meshPool, u, uReconstructX, uReconstructY, uRecon
     !   temporary arrays needed in the compute procedure
     logical :: includeHalosLocal
     integer, pointer :: nCells_ptr, nVertLevels_ptr
-    integer          :: nCells, nVertLevels
+    integer :: nCells, nVertLevels
     integer, dimension(:,:), pointer :: edgesOnCell
     integer, dimension(:), pointer :: nEdgesOnCell
     integer :: iCell,iEdge, i, k
@@ -249,7 +249,7 @@ subroutine mpas_reconstruct_2d(meshPool, u, uReconstructX, uReconstructY, uRecon
        call mpas_pool_get_dimension(meshPool, 'nCellsSolve', nCells_ptr)
     end if
     call mpas_pool_get_dimension(meshPool, 'nVertLevels', nVertLevels_ptr)
-    ! Dereference pointers so that OpenACC copies them correctly
+    ! Dereference pointers to make sure OpenACC copies the value pointed to implicitly
     nCells = nCells_ptr
     nVertLevels = nVertLevels_ptr
 
@@ -275,7 +275,7 @@ subroutine mpas_reconstruct_2d(meshPool, u, uReconstructX, uReconstructY, uRecon
     do iCell = 1, nCells
       ! initialize the reconstructed vectors
       !$acc loop vector
-      do k=1,nVertLevels
+      do k = 1,nVertLevels
         uReconstructX(k,iCell) = 0.0
         uReconstructY(k,iCell) = 0.0
         uReconstructZ(k,iCell) = 0.0
@@ -284,10 +284,10 @@ subroutine mpas_reconstruct_2d(meshPool, u, uReconstructX, uReconstructY, uRecon
       ! a more efficient reconstruction where rbf_values*matrix_reconstruct
       ! has been precomputed in coeffs_reconstruct
       !$acc loop seq
-      do i=1,nEdgesOnCell(iCell)
+      do i = 1,nEdgesOnCell(iCell)
         iEdge = edgesOnCell(i,iCell)
         !$acc loop vector
-        do k=1,nVertLevels
+        do k = 1,nVertLevels
           uReconstructX(k,iCell) = uReconstructX(k,iCell) &
             + coeffs_reconstruct(1,i,iCell) * u(k,iEdge)
           uReconstructY(k,iCell) = uReconstructY(k,iCell) &
@@ -313,7 +313,7 @@ subroutine mpas_reconstruct_2d(meshPool, u, uReconstructX, uReconstructY, uRecon
         clon = cos(lonCell(iCell))
         slon = sin(lonCell(iCell))
         !$acc loop vector
-        do k=1,nVertLevels
+        do k = 1,nVertLevels
           uReconstructZonal(k,iCell) = -uReconstructX(k,iCell)*slon + &
                                         uReconstructY(k,iCell)*clon
           uReconstructMeridional(k,iCell) = -(uReconstructX(k,iCell)*clon       &
@@ -328,7 +328,7 @@ subroutine mpas_reconstruct_2d(meshPool, u, uReconstructX, uReconstructY, uRecon
       !$acc parallel default(present)
       !$acc loop gang vector collapse(2)
       do iCell = 1, nCells
-        do k=1,nVertLevels
+        do k = 1,nVertLevels
           uReconstructZonal     (k,iCell) = uReconstructX(k,iCell)
           uReconstructMeridional(k,iCell) = uReconstructY(k,iCell)
         end do