change_pr_units only works for canonical units

The change_pr_units function assumes incoming data in canonical units. Also, it seems to always use 1s as the integration time and the density of water as density when changing between amount based and rate based units, or mass based and thickness based units, respectively.