[Get max proper order and reflection for an axis]

    .. todo:: Complete geom_parse_axis docstring

def geom_check_axis(g, atwts, ax,
        nmax=_DEF.SYMM_MATCH_NMAX,
        tol=_DEF.SYMM_MATCH_TOL):
    """ [Get max proper order and reflection for an axis]

    .. todo:: Complete geom_parse_axis docstring

    """

    # Imports
    import numpy as np

    # Store the max found rotation order of the geometry.
    order = geom_find_rotsymm(g, atwts, ax, \
                                        False, nmax, tol)[0]

    # Store the presence/absence of a reflection plane.
    refl = geom_symm_match(g, atwts, ax, 0, True) < tol

    # Return the pair of values for outside handling
    return order, refl