$\mathrm{k1\_}\mathrm{Ca}+\mathrm{k2\_}\mathrm{Ci}+D-\mathrm{k1}\mathrm{Ra}P+\mathrm{k2}\mathrm{Ri}P+\mathrm{kcl}P$ $\mathrm{k1\_}\mathrm{Ca}+\mathrm{k3}\mathrm{Ri}-\mathrm{k1}\mathrm{Ra}P+\mathrm{k3\_}\mathrm{Ra}$ $\mathrm{k2\_}\mathrm{Ci}+\mathrm{k3\_}\mathrm{Ra}-\mathrm{k2}\mathrm{Ri}P+\mathrm{k3}\mathrm{Ri}$ $\mathrm{k1}\mathrm{Ra}P+\mathrm{k4}\mathrm{Ci}-\mathrm{k1\_}\mathrm{Ca}+\mathrm{k4\_}\mathrm{Ca}$ $\mathrm{k2}\mathrm{Ri}P+\mathrm{k4\_}\mathrm{Ca}-\mathrm{k2\_}\mathrm{Ci}+\mathrm{k4}\mathrm{Ci}$ $\frac{\mathrm{Ra}+\mathrm{Ca}}{\mathrm{Ra}+\mathrm{Ca}+\mathrm{Ri}+\mathrm{Ci}}$ $\mathrm{De}+\mathrm{Dd}$ $\begin{cases}\mathrm{dmin} & \text{if $\frac{\mathrm{time\_hour}}{\mathrm{cycle\_length}}\ge j$}\\ \mathrm{dmax} & \text{if $(\mathrm{time\_hour}-\lfloor \frac{\mathrm{time\_hour}}{\mathrm{cycle\_length}}\rfloor \mathrm{cycle\_length}\ge 0)\land (\mathrm{time\_hour}-\lfloor \frac{\mathrm{time\_hour}}{\mathrm{cycle\_length}}\rfloor \mathrm{cycle\_length}< \mathrm{tau\_on})$}\\ \mathrm{dmin} & \text{if $(\mathrm{time\_hour}-\lfloor \frac{\mathrm{time\_hour}}{\mathrm{cycle\_length}}\rfloor \mathrm{cycle\_length}\ge \mathrm{tau\_on})\land (\mathrm{time\_hour}-\lfloor \frac{\mathrm{time\_hour}}{\mathrm{cycle\_length}}\rfloor \mathrm{cycle\_length}< \mathrm{cycle\_length})$}\end{cases}$ $\mathrm{tau\_on}+\mathrm{tau\_off}$ $\frac{\mathrm{time\_}\times 1}{3600}$