this setting is not constistent with the original OGS-PQC setting. correct it.
we could have _J_tot_kin > 0, but _n_xi_kin = 0 case. Therefore we should use _J_tot_kin > 0 here.
same reason as the previous commit.
switch off the activity model. Do not complicate the problem for now.
print out all time steps to see all details.
protection to remove negative values.
1e-9 for nonlinear iterations is lose enough.