| 62 | 
  | 
  args_info->theta_omega_given = 0 ; | 
| 63 | 
  | 
  args_info->gxyz_given = 0 ; | 
| 64 | 
  | 
  args_info->p2_given = 0 ; | 
| 65 | 
+ | 
  args_info->rp2_given = 0 ; | 
| 66 | 
  | 
  args_info->scd_given = 0 ; | 
| 67 | 
  | 
  args_info->density_given = 0 ; | 
| 68 | 
  | 
  args_info->slab_density_given = 0 ; | 
| 121 | 
  | 
  printf("%s\n","  -i, --input=filename          input dump file"); | 
| 122 | 
  | 
  printf("%s\n","  -o, --output=filename         output file name"); | 
| 123 | 
  | 
  printf("%s\n","  -n, --step=INT                process every n frame  (default=`1')"); | 
| 124 | 
< | 
  printf("%s\n","  -x,  --nbins_x=INT              number of bins in x axis  (default=`100')"); | 
| 125 | 
< | 
  printf("%s\n","  -y , --nbins_y=INT              number of bins in y axis  (default=`100')"); | 
| 126 | 
< | 
   printf("%s\n","  -r, --nrbins=INT              number of bins for distance  (default=`100')"); | 
| 124 | 
> | 
  printf("%s\n","  -r, --nrbins=INT              number of bins for distance  (default=`100')"); | 
| 125 | 
> | 
  printf("%s\n","  -x, --nbins_x=INT             number of bins in x axis  (default=`100')"); | 
| 126 | 
> | 
  printf("%s\n","  -y, --nbins_y=INT             number of bins in y axis  (default=`100')"); | 
| 127 | 
  | 
  printf("%s\n","  -a, --nanglebins=INT          number of bins for cos(angle)  (default=`50')"); | 
| 128 | 
  | 
  printf("%s\n","  -l, --length=DOUBLE           maximum length (Defaults to 1/2 smallest length \n                                  of first frame)"); | 
| 129 | 
  | 
  printf("%s\n","  -z, --zoffset=DOUBLE          Where to set the zero for the slab_density \n                                  calculation  (default=`0')"); | 
| 141 | 
  | 
  printf("%s\n","      --theta_omega             g(cos(theta), cos(omega))"); | 
| 142 | 
  | 
  printf("%s\n","      --gxyz                    g(x, y, z)"); | 
| 143 | 
  | 
  printf("%s\n","      --p2                      p2 order parameter (--sele1 and --sele2 must be \n                                  specified)"); | 
| 144 | 
+ | 
  printf("%s\n","      --rp2                     rp2 order parameter (--sele1 and --sele2 must \n                                  be specified)"); | 
| 145 | 
  | 
  printf("%s\n","      --scd                     scd order parameter(either --sele1, --sele2, \n                                  --sele3 are specified or --molname, --begin, \n                                  --end are specified)"); | 
| 146 | 
  | 
  printf("%s\n","      --density                 density plot (--sele1 must be specified)"); | 
| 147 | 
  | 
  printf("%s\n","      --slab_density            slab density (--sele1 must be specified)"); | 
| 148 | 
< | 
  printf("%s\n","      --hxy            hxy (--sele1 must be specified)"); | 
| 148 | 
> | 
  printf("%s\n","      --hxy                     hxy (--sele1 must be specified)"); | 
| 149 | 
  | 
   | 
| 150 | 
  | 
} | 
| 151 | 
  | 
 | 
| 190 | 
  | 
      free (args_info->nrbins_orig); /* free previous argument */ | 
| 191 | 
  | 
      args_info->nrbins_orig = 0; | 
| 192 | 
  | 
    } | 
| 193 | 
< | 
    if (args_info->nbins_x_orig) | 
| 193 | 
> | 
  if (args_info->nbins_x_orig) | 
| 194 | 
  | 
    { | 
| 195 | 
  | 
      free (args_info->nbins_x_orig); /* free previous argument */ | 
| 196 | 
  | 
      args_info->nbins_x_orig = 0; | 
| 200 | 
  | 
      free (args_info->nbins_y_orig); /* free previous argument */ | 
| 201 | 
  | 
      args_info->nbins_y_orig = 0; | 
| 202 | 
  | 
    } | 
| 201 | 
– | 
         | 
| 203 | 
  | 
  if (args_info->nanglebins_orig) | 
| 204 | 
  | 
    { | 
| 205 | 
  | 
      free (args_info->nanglebins_orig); /* free previous argument */ | 
| 428 | 
  | 
  } | 
| 429 | 
  | 
  if (args_info->p2_given) { | 
| 430 | 
  | 
    fprintf(outfile, "%s\n", "p2"); | 
| 431 | 
+ | 
  } | 
| 432 | 
+ | 
  if (args_info->rp2_given) { | 
| 433 | 
+ | 
    fprintf(outfile, "%s\n", "rp2"); | 
| 434 | 
  | 
  } | 
| 435 | 
  | 
  if (args_info->scd_given) { | 
| 436 | 
  | 
    fprintf(outfile, "%s\n", "scd"); | 
| 444 | 
  | 
  if (args_info->hxy_given) { | 
| 445 | 
  | 
    fprintf(outfile, "%s\n", "hxy"); | 
| 446 | 
  | 
  } | 
| 447 | 
+ | 
   | 
| 448 | 
  | 
  fclose (outfile); | 
| 449 | 
  | 
 | 
| 450 | 
  | 
  i = EXIT_SUCCESS; | 
| 489 | 
  | 
  args_info->theta_omega_given = 0 ; | 
| 490 | 
  | 
  args_info->gxyz_given = 0 ; | 
| 491 | 
  | 
  args_info->p2_given = 0 ; | 
| 492 | 
+ | 
  args_info->rp2_given = 0 ; | 
| 493 | 
  | 
  args_info->scd_given = 0 ; | 
| 494 | 
  | 
  args_info->density_given = 0 ; | 
| 495 | 
  | 
  args_info->slab_density_given = 0 ; | 
| 496 | 
< | 
  args_info->hxy_given= 0 ; | 
| 496 | 
> | 
  args_info->hxy_given = 0 ; | 
| 497 | 
  | 
 | 
| 498 | 
  | 
  args_info->staticProps_group_counter = 0; | 
| 499 | 
  | 
} | 
| 587 | 
  | 
        { "output",     1, NULL, 'o' }, | 
| 588 | 
  | 
        { "step",       1, NULL, 'n' }, | 
| 589 | 
  | 
        { "nrbins",     1, NULL, 'r' }, | 
| 590 | 
+ | 
        { "nbins_x",    1, NULL, 'x' }, | 
| 591 | 
+ | 
        { "nbins_y",    1, NULL, 'y' }, | 
| 592 | 
  | 
        { "nanglebins", 1, NULL, 'a' }, | 
| 593 | 
  | 
        { "length",     1, NULL, 'l' }, | 
| 594 | 
  | 
        { "zoffset",    1, NULL, 'z' }, | 
| 605 | 
  | 
        { "theta_omega",        0, NULL, 0 }, | 
| 606 | 
  | 
        { "gxyz",       0, NULL, 0 }, | 
| 607 | 
  | 
        { "p2", 0, NULL, 0 }, | 
| 608 | 
+ | 
        { "rp2",        0, NULL, 0 }, | 
| 609 | 
  | 
        { "scd",        0, NULL, 0 }, | 
| 610 | 
  | 
        { "density",    0, NULL, 0 }, | 
| 611 | 
  | 
        { "slab_density",       0, NULL, 0 }, | 
| 614 | 
  | 
      }; | 
| 615 | 
  | 
 | 
| 616 | 
  | 
      stop_char = 0; | 
| 617 | 
< | 
      c = getopt_long (argc, argv, "hVi:o:n:r:a:l:z:", long_options, &option_index); | 
| 617 | 
> | 
      c = getopt_long (argc, argv, "hVi:o:n:r:x:y:a:l:z:", long_options, &option_index); | 
| 618 | 
  | 
 | 
| 619 | 
  | 
      if (c == -1) break;       /* Exit from `while (1)' loop.  */ | 
| 620 | 
  | 
 | 
| 709 | 
  | 
        case 'x':       /* number of bins in x axis.  */ | 
| 710 | 
  | 
          if (local_args_info.nbins_x_given) | 
| 711 | 
  | 
            { | 
| 712 | 
< | 
              fprintf (stderr, "%s: `--nbins_x' (`-r') option given more than once%s\n", argv[0], (additional_error ? additional_error : "")); | 
| 712 | 
> | 
              fprintf (stderr, "%s: `--nbins_x' (`-x') option given more than once%s\n", argv[0], (additional_error ? additional_error : "")); | 
| 713 | 
  | 
              goto failure; | 
| 714 | 
  | 
            } | 
| 715 | 
  | 
          if (args_info->nbins_x_given && ! override) | 
| 729 | 
  | 
        case 'y':       /* number of bins in y axis.  */ | 
| 730 | 
  | 
          if (local_args_info.nbins_y_given) | 
| 731 | 
  | 
            { | 
| 732 | 
< | 
              fprintf (stderr, "%s: `--nbins_y' (`-r') option given more than once%s\n", argv[0], (additional_error ? additional_error : "")); | 
| 732 | 
> | 
              fprintf (stderr, "%s: `--nbins_y' (`-y') option given more than once%s\n", argv[0], (additional_error ? additional_error : "")); | 
| 733 | 
  | 
              goto failure; | 
| 734 | 
  | 
            } | 
| 735 | 
  | 
          if (args_info->nbins_y_given && ! override) | 
| 1047 | 
  | 
            args_info->staticProps_group_counter += 1; | 
| 1048 | 
  | 
            break; | 
| 1049 | 
  | 
          } | 
| 1050 | 
+ | 
          /* rp2 order parameter (--sele1 and --sele2 must be specified).  */ | 
| 1051 | 
+ | 
          else if (strcmp (long_options[option_index].name, "rp2") == 0) | 
| 1052 | 
+ | 
          { | 
| 1053 | 
+ | 
            if (local_args_info.rp2_given) | 
| 1054 | 
+ | 
              { | 
| 1055 | 
+ | 
                fprintf (stderr, "%s: `--rp2' option given more than once%s\n", argv[0], (additional_error ? additional_error : "")); | 
| 1056 | 
+ | 
                goto failure; | 
| 1057 | 
+ | 
              } | 
| 1058 | 
+ | 
            if (args_info->rp2_given && ! override) | 
| 1059 | 
+ | 
              continue; | 
| 1060 | 
+ | 
            local_args_info.rp2_given = 1; | 
| 1061 | 
+ | 
            args_info->rp2_given = 1; | 
| 1062 | 
+ | 
            if (args_info->staticProps_group_counter && override) | 
| 1063 | 
+ | 
              reset_group_staticProps (args_info); | 
| 1064 | 
+ | 
            args_info->staticProps_group_counter += 1; | 
| 1065 | 
+ | 
            break; | 
| 1066 | 
+ | 
          } | 
| 1067 | 
  | 
          /* scd order parameter(either --sele1, --sele2, --sele3 are specified or --molname, --begin, --end are specified).  */ | 
| 1068 | 
  | 
          else if (strcmp (long_options[option_index].name, "scd") == 0) | 
| 1069 | 
  | 
          { |