def METHOD_NAME(Nz, dz, cdt, n_order, v_comoving, use_galilean):
    """
    Return the stencil reach (in spatial space) for a given finite order
    stencil and for a given simulation setup with Nz cells and spacing dz.
    The stencil reach depends only slightly on the transverse k (kr) and will
    be calculated at a fixed kperp of 0.5 . The stencil reach is needed to
    define the number of guard cells between domains when running in parallel.
    Parameters:
    ----------
    Nz: int
        The number of cells in the logintudinal direction
    dz: float (microns)
        The cell spacing in the longitudinal direction
    cdt: float
        Timestep times speed of light
    n_order: int (multiple of 2)
        The (finite) order of the arbitrary order spectral,
        Maxwell (PSATD) solver, which defines the stencil reach.
    v_comoving: float or None, optional
        If this variable is None, the standard PSATD is used (default).
        Otherwise, the current is assumed to be "comoving",
        i.e. constant with respect to (z - v_comoving * t).
        This can be done in two ways: either by
        - Using a PSATD scheme that takes this hypothesis into account
        - Solving the PSATD scheme in a Galilean frame
    use_galilean: bool, optional
        Determines which one of the two above schemes is used
        When use_galilean is true, the whole grid moves
        with a speed v_comoving
    Returns:
    -------
    Number of cells needed for the stencil to decrease to
    machine precision at kperp = 0.5
    """
    # Calculate the real kz for the given grid (Nz)
    real_kz = 2 * np.pi * np.fft.fftfreq(Nz, d=dz)
    # Get the modified, finite order kz
    kz = get_modified_k(real_kz, n_order, dz=dz)
    # Calculate the stencil reach at an arbitrary kperp = 0.5
    # (Note: The stencil reach depends only weakly on kperp)
    return stencil_reach(kz, 0.5, cdt, v_comoving, use_galilean)