Added unit system definition. Default is kPa and kJ for backward compatibility.

SI can be set with the following command:
set_standard_unit_system(UNIT_SYSTEM_SI)

This should be implemented at a higher level when brines and other fluids take SI for internal calculation.