def METHOD_NAME(mesh, dDens, gravPoints, dz):
    dpi=92
    orientation = 'horizontal'
    fig = plt.figure(facecolor='white', figsize=(2*800/dpi, 2*490/dpi), dpi=dpi)
    axDDe = fig.add_subplot(3,1,1)
    axGra = fig.add_subplot(3,1,2)
    #axGra = fig.fig.add_subplot(1,3,2)
        # ** Density **
    gciDDe = pg.viewer.mpl.drawModel(axDDe, mesh, data=dDens[1],
                                    cMin=0, cMax=20,
                                    )
    cbar = createColorbar(gciDDe, orientation=orientation,
                          label='Delta density in kg/m$^3$')
    def ani(i):
        axGra.clear()
        axGra.plot(pg.x(gravPoints), dz[i])
        axGra.plot(pg.x(gravPoints), pg.y(gravPoints), 'v', color='black')
        axGra.set_ylabel('Grav in mGal')
        axGra.set_xlim((-20, 20))
        axGra.set_ylim((0, 0.001))
        axGra.grid()
        pg.viewer.mpl.setMappableData(gciDDe, abs(dDens[i]),
                                     cMin=0, cMax=20,
                                     logScale=False)
    for i in range(len(dz)):
        ani(i)
        plt.pause(0.001)