| 12 |
|
#include "directorHead.h" |
| 13 |
|
#include "directorWhole.h" |
| 14 |
|
#include "rmsd.h" |
| 15 |
+ |
#include "gofz.h" |
| 16 |
|
|
| 17 |
|
|
| 18 |
|
#define VERSION_MAJOR 0 |
| 49 |
|
int scdCorr; |
| 50 |
|
double startTime; |
| 51 |
|
double maxLength; |
| 52 |
< |
int directorHead, directorWhole, doRMSD; |
| 52 |
> |
int directorHead, directorWhole, doRMSD, doGofz; |
| 53 |
|
enum atomNames rmsdType; |
| 54 |
|
|
| 55 |
|
// system initialization |
| 71 |
|
directorHead = 0; |
| 72 |
|
directorWhole = 0; |
| 73 |
|
doRMSD = 0; |
| 74 |
+ |
doGofz = 0; |
| 75 |
|
|
| 76 |
|
|
| 77 |
|
// parse the command line |
| 320 |
|
break; |
| 321 |
|
|
| 322 |
|
case 'w': |
| 323 |
< |
// -h turn on director head |
| 323 |
> |
// -w turn on director head |
| 324 |
|
|
| 325 |
|
directorWhole = 1; |
| 326 |
+ |
break; |
| 327 |
+ |
|
| 328 |
+ |
case 'g': |
| 329 |
+ |
// -g turn on gofZ |
| 330 |
+ |
|
| 331 |
+ |
doGofz = 1; |
| 332 |
|
break; |
| 333 |
|
|
| 334 |
|
|
| 502 |
|
|
| 503 |
|
rmsd( rmsdType, startTime, outPrefix ); |
| 504 |
|
} |
| 505 |
+ |
|
| 506 |
+ |
if(doGofz){ |
| 507 |
+ |
|
| 508 |
+ |
fprintf ( stdout, |
| 509 |
+ |
"Calculating the gofZ\n" ); |
| 510 |
+ |
fflush( stdout ); |
| 511 |
+ |
|
| 512 |
+ |
calcGofz( startTime, atoms, outPrefix ); |
| 513 |
+ |
} |
| 514 |
|
|
| 515 |
|
|
| 516 |
|
closeFile(); |
| 537 |
|
" *Defaults to 1/2 smallest length of first frame.\n" |
| 538 |
|
" -s Calculate the Scd chain correlation.\n" |
| 539 |
|
" -h Calculate the directors for the head groups\n" |
| 540 |
< |
" -w Calculate the director fro the bilayers\n" |
| 540 |
> |
" -w Calculate the director from the bilayers\n" |
| 541 |
> |
" -g Calculate the gofz profile\n" |
| 542 |
|
"\n" |
| 543 |
|
" long:\n" |
| 544 |
|
" -----\n" |