rate8(u,p,t) = (1-(mod(t, 3600) > 1800)) * u[7] * p["H_tet"]
function repressor_prod!(integrator)
    integrator.u[5] += 1
end

rate9(u,p,t) = p["P_lon"] u[5] / (p["K_lon"] * (1 + u[4]/p["K_lon"] +u[5]/p["K_lon"]))
rate5_delt(u,p,t) = p["delt"] * u[5]
function repressor_deg!(integrator)
    integrator.u[5] -= 1
end