def METHOD_NAME(mesh, nodes=False, seed=None, **kwargs):
    """Generate geostatistical model (cell or node) for given mesh.
    Parameters
    ----------
    mesh : gimliapi:`GIMLI::Mesh`
        Mesh
    nodes : bool [False]
        use node positions, otherwise (default) cell centers are used
    seed : int, array_like[ints], SeedSequence, BitGenerator, Generator}, optional
        A seed to initialize the BitGenerator. If None, then fresh, unpredictable 
        entropy will be used. The `seed` variable is passed to :func:`numpy.random.default_rng`
    **kwargs
        I : float or list of floats
            correlation lengths (range) in individual directions
        dip : float
            dip angle of major axis (I[0])
        strike : float
            strike angle (for 3D)
    Returns
    -------
    res : np.array of size cellCount or nodeCount (nodes=True)
    """
    rng = np.random.default_rng(seed=seed)
    return rng.multivariate_normal(np.ones(mesh.cellCount()),
                                   covarianceMatrix(mesh, nodes=nodes, **kwargs))