66 |
|
newMdFile.open(newMdFileName.c_str()); |
67 |
|
oldMdFile.getline(buffer, MAXLEN); |
68 |
|
|
69 |
– |
unsigned int i = 0; |
69 |
|
while (!oldMdFile.eof()) { |
70 |
|
|
71 |
|
//correct molecule number |
166 |
|
CurlingStoneDecahedron* csd = new CurlingStoneDecahedron(columnAtoms, nShells, twinAtoms, truncatedPlanes); |
167 |
|
Points = csd->getPoints(); |
168 |
|
} |
170 |
– |
|
171 |
– |
|
169 |
|
|
170 |
|
outputFileName = args_info.output_arg; |
171 |
|
|
172 |
< |
//creat new .md file on fly which corrects the number of molecule |
172 |
> |
// create a new .md file on fly which corrects the number of |
173 |
> |
// molecules |
174 |
|
|
175 |
|
createMdFile(inputFileName, outputFileName, Points.size()); |
176 |
|
|
192 |
|
Vector3d boxMax; |
193 |
|
Vector3d boxMin; |
194 |
|
|
195 |
< |
for (int n = 0; n < Points.size(); n++) { |
195 |
> |
for (unsigned int n = 0; n < Points.size(); n++) { |
196 |
|
mol = NewInfo->getMoleculeByGlobalIndex(l); |
197 |
|
Vector3d location = Points[n] * latticeConstant; |
198 |
|
Vector3d orientation = Vector3d(0, 0, 1.0); |
199 |
< |
|
199 |
> |
|
200 |
|
if (n == 0) { |
201 |
|
boxMax = location; |
202 |
|
boxMin = location; |
234 |
|
|
235 |
|
delete writer; |
236 |
|
|
237 |
< |
// cleanup a by calling sim error..... |
237 |
> |
// clean up by calling simError..... |
238 |
|
sprintf(painCave.errMsg, "A new OpenMD file called \"%s\" has been " |
239 |
|
"generated.\n", outputFileName.c_str()); |
240 |
|
painCave.isFatal = 0; |