| 101 |  | DefineOptionalParameter(HydroPropFile, "HydroPropFile"); | 
| 102 |  | DefineOptionalParameter(Viscosity, "viscosity"); | 
| 103 |  | DefineOptionalParameter(BeadSize, "beadSize"); | 
| 104 | + | DefineOptionalParameter(FrozenBufferRadius, "frozenBufferRadius"); | 
| 105 | + | DefineOptionalParameter(LangevinBufferRadius, "langevinBufferRadius"); | 
| 106 |  |  | 
| 107 |  | DefineOptionalParameterWithDefaultValue(UsePeriodicBoundaryConditions, "usePeriodicBoundaryConditions", true); | 
| 108 |  | DefineOptionalParameterWithDefaultValue(UseInitalTime, "useInitialTime", false); | 
| 118 |  | DefineOptionalParameterWithDefaultValue(OutputForceVector, "outputForceVector", 0); | 
| 119 |  | DefineOptionalParameterWithDefaultValue(SkinThickness, "skinThickness", 1.0); | 
| 120 |  | DefineOptionalParameterWithDefaultValue(StatFileFormat, "statFileFormat", "TIME|TOTAL_ENERGY|POTENTIAL_ENERGY|KINETIC_ENERGY|TEMPERATURE|PRESSURE|VOLUME|CONSERVED_QUANTITY"); | 
| 121 | + | DefineOptionalParameterWithDefaultValue(UseSphericalBoundaryConditions, "useSphericalBoundaryConditions", false); | 
| 122 |  |  | 
| 123 |  |  | 
| 124 |  | deprecatedKeywords_.insert("nComponents"); | 
| 136 |  |  | 
| 137 |  | CheckParameter(ForceField, isNotEmpty()); | 
| 138 |  | CheckParameter(TargetTemp, isPositive()); | 
| 139 | < | CheckParameter(Ensemble, isEqualIgnoreCase("NVE") || isEqualIgnoreCase("NVT") || isEqualIgnoreCase("NPTi") || isEqualIgnoreCase("NPTf") || isEqualIgnoreCase("NPTxyz") || isEqualIgnoreCase("NPAT")  || isEqualIgnoreCase("LANGEVINDYNAMICS") || isEqualIgnoreCase("NPRT")); | 
| 139 | > | CheckParameter(Ensemble, isEqualIgnoreCase("NVE") || isEqualIgnoreCase("NVT") || isEqualIgnoreCase("NPTi") || isEqualIgnoreCase("NPTf") || isEqualIgnoreCase("NPTxyz") || isEqualIgnoreCase("NPAT")  || isEqualIgnoreCase("LANGEVINDYNAMICS") || isEqualIgnoreCase("LD") || isEqualIgnoreCase("NPRT") || isEqualIgnoreCase("NPGT")); | 
| 140 |  | CheckParameter(Dt, isPositive()); | 
| 141 |  | CheckParameter(RunTime, isPositive()); | 
| 142 |  | CheckParameter(InitialConfig, isNotEmpty()); | 
| 183 |  | CheckParameter(ThermIntOmegaSpringConst, isPositive()); | 
| 184 |  | CheckParameter(DampingAlpha,isNonNegative()); | 
| 185 |  | CheckParameter(SkinThickness, isPositive()); | 
| 186 | < | CheckParameter(Viscosity,isNonNegative()); | 
| 187 | < | CheckParameter(BeadSize,isPositive()); | 
| 186 | > | CheckParameter(Viscosity, isNonNegative()); | 
| 187 | > | CheckParameter(BeadSize, isPositive()); | 
| 188 | > | CheckParameter(FrozenBufferRadius, isPositive()); | 
| 189 | > | CheckParameter(LangevinBufferRadius, isPositive()); | 
| 190 |  | for(std::vector<Component*>::iterator i = components_.begin(); i != components_.end(); ++i) { | 
| 191 |  | if (!(*i)->findMoleculeStamp(moleculeStamps_)) { | 
| 192 |  | std::ostringstream oss; |