| 228 |  | catch (OpenMDException& e) { | 
| 229 |  | sprintf(painCave.errMsg, | 
| 230 |  | "%s\n", | 
| 231 | < | e.getMessage().c_str()); | 
| 231 | > | e.what()); | 
| 232 |  | painCave.isFatal = 1; | 
| 233 |  | simError(); | 
| 234 |  | } | 
| 253 |  | std::string mdRawData; | 
| 254 |  | int metaDataBlockStart = -1; | 
| 255 |  | int metaDataBlockEnd = -1; | 
| 256 | – | int i, j; | 
| 256 |  | streamoff mdOffset; | 
| 257 |  | int mdFileVersion; | 
| 258 |  |  | 
| 293 |  | mdFile_.getline(buffer, bufferSize); | 
| 294 |  | ++lineNo; | 
| 295 |  | std::string line = trimLeftCopy(buffer); | 
| 296 | < | i = CaseInsensitiveFind(line, "<OpenMD"); | 
| 297 | < | if (static_cast<size_t>(i) == string::npos) { | 
| 296 | > | std::size_t i = CaseInsensitiveFind(line, "<OpenMD"); | 
| 297 | > | if (i == string::npos) { | 
| 298 |  | // try the older file strings to see if that works: | 
| 299 |  | i = CaseInsensitiveFind(line, "<OOPSE"); | 
| 300 |  | } | 
| 301 |  |  | 
| 302 | < | if (static_cast<size_t>(i) == string::npos) { | 
| 302 | > | if (i == string::npos) { | 
| 303 |  | // still no luck! | 
| 304 |  | sprintf(painCave.errMsg, | 
| 305 |  | "SimCreator: File: %s is not a valid OpenMD file!\n", | 
| 334 |  |  | 
| 335 |  | std::string line = trimLeftCopy(buffer); | 
| 336 |  | if (metaDataBlockStart == -1) { | 
| 337 | < | i = CaseInsensitiveFind(line, "<MetaData>"); | 
| 337 | > | std::size_t i = CaseInsensitiveFind(line, "<MetaData>"); | 
| 338 |  | if (i != string::npos) { | 
| 339 |  | metaDataBlockStart = lineNo; | 
| 340 |  | mdOffset = mdFile_.tellg(); | 
| 341 |  | } | 
| 342 |  | } else { | 
| 343 | < | i = CaseInsensitiveFind(line, "</MetaData>"); | 
| 343 | > | std::size_t i = CaseInsensitiveFind(line, "</MetaData>"); | 
| 344 |  | if (i != string::npos) { | 
| 345 |  | metaDataBlockEnd = lineNo; | 
| 346 |  | } | 
| 373 |  | for (int i = 0; i < metaDataBlockEnd - metaDataBlockStart - 1; ++i) { | 
| 374 |  | mdFile_.getline(buffer, bufferSize); | 
| 375 |  | std::string line = trimLeftCopy(buffer); | 
| 376 | < | j = CaseInsensitiveFind(line, "## Last run using OpenMD Version"); | 
| 377 | < | if (static_cast<size_t>(j) != string::npos) { | 
| 376 | > | std::size_t j = CaseInsensitiveFind(line, | 
| 377 | > | "## Last run using OpenMD Version"); | 
| 378 | > | if (j != string::npos) { | 
| 379 |  | foundVersion = true; | 
| 380 |  | mdRawData += version; | 
| 381 |  | } else { | 
| 785 |  | } | 
| 786 |  | } | 
| 787 |  |  | 
| 788 | < | if (simParams->getOutputElectricField() | simParams->haveElectricField()) { | 
| 788 | > | if (simParams->getOutputElectricField() | | 
| 789 | > | simParams->haveElectricField() | simParams->haveUniformField() | | 
| 790 | > | simParams->haveUniformGradientStrength() | | 
| 791 | > | simParams->haveUniformGradientDirection1() | | 
| 792 | > | simParams->haveUniformGradientDirection2() ) { | 
| 793 |  | storageLayout |= DataStorage::dslElectricField; | 
| 794 |  | } | 
| 795 |  |  | 
| 796 | + | if (simParams->getOutputSitePotential() ) { | 
| 797 | + | storageLayout |= DataStorage::dslSitePotential; | 
| 798 | + | } | 
| 799 | + |  | 
| 800 |  | if (simParams->getOutputFluctuatingCharges()) { | 
| 801 |  | storageLayout |= DataStorage::dslFlucQPosition; | 
| 802 |  | storageLayout |= DataStorage::dslFlucQVelocity; | 
| 833 |  | int beginBendIndex; | 
| 834 |  | int beginTorsionIndex; | 
| 835 |  | int beginInversionIndex; | 
| 836 | + | #ifdef IS_MPI | 
| 837 |  | int nGlobalAtoms = info->getNGlobalAtoms(); | 
| 838 |  | int nGlobalRigidBodies = info->getNGlobalRigidBodies(); | 
| 839 | + | #endif | 
| 840 |  |  | 
| 841 |  | beginAtomIndex = 0; | 
| 842 |  | // The rigid body indices begin immediately after the atom indices: |