[6610] | 1 | |
---|
| 2 | #include <inca_define.h> |
---|
| 3 | |
---|
| 4 | SUBROUTINE REAC_FLX(base_sol, reaction_rates, invariants, reacflux) |
---|
| 5 | ! Sophie SZOPA, LSCE, 2016. |
---|
| 6 | |
---|
| 7 | USE INCA_DIM |
---|
| 8 | USE SPECIES_NAMES |
---|
| 9 | USE RATE_INDEX_MOD |
---|
| 10 | IMPLICIT NONE |
---|
| 11 | |
---|
| 12 | !-------------------------------------------------------------------- |
---|
| 13 | ! ... Dummy args |
---|
| 14 | !-------------------------------------------------------------------- |
---|
| 15 | REAL, DIMENSION(PLNPLV,RXNCNT), INTENT(out) :: reacflux |
---|
| 16 | REAL, DIMENSION(PLNPLV,RXNCNT), INTENT(in) :: reaction_rates |
---|
| 17 | REAL, DIMENSION(PLNPLV,NFS), INTENT(in) :: invariants |
---|
| 18 | REAL, DIMENSION(PLNPLV,PCNST), INTENT(in) :: base_sol |
---|
| 19 | |
---|
| 20 | |
---|
| 21 | |
---|
| 22 | reacflux(:,jN2O) = invariants(:,INDEXM) * & |
---|
| 23 | reaction_rates(:,jN2O) & |
---|
| 24 | * base_sol(:,id_N2O) |
---|
| 25 | |
---|
| 26 | reacflux(:,jMCF) = invariants(:,INDEXM) * & |
---|
| 27 | reaction_rates(:,jMCF) & |
---|
| 28 | * base_sol(:,id_MCF) |
---|
| 29 | |
---|
| 30 | reacflux(:,kCOOH) = invariants(:,INDEXM) * & |
---|
| 31 | reaction_rates(:,kCOOH) & |
---|
| 32 | * base_sol(:,id_CO) & |
---|
| 33 | * invariants(:,inv_OH) |
---|
| 34 | |
---|
| 35 | reacflux(:,kRn222) = invariants(:,INDEXM) * & |
---|
| 36 | reaction_rates(:,kRn222) & |
---|
| 37 | * base_sol(:,id_Rn222) |
---|
| 38 | |
---|
| 39 | reacflux(:,kMCFOH) = invariants(:,INDEXM) * & |
---|
| 40 | reaction_rates(:,kMCFOH) & |
---|
| 41 | * base_sol(:,id_MCF) & |
---|
| 42 | * invariants(:,inv_OH) |
---|
| 43 | |
---|
| 44 | reacflux(:,kBe7) = invariants(:,INDEXM) * & |
---|
| 45 | reaction_rates(:,kBe7) & |
---|
| 46 | * base_sol(:,id_Be7) |
---|
| 47 | |
---|
| 48 | reacflux(:,kN2OO1D2NO) = invariants(:,INDEXM) * & |
---|
| 49 | reaction_rates(:,kN2OO1D2NO) & |
---|
| 50 | * base_sol(:,id_N2O) & |
---|
| 51 | * invariants(:,inv_O1D) |
---|
| 52 | |
---|
| 53 | reacflux(:,kN2OO1DN2) = invariants(:,INDEXM) * & |
---|
| 54 | reaction_rates(:,kN2OO1DN2) & |
---|
| 55 | * base_sol(:,id_N2O) & |
---|
| 56 | * invariants(:,inv_O1D) |
---|
| 57 | |
---|
| 58 | reacflux(:,kCH4OH) = invariants(:,INDEXM) * & |
---|
| 59 | reaction_rates(:,kCH4OH) & |
---|
| 60 | * base_sol(:,id_CH4) & |
---|
| 61 | * invariants(:,inv_OH) |
---|
| 62 | |
---|
| 63 | reacflux(:,kCH4O1DCH3O2) = invariants(:,INDEXM) * & |
---|
| 64 | reaction_rates(:,kCH4O1DCH3O2) & |
---|
| 65 | * base_sol(:,id_CH4) & |
---|
| 66 | * invariants(:,inv_O1D) |
---|
| 67 | |
---|
| 68 | reacflux(:,kCH4O1DH2) = invariants(:,INDEXM) * & |
---|
| 69 | reaction_rates(:,kCH4O1DH2) & |
---|
| 70 | * base_sol(:,id_CH4) & |
---|
| 71 | * invariants(:,inv_O1D) |
---|
| 72 | |
---|
| 73 | reacflux(:,kCH4O1DCH3OH) = invariants(:,INDEXM) * & |
---|
| 74 | reaction_rates(:,kCH4O1DCH3OH) & |
---|
| 75 | * base_sol(:,id_CH4) & |
---|
| 76 | * invariants(:,inv_O1D) |
---|
| 77 | |
---|
| 78 | |
---|
| 79 | END SUBROUTINE REAC_FLX |
---|