rrates_init.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: RRATES_INIT(IER)                                       C
!  Purpose: Initialize reaction rate arrays                            C
!                                                                      C
!  Author:                                            Date:            C
!  Reviewer:                                          Date:            C
!                                                                      C
!  Revision Number:                                                    C
!  Purpose:                                                            C
!  Author:                                            Date: dd-mmm-yy  C
!  Reviewer:                                          Date: dd-mmm-yy  C
!                                                                      C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
      SUBROUTINE rrates_init()

! Global Variables:
!---------------------------------------------------------------------//
      use energy,     only : hor_g, hor_s
      use param1,     only : zero
      use rxns,       only : r_gp, rox_gc, sum_r_g
      use rxns,       only : r_sp, rox_sc, sum_r_s
      use rxns,       only : r_phase, no_of_rxns, use_rrates, rrate
      use stiff_chem, only : stiff_chemistry

      implicit none

! Local Variables:
!----------------------------------------------------------------------!
! NONE

! Flag for invoking reaction rate calculations.
      rrate = .false.
      IF(no_of_rxns > 0)  rrate = .true.   ! automated mass balance
      IF(use_rrates)      rrate = .true.   ! legacy hook
      IF(stiff_chemistry) rrate = .false.  ! stiff chemistry solver

! Gas phase source terms:
      r_gp    = zero  ! Rate of species formation
      rox_gc  = zero  ! Rate of species consumption (divided X_g)
      sum_r_g = zero  ! Net rate of gas formation/consumption
      hor_g   = zero  ! Heat of reaction

! Solids phase source terms:
      r_sp    = zero  ! Rate of species formation
      rox_sc  = zero  ! Rate of species consumption (divided X_s)
      sum_r_s = zero  ! Net rate of solids formation/consumption
      hor_s   = zero  ! Heat of reaction

! Interphase mass transfer.
      r_phase = zero

      RETURN
      END