Commit 45c7b31e by Francois Gygi

added stepper kinetic energy to econst

git-svn-id: http://qboxcode.org/svn/qb/trunk@1552 cba15fb0-1239-40c8-b417-11db7ca47a34
parent 860a212d
......@@ -244,8 +244,10 @@ void CPSampleStepper::step(int niter)
cout << " <temp_ion> " << mdionic_stepper->temp() << " </temp_ion>\n";
cout << " <eta_ion> " << mdionic_stepper->eta() << " </eta_ion>\n";
}
cout << " <econst> " << energy+ekin_ion+ekin_e << " </econst>\n";
cout << " <ekin_ec> " << energy+ekin_ion+2*ekin_e << " </ekin_ec>\n";
double econst = energy + ekin_ion + ekin_e +
mdionic_stepper->ekin_stepper();
cout << " <econst> " << econst << " </econst>\n";
cout << " <ekin_ec> " << econst + ekin_e << " </ekin_ec>\n";
}
if ( compute_stress )
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment