def METHOD_NAME(
    network,
    throat_diameter='throat.diameter',
    throat_length='throat.length',
):
    r"""
    Calculate throat volume assuing a square cross-section
    Parameters
    ----------
    %(network)s
    %(Dt)s
    %(Lt)s
    Returns
    -------
    Notes
    -----
    At present this models does NOT account for the volume reprsented by the
    intersection of the throat with a spherical pore body.
    """
    leng = network[throat_length]
    diam = network[throat_diameter]
    value = leng*diam**2
    return value