| 44 | 
  | 
  char* conversionCheck; | 
| 45 | 
  | 
  int conversionError; | 
| 46 | 
  | 
  int optionError; | 
| 47 | 
+ | 
  char* rmsdName; | 
| 48 | 
  | 
  char* pair1; | 
| 49 | 
  | 
  char* pair2; | 
| 50 | 
  | 
  int scdCorr; | 
| 51 | 
  | 
  double startTime; | 
| 52 | 
+ | 
  double avgXY; | 
| 53 | 
  | 
  double maxLength; | 
| 54 | 
< | 
  int directorHead, directorWhole, doRMSD, doGofz; | 
| 54 | 
> | 
 | 
| 55 | 
> | 
  int directorHead, directorWhole, doRMSD, doGofz, doAvgXY; | 
| 56 | 
  | 
  enum atomNames rmsdType; | 
| 57 | 
  | 
 | 
| 58 | 
  | 
  // system initialization | 
| 75 | 
  | 
  directorWhole = 0; | 
| 76 | 
  | 
  doRMSD = 0; | 
| 77 | 
  | 
  doGofz = 0; | 
| 78 | 
+ | 
  doAvgXY = 0; | 
| 79 | 
  | 
   | 
| 80 | 
  | 
 | 
| 81 | 
  | 
  // parse the command line | 
| 103 | 
  | 
            usage(); | 
| 104 | 
  | 
            exit(0); | 
| 105 | 
  | 
          }        | 
| 106 | 
< | 
          pair1 = argV[i]; | 
| 106 | 
> | 
          rmsdName = argV[i]; | 
| 107 | 
  | 
 | 
| 108 | 
< | 
          if( !strcmp(pair1, "HEAD") )    rmsdType = HEAD; | 
| 109 | 
< | 
          else if( !strcmp(pair1, "CH") ) rmsdType = CH; | 
| 110 | 
< | 
          else if( !strcmp(pair1, "CH2") )rmsdType = CH2; | 
| 111 | 
< | 
          else if( !strcmp(pair1, "CH3") )rmsdType = CH3; | 
| 112 | 
< | 
          else if( !strcmp(pair1, "SSD") )rmsdType = SSD; | 
| 113 | 
< | 
          else if( !strcmp(pair1, "COM") )rmsdType = COM; | 
| 108 | 
> | 
          if( !strcmp(rmsdName, "HEAD") )    rmsdType = HEAD; | 
| 109 | 
> | 
          else if( !strcmp(rmsdName, "CH") ) rmsdType = CH; | 
| 110 | 
> | 
          else if( !strcmp(rmsdName, "CH2") )rmsdType = CH2; | 
| 111 | 
> | 
          else if( !strcmp(rmsdName, "CH3") )rmsdType = CH3; | 
| 112 | 
> | 
          else if( !strcmp(rmsdName, "SSD") )rmsdType = SSD; | 
| 113 | 
> | 
          else if( !strcmp(rmsdName, "COM") )rmsdType = COM; | 
| 114 | 
  | 
          else{ | 
| 115 | 
  | 
            fprintf( stderr, | 
| 116 | 
  | 
                     "Unrecognized rmsd atom type \"%s\"\n", | 
| 117 | 
< | 
                     pair1 ); | 
| 117 | 
> | 
                     rmsdName ); | 
| 118 | 
  | 
            exit(0); | 
| 119 | 
  | 
          } | 
| 120 | 
  | 
 | 
| 335 | 
  | 
            doGofz = 1; | 
| 336 | 
  | 
            break; | 
| 337 | 
  | 
 | 
| 338 | 
+ | 
          case 'v': | 
| 339 | 
+ | 
            // -v turn on average XY | 
| 340 | 
  | 
 | 
| 341 | 
+ | 
            doAvgXY = 1; | 
| 342 | 
+ | 
            break; | 
| 343 | 
+ | 
 | 
| 344 | 
+ | 
 | 
| 345 | 
  | 
          default: | 
| 346 | 
  | 
 | 
| 347 | 
  | 
            fprintf( stderr,  | 
| 510 | 
  | 
              "Calculating the RMSD\n" ); | 
| 511 | 
  | 
    fflush( stdout ); | 
| 512 | 
  | 
     | 
| 513 | 
< | 
    rmsd( rmsdType, startTime, outPrefix ); | 
| 513 | 
> | 
    rmsd( rmsdType, startTime, outPrefix, rmsdName ); | 
| 514 | 
  | 
  } | 
| 515 | 
  | 
 | 
| 516 | 
  | 
  if(doGofz){ | 
| 521 | 
  | 
     | 
| 522 | 
  | 
    calcGofz( startTime, atoms, outPrefix ); | 
| 523 | 
  | 
  } | 
| 524 | 
+ | 
 | 
| 525 | 
+ | 
  if(doAvgXY){ | 
| 526 | 
+ | 
 | 
| 527 | 
+ | 
    fprintf ( stdout, | 
| 528 | 
+ | 
              "Calculating the average XY\n" ); | 
| 529 | 
+ | 
    fflush( stdout ); | 
| 530 | 
+ | 
     | 
| 531 | 
+ | 
    avgXY = calcAverageXY( startTime ); | 
| 532 | 
+ | 
     | 
| 533 | 
+ | 
    fprintf ( stdout, | 
| 534 | 
+ | 
              "Done, average XY = %6G\n", avgXY ); | 
| 535 | 
+ | 
    fflush( stdout ); | 
| 536 | 
+ | 
     | 
| 537 | 
+ | 
  } | 
| 538 | 
  | 
   | 
| 539 | 
  | 
 | 
| 540 | 
  | 
  closeFile(); |