#: These guys have to be zero at the steady-state (steady).
zero0 = dict.fromkeys(mod.t * mod.ncstr)
for key in zero0.keys():
zero0[key] = 0.0
mod.steady = False
if mod.steady:
mod.Cadot = zero0
mod.Tdot = zero0
mod.Tjdot = zero0
else:
mod.Cadot = DerivativeVar(mod.Ca, initialize=-3.58709135E+01)
mod.Tdot = DerivativeVar(mod.T, initialize=5.19191848E+03)
mod.Tjdot = DerivativeVar(mod.Tj, initialize=-9.70467399E+02)
#: These guys as well (steady).
mod.Ca_ic = Param(mod.ncstr, default=1.9193793974995963E-02, mutable=True)
mod.T_ic = Param(mod.ncstr, default=3.8400724261199036E+02, mutable=True)
mod.Tj_ic = Param(mod.ncstr, default=3.7127352272578315E+02, mutable=True)
# m.Ca_ic = Param(m.ncstr, default=1.9193793974995963E-02)
# m.T_ic = Param(m.ncstr, default=3.8400724261199036E+02)
# m.Tj_ic = Param(m.ncstr, default=3.7127352272578315E+02)
mod.de_ca = Constraint(mod.t, mod.ncstr)
mod.de_T = Constraint(mod.t, mod.ncstr)
mod.de_Tj = Constraint(mod.t, mod.ncstr)
#: No need of these guys at steady.
if mod.steady:
mod.Ca_icc = None
mod.T_icc = None
