zero_norm_vel.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Subroutine: ZERO_NORM_VEL                                           C
!  Purpose: Set the velocity component normal to a wall to zero        C
!                                                                      C
!  Author: M. Syamlal                                 Date: 14-MAY-92  C
!  Reviewer: S. Venkatesan                            Date: 11-DEC-92  C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced:                                               C
!  Variables modified:                                                 C
!  Local variables:                                                    C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C

      SUBROUTINE zero_norm_vel

!-----------------------------------------------
! Modules
!-----------------------------------------------
      USE param
      USE param1
      USE parallel
      USE geometry
      USE physprop
      USE fldvar
      USE indices
      USE is
      USE compar
      USE discretelement
      USE mfix_pic
      USE functions

      IMPLICIT NONE
!-----------------------------------------------
! Local variables
!-----------------------------------------------
!
      INTEGER :: ISV
! Indicies
      INTEGER :: IJK, IMJK, IJMK, IJKM
      INTEGER :: M

!!$omp  parallel do private( IMJK, IJMK, IJKM)

      DO ijk = ijkstart3, ijkend3

         IF (.NOT.wall_at(ijk)) THEN
            IF (ip_at_e(ijk)) u_g(ijk) = zero
            IF (ip_at_n(ijk)) v_g(ijk) = zero
            IF (ip_at_t(ijk)) w_g(ijk) = zero
         ELSE
            imjk = im_of(ijk)
            ijmk = jm_of(ijk)
            ijkm = km_of(ijk)
            u_g(ijk) = zero
            v_g(ijk) = zero
            w_g(ijk) = zero
            IF (.NOT.(cyclic_at(ijk) .AND. (i_of(ijk)==imax2 .OR. &
                i_of(ijk)==imax3))) u_g(imjk) = zero
            IF (.NOT.(cyclic_at(ijk) .AND. (j_of(ijk)==jmax2 .OR. &
                j_of(ijk)==jmax3))) v_g(ijmk) = zero
            IF (.NOT.(cyclic_at(ijk) .AND. (k_of(ijk)==kmax2 .OR. &
                k_of(ijk)==kmax3))) w_g(ijkm) = zero
         ENDIF
      ENDDO   ! end do (ijk=ijkstart3,ijkend3)


      IF (.NOT.discrete_element .OR. des_continuum_hybrid .OR.mppic) THEN
         DO m = 1, mmax
!!$omp  parallel do private( ISV,  IMJK, IJMK, IJKM)
            DO ijk = ijkstart3, ijkend3

               IF (.NOT.wall_at(ijk)) THEN
                  IF (ip_at_e(ijk)) THEN
                     u_s(ijk,m) = zero
                  ELSEIF (sip_at_e(ijk)) THEN
                     isv = flag_e(ijk) - 1000
                     u_s(ijk,m) = is_vel_s(isv,m)
                  ENDIF
                  IF (ip_at_n(ijk)) THEN
                     v_s(ijk,m) = zero
                  ELSEIF (sip_at_n(ijk)) THEN
                     isv = flag_n(ijk) - 1000
                     v_s(ijk,m) = is_vel_s(isv,m)
                  ENDIF
                  IF (ip_at_t(ijk)) THEN
                     w_s(ijk,m) = zero
                  ELSEIF (sip_at_t(ijk)) THEN
                     isv = flag_t(ijk) - 1000
                     w_s(ijk,m) = is_vel_s(isv,m)
                  ENDIF
               ELSE
                  imjk = im_of(ijk)
                  ijmk = jm_of(ijk)
                  ijkm = km_of(ijk)
                  u_s(ijk,m) = zero
                  v_s(ijk,m) = zero
                  w_s(ijk,m) = zero
                  IF (.NOT.(cyclic_at(ijk) .AND. (i_of(ijk)==imax2 .OR. &
                      i_of(ijk)==imax3))) u_s(imjk,m) = zero
                  IF (.NOT.(cyclic_at(ijk) .AND. (j_of(ijk)==jmax2 .OR. &
                      j_of(ijk)==jmax3))) v_s(ijmk,m) = zero
                  IF (.NOT.(cyclic_at(ijk) .AND. (k_of(ijk)==kmax2 .OR. &
                      k_of(ijk)==kmax3))) w_s(ijkm,m) = zero
               ENDIF
            ENDDO   ! end do (ijk=ijkstart3,ijkend3)
         ENDDO   ! end do (m=1,mmax)
      ENDIF   ! endif (.not.discrete_element)

      RETURN

      END SUBROUTINE zero_norm_vel