set_ro_g.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Subroutine: SET_RO_g                                                C
!  Purpose: Initialize the gas densities                               C
!                                                                      C
!  Author: M. Syamlal                                 Date: 21-JAN-92  C
!  Reviewer:M. Syamlal, S. Venkatesan, P. Nicoletti,  Date: 29-JAN-92  C
!           W. Rogers                                                  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
!  Variables referenced: MW_MIX_g, P_g, T_g, EP_g, RO_g0               C
!  Variables modified: RO_g, ROP_g                                     C
!  Local variables: IJK                                                C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C

      SUBROUTINE set_ro_g

!-----------------------------------------------
! Modules
!-----------------------------------------------
      USE compar
      USE constant
      USE eos, only: eosg
      USE fldvar
      USE functions
      USE geometry
      USE indices
      USE parallel
      USE param
      USE param1
      USE physprop
      IMPLICIT NONE
!-----------------------------------------------
! Local variables
!-----------------------------------------------
      INTEGER :: IJK
!-----------------------------------------------
      IF (ro_g0 == undefined) THEN   ! compressible case

!!$omp parallel do private(IJK)
         DO ijk = ijkstart3, ijkend3
! calculate ro_g and rop_g in all fluid and flow boundary cells
            IF (.NOT.wall_at(ijk)) THEN
! set_bc0 will have already defined ro_g and rop_g in MI and PI
! boundary cells (redundant-remove in set_bc0?)
               ro_g(ijk) = eosg(mw_mix_g(ijk),p_g(ijk),t_g(ijk))
               rop_g(ijk) = ep_g(ijk)*ro_g(ijk)
            ENDIF
         ENDDO

      ELSE   ! incompressible case

!!$omp   parallel do private(ijk)
         DO ijk = ijkstart3, ijkend3
            IF (.NOT.wall_at(ijk)) THEN
! assign ro_g and calculate rop_g in all fluid and flow boundary cells
! set_constprop will have already defined ro_g in fluid and flow
! boundary cells (redundant- remove here?)
               ro_g(ijk) = ro_g0
! set_bc0 will have already defined rop_g in MI and PI boundary cells
! (redundant-remove in set_bc0?)
               rop_g(ijk) = ep_g(ijk)*ro_g(ijk)
            ENDIF
         ENDDO
      ENDIF

      RETURN
      END SUBROUTINE set_ro_g