gcmc functions now take in pressure, in units of bar. the pressure is converted to fugacity (in units Pa) internally using either the ideal gas or @huynmela's PREOS code. So except for internally, everything is in bar now. Even querying results["fugacity (bar)"].
pass eos=:ideal or eos=:PengRobinson to convert pressure to fugacity.
results["fugacity (bar)"].eos=:idealoreos=:PengRobinsonto convert pressure to fugacity.