def METHOD_NAME(T):
    # Initial values of the carbonates and mass of water with nh4-accetate
    n0Calcite = 10.0
    water_kg = 1.0
    # Define initial equilibrium state corresponding to the seawater
    state_sw = ChemicalState(system)
    state_sw.temperature(T, "celsius")
    state_sw.pressure(1.0, "atm")
    state_sw.set("H2O", water_kg, "kg")
    state_sw.set("NH4+", 1, "mol")
    # Calculate chemical state corresponding to the seawater
    res = solver.solve(state_sw)
    # Update aqueous properties
    aqprops.update(state_sw)
    print("pH (before adding Calcite) = ", aqprops.pH())
    # Throw exception if the equilibrium couldn't be found
    assert res.succeeded(), "The equilibrium calculation did not succeed!"
    # Add carbonates
    state_sw.set("Calcite", n0Calcite, "mol")
    # Equilibrate the seawater with carbonates
    solver.solve(state_sw)
    # Updated aqueous properties
    aqprops.update(state_sw)
    print("pH (after adding Calcite)  = ", aqprops.pH())
    # Fetch values of the specified species
    nCalcite = state_sw.speciesAmount("Calcite")
    nCa2 = state_sw.speciesAmount("Ca+2")
    pH = aqprops.pH()
    return (pH[0], nCalcite[0], nCa2[0])