R2d_mod.F90

Go to the documentation of this file.

module r2d_mod

  use real_type, only: rp
  use basic_tools
  !$ use OMP_LIB

  implicit none
  private

  !       %----------------------------------------%
  !       |                                        |
  !       |          PUBLIC DATA                   |
  !       |                                        |
  !       %----------------------------------------%
  public :: copy 

  !       %----------------------------------------%
  !       |                                        |
  !       |       GENERIc SUBROUTINES              |
  !       |                                        |
  !       %----------------------------------------%

  interface copy
     module procedure r2d_copy    !! TESTED
  end interface copy


contains


  subroutine r2d_copy(y, x)
    real(RP), dimension(:,:), intent(out) :: y
    real(RP), dimension(:,:), intent(in)  :: x

    integer :: ii, n1, n2

#if DBG>0
    if ( size(x,1) < size(y,1)) call quit(&
         &"R2d_mod: R2d_copy: 1" )
    if ( size(x,2) < size(y,2)) call quit(&
         &"R2d_mod: R2d_copy: 2" )
#endif

    n1 = size(y,1)
    n2 = size(y,2)
    !$OMP PARALLEL 
    !$OMP DO
    do ii=1, n2
       y(1:n1, ii) = x(1:n1, ii)
    end do
    !$OMP END DO
    !$OMP END PARALLEL
    
  end subroutine r2d_copy


end module r2d_mod