function rho = rhol(d)
rho= 0.64582/d.MW*1e6; % mol/m^3
