function [xy_local_box,z_box] = extrapolate_local_new(xy_local,z,x_res,y_res,n_pixels_boundary,R) % [xy_local_box,z_box] = extrapolate_local(xy_local,z,x_res,y_res) % x_res Output resolution in the ...
%% for now, just use nearest neighbor extrapolation. for i = SIM.mbc+MESH.NX(1):MESH.NX(1)+2*SIM.mbc uf{1}(i,:,:) = uf{1}(SIM.mbc+MESH.NX(1),:,:); uf{2}(i,:,:) = uf{2 ...