function vq = upsampling(globalx,scalar)
    [nely,nelx] = size(globalx);
    F = griddedInterpolant(globalx);
    xq = (0 : 1/scalar : nely)';
    yq = (0 : 1/scalar : nelx)';
    vq = F({xq,yq});
end