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() 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 ! Local Variables: 30 !----------------------------------------------------------------------! 31 ! NONE 32 33 ! Flag for invoking reaction rate calculations. 34 RRATE = .FALSE. 35 IF(NO_OF_RXNS > 0) RRATE = .TRUE. ! automated mass balance 36 IF(USE_RRATES) RRATE = .TRUE. ! legacy hook 37 IF(STIFF_CHEMISTRY) RRATE = .FALSE. ! stiff chemistry solver 38 39 ! Gas phase source terms: 40 R_gp = ZERO ! Rate of species formation 41 RoX_gc = ZERO ! Rate of species consumption (divided X_g) 42 SUM_R_G = ZERO ! Net rate of gas formation/consumption 43 HOR_G = ZERO ! Heat of reaction 44 45 ! Solids phase source terms: 46 R_sp = ZERO ! Rate of species formation 47 RoX_sc = ZERO ! Rate of species consumption (divided X_s) 48 SUM_R_S = ZERO ! Net rate of solids formation/consumption 49 HOR_S = ZERO ! Heat of reaction 50 51 ! Interphase mass transfer. 52 R_PHASE = ZERO 53 54 RETURN 55 END 56