Recipe: Matlab program fills the seam

function [lons1 lats1 data1]=exendFace1(lons,lats,cubed_data)
%face 1 includes extra east and west neighbors
     N=size(lons,1);

     lats1(N+2,N)=0.0;
     lons1(N+2,N)=0.0;
     data1(N+2,N)=0.0;

     lats1(2:N+1,:)=lats(:,:,1);
     lons1(2:N+1,:)=lons(:,:,1);
     data1(2:N+1,:)=cubed_data(:,:,1);
     lats1(1,:) = flipud(lats(:,N,5));
     lons1(1,:) = flipud(lons(:,N,5));
     data1(1,:) = flipud(cubed_data(:,N,5));
     lats1(N+2,:)=lats(1,:,2);
     lons1(N+2,:)=lons(1,:,2);
     data1(N+2,:)=cubed_data(1,:,2);