include enclosed_gas_mass.h routine to calculing enclosed gas mass