56 |
|
#include "brains/SimCreator.hpp" |
57 |
|
#include "brains/SimSnapshotManager.hpp" |
58 |
|
#include "io/DumpReader.hpp" |
59 |
< |
#include "UseTheForce/ForceFieldFactory.hpp" |
59 |
> |
#include "brains/ForceField.hpp" |
60 |
|
#include "utils/simError.h" |
61 |
|
#include "utils/StringUtils.hpp" |
62 |
|
#include "math/SeqRandNumGen.hpp" |
371 |
|
metaDataBlockStart + 1); |
372 |
|
|
373 |
|
//create the force field |
374 |
< |
ForceField * ff = ForceFieldFactory::getInstance()->createForceField(simParams->getForceField()); |
374 |
> |
ForceField * ff = new ForceField(simParams->getForceField()); |
375 |
|
|
376 |
|
if (ff == NULL) { |
377 |
|
sprintf(painCave.errMsg, |
751 |
|
if (simParams->getOutputParticlePotential()) { |
752 |
|
storageLayout |= DataStorage::dslParticlePot; |
753 |
|
} |
754 |
+ |
|
755 |
+ |
if (simParams->havePrintHeatFlux()) { |
756 |
+ |
if (simParams->getPrintHeatFlux()) { |
757 |
+ |
storageLayout |= DataStorage::dslParticlePot; |
758 |
+ |
} |
759 |
+ |
} |
760 |
+ |
|
761 |
|
if (simParams->getOutputElectricField()) { |
762 |
|
storageLayout |= DataStorage::dslElectricField; |
763 |
|
} |