This paper describes a local method for the construction of differenti
able functions interpolating a set of data connected in a general tens
er-product topology, These functions are shape-preserving in the sense
that they locally maintain the monotonicity of the data along curves
''parallel'' to the related coordinates lines. Error estimates and gra
phical examples are provided.