Get volume

import molsysmt as msm

molsys = msm.convert(msm.systems['T4 lysozyme L99A']['181l.h5msm'], selection='molecule_type=="protein"')

volume = msm.physchem.get_volume(molsys, definition='grantham')

view = msm.view(molsys)
msm.thirds.nglview.set_color_by_value(view, values=volume, min_value=0.0, cmap='Purples')

view