max/iphreeqc
1
0
Fork 0
mirror of https://git.gfz-potsdam.de/naaice/iphreeqc.git synced 2025-12-16 16:44:49 +01:00
Code Issues Packages Projects Releases Wiki Activity

Tony trying to fix v_m in model.cpp

Browse Source 
git-svn-id: svn://136.177.114.72/svn_GW/phreeqc3/trunk@7957 1feff8c3-07ed-0310-ac33-dd36852eb9cd
This commit is contained in:
David L Parkhurst 2013-08-23 19:46:12 +00:00
parent ec4fa77cba
commit 695d6d4e2d
1 changed files with 11 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
model.cpp
View File

@ -2714,22 +2714,21 @@ calc_gas_pressures(void)
V_m = gas_phase_ptr->Get_volume() / gas_phase_ptr->Get_total_moles();
if (V_m < 0.016)
{
V_m = 0.016;
V_m = 0.016;
} else if (V_m > 1e4)
{
V_m = 1e4;
}
if (gas_phase_ptr->Get_v_m() > 0.016)
{
if (V_m < 0.04)
V_m = (9. * gas_phase_ptr->Get_v_m() + V_m) / 10;
else if (V_m < 0.045)
V_m = (6. * gas_phase_ptr->Get_v_m() + V_m) / 7;
else if (V_m < 0.07)
V_m = (3. * gas_phase_ptr->Get_v_m() + V_m) / 4;
else
V_m = (1. * gas_phase_ptr->Get_v_m() + V_m) / 2;
}
if (V_m < 0.02)
V_m = (8. * gas_phase_ptr->Get_v_m() + V_m) / 9;
else if (V_m < 0.03)
V_m = (6. * gas_phase_ptr->Get_v_m() + V_m) / 7;
else if (V_m < 0.05)
V_m = (4. * gas_phase_ptr->Get_v_m() + V_m) / 5;
else if (V_m < 0.07)
V_m = (2. * gas_phase_ptr->Get_v_m() + V_m) / 3;
else
V_m = (1. * gas_phase_ptr->Get_v_m() + V_m) / 2;
if (iterations > 99)
{
//V_m *= 1; /* debug */