# Load the parameters database
param_database = CosmoRSParameterSets()

# Print the names of the available parameter sets
parameter_names = param_database.parameterSets()
nlprint('Available parameter set names:')
for name in parameter_names:
    nlprint(name)
nlprint('')

# Load a specific parameter set
parameters = param_database.loadParameterSet('klamt_1998')

# Calculate the chemical potential of water with the new parameter set
molecule_database = CosmoRSSpeciesDatabase()
water = molecule_database.exportSpecies('water')
solvent = CosmoRealSolvent(298*Kelvin, water, parameters)
chemical_potential = solvent.chemicalPotential(water)
nlprint(f'The estimated chemical potential of water is {chemical_potential}')