/* There are a total of 9 entries in the algebraic variable array. There are a total of 2 entries in each of the rate and state variable arrays. There are a total of 13 entries in the constant variable array. */ /* * VOI is t in component main (second). * CONSTANTS[0] is V_m in component main (J_per_C). * CONSTANTS[1] is RTF in component main (J_per_C). * CONSTANTS[10] is epsilon in component main (dimensionless). * STATES[0] is TA_i in component main (mol_per_m3). * ALGEBRAIC[0] is CO2_o in component main (mol_per_m3). * ALGEBRAIC[4] is CO2_i in component main (mol_per_m3). * ALGEBRAIC[5] is HCO3_i in component main (mol_per_m3). * ALGEBRAIC[1] is HCO3_o in component main (mol_per_m3). * STATES[1] is H_i in component main (mol_per_m3). * CONSTANTS[2] is H_o in component main (mol_per_m3). * CONSTANTS[3] is H_Lim in component main (mol_per_m3). * ALGEBRAIC[6] is M_CO2 in component main (mol_per_m2_s). * ALGEBRAIC[7] is M_HCO3 in component main (mol_per_m2_s). * ALGEBRAIC[8] is M_H in component main (mol_per_m2_s). * ALGEBRAIC[2] is pH_i in component main (dimensionless). * CONSTANTS[11] is pH_o in component main (dimensionless). * CONSTANTS[12] is pH_Lim in component main (dimensionless). * CONSTANTS[4] is P_CO2 in component main (m_per_s). * CONSTANTS[5] is P_HCO3 in component main (m_per_s). * CONSTANTS[6] is k in component main (m_per_s). * CONSTANTS[7] is K_A in component main (mol_per_m3). * CONSTANTS[8] is rho in component main (per_m). * ALGEBRAIC[3] is alpha_i in component main (dimensionless). * CONSTANTS[9] is beta in component main (mol_per_m3). * RATES[0] is d/dt TA_i in component main (mol_per_m3). * RATES[1] is d/dt H_i in component main (mol_per_m3). * There are a total of 3 condition variables. */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { CONSTANTS[0] = -0.057; CONSTANTS[1] = 0.0256796; STATES[0] = 0.0; STATES[1] = 3.981071705534970e-05; CONSTANTS[2] = 1.995262314968879e-05; CONSTANTS[3] = 3.981071705534970e-05; CONSTANTS[4] = 6e-5; CONSTANTS[5] = 5e-9; CONSTANTS[6] = 0.0375; CONSTANTS[7] = 1e-3; CONSTANTS[8] = 8000; CONSTANTS[9] = -26; CONSTANTS[10] = exp(- CONSTANTS[0]/CONSTANTS[1]); CONSTANTS[11] = - arbitrary_log( 0.00100000*CONSTANTS[2], 10); CONSTANTS[12] = - arbitrary_log( 0.00100000*CONSTANTS[3], 10); RATES[0] = 0.1001; RATES[1] = 0.1001; } void computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { resid[0] = RATES[0] - CONSTANTS[8]*(ALGEBRAIC[6]+ALGEBRAIC[7]); resid[1] = RATES[1] - (( - 2.30300*STATES[1])/CONSTANTS[9])*CONSTANTS[8]*(( (1.00000 - ALGEBRAIC[3])*ALGEBRAIC[6] - ALGEBRAIC[3]*ALGEBRAIC[7]) - ALGEBRAIC[8]); } void computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { } void computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[3] = STATES[1]/(STATES[1]+CONSTANTS[7]); ALGEBRAIC[0] = (CONDVAR[0]<0.00000 ? 0.00000 : CONDVAR[1]<0.00000 ? 1.18770 : 0.00000); ALGEBRAIC[4] = ALGEBRAIC[3]*STATES[0]; ALGEBRAIC[6] = CONSTANTS[4]*(ALGEBRAIC[0] - ALGEBRAIC[4]); ALGEBRAIC[5] = (1.00000 - ALGEBRAIC[3])*STATES[0]; ALGEBRAIC[1] = ( CONSTANTS[7]*ALGEBRAIC[0])/CONSTANTS[2]; ALGEBRAIC[7] = ( (( CONSTANTS[5]*CONSTANTS[0])/CONSTANTS[1])*(ALGEBRAIC[1] - ALGEBRAIC[5]*CONSTANTS[10]))/(1.00000 - CONSTANTS[10]); ALGEBRAIC[2] = - arbitrary_log( 0.00100000*STATES[1], 10); ALGEBRAIC[8] = (CONDVAR[2]<0.00000 ? CONSTANTS[6]*(STATES[1] - CONSTANTS[3]) : 0.00000); } void getStateInformation(double* SI) { SI[0] = 1.0; SI[1] = 1.0; } void computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { CONDVAR[0] = VOI - 100.000; CONDVAR[1] = VOI - 2800.00; CONDVAR[2] = ALGEBRAIC[2] - CONSTANTS[12]; }