On The Modeling Of Multiple Functional Outcomes With Spatially Heterogeneous Shape Characteristics