| 43 | 
  | 
  args_info->input_given = 0 ; | 
| 44 | 
  | 
  args_info->output_given = 0 ; | 
| 45 | 
  | 
  args_info->step_given = 0 ; | 
| 46 | 
< | 
  args_info->nrbins_given = 0 ; | 
| 46 | 
> | 
  args_info->nbins_given = 0 ; | 
| 47 | 
  | 
  args_info->nbins_x_given = 0 ; | 
| 48 | 
  | 
  args_info->nbins_y_given = 0 ; | 
| 49 | 
  | 
  args_info->nanglebins_given = 0 ; | 
| 50 | 
  | 
  args_info->length_given = 0 ; | 
| 51 | 
– | 
  args_info->LegendreL_given = 0 ; | 
| 51 | 
  | 
  args_info->rcut_given = 0 ; | 
| 52 | 
  | 
  args_info->zoffset_given = 0 ; | 
| 53 | 
  | 
  args_info->sele1_given = 0 ; | 
| 81 | 
  | 
  args_info->output_orig = NULL; | 
| 82 | 
  | 
  args_info->step_arg = 1; | 
| 83 | 
  | 
  args_info->step_orig = NULL; | 
| 84 | 
< | 
  args_info->nrbins_arg = 100; | 
| 85 | 
< | 
  args_info->nrbins_orig = NULL; | 
| 84 | 
> | 
  args_info->nbins_arg = 100; | 
| 85 | 
> | 
  args_info->nbins_orig = NULL; | 
| 86 | 
  | 
  args_info->nbins_x_arg = 100; | 
| 87 | 
  | 
  args_info->nbins_x_orig = NULL; | 
| 88 | 
  | 
  args_info->nbins_y_arg = 100; | 
| 90 | 
  | 
  args_info->nanglebins_arg = 50; | 
| 91 | 
  | 
  args_info->nanglebins_orig = NULL; | 
| 92 | 
  | 
  args_info->length_orig = NULL; | 
| 94 | 
– | 
  args_info->LegendreL_orig = NULL; | 
| 93 | 
  | 
  args_info->rcut_orig = NULL; | 
| 94 | 
  | 
  args_info->zoffset_arg = 0; | 
| 95 | 
  | 
  args_info->zoffset_orig = NULL; | 
| 124 | 
  | 
  printf("%s\n","  -i, --input=filename          input dump file"); | 
| 125 | 
  | 
  printf("%s\n","  -o, --output=filename         output file name"); | 
| 126 | 
  | 
  printf("%s\n","  -n, --step=INT                process every n frame  (default=`1')"); | 
| 127 | 
< | 
  printf("%s\n","  -r, --nrbins=INT              number of bins for distance  (default=`100')"); | 
| 127 | 
> | 
  printf("%s\n","  -b, --nbins=INT               number of bins (general purpose)  \n                                  (default=`100')"); | 
| 128 | 
  | 
  printf("%s\n","  -x, --nbins_x=INT             number of bins in x axis  (default=`100')"); | 
| 129 | 
  | 
  printf("%s\n","  -y, --nbins_y=INT             number of bins in y axis  (default=`100')"); | 
| 130 | 
  | 
  printf("%s\n","  -a, --nanglebins=INT          number of bins for cos(angle)  (default=`50')"); | 
| 131 | 
  | 
  printf("%s\n","      --length=DOUBLE           maximum length (Defaults to 1/2 smallest length \n                                  of first frame)"); | 
| 134 | 
– | 
  printf("%s\n","  -l, --LegendreL=INT           Order of Legendre Polynomial (used for Bond \n                                  Order calculations)"); | 
| 132 | 
  | 
  printf("%s\n","  -c, --rcut=DOUBLE             cutoff radius (rcut)"); | 
| 133 | 
  | 
  printf("%s\n","  -z, --zoffset=DOUBLE          Where to set the zero for the slab_density \n                                  calculation  (default=`0')"); | 
| 134 | 
  | 
  printf("%s\n","      --sele1=selection script  select the first stuntdouble set"); | 
| 139 | 
  | 
  printf("%s\n","      --begin=INT               begin internal index"); | 
| 140 | 
  | 
  printf("%s\n","      --end=INT                 end internal index"); | 
| 141 | 
  | 
  printf("%s\n","\n Group: staticProps\n   an option of this group is required"); | 
| 142 | 
< | 
  printf("%s\n","  -b, --bo                      bond order parameter (--rcut and --LegendreL \n                                  must be specified"); | 
| 142 | 
> | 
  printf("%s\n","      --bo                      bond order parameter (--rcut must be specified"); | 
| 143 | 
  | 
  printf("%s\n","  -g, --gofr                    g(r)"); | 
| 144 | 
  | 
  printf("%s\n","      --r_theta                 g(r, cos(theta))"); | 
| 145 | 
  | 
  printf("%s\n","      --r_omega                 g(r, cos(omega))"); | 
| 190 | 
  | 
      free (args_info->step_orig); /* free previous argument */ | 
| 191 | 
  | 
      args_info->step_orig = 0; | 
| 192 | 
  | 
    } | 
| 193 | 
< | 
  if (args_info->nrbins_orig) | 
| 193 | 
> | 
  if (args_info->nbins_orig) | 
| 194 | 
  | 
    { | 
| 195 | 
< | 
      free (args_info->nrbins_orig); /* free previous argument */ | 
| 196 | 
< | 
      args_info->nrbins_orig = 0; | 
| 195 | 
> | 
      free (args_info->nbins_orig); /* free previous argument */ | 
| 196 | 
> | 
      args_info->nbins_orig = 0; | 
| 197 | 
  | 
    } | 
| 198 | 
  | 
  if (args_info->nbins_x_orig) | 
| 199 | 
  | 
    { | 
| 215 | 
  | 
      free (args_info->length_orig); /* free previous argument */ | 
| 216 | 
  | 
      args_info->length_orig = 0; | 
| 217 | 
  | 
    } | 
| 221 | 
– | 
  if (args_info->LegendreL_orig) | 
| 222 | 
– | 
    { | 
| 223 | 
– | 
      free (args_info->LegendreL_orig); /* free previous argument */ | 
| 224 | 
– | 
      args_info->LegendreL_orig = 0; | 
| 225 | 
– | 
    } | 
| 218 | 
  | 
  if (args_info->rcut_orig) | 
| 219 | 
  | 
    { | 
| 220 | 
  | 
      free (args_info->rcut_orig); /* free previous argument */ | 
| 330 | 
  | 
      fprintf(outfile, "%s\n", "step"); | 
| 331 | 
  | 
    } | 
| 332 | 
  | 
  } | 
| 333 | 
< | 
  if (args_info->nrbins_given) { | 
| 334 | 
< | 
    if (args_info->nrbins_orig) { | 
| 335 | 
< | 
      fprintf(outfile, "%s=\"%s\"\n", "nrbins", args_info->nrbins_orig); | 
| 333 | 
> | 
  if (args_info->nbins_given) { | 
| 334 | 
> | 
    if (args_info->nbins_orig) { | 
| 335 | 
> | 
      fprintf(outfile, "%s=\"%s\"\n", "nbins", args_info->nbins_orig); | 
| 336 | 
  | 
    } else { | 
| 337 | 
< | 
      fprintf(outfile, "%s\n", "nrbins"); | 
| 337 | 
> | 
      fprintf(outfile, "%s\n", "nbins"); | 
| 338 | 
  | 
    } | 
| 339 | 
  | 
  } | 
| 340 | 
  | 
  if (args_info->nbins_x_given) { | 
| 365 | 
  | 
      fprintf(outfile, "%s\n", "length"); | 
| 366 | 
  | 
    } | 
| 367 | 
  | 
  } | 
| 376 | 
– | 
  if (args_info->LegendreL_given) { | 
| 377 | 
– | 
    if (args_info->LegendreL_orig) { | 
| 378 | 
– | 
      fprintf(outfile, "%s=\"%s\"\n", "LegendreL", args_info->LegendreL_orig); | 
| 379 | 
– | 
    } else { | 
| 380 | 
– | 
      fprintf(outfile, "%s\n", "LegendreL"); | 
| 381 | 
– | 
    } | 
| 382 | 
– | 
  } | 
| 368 | 
  | 
  if (args_info->rcut_given) { | 
| 369 | 
  | 
    if (args_info->rcut_orig) { | 
| 370 | 
  | 
      fprintf(outfile, "%s=\"%s\"\n", "rcut", args_info->rcut_orig); | 
| 607 | 
  | 
        { "input",      1, NULL, 'i' }, | 
| 608 | 
  | 
        { "output",     1, NULL, 'o' }, | 
| 609 | 
  | 
        { "step",       1, NULL, 'n' }, | 
| 610 | 
< | 
        { "nrbins",     1, NULL, 'r' }, | 
| 610 | 
> | 
        { "nbins",      1, NULL, 'b' }, | 
| 611 | 
  | 
        { "nbins_x",    1, NULL, 'x' }, | 
| 612 | 
  | 
        { "nbins_y",    1, NULL, 'y' }, | 
| 613 | 
  | 
        { "nanglebins", 1, NULL, 'a' }, | 
| 614 | 
  | 
        { "length",     1, NULL, 0 }, | 
| 630 | 
– | 
        { "LegendreL",  1, NULL, 'l' }, | 
| 615 | 
  | 
        { "rcut",       1, NULL, 'c' }, | 
| 616 | 
  | 
        { "zoffset",    1, NULL, 'z' }, | 
| 617 | 
  | 
        { "sele1",      1, NULL, 0 }, | 
| 621 | 
  | 
        { "molname",    1, NULL, 0 }, | 
| 622 | 
  | 
        { "begin",      1, NULL, 0 }, | 
| 623 | 
  | 
        { "end",        1, NULL, 0 }, | 
| 624 | 
< | 
        { "bo", 0, NULL, 'b' }, | 
| 624 | 
> | 
        { "bo", 0, NULL, 0 }, | 
| 625 | 
  | 
        { "gofr",       0, NULL, 'g' }, | 
| 626 | 
  | 
        { "r_theta",    0, NULL, 0 }, | 
| 627 | 
  | 
        { "r_omega",    0, NULL, 0 }, | 
| 637 | 
  | 
      }; | 
| 638 | 
  | 
 | 
| 639 | 
  | 
      stop_char = 0; | 
| 640 | 
< | 
      c = getopt_long (argc, argv, "hVi:o:n:r:x:y:a:l:c:z:bgpsd", long_options, &option_index); | 
| 640 | 
> | 
      c = getopt_long (argc, argv, "hVi:o:n:b:x:y:a:c:z:gpsd", long_options, &option_index); | 
| 641 | 
  | 
 | 
| 642 | 
  | 
      if (c == -1) break;       /* Exit from `while (1)' loop.  */ | 
| 643 | 
  | 
 | 
| 709 | 
  | 
          args_info->step_orig = gengetopt_strdup (optarg); | 
| 710 | 
  | 
          break; | 
| 711 | 
  | 
 | 
| 712 | 
< | 
        case 'r':       /* number of bins for distance.  */ | 
| 713 | 
< | 
          if (local_args_info.nrbins_given) | 
| 712 | 
> | 
        case 'b':       /* number of bins (general purpose).  */ | 
| 713 | 
> | 
          if (local_args_info.nbins_given) | 
| 714 | 
  | 
            { | 
| 715 | 
< | 
              fprintf (stderr, "%s: `--nrbins' (`-r') option given more than once%s\n", argv[0], (additional_error ? additional_error : "")); | 
| 715 | 
> | 
              fprintf (stderr, "%s: `--nbins' (`-b') option given more than once%s\n", argv[0], (additional_error ? additional_error : "")); | 
| 716 | 
  | 
              goto failure; | 
| 717 | 
  | 
            } | 
| 718 | 
< | 
          if (args_info->nrbins_given && ! override) | 
| 718 | 
> | 
          if (args_info->nbins_given && ! override) | 
| 719 | 
  | 
            continue; | 
| 720 | 
< | 
          local_args_info.nrbins_given = 1; | 
| 721 | 
< | 
          args_info->nrbins_given = 1; | 
| 722 | 
< | 
          args_info->nrbins_arg = strtol (optarg, &stop_char, 0); | 
| 720 | 
> | 
          local_args_info.nbins_given = 1; | 
| 721 | 
> | 
          args_info->nbins_given = 1; | 
| 722 | 
> | 
          args_info->nbins_arg = strtol (optarg, &stop_char, 0); | 
| 723 | 
  | 
          if (!(stop_char && *stop_char == '\0')) { | 
| 724 | 
  | 
            fprintf(stderr, "%s: invalid numeric value: %s\n", argv[0], optarg); | 
| 725 | 
  | 
            goto failure; | 
| 726 | 
  | 
          } | 
| 727 | 
< | 
          if (args_info->nrbins_orig) | 
| 728 | 
< | 
            free (args_info->nrbins_orig); /* free previous string */ | 
| 729 | 
< | 
          args_info->nrbins_orig = gengetopt_strdup (optarg); | 
| 727 | 
> | 
          if (args_info->nbins_orig) | 
| 728 | 
> | 
            free (args_info->nbins_orig); /* free previous string */ | 
| 729 | 
> | 
          args_info->nbins_orig = gengetopt_strdup (optarg); | 
| 730 | 
  | 
          break; | 
| 731 | 
  | 
 | 
| 732 | 
  | 
        case 'x':       /* number of bins in x axis.  */ | 
| 789 | 
  | 
          args_info->nanglebins_orig = gengetopt_strdup (optarg); | 
| 790 | 
  | 
          break; | 
| 791 | 
  | 
 | 
| 808 | 
– | 
        case 'l':       /* Order of Legendre Polynomial (used for Bond Order calculations).  */ | 
| 809 | 
– | 
          if (local_args_info.LegendreL_given) | 
| 810 | 
– | 
            { | 
| 811 | 
– | 
              fprintf (stderr, "%s: `--LegendreL' (`-l') option given more than once%s\n", argv[0], (additional_error ? additional_error : "")); | 
| 812 | 
– | 
              goto failure; | 
| 813 | 
– | 
            } | 
| 814 | 
– | 
          if (args_info->LegendreL_given && ! override) | 
| 815 | 
– | 
            continue; | 
| 816 | 
– | 
          local_args_info.LegendreL_given = 1; | 
| 817 | 
– | 
          args_info->LegendreL_given = 1; | 
| 818 | 
– | 
          args_info->LegendreL_arg = strtol (optarg, &stop_char, 0); | 
| 819 | 
– | 
          if (!(stop_char && *stop_char == '\0')) { | 
| 820 | 
– | 
            fprintf(stderr, "%s: invalid numeric value: %s\n", argv[0], optarg); | 
| 821 | 
– | 
            goto failure; | 
| 822 | 
– | 
          } | 
| 823 | 
– | 
          if (args_info->LegendreL_orig) | 
| 824 | 
– | 
            free (args_info->LegendreL_orig); /* free previous string */ | 
| 825 | 
– | 
          args_info->LegendreL_orig = gengetopt_strdup (optarg); | 
| 826 | 
– | 
          break; | 
| 827 | 
– | 
 | 
| 792 | 
  | 
        case 'c':       /* cutoff radius (rcut).  */ | 
| 793 | 
  | 
          if (local_args_info.rcut_given) | 
| 794 | 
  | 
            { | 
| 827 | 
  | 
          if (args_info->zoffset_orig) | 
| 828 | 
  | 
            free (args_info->zoffset_orig); /* free previous string */ | 
| 829 | 
  | 
          args_info->zoffset_orig = gengetopt_strdup (optarg); | 
| 866 | 
– | 
          break; | 
| 867 | 
– | 
 | 
| 868 | 
– | 
        case 'b':       /* bond order parameter (--rcut and --LegendreL must be specified.  */ | 
| 869 | 
– | 
          if (local_args_info.bo_given) | 
| 870 | 
– | 
            { | 
| 871 | 
– | 
              fprintf (stderr, "%s: `--bo' (`-b') option given more than once%s\n", argv[0], (additional_error ? additional_error : "")); | 
| 872 | 
– | 
              goto failure; | 
| 873 | 
– | 
            } | 
| 874 | 
– | 
          if (args_info->bo_given && ! override) | 
| 875 | 
– | 
            continue; | 
| 876 | 
– | 
          local_args_info.bo_given = 1; | 
| 877 | 
– | 
          args_info->bo_given = 1; | 
| 878 | 
– | 
          if (args_info->staticProps_group_counter && override) | 
| 879 | 
– | 
            reset_group_staticProps (args_info); | 
| 880 | 
– | 
          args_info->staticProps_group_counter += 1; | 
| 830 | 
  | 
          break; | 
| 831 | 
  | 
 | 
| 832 | 
  | 
        case 'g':       /* g(r).  */ | 
| 1048 | 
  | 
            if (args_info->end_orig) | 
| 1049 | 
  | 
              free (args_info->end_orig); /* free previous string */ | 
| 1050 | 
  | 
            args_info->end_orig = gengetopt_strdup (optarg); | 
| 1051 | 
+ | 
          } | 
| 1052 | 
+ | 
          /* bond order parameter (--rcut must be specified.  */ | 
| 1053 | 
+ | 
          else if (strcmp (long_options[option_index].name, "bo") == 0) | 
| 1054 | 
+ | 
          { | 
| 1055 | 
+ | 
            if (local_args_info.bo_given) | 
| 1056 | 
+ | 
              { | 
| 1057 | 
+ | 
                fprintf (stderr, "%s: `--bo' option given more than once%s\n", argv[0], (additional_error ? additional_error : "")); | 
| 1058 | 
+ | 
                goto failure; | 
| 1059 | 
+ | 
              } | 
| 1060 | 
+ | 
            if (args_info->bo_given && ! override) | 
| 1061 | 
+ | 
              continue; | 
| 1062 | 
+ | 
            local_args_info.bo_given = 1; | 
| 1063 | 
+ | 
            args_info->bo_given = 1; | 
| 1064 | 
+ | 
            if (args_info->staticProps_group_counter && override) | 
| 1065 | 
+ | 
              reset_group_staticProps (args_info); | 
| 1066 | 
+ | 
            args_info->staticProps_group_counter += 1; | 
| 1067 | 
+ | 
            break; | 
| 1068 | 
  | 
          } | 
| 1069 | 
  | 
          /* g(r, cos(theta)).  */ | 
| 1070 | 
  | 
          else if (strcmp (long_options[option_index].name, "r_theta") == 0) |