Rate equation modelling of molecular multiphoton ionization dynamics