def default_qp(depth, qs):
    """Default rule for Qp as a function of Qs.

    Args:
        depth (float)
            Depth of location in m.
        qs (float)
            Quality factor for S wave.

    Returns:
        Quality factor for P wave.
    """
    return 2.0*qs