.. function:: py__init_n_domains_p_pdaf

    Get number of analysis domains.

    Parameters
    ----------
    step: int
        Current step
    n_domains_p: int
        Number of analysis domains.

    Returns
    -------
    n_domains_p: int
        Number of analysis domains.