def METHOD_NAME(poro, densMatrix=2510, densFluid=1000, satur=1,
            mesh=None):
    r"""
        densMatrix, densFluid in kg/m^3
    """
    poro = parseArgToArray(poro, mesh.cellCount(), mesh)
    densMatrix = parseArgToArray(densMatrix, mesh.cellCount(), mesh)
    densFluid = parseArgToArray(densFluid, mesh.cellCount(), mesh)
    satur = parseArgToArray(satur, mesh.cellCount(), mesh)
    dens = np.array(densMatrix * (1.-poro)) + densFluid * poro * satur
    return dens