Towards Adaptive Refinement for Multivariate Functional Approximation of Scientific Data