Pupa first undergo one extra day of survival, calculated as $$\overline{P_{[t-1]}} * (1-\mu_{ad})$$. This is an artifact of the conversion from continuous to discrete time (as mentioned in the original Hancock paper this model is derived from).
Then, pupation into adult males is calculated as $$(1-\overline{\phi}) * \overline{P_{[t]}}$$ and into adult females as $$\overline{\phi} * \overline{P_{[t]}}$$

oneDay_pupation_deterministic_Patch()