#include <inca_define.h>

SUBROUTINE REAC_FLX(base_sol, reaction_rates, invariants, reacflux)
  ! Sophie SZOPA, LSCE, 2016.

  USE INCA_DIM
  USE SPECIES_NAMES
  USE RATE_INDEX_MOD
  IMPLICIT NONE

  !--------------------------------------------------------------------
  !     ... Dummy args                                                          
  !--------------------------------------------------------------------
  REAL, DIMENSION(PLNPLV,RXNCNT), INTENT(out) :: reacflux
  REAL, DIMENSION(PLNPLV,RXNCNT), INTENT(in)  ::  reaction_rates
  REAL, DIMENSION(PLNPLV,NFS), INTENT(in)     :: invariants
  REAL, DIMENSION(PLNPLV,PCNST), INTENT(in)   ::  base_sol



  reacflux(:,jN2O) = invariants(:,INDEXM) *  &
       reaction_rates(:,jN2O)    &
       * base_sol(:,id_N2O)

  reacflux(:,jMCF) = invariants(:,INDEXM) *  &
       reaction_rates(:,jMCF)    &
       * base_sol(:,id_MCF)

  reacflux(:,kCOOH) = invariants(:,INDEXM) *  &
       reaction_rates(:,kCOOH)    &
       * base_sol(:,id_CO)   &
       * invariants(:,inv_OH)

  reacflux(:,kRn222) = invariants(:,INDEXM) *  &
       reaction_rates(:,kRn222)    &
       * base_sol(:,id_Rn222)

  reacflux(:,kMCFOH) = invariants(:,INDEXM) *  &
       reaction_rates(:,kMCFOH)    &
       * base_sol(:,id_MCF)   &
       * invariants(:,inv_OH)

  reacflux(:,kBe7) = invariants(:,INDEXM) *  &
       reaction_rates(:,kBe7)    &
       * base_sol(:,id_Be7)

  reacflux(:,kN2OO1D2NO) = invariants(:,INDEXM) *  &
       reaction_rates(:,kN2OO1D2NO)    &
       * base_sol(:,id_N2O)   &
       * invariants(:,inv_O1D)

  reacflux(:,kN2OO1DN2) = invariants(:,INDEXM) *  &
       reaction_rates(:,kN2OO1DN2)    &
       * base_sol(:,id_N2O)   &
       * invariants(:,inv_O1D)

  reacflux(:,kCH4OH) = invariants(:,INDEXM) *  &
       reaction_rates(:,kCH4OH)    &
       * base_sol(:,id_CH4)   &
       * invariants(:,inv_OH)

  reacflux(:,kCH4O1DCH3O2) = invariants(:,INDEXM) *  &
       reaction_rates(:,kCH4O1DCH3O2)    &
       * base_sol(:,id_CH4)   &
       * invariants(:,inv_O1D)

  reacflux(:,kCH4O1DH2) = invariants(:,INDEXM) *  &
       reaction_rates(:,kCH4O1DH2)    &
       * base_sol(:,id_CH4)   &
       * invariants(:,inv_O1D)

  reacflux(:,kCH4O1DCH3OH) = invariants(:,INDEXM) *  &
       reaction_rates(:,kCH4O1DCH3OH)    &
       * base_sol(:,id_CH4)   &
       * invariants(:,inv_O1D)


END SUBROUTINE REAC_FLX