Partially monotone tensor spline estimation of the joint distribution