in calc_dRRH (H4), rho 1025.. surely it should be 1000 like in the other scripts and the rules?