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 169 of file scale_filter.F90.

References scale_prc::prc_abort().

    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

    do k = ks, ke

       !$omp parallel do
       do j = js, je
       do i = is, ie
          work_data(i,j) = data(k,i,j)
       end do
       end do
       if ( flag ) then
          !$omp parallel do
          do j = js, je
          do i = is, ie
             work_sign(i,j) = limiter_sign(k,i,j)
          end do
          end do

          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
       do j = 1, ja
       do i = 1, ia
          data(k,i,j) = work_data(i,j)
       end do
       end do

    end do


    return

Here is the call graph for this function: