import chemics as cm

formula = input("Chemical Formula e.g CH4, H3N, ClH = ")
press = float(input("Process pressure kPaA = "))
temp = float(input("Process temperature degC = "))

mw = cm.mw(formula)    # molecular weight of N2 gas [g/mol]
p = press*1000      # pressure [Pa]
tk = temp+273.15        # temperature of gas [K]

rhog = cm.rhog(mw, p, tk)   # density of nitrogen gas [kg/m^3]

print(f'Density of {formula} gas = {rhog:.4f} kg/m^3')

#hc = cm.cp_gas(formula, tk, full=True)
#print(f'Heat capacity of {formula} gas = {hc[0]:.4f} J/(mol.K)')

results = cm.mu_gas(formula, tk, full=True)
tocP = results[0] * 0.0001
toPas= results[0] * 0.0000001
print(f'Absolute / Dynamic Viscosity of {formula} gas =\n{results[0]:.4f} µP or {tocP:.4f} cP or {toPas:.4g} Pa.s')