function field_large_out = merge_small_large (field_large_in,field_small,...
    Nbcs,Node_in,Node,indices_in)

%-------------------------------------------------------------------------------
% Reshape large domain.
%-------------------------------------------------------------------------------

field_large_2D = field_large_in;
field_large_2D(1:Nbcs) = NaN;
field_large_2D = field_large_2D(Node_in);

%-------------------------------------------------------------------------------
% Reshape small domain.
%-------------------------------------------------------------------------------

field_small_2D = field_small;
field_small_2D(1:Nbcs) = NaN;
field_small_2D = field_small_2D(Node);

%-------------------------------------------------------------------------------
% Put small domain data in large domain.
%-------------------------------------------------------------------------------

[Ny, Nx] = size(field_small_2D);

field_large_2D(1:Ny,1:Nx) = field_small_2D;

field_large_out                              = field_large_in;
field_large_out(Nbcs+(1:length(indices_in))) = field_large_2D(indices_in);

