added subroutine to calculate chemical equilibrium

need to check if mu+phi_e+phi_g=const
added pressure scale height to chem_equil routine
added lib directory, created module for all energy density components