bodyforce_mod.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module:  BODYFORCE                                                  C
!  Purpose: Include file for all body force statement functions        C
!                                                                      C
!  Author: M. Syamlal                                 Date:  6-MAR-92  C
!                                                                      C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C

      MODULE bodyforce

      USE constant, only: gravity_x, gravity_y, gravity_z
      IMPLICIT NONE

      CONTAINS

! Body force on gas at i+1/2, j, k
      DOUBLE PRECISION FUNCTION bfx_g(IJK)
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: ijk
      bfx_g = gravity_x
      END FUNCTION bfx_g

! Body force on gas at i, j+1/2, k
      DOUBLE PRECISION FUNCTION bfy_g(IJK)
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: ijk
      bfy_g = gravity_y
      END FUNCTION bfy_g

! Body force on gas at i, j, k+1/2
      DOUBLE PRECISION FUNCTION bfz_g(IJK)
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: ijk
      bfz_g = gravity_z
      END FUNCTION bfz_g

! Body force on solids m at i+1/2, j, k
      DOUBLE PRECISION FUNCTION bfx_s(IJK,M)
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: ijk,m
      bfx_s = gravity_x
      END FUNCTION bfx_s

! Body force on solids m at i, j+1/2, k
      DOUBLE PRECISION FUNCTION bfy_s(IJK,M)
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: ijk,m
      bfy_s = gravity_y
      END FUNCTION bfy_s

! Body force on solids m at i, j, k+1/2
      DOUBLE PRECISION FUNCTION bfz_s(IJK,M)
      IMPLICIT NONE
      INTEGER, INTENT(IN) :: ijk,m
      bfz_s = gravity_z
      END FUNCTION bfz_s

      END MODULE bodyforce