def lTdot(m, i, j, t):
if j > 0:
return m.Tdot[i, j, t] * m.hi_t[i] ==\
-(m.pm[i, j, t] - m.pn[i, j, t]) * m.dx_dt[i, j, t] / \
(m.x[i, j, t] *
exp(m.CapAm - m.CapBm/(m.T[i, j, t] + m.CapCm)) * m.CapBm/(m.T[i, j, t] + m.CapCm)**2 +
(1 - m.x[i, j, t]) *
exp(m.CapAn - m.CapBn/(m.T[i, j, t] + m.CapCn)) * m.CapBn/(m.T[i, j, t] + m.CapCn)**2)
else:
return Constraint.Skip
def lpn(m, i, j, t):
if j > 0:
return m.pn[i, j,
t] == exp(m.CapAn - m.CapBn / (m.T[i, j, t] + m.CapCn))
else:
return Constraint.Skip
def _lpm(m, i, j, t):
if j > 0:
return m.pm[i, j,
t] == exp(m.CapAm - m.CapBm / (m.T[i, j, t] + m.CapCm))
else:
return Constraint.Skip