equilibrium_data.cpp
36b376748b48reactivetransport

equilibrium_data.cpp
View Options

			/*-------------------------------------------------------

			- Module : specmicp
			- File : equilibrium_data.cpp
			- Author : Fabien Georget

			Copyright (c) 2014, Fabien Georget, Princeton University

			---------------------------------------------------------*/

			#include "equilibrium_data.hpp"

			namespace specmicp {


			//! return log(IAP) for mineral m
			double EquilibriumState::logIAP(int m)
			{
			double logiap = 0.0;
			for (int i=1; i<m_data->nb_component; ++i) logiap += m_main_variables(i)*m_data->nu_mineral(m, i);
			return logiap;
			}

			//! return log(IAP) for mineral m
			double EquilibriumState::logIAP_kinetic(int m)
			{
			double logiap = 0.0;
			for (int i=1; i<m_data->nb_component; ++i) logiap += m_main_variables(i)*m_data->nu_mineral_kinetic(m, i);
			return logiap;
			}

			} // end namespace specmicp