[PD] pmpd simulation and numerical output