*log << "SIM: " << nequileft << " of " << nequi << " equilibration steps left after virtual equilibration." << endl << flush;
progress equiProg(nequileft);
for( int iequi=0; iequi<nequileft; iequi++ ) {
*log << "SIM: " << nequileft << " of " << nequi << " equilibration steps left after virtual equilibration." << endl << flush;
progress equiProg(nequileft);
for( int iequi=0; iequi<nequileft; iequi++ ) {