Recipe: Matlab program fills the seam4

function [lons4 lats4 data4]=exendFace4(lons,lats,cubed_data)
% Face 4 includes extra east and west neighbors
     N=size(lons,1);

     lats4(N,N+2)=0.0;
     lons4(N,N+2)=0.0;
     data4(N,N+2)=0.0;

     lats4(:,2:N+1)=lats(:,:,4);
     lons4(:,2:N+1)=lons(:,:,4);
     data4(:,2:N+1)=cubed_data(:,:,4);

     lats4(:,N+2)=lats(:,1,5);
     lons4(:,N+2)=lons(:,1,5);
     data4(:,N+2)=cubed_data(:,1,5);
     lats4(:,1)=fliplr(lats(N,:,2));
     lons4(:,1)=fliplr(lons(N,:,2));
     data4(:,1)=fliplr(cubed_data(N,:,2));