def nbi_growth(nbi, nvi, nb, nb0, nv, C, f, R, g, c0, alpha, e, pv, eta):
    """
    Single clone growth rate for bacteria.
    Inputs:
    nbi : 1D vector in time of bacteria clone sizes
    nvi : 1D vector in time of corresponding phage clone sizes
    nb : 1D vector in time of total bacteria population size
    nb0 : 1D vector in time of number of bacteria without spacers
    nv : 1D vector in time of total phage population size
    C : 1D vector in time of total nutrients
    f, g, c0, alpha, e, pv, eta : simulation parameters

    Output: 
    s : 1D vector of phage clone growth rate per bacterial generation (not multiplied by population size)
    """
    F = f*g*c0
    r = R*g*c0

    s = g*C - F - r - alpha*pv*(nv - e*nvi) + alpha*eta*nb0*nvi*(1-pv)/nbi
        
    return s/(g*c0)