qmomk_set_bc.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: QMOMK_SET_BC                                           C
!  Author: Alberto Passalacqua                        Date:            C
!  Reviewer:                                          Date:            C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C

SUBROUTINE qmomk_set_bc

!-----------------------------------------------
!     M o d u l e s
!-----------------------------------------------
  USE param
  USE param1
  USE constant
  USE physprop
  USE fldvar
  USE geometry
  USE compar
  USE indices
  USE bc
  USE qmom_kinetic_equation
  USE qmomk_quadrature
  USE qmomk_bc
  USE functions

  IMPLICIT NONE

  INTEGER :: L
  INTEGER :: I1, I2, J1, J2, K1, K2

  DO l = 1, dimension_bc
      IF (bc_defined(l)) THEN

        i1 = bc_i_w(l)
        i2 = bc_i_e(l)
        j1 = bc_j_s(l)
        j2 = bc_j_n(l)
        k1 = bc_k_b(l)
        k2 = bc_k_t(l)

        ! Wall BC's
        IF (bc_type_enum(l)==free_slip_wall .OR. bc_type_enum(l)==no_slip_wall&
            .OR. bc_type_enum(l)==par_slip_wall) THEN
          IF (qmomk_wall_bc_type == 'SPECULAR_REFLECTIVE') THEN
            CALL qmomk_reflective_wall_bc(l, i1, i2, j1, j2, k1, k2, .false.)
          ENDIF
        ! Outlet BC's
        ELSEIF (bc_type_enum(l) == mass_outflow .OR. bc_type_enum(l)==p_outflow&
            .OR. bc_type_enum(l)==outflow) THEN
          CALL qmomk_outlet_bc(l, i1, i2, j1, j2, k1, k2, .false.)
        ! Inlet BC's
        ELSEIF (bc_type_enum(l) == mass_inflow) THEN
          CALL qmomk_inlet_bc(l, .false.)
        ! Cyclic BC's
        ELSEIF (cyclic) THEN
          CALL qmomk_cyclic_bc(.false.)
        END IF
      END IF
  END DO

END SUBROUTINE qmomk_set_bc