mod.Cp = Param(initialize=4.2)
mod.Vw = Param(initialize=10)
mod.rhow = Param(initialize=1000)
mod.Cpw = Param(initialize=4.2)
mod.k0 = Param(initialize=4.11e13)
mod.E = Param(initialize=76534.704)
mod.R = Param(initialize=8.314472)
mod.Er = Param(initialize=lambda m: (value(mod.E) / value(mod.R)))
mod.dH = Param(initialize=596619.)
mod.F = Param(mod.t, mutable=True, default=1.2000000000000000E+02)
mod.Fw = Param(mod.t, mutable=True, default=3.0000000000000000E+01)
# States
mod.Ca = Var(mod.t,
mod.ncstr,
initialize=1.60659680385930765667001907104350E-02)
mod.T = Var(mod.t,
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
