Recipe: Matlab program fills the seam3

function [lons3 lats3 data3]=exendFace3(lons,lats,cubed_data)
% This side faces north, it includes extra four southern neighbors
     N=size(lons,1);

     lats3(N+2,N+2)=0;
     lons3(N+2,N+2)=0;
     data3(N+2,N+2)=0;

     lats3(2:N+1,2:N+1)=lats(:,:,3);
     lons3(2:N+1,2:N+1)=lons(:,:,3);
     data3(2:N+1,2:N+1)=cubed_data(:,:,3);

     lats3(1,2:N+1)=flipud(lats(:,N,1)) ;
     lons3(1,2:N+1)=flipud(lons(:,N,1)) ;
     data3(1,2:N+1)=flipud(cubed_data(:,N,1)) ;
     lats3(N+2,2:N+1)=lats(1,:,4) ;
     lons3(N+2,2:N+1)=lons(1,:,4) ;
     data3(N+2,2:N+1)=cubed_data(1,:,4);

     lats3(2:N+1,1) = lats(:,N,2);
     lons3(2:N+1,1) = lons(:,N,2);
     data3(2:N+1,1) = cubed_data(:,N,2);

     lats3(2:N+1,N+2) = fliplr(lats(1,:,5));
     lons3(2:N+1,N+2) = fliplr(lons(1,:,5));
     data3(2:N+1,N+2) = fliplr(cubed_data(1,:,5));

     lats3(1,1) =lats(N,N,1);
     lons3(1,1) =lons(N,N,1);
     data3(1,1) =cubed_data(N,N,1);

     lats3(1,N+2) =lats(1,N,1);
     lons3(1,N+2) =lons(1,N,1);
     data3(1,N+2) =cubed_data(1,N,1);

     lats3(N+2,1) =lats(1,1,4);
     lons3(N+2,1) =lons(1,1,4);
     data3(N+2,1) =cubed_data(1,1,4);

     lats3(N+2,N+2) =lats(1,1,5);
     lons3(N+2,N+2) =lons(1,1,5);
     data3(N+2,N+2) =cubed_data(1,1,5);