Changing:
function prob_func(prob, i, repeat)
rlon = deg2rad(puffs[i][1])
rlat = deg2rad(puffs[i][2])
rlev = puffs[i][3]
u0 = [
sys.Puff₊lon => rlon,
sys.Puff₊lat => rlat,
sys.Puff₊lev => rlev,
sys.ElementalCarbon₊EC => ec_0
]
p = [
sys.GaussianKC₊Δz => 100,
]
ts = (tspan[1] + floor((i-1) / puff_release_per_interval) * puff_release_interval, tspan[2])
remake(prob; u0 = u0, tspan = ts, p = p)
end
to:
set_params = setp(prob,[sys.GaussianKC₊Δz])
set_u0 = setsym(prob, [sys.Puff₊lon, sys.Puff₊lat, sys.Puff₊lev, sys.ElementalCarbon₊EC,])
function prob_func(prob, i, repeat)
rlon = deg2rad(puffs[i][1])
rlat = deg2rad(puffs[i][2])
rlev = puffs[i][3]
ts = (tspan[1] + floor((i-1) / puff_release_per_interval) * puff_release_interval, tspan[2])
set_params(prob, (100))
set_u0(prob, (rlon, rlat, rlev, ec_0))
prob, ts
end
results in this error:
ERROR: UndefVarError: setp not defined in Main
Suggestion: check for spelling errors or missing imports.
Hint: a global variable of this name may be made accessible by importing SymbolicIndexingInterface in the current active module Main
Stacktrace:
[1] top-level scope
Changing:
function prob_func(prob, i, repeat)
rlon = deg2rad(puffs[i][1])
rlat = deg2rad(puffs[i][2])
rlev = puffs[i][3]
end
to:
set_params = setp(prob,[sys.GaussianKC₊Δz])
set_u0 = setsym(prob, [sys.Puff₊lon, sys.Puff₊lat, sys.Puff₊lev, sys.ElementalCarbon₊EC,])
function prob_func(prob, i, repeat)
rlon = deg2rad(puffs[i][1])
rlat = deg2rad(puffs[i][2])
rlev = puffs[i][3]
end
results in this error:
ERROR: UndefVarError:
setpnot defined inMainSuggestion: check for spelling errors or missing imports.
Hint: a global variable of this name may be made accessible by importing SymbolicIndexingInterface in the current active module Main
Stacktrace:
[1] top-level scope