mod.ncstr,
initialize=3.92336059452774350120307644829154E+02)
mod.Tj = Var(mod.t,
mod.ncstr,
initialize=3.77995395658401662331016268581152E+02)
mod.k = Var(mod.t, mod.ncstr, initialize=4.70706140E+02)
mod.kdef = Constraint(mod.t, mod.ncstr)
#: 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)
