function[Data]=mapGlobal2Local(Data)

% mapGlobal2Local: Function which maps the global degrees of freedom to the
% local ones for the INTERNODES matrix
% INPUT:
% Data:     Structure containing all the data parameters
% OUTPUT:
% Data:     Updated structure with the local degrees of freedom

Nf=Data.Nf;
[~, Data.body{1}.indx_local, ~]=intersect(Nf, Data.body{1}.interface_dof);
[~, Data.body{2}.indx_local, ~]=intersect(Nf, Data.body{2}.interface_dof);
end