425 |
|
//find the storage layout |
426 |
|
|
427 |
|
int storageLayout = computeStorageLayout(info); |
428 |
– |
|
429 |
– |
cerr << "computed Storage Layout = " << storageLayout << "\n"; |
428 |
|
|
429 |
|
//allocate memory for DataStorage(circular reference, need to |
430 |
|
//break it) |
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 |
|
} |