scale_filter Module Reference

module FILTER More...

Functions/Subroutines
subroutine	filter_hyperdiff_3d (KA, KS, KE, IA, IS, IE, JA, JS, JE, data, order, nite, limiter_sign)
	Hyper diffusion filter 3D. More...

Detailed Description

module FILTER

Description

Horizontal filter

Author

Team SCALE

Function/Subroutine Documentation

filter_hyperdiff_3d()

subroutine scale_filter::filter_hyperdiff_3d	(	integer, intent(in)	KA,
		integer, intent(in)	KS,
		integer, intent(in)	KE,
		integer, intent(in)	IA,
		integer, intent(in)	IS,
		integer, intent(in)	IE,
		integer, intent(in)	JA,
		integer, intent(in)	JS,
		integer, intent(in)	JE,
		real(rp), dimension(ka,ia,ja), intent(inout)	data,
		integer, intent(in)	order,
		integer, intent(in)	nite,
		real(rp), dimension(ka,ia,ja), intent(in), optional	limiter_sign
	)

Hyper diffusion filter 3D.

Definition at line 252 of file scale_filter.F90.

    use scale_comm_cartesc, only: &
       comm_vars8, &
       comm_wait
    use scale_prc, only: &
       prc_abort
    integer, intent(in) :: KA, KS, KE
    integer, intent(in) :: IA, IS, IE
    integer, intent(in) :: JA, JS, JE
 
    real(RP), intent(inout) :: data(KA,IA,JA)
 
    integer,  intent(in)    :: order
    integer,  intent(in)    :: nite
 
    real(RP), intent(in), optional :: limiter_sign(KA,IA,JA)
 
    real(RP) :: work_data(IA,JA)
    real(RP) :: work_sign(IA,JA)
 
    logical :: flag
 
    integer :: k, i, j
 
    if ( present(limiter_sign) ) then
       flag = .true.
    else
       flag = .false.
    end if
 
    !$acc data copy(data) create(work_data)
    !$acc data copyin(limiter_sign) create(work_sign) if(flag)
 
    do k = ks, ke
 
       !$omp parallel do
       !$acc kernels
       do j = js, je
       do i = is, ie
          work_data(i,j) = data(k,i,j)
       end do
       end do
       !$acc end kernels
       if ( flag ) then
          !$omp parallel do
          !$acc kernels
          do j = js, je
          do i = is, ie
             work_sign(i,j) = limiter_sign(k,i,j)
          end do
          end do
          !$acc end kernels
 
          call filter_hyperdiff_2d( ia, is, ie, ja, js, je, &
                                    work_data(:,:), order, nite, &
                                    limiter_sign = work_sign )
 
       else
 
          call filter_hyperdiff_2d( ia, is, ie, ja, js, je, &
                                    work_data(:,:), order, nite )
       end if
 
 
       !$omp parallel do
       !$acc kernels
       do j = 1, ja
       do i = 1, ia
          data(k,i,j) = work_data(i,j)
       end do
       end do
       !$acc end kernels
 
    end do
 
    !$acc end data
    !$acc end data
 
    return

References scale_prc::prc_abort().

Here is the call graph for this function: