Hi Paz,
That was a long time ago, but there is a simple example in the PHREEQC users' forum:
iphreeqc = actxserver('IPhreeqcCOM.Object'); iphreeqc.LoadDatabase(['whatever path' '\USGS\IPhreeqcCOM 3.3.7-11094\database\phreeqc.dat']);
iphreeqc.ClearAccumulatedLines;
iphreeqc.AccumulateLine ('USER_PUNCH') iphreeqc.AccumulateLine ('-head Spec_Cond'); iphreeqc.AccumulateLine ('PUNCH SC');
iphreeqc.AccumulateLine ('SELECTED_OUTPUT'); iphreeqc.AccumulateLine ('-activities Ca+2 CO3-2'); iphreeqc.AccumulateLine ('-ionic_strength'); iphreeqc.AccumulateLine ('-pH');
iphreeqc.AccumulateLine ('SOLUTION 1 '); iphreeqc.AccumulateLine ('-pH 7 charge'); iphreeqc.AccumulateLine ('-temp 25'); iphreeqc.AccumulateLine ('-units mol/L'); iphreeqc.AccumulateLine (['C(4) ' num2str(c(1))]); iphreeqc.AccumulateLine (['Ca ' num2str(c(2))]); iphreeqc.AccumulateLine (['Cl ' num2str(c(3))]); iphreeqc.AccumulateLine (['N(-3) ' num2str(c(4))]);
try iphreeqc.RunAccumulated; catch break; end
out_PHREEQC = iphreeqc.GetSelectedOutputArray;
The original thread is at the address below.
Best regards,
Leo.
