def get_points_data(dryspells, var1, var2, weights, points):
    """
    Extract dryspells, full time-series of var1 and var2, and averaging
    weights for a list of land points from input data on the full domain.

    The variables are typically near-surface air temperature and land
    surface temperature, and the weights are 1 km MODIS pixel counts.
    The full domain is typically global.

    Parameters
    ----------
    dryspells : list of list of <event.Event> instances.
        Descriptions of all dry spells for each land point on the full
        domain.
    var1: MaskedArray, shape(time, land)
        First input variable.
    var2: MaskedArray, shape(time, land)
        Second input variable.
    weights : ndarray, shape(time, land)
        Weights used in calculating the composite means.
    points : list of ints
       Zero-based list of land points to be extracted.

    Returns
    -------
    dryspells_p : list of list of <event.Event> instances.
        Descriptions of all dry spells for each requested land point.
    var1_p : list of MaskedArrays
        Full time series of variable 1 for each requested land point.
    var2_p : list of MaskedArrays
        Full time series of variable 2 for each requested land point.
    weights_p : list of ndarrays
        Full time series of weights for each requested land point.
    """

    dryspells_p = []
    var1_p = []
    var2_p = []
    weights_p = []

    for p in points:
        dryspells_p.append(dryspells[p])
        var1_p.append(var1[:, p])
        var2_p.append(var2[:, p])
        weights_p.append(weights[:, p])

    return dryspells_p, var1_p, var2_p, weights_p