def METHOD_NAME(vmin, vmax, num_per_mag=10):
    """Initiate pseudo-log discrete colormap levels
    Note
    ----
        This is a beta version and aims to
    Parameters
    ----------
    vmin : float
        lower end of colormap (e.g. minimum value of data)
    vmax : float
        upper value of colormap (e.g. maximum value of data)
    """
    high = float(exponent(vmax))
    low = -3.0 if vmin == 0 else float(exponent(vmin))
    lvls = []
    if 1 % vmax * 10 ** (-high) == 0:
        low += 1
        high -= 1
    for mag in range(int(low), int(high)):
        lvls.extend(np.linspace(1, 10, num_per_mag - 1, endpoint=0) * 10 ** (mag))
    lvls.extend(np.linspace(1, vmax * 10 ** (-high), num_per_mag, endpoint=1) * 10 ** (high))
    return lvls