Обратные вызовы в Джулии, когда я использую реакционную сеть

Я новичок в использовании Джулии и пытаюсь смоделировать сеть реакций, в то время как возмущение меняет концентрацию одного вида с течением времени. Я узнал, что могу использовать обратный вызов, чтобы изменить ценность определенного вида в определенный момент времени, чтобы имитировать это внешнее возмущение. Хотя я прочитал документацию, но не смог найти, как это сделать. Вот код, в котором я попытался реализовать эту симуляцию. Сначала я определил реакционную сеть и ее параметры:

using Catalyst
using DifferentialEquations

rn = @reaction_network begin
α, A + B --> 2B
β, B --> A
end α,β

p     = [:α => 1, :β => 2]
tspan = (0.0,20.0)
u0    = [:A => 5.0, :B=> 5.0]
op    = ODEProblem(rn, u0, tspan, p)

Далее, для обратного вызова, который срабатывает в момент t=10

condition(u, t, integrator) = t == 10

увеличить ценность вида А в два раза

affect!(integrator) = integrator.??? += 2

Но здесь я не уверен, какие параметры использовать для изменения значения A.

И остальная часть моего кода:

cb = DiscreteCallback(condition, affect!)

sol = solve(prob, Tsit5(), callback = cb, tstops = [10.0])

using Plots
plot(sol) 

Буду признателен за любую помощь в улучшении моего понимания в этой области.