51 |
|
double diagonal; /* the diagonal length of the sim box */ |
52 |
|
|
53 |
|
unsigned int n_atoms; /*the number of atoms in each time step */ |
54 |
< |
char read_buffer[120]; /*the line buffer for reading */ |
54 |
> |
char read_buffer[2000]; /*the line buffer for reading */ |
55 |
|
char *eof_test; /*ptr to see when we reach the end of the file */ |
56 |
|
char *foo; /*the pointer to the current string token */ |
57 |
|
FILE *in_file; /* the input file */ |
271 |
|
out_prefix = strtok(in_name, "."); |
272 |
|
} |
273 |
|
|
274 |
< |
sprintf( out_format, "%s%s%%0%dd.pov", pov_dir, out_prefix, nZeroes ); |
274 |
> |
sprintf( out_format, "%s%s-%%0%dd.pov", pov_dir, out_prefix, nZeroes ); |
275 |
|
|
276 |
|
// start reading the first frame |
277 |
|
|
303 |
|
if( draw_box ){ |
304 |
|
foo = strtok(read_buffer, " ,;\t"); |
305 |
|
if(foo == NULL){ |
306 |
< |
printf("error in reading file\n"); |
306 |
> |
printf("error in reading file time\n"); |
307 |
|
exit(8); |
308 |
|
} |
309 |
|
|
310 |
|
foo = strtok(NULL, " ,;\t"); |
311 |
|
if(foo == NULL){ |
312 |
< |
printf("error in reading file\n"); |
312 |
> |
printf("error in reading file h00\n"); |
313 |
|
exit(8); |
314 |
|
} |
315 |
|
current_frame->Hmat[0][0] = atof( foo ); |
316 |
|
|
317 |
|
foo = strtok(NULL, " ,;\t"); |
318 |
|
if(foo == NULL){ |
319 |
< |
printf("error in reading file\n"); |
319 |
> |
printf("error in reading file h10\n"); |
320 |
|
exit(8); |
321 |
|
} |
322 |
|
current_frame->Hmat[1][0] = atof( foo ); |
323 |
|
|
324 |
|
foo = strtok(NULL, " ,;\t"); |
325 |
|
if(foo == NULL){ |
326 |
< |
printf("error in reading file\n"); |
326 |
> |
printf("error in reading file h20\n"); |
327 |
|
exit(8); |
328 |
|
} |
329 |
|
current_frame->Hmat[2][0] = atof( foo ); |
330 |
|
|
331 |
|
foo = strtok(NULL, " ,;\t"); |
332 |
|
if(foo == NULL){ |
333 |
< |
printf("error in reading file\n"); |
333 |
> |
printf("error in reading file h01\n"); |
334 |
|
exit(8); |
335 |
|
} |
336 |
|
current_frame->Hmat[0][1] = atof( foo ); |
337 |
|
|
338 |
|
foo = strtok(NULL, " ,;\t"); |
339 |
|
if(foo == NULL){ |
340 |
< |
printf("error in reading file\n"); |
340 |
> |
printf("error in reading file h11\n"); |
341 |
|
exit(8); |
342 |
|
} |
343 |
|
current_frame->Hmat[1][1] = atof( foo ); |
344 |
|
|
345 |
|
foo = strtok(NULL, " ,;\t"); |
346 |
|
if(foo == NULL){ |
347 |
< |
printf("error in reading file\n"); |
347 |
> |
printf("error in reading file h21\n"); |
348 |
|
exit(8); |
349 |
|
} |
350 |
|
current_frame->Hmat[2][1] = atof( foo ); |
351 |
|
|
352 |
|
foo = strtok(NULL, " ,;\t"); |
353 |
|
if(foo == NULL){ |
354 |
< |
printf("error in reading file\n"); |
354 |
> |
printf("error in reading file h02\n"); |
355 |
|
exit(8); |
356 |
|
} |
357 |
|
current_frame->Hmat[0][2] = atof( foo ); |
358 |
|
|
359 |
|
foo = strtok(NULL, " ,;\t"); |
360 |
|
if(foo == NULL){ |
361 |
< |
printf("error in reading file\n"); |
361 |
> |
printf("error in reading file h12\n"); |
362 |
|
exit(8); |
363 |
|
} |
364 |
|
current_frame->Hmat[1][2] = atof( foo ); |
365 |
|
|
366 |
|
foo = strtok(NULL, " ,;\t"); |
367 |
|
if(foo == NULL){ |
368 |
< |
printf("error in reading file\n"); |
368 |
> |
printf("error in reading file h22\n"); |
369 |
|
exit(8); |
370 |
|
} |
371 |
|
current_frame->Hmat[2][2] = atof( foo ); |
376 |
|
|
377 |
|
eof_test = fgets(read_buffer, sizeof(read_buffer), in_file); |
378 |
|
if(eof_test == NULL){ |
379 |
< |
printf("error in reading file\n"); |
379 |
> |
printf("error in reading file line at atom %d\n", i); |
380 |
|
exit(8); |
381 |
|
} |
382 |
|
|
383 |
|
foo = strtok(read_buffer, " ,;\t"); |
384 |
|
if(foo == NULL){ |
385 |
< |
printf("error in reading file\n"); |
385 |
> |
printf("error in reading atom %d name\n", i); |
386 |
|
exit(8); |
387 |
|
} |
388 |
|
(void)strcpy(current_frame->r[i].name, foo); /*copy the atom name */ |
391 |
|
|
392 |
|
foo = strtok(NULL, " ,;\t"); |
393 |
|
if(foo == NULL){ |
394 |
< |
printf("error in reading file\n"); |
394 |
> |
printf("error in reading atom %d position x\n", i); |
395 |
|
exit(8); |
396 |
|
} |
397 |
|
(void)sscanf(foo, "%lf",¤t_frame->r[i].x); |
401 |
|
|
402 |
|
foo = strtok(NULL, " ,;\t"); |
403 |
|
if(foo == NULL){ |
404 |
< |
printf("error in reading file\n"); |
404 |
> |
printf("error in reading atom %d position y\n", i); |
405 |
|
exit(8); |
406 |
|
} |
407 |
|
(void)sscanf(foo, "%lf", ¤t_frame->r[i].y); |
411 |
|
|
412 |
|
foo = strtok(NULL, " ,;\t"); |
413 |
|
if(foo == NULL){ |
414 |
< |
printf("error in reading file\n"); |
414 |
> |
printf("error in reading atom %d position z\n", i); |
415 |
|
exit(8); |
416 |
|
} |
417 |
|
(void)sscanf(foo, "%lf", ¤t_frame->r[i].z); |