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 |
|
{ |