use subdivision_rate (P_rate) instead of constant 2