In this paper we present a simple automatic algorithm for constructing
C2 continuous planar parametric interpolants, which preserve the loca
l-convexity information contained in the interpolation point set and t
he imposed boundary conditions. The algorithm exploits the shape-prese
rving capabilities of a special family of polynomial splines of non-un
iform degree, introduced, for the case of functional data, in (Kaklis
and Pandelis, 1990). The numerical performance of the algorithm is dis
cussed for three data sets and various parameterizations.