calculate Omega calcite for comparison Omegacal=Ca\_chamber\*CO3\_chamber/Kspcal print 'Omega calcite=', Omegacal Omegaara=Ca\_chamber\*CO3\_chamber/Ksp print 'Omega aragonite=', Omegaara
