1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! Module name: RRATES_INIT(IER) C 4 ! Purpose: Initialize reaction rate arrays C 5 ! C 6 ! Author: Date: C 7 ! Reviewer: Date: C 8 ! C 9 ! Revision Number: C 10 ! Purpose: C 11 ! Author: Date: dd-mmm-yy C 12 ! Reviewer: Date: dd-mmm-yy C 13 ! C 14 ! C 15 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 16 SUBROUTINE RRATES_INIT(IER) 17 18 ! Global Variables: 19 !---------------------------------------------------------------------// 20 use energy, only : HOR_g, HOR_s 21 use param1, only : ZERO 22 use rxns, only : R_gp, RoX_gc, SUM_R_g 23 use rxns, only : R_sp, RoX_sc, SUM_R_s 24 use rxns, only : R_PHASE, NO_OF_RXNS, USE_RRATES, RRATE 25 use stiff_chem, only : STIFF_CHEMISTRY 26 27 implicit none 28 29 ! Passed Variables: 30 !----------------------------------------------------------------------! 31 ! Error index 32 INTEGER, intent(INOUT) :: IER 33 34 ! Local Variables: 35 !----------------------------------------------------------------------! 36 ! NONE 37 38 ! Flag for invoking reaction rate calculations. 39 RRATE = .FALSE. 40 IF(NO_OF_RXNS > 0) RRATE = .TRUE. ! automated mass balance 41 IF(USE_RRATES) RRATE = .TRUE. ! legacy hook 42 IF(STIFF_CHEMISTRY) RRATE = .FALSE. ! stiff chemistry solver 43 44 ! Gas phase source terms: 45 R_gp = ZERO ! Rate of species formation 46 RoX_gc = ZERO ! Rate of species consumption (divided X_g) 47 SUM_R_G = ZERO ! Net rate of gas formation/consumption 48 HOR_G = ZERO ! Heat of reaction 49 50 ! Solids phase source terms: 51 R_sp = ZERO ! Rate of species formation 52 RoX_sc = ZERO ! Rate of species consumption (divided X_s) 53 SUM_R_S = ZERO ! Net rate of solids formation/consumption 54 HOR_S = ZERO ! Heat of reaction 55 56 ! Interphase mass transfer. 57 R_PHASE = ZERO 58 59 RETURN 60 END 61