| 36 |  | "  -b, --nbins=INT               number of bins (general purpose)  \n                                  (default=`100')", | 
| 37 |  | "  -x, --nbins_x=INT             number of bins in x axis  (default=`100')", | 
| 38 |  | "  -y, --nbins_y=INT             number of bins in y axis  (default=`100')", | 
| 39 | + | "      --nbins_z=INT             number of bins in z axis  (default=`100')", | 
| 40 |  | "  -a, --nanglebins=INT          number of bins for cos(angle)  (default=`50')", | 
| 41 |  | "      --length=DOUBLE           maximum length (Defaults to 1/2 smallest length \n                                  of first frame)", | 
| 42 |  | "  -c, --rcut=DOUBLE             cutoff radius (rcut)", | 
| 58 |  | "      --gofz                    g(z)", | 
| 59 |  | "      --r_theta                 g(r, cos(theta))", | 
| 60 |  | "      --r_omega                 g(r, cos(omega))", | 
| 61 | + | "      --r_z                     g(r, z)", | 
| 62 |  | "      --theta_omega             g(cos(theta), cos(omega))", | 
| 63 |  | "      --gxyz                    g(x, y, z)", | 
| 64 |  | "  -p, --p2                      p2 order parameter (--sele1 and --sele2 must be \n                                  specified)", | 
| 105 |  | args_info->nbins_given = 0 ; | 
| 106 |  | args_info->nbins_x_given = 0 ; | 
| 107 |  | args_info->nbins_y_given = 0 ; | 
| 108 | + | args_info->nbins_z_given = 0 ; | 
| 109 |  | args_info->nanglebins_given = 0 ; | 
| 110 |  | args_info->length_given = 0 ; | 
| 111 |  | args_info->rcut_given = 0 ; | 
| 126 |  | args_info->gofz_given = 0 ; | 
| 127 |  | args_info->r_theta_given = 0 ; | 
| 128 |  | args_info->r_omega_given = 0 ; | 
| 129 | + | args_info->r_z_given = 0 ; | 
| 130 |  | args_info->theta_omega_given = 0 ; | 
| 131 |  | args_info->gxyz_given = 0 ; | 
| 132 |  | args_info->p2_given = 0 ; | 
| 156 |  | args_info->nbins_x_orig = NULL; | 
| 157 |  | args_info->nbins_y_arg = 100; | 
| 158 |  | args_info->nbins_y_orig = NULL; | 
| 159 | + | args_info->nbins_z_arg = 100; | 
| 160 | + | args_info->nbins_z_orig = NULL; | 
| 161 |  | args_info->nanglebins_arg = 50; | 
| 162 |  | args_info->nanglebins_orig = NULL; | 
| 163 |  | args_info->length_orig = NULL; | 
| 195 |  | args_info->nbins_help = gengetopt_args_info_help[5] ; | 
| 196 |  | args_info->nbins_x_help = gengetopt_args_info_help[6] ; | 
| 197 |  | args_info->nbins_y_help = gengetopt_args_info_help[7] ; | 
| 198 | < | args_info->nanglebins_help = gengetopt_args_info_help[8] ; | 
| 199 | < | args_info->length_help = gengetopt_args_info_help[9] ; | 
| 200 | < | args_info->rcut_help = gengetopt_args_info_help[10] ; | 
| 201 | < | args_info->zoffset_help = gengetopt_args_info_help[11] ; | 
| 202 | < | args_info->sele1_help = gengetopt_args_info_help[12] ; | 
| 203 | < | args_info->sele2_help = gengetopt_args_info_help[13] ; | 
| 204 | < | args_info->sele3_help = gengetopt_args_info_help[14] ; | 
| 205 | < | args_info->refsele_help = gengetopt_args_info_help[15] ; | 
| 206 | < | args_info->comsele_help = gengetopt_args_info_help[16] ; | 
| 207 | < | args_info->molname_help = gengetopt_args_info_help[17] ; | 
| 208 | < | args_info->begin_help = gengetopt_args_info_help[18] ; | 
| 209 | < | args_info->end_help = gengetopt_args_info_help[19] ; | 
| 210 | < | args_info->radius_help = gengetopt_args_info_help[20] ; | 
| 211 | < | args_info->bo_help = gengetopt_args_info_help[22] ; | 
| 212 | < | args_info->bor_help = gengetopt_args_info_help[23] ; | 
| 213 | < | args_info->bad_help = gengetopt_args_info_help[24] ; | 
| 214 | < | args_info->gofr_help = gengetopt_args_info_help[25] ; | 
| 215 | < | args_info->gofz_help = gengetopt_args_info_help[26] ; | 
| 216 | < | args_info->r_theta_help = gengetopt_args_info_help[27] ; | 
| 217 | < | args_info->r_omega_help = gengetopt_args_info_help[28] ; | 
| 218 | < | args_info->theta_omega_help = gengetopt_args_info_help[29] ; | 
| 219 | < | args_info->gxyz_help = gengetopt_args_info_help[30] ; | 
| 220 | < | args_info->p2_help = gengetopt_args_info_help[31] ; | 
| 221 | < | args_info->rp2_help = gengetopt_args_info_help[32] ; | 
| 222 | < | args_info->scd_help = gengetopt_args_info_help[33] ; | 
| 223 | < | args_info->density_help = gengetopt_args_info_help[34] ; | 
| 224 | < | args_info->slab_density_help = gengetopt_args_info_help[35] ; | 
| 225 | < | args_info->p_angle_help = gengetopt_args_info_help[36] ; | 
| 226 | < | args_info->hxy_help = gengetopt_args_info_help[37] ; | 
| 227 | < | args_info->rho_r_help = gengetopt_args_info_help[38] ; | 
| 228 | < | args_info->hullvol_help = gengetopt_args_info_help[39] ; | 
| 198 | > | args_info->nbins_z_help = gengetopt_args_info_help[8] ; | 
| 199 | > | args_info->nanglebins_help = gengetopt_args_info_help[9] ; | 
| 200 | > | args_info->length_help = gengetopt_args_info_help[10] ; | 
| 201 | > | args_info->rcut_help = gengetopt_args_info_help[11] ; | 
| 202 | > | args_info->zoffset_help = gengetopt_args_info_help[12] ; | 
| 203 | > | args_info->sele1_help = gengetopt_args_info_help[13] ; | 
| 204 | > | args_info->sele2_help = gengetopt_args_info_help[14] ; | 
| 205 | > | args_info->sele3_help = gengetopt_args_info_help[15] ; | 
| 206 | > | args_info->refsele_help = gengetopt_args_info_help[16] ; | 
| 207 | > | args_info->comsele_help = gengetopt_args_info_help[17] ; | 
| 208 | > | args_info->molname_help = gengetopt_args_info_help[18] ; | 
| 209 | > | args_info->begin_help = gengetopt_args_info_help[19] ; | 
| 210 | > | args_info->end_help = gengetopt_args_info_help[20] ; | 
| 211 | > | args_info->radius_help = gengetopt_args_info_help[21] ; | 
| 212 | > | args_info->bo_help = gengetopt_args_info_help[23] ; | 
| 213 | > | args_info->bor_help = gengetopt_args_info_help[24] ; | 
| 214 | > | args_info->bad_help = gengetopt_args_info_help[25] ; | 
| 215 | > | args_info->gofr_help = gengetopt_args_info_help[26] ; | 
| 216 | > | args_info->gofz_help = gengetopt_args_info_help[27] ; | 
| 217 | > | args_info->r_theta_help = gengetopt_args_info_help[28] ; | 
| 218 | > | args_info->r_omega_help = gengetopt_args_info_help[29] ; | 
| 219 | > | args_info->r_z_help = gengetopt_args_info_help[30] ; | 
| 220 | > | args_info->theta_omega_help = gengetopt_args_info_help[31] ; | 
| 221 | > | args_info->gxyz_help = gengetopt_args_info_help[32] ; | 
| 222 | > | args_info->p2_help = gengetopt_args_info_help[33] ; | 
| 223 | > | args_info->rp2_help = gengetopt_args_info_help[34] ; | 
| 224 | > | args_info->scd_help = gengetopt_args_info_help[35] ; | 
| 225 | > | args_info->density_help = gengetopt_args_info_help[36] ; | 
| 226 | > | args_info->slab_density_help = gengetopt_args_info_help[37] ; | 
| 227 | > | args_info->p_angle_help = gengetopt_args_info_help[38] ; | 
| 228 | > | args_info->hxy_help = gengetopt_args_info_help[39] ; | 
| 229 | > | args_info->rho_r_help = gengetopt_args_info_help[40] ; | 
| 230 | > | args_info->hullvol_help = gengetopt_args_info_help[41] ; | 
| 231 |  |  | 
| 232 |  | } | 
| 233 |  |  | 
| 317 |  | free_string_field (&(args_info->nbins_orig)); | 
| 318 |  | free_string_field (&(args_info->nbins_x_orig)); | 
| 319 |  | free_string_field (&(args_info->nbins_y_orig)); | 
| 320 | + | free_string_field (&(args_info->nbins_z_orig)); | 
| 321 |  | free_string_field (&(args_info->nanglebins_orig)); | 
| 322 |  | free_string_field (&(args_info->length_orig)); | 
| 323 |  | free_string_field (&(args_info->rcut_orig)); | 
| 387 |  | write_into_file(outfile, "nbins_x", args_info->nbins_x_orig, 0); | 
| 388 |  | if (args_info->nbins_y_given) | 
| 389 |  | write_into_file(outfile, "nbins_y", args_info->nbins_y_orig, 0); | 
| 390 | + | if (args_info->nbins_z_given) | 
| 391 | + | write_into_file(outfile, "nbins_z", args_info->nbins_z_orig, 0); | 
| 392 |  | if (args_info->nanglebins_given) | 
| 393 |  | write_into_file(outfile, "nanglebins", args_info->nanglebins_orig, 0); | 
| 394 |  | if (args_info->length_given) | 
| 429 |  | write_into_file(outfile, "r_theta", 0, 0 ); | 
| 430 |  | if (args_info->r_omega_given) | 
| 431 |  | write_into_file(outfile, "r_omega", 0, 0 ); | 
| 432 | + | if (args_info->r_z_given) | 
| 433 | + | write_into_file(outfile, "r_z", 0, 0 ); | 
| 434 |  | if (args_info->theta_omega_given) | 
| 435 |  | write_into_file(outfile, "theta_omega", 0, 0 ); | 
| 436 |  | if (args_info->gxyz_given) | 
| 513 |  | args_info->gofz_given = 0 ; | 
| 514 |  | args_info->r_theta_given = 0 ; | 
| 515 |  | args_info->r_omega_given = 0 ; | 
| 516 | + | args_info->r_z_given = 0 ; | 
| 517 |  | args_info->theta_omega_given = 0 ; | 
| 518 |  | args_info->gxyz_given = 0 ; | 
| 519 |  | args_info->p2_given = 0 ; | 
| 773 |  | { "nbins",      1, NULL, 'b' }, | 
| 774 |  | { "nbins_x",    1, NULL, 'x' }, | 
| 775 |  | { "nbins_y",    1, NULL, 'y' }, | 
| 776 | + | { "nbins_z",    1, NULL, 0 }, | 
| 777 |  | { "nanglebins", 1, NULL, 'a' }, | 
| 778 |  | { "length",     1, NULL, 0 }, | 
| 779 |  | { "rcut",       1, NULL, 'c' }, | 
| 794 |  | { "gofz",       0, NULL, 0 }, | 
| 795 |  | { "r_theta",    0, NULL, 0 }, | 
| 796 |  | { "r_omega",    0, NULL, 0 }, | 
| 797 | + | { "r_z",        0, NULL, 0 }, | 
| 798 |  | { "theta_omega",        0, NULL, 0 }, | 
| 799 |  | { "gxyz",       0, NULL, 0 }, | 
| 800 |  | { "p2", 0, NULL, 'p' }, | 
| 995 |  | break; | 
| 996 |  |  | 
| 997 |  | case 0: /* Long option with no short option */ | 
| 998 | + | /* number of bins in z axis.  */ | 
| 999 | + | if (strcmp (long_options[option_index].name, "nbins_z") == 0) | 
| 1000 | + | { | 
| 1001 | + |  | 
| 1002 | + |  | 
| 1003 | + | if (update_arg( (void *)&(args_info->nbins_z_arg), | 
| 1004 | + | &(args_info->nbins_z_orig), &(args_info->nbins_z_given), | 
| 1005 | + | &(local_args_info.nbins_z_given), optarg, 0, "100", ARG_INT, | 
| 1006 | + | check_ambiguity, override, 0, 0, | 
| 1007 | + | "nbins_z", '-', | 
| 1008 | + | additional_error)) | 
| 1009 | + | goto failure; | 
| 1010 | + |  | 
| 1011 | + | } | 
| 1012 |  | /* maximum length (Defaults to 1/2 smallest length of first frame).  */ | 
| 1013 | < | if (strcmp (long_options[option_index].name, "length") == 0) | 
| 1013 | > | else if (strcmp (long_options[option_index].name, "length") == 0) | 
| 1014 |  | { | 
| 1015 |  |  | 
| 1016 |  |  | 
| 1251 |  | goto failure; | 
| 1252 |  |  | 
| 1253 |  | } | 
| 1254 | + | /* g(r, z).  */ | 
| 1255 | + | else if (strcmp (long_options[option_index].name, "r_z") == 0) | 
| 1256 | + | { | 
| 1257 | + |  | 
| 1258 | + | if (args_info->staticProps_group_counter && override) | 
| 1259 | + | reset_group_staticProps (args_info); | 
| 1260 | + | args_info->staticProps_group_counter += 1; | 
| 1261 | + |  | 
| 1262 | + | if (update_arg( 0 , | 
| 1263 | + | 0 , &(args_info->r_z_given), | 
| 1264 | + | &(local_args_info.r_z_given), optarg, 0, 0, ARG_NO, | 
| 1265 | + | check_ambiguity, override, 0, 0, | 
| 1266 | + | "r_z", '-', | 
| 1267 | + | additional_error)) | 
| 1268 | + | goto failure; | 
| 1269 | + |  | 
| 1270 | + | } | 
| 1271 |  | /* g(cos(theta), cos(omega)).  */ | 
| 1272 |  | else if (strcmp (long_options[option_index].name, "theta_omega") == 0) | 
| 1273 |  | { |