| 57 |  | DefineOptionalParameter(Ensemble, "ensemble"); | 
| 58 |  | DefineOptionalParameter(Dt, "dt"); | 
| 59 |  | DefineOptionalParameter(RunTime, "runTime"); | 
| 60 | < | DefineOptionalParameter(InitialConfig, "initialConfig"); | 
| 60 | > | //DefineOptionalParameter(InitialConfig, "initialConfig"); | 
| 61 |  | DefineOptionalParameter(FinalConfig, "finalConfig"); | 
| 62 |  | DefineOptionalParameter(SampleTime, "sampleTime"); | 
| 63 |  | DefineOptionalParameter(ResetTime, "resetTime"); | 
| 99 |  | DefineOptionalParameter(CutoffPolicy, "cutoffPolicy"); | 
| 100 |  | DefineOptionalParameter(SwitchingFunctionType, "switchingFunctionType"); | 
| 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); | 
| 109 |  | DefineOptionalParameterWithDefaultValue(UseIntialExtendedSystemState, "useInitialExtendedSystemState", 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 | + | DefineOptionalParameterWithDefaultValue(AccumulateBoxDipole, "accumulateBoxDipole", false); | 
| 123 |  |  | 
| 124 |  |  | 
| 125 |  | deprecatedKeywords_.insert("nComponents"); | 
| 126 |  | deprecatedKeywords_.insert("nZconstraints"); | 
| 127 | + | deprecatedKeywords_.insert("initialConfig"); | 
| 128 |  |  | 
| 129 |  | } | 
| 130 |  |  | 
| 138 |  |  | 
| 139 |  | CheckParameter(ForceField, isNotEmpty()); | 
| 140 |  | CheckParameter(TargetTemp, isPositive()); | 
| 141 | < | CheckParameter(Ensemble, isEqualIgnoreCase("NVE") || isEqualIgnoreCase("NVT") || isEqualIgnoreCase("NPTi") || isEqualIgnoreCase("NPTf") || isEqualIgnoreCase("NPTxyz") || isEqualIgnoreCase("NPAT")  || isEqualIgnoreCase("LANGEVINDYNAMICS") ); | 
| 141 | > | CheckParameter(Ensemble, isEqualIgnoreCase("NVE") || isEqualIgnoreCase("NVT") || isEqualIgnoreCase("NPTi") || isEqualIgnoreCase("NPTf") || isEqualIgnoreCase("NPTxyz") || isEqualIgnoreCase("NPAT")  || isEqualIgnoreCase("LANGEVINDYNAMICS") || isEqualIgnoreCase("LD") || isEqualIgnoreCase("NPRT") || isEqualIgnoreCase("NPGT")); | 
| 142 |  | CheckParameter(Dt, isPositive()); | 
| 143 |  | CheckParameter(RunTime, isPositive()); | 
| 144 | < | CheckParameter(InitialConfig, isNotEmpty()); | 
| 144 | > | //CheckParameter(InitialConfig, isNotEmpty()); | 
| 145 |  | CheckParameter(FinalConfig, isNotEmpty()); | 
| 146 |  | CheckParameter(SampleTime, isNonNegative()); | 
| 147 |  | CheckParameter(ResetTime, isNonNegative()); | 
| 150 |  | CheckParameter(SwitchingRadius, isNonNegative()); | 
| 151 |  | CheckParameter(Dielectric, isPositive()); | 
| 152 |  | CheckParameter(ThermalTime,  isNonNegative()); | 
| 153 | < | CheckParameter(TargetPressure,  isPositive()); | 
| 153 | > | //  CheckParameter(TargetPressure,  isPositive()); | 
| 154 |  | CheckParameter(TauThermostat, isPositive()); | 
| 155 |  | CheckParameter(TauBarostat, isPositive()); | 
| 156 |  | CheckParameter(ZconsTime, isPositive()); | 
| 166 |  | CheckParameter(MinimizerLSMaxIter, isPositive()); | 
| 167 |  | CheckParameter(ZconsGap, isPositive()); | 
| 168 |  | CheckParameter(ZconsFixtime, isPositive()); | 
| 169 | < | CheckParameter(ThermodynamicIntegrationLambda, isPositive()); | 
| 169 | > | CheckParameter(ThermodynamicIntegrationLambda, isNonNegative()); | 
| 170 |  | CheckParameter(ThermodynamicIntegrationK, isPositive()); | 
| 171 |  | CheckParameter(ForceFieldVariant, isNotEmpty()); | 
| 172 |  | CheckParameter(ForceFieldFileName, isNotEmpty()); | 
| 173 |  | CheckParameter(ThermIntDistSpringConst, isPositive()); | 
| 174 |  | CheckParameter(ThermIntThetaSpringConst, isPositive()); | 
| 175 |  | CheckParameter(ThermIntOmegaSpringConst, isPositive()); | 
| 176 | < | CheckParameter(SurfaceTension, isPositive()); | 
| 176 | > | //  CheckParameter(SurfaceTension, isNonNegative()); | 
| 177 |  | CheckParameter(ElectrostaticSummationMethod, isEqualIgnoreCase("NONE") || isEqualIgnoreCase("SHIFTED_POTENTIAL") || isEqualIgnoreCase("SHIFTED_FORCE") || isEqualIgnoreCase("REACTION_FIELD")); | 
| 178 |  | CheckParameter(ElectrostaticScreeningMethod, isEqualIgnoreCase("UNDAMPED") || isEqualIgnoreCase("DAMPED")); | 
| 179 |  | CheckParameter(CutoffPolicy, isEqualIgnoreCase("MIX") || isEqualIgnoreCase("MAX") || isEqualIgnoreCase("TRADITIONAL")); | 
| 185 |  | CheckParameter(ThermIntOmegaSpringConst, isPositive()); | 
| 186 |  | CheckParameter(DampingAlpha,isNonNegative()); | 
| 187 |  | CheckParameter(SkinThickness, isPositive()); | 
| 188 | < |  | 
| 188 | > | CheckParameter(Viscosity, isNonNegative()); | 
| 189 | > | CheckParameter(BeadSize, isPositive()); | 
| 190 | > | CheckParameter(FrozenBufferRadius, isPositive()); | 
| 191 | > | CheckParameter(LangevinBufferRadius, isPositive()); | 
| 192 |  | for(std::vector<Component*>::iterator i = components_.begin(); i != components_.end(); ++i) { | 
| 193 |  | if (!(*i)->findMoleculeStamp(moleculeStamps_)) { | 
| 194 |  | std::ostringstream oss; |