rewrite iapws95 calculate procedure