def METHOD_NAME(out, stream=False):
    mesh = pg.load(out + '.bms')
    meshC = pg.load(out + 'C.bms')
    densMatrix = pg.load(out + 'density.bmat')
    vels = np.load(out+'velo.bmat.npy')
    cMin= 9e99
    cMax= -9e99
    for d in densMatrix:
        cMin = min(cMin, pg.min(d))
        cMax = max(cMax, pg.max(d))
    if stream:
        pg.viewer.mpl.saveAnimation(mesh, densMatrix, out + '-stream',
                               plc=None, label='',
                               cMin=cMin, cMax=cMax, logScale=False, cmap=None,
                               vData=vels, coarseMesh=meshC, color='white'
                               )
    else:
        pg.viewer.mpl.saveAnimation(mesh, densMatrix, out,
                               plc=None, label='',
                               cMin=cMin, cMax=cMax, logScale=False, cmap=None)