Calculation of simultaneous isotopic and trace element variations during water-rock interaction with applications to carbonate diagenesis