var/da/da_par_util/da_transpose_x2y.inc

   1 subroutine da_transpose_x2y (grid)
   2
   3    implicit none
   4
   5    type(domain), intent(inout)               :: grid
   6    integer                                   :: ij, i, j, k
   7
   8    if (trace_use_dull) call da_trace_entry("da_transpose_x2y")
   9
  10 #ifdef DM_PARALLEL
  11 #include "XPOSE_V1_x2y.inc"
  12 #else
  13    !$OMP PARALLEL DO &
  14    !$OMP PRIVATE ( ij, i, j, k )
  15    do ij = 1 , grid%num_tiles
  16       do k = grid%xp%kds, grid%xp%kde
  17          do j = grid%j_start(ij), grid%j_end(ij)
  18             do i = grid%xp%ids, grid%xp%ide
  19                grid%xp % v1y(i,j,k) = grid%xp % v1x(i,j,k)
  20             end do
  21          end do
  22       end do
  23    end do
  24    !$OMP END PARALLEL DO
  25 #endif
  26
  27    if (trace_use_dull) call da_trace_exit("da_transpose_x2y")
  28
  29 end subroutine da_transpose_x2y
  30
  31