always reset gases

This commit is contained in:
David Parkhurst 2021-07-14 13:53:28 -06:00
parent 13ec2fcd3e
commit 20281a0e7a

View File

@ -2135,17 +2135,20 @@ Restart:
if (gas_in == FALSE)
continue;
x[i]->moles -= d2;
*use.Get_gas_phase_ptr() = base_gas_phase;
for (size_t g = 0; g < base_phases.size(); g++)
{
*phase_ptrs[g] = base_phases[g];
}
break;
case SS_MOLES:
delta[i] = -d2;
reset();
break;
}
if (use.Get_gas_phase_ptr() != NULL)
{
*use.Get_gas_phase_ptr() = base_gas_phase;
for (size_t g = 0; g < base_phases.size(); g++)
{
*phase_ptrs[g] = base_phases[g];
}
}
}
molalities(TRUE);
if (full_pitzer == TRUE)