check_plane.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: CHECK_PLANE                                            C
!  Purpose: make sure the flow boundary condition or internal surface  C
!           is a plane                                                 C
!                                                                      C
!  Author: P. Nicoletti                               Date: 10-DEC-91  C
!  Reviewer: M.SYAMLAL, W.ROGERS, P.NICOLETTI         Date: 24-JAN-92  C
!                                                                      C
!  Revision Number:                                                    C
!  Purpose:                                                            C
!  Author:                                            Date: dd-mmm-yy  C
!  Reviewer:                                          Date: dd-mmm-yy  C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C

      SUBROUTINE check_plane(X_CONSTANT, Y_CONSTANT, Z_CONSTANT, BC, NAME)

!-----------------------------------------------
! Modules
!-----------------------------------------------
         USE compar
         USE exit, only: mfix_exit
         USE funits

      IMPLICIT NONE

!-----------------------------------------------
! Dummy arguments
!-----------------------------------------------
! surface indicators
      LOGICAL, INTENT(IN) :: X_CONSTANT,Y_CONSTANT,Z_CONSTANT
! boundary condition or internal surface index
      INTEGER, INTENT(IN) ::  BC
! BC or IS
      CHARACTER(LEN=2), INTENT(IN) :: NAME
!-----------------------------------------------
! Local variables
!-----------------------------------------------
! number of directions that are not constant (must equal 2)
      INTEGER :: N
!-----------------------------------------------


! number of directions that are not constant (must equal 2)
      n = 3
      IF (x_constant) n = n - 1
      IF (y_constant) n = n - 1
      IF (z_constant) n = n - 1

      IF (n /= 2) THEN
         IF(dmp_log)WRITE (unit_log, 1000) name, bc
         call mfix_exit(mype)
      ENDIF

      RETURN

 1000 FORMAT(/70('*')//' From: CHECK_PLANE',/'Message: ',a,' No ',i3,&
         ' is not a plane',/70('*')/)
      END SUBROUTINE check_plane