A specie activited on one proc and not on another is properly transported thanks to the reduce function in YEqn.H as added in reactingFoam solver. However, now the species is well activited in every proc so that it is also written and identified as activated
code added in the inEOA function to give the proportion of the distance to ellispoid center due to temperature and pressure