function [figure1] = Plot_TrendMaps(T,TACC,TD,TACCD,Lcg)

lon = [-179.5:179.5]';
lat = [-89.5:89.5]';
[lat,lon] = meshgrid(lat,lon);
L1(:,1) = reshape(lon,size(lon,1)*size(lon,2),1);
L1(:,2) = reshape(lat,size(lon,1)*size(lon,2),1);
idx = knnsearch(Lcg,L1);

T1 = NaN(size(L1,1),1);
T1 = T(idx,1)*1000;
for i = 1:numel(idx);
    [d1km(i,1) d2km(i,1)]=lldistkm(L1(i,:),Lcg(idx(i),:));
end
sdis = find(d1km>120);
T1(sdis) = NaN;

T1 = reshape(T1,size(lon));

figure1 = figure;

axes1 = axes('Parent',figure1,...
    'Position',[0.1 0.64 0.39 0.27]);
m_proj('robinson','lon',[-180 180]);
m_pcolor(lon,lat,T1);shading flat
m_coast('patch',[1 1 1],'edgecolor','k');
m_grid('XTick',[],'YTick',[]);
caxis([-7 7])
load('roma.mat')
colormap(roma);
title('Trend [mmyr^{-1}]','Fontsize',10)
ylabel('Altimetry','Fontsize',10)

T2 = NaN(size(L1,1),1);
T2 = T(idx,2)*1000;
T2(sdis) = NaN;
T2 = reshape(T2,size(lon));

axes2 = axes('Parent',figure1,...
    'Position',[0.1 0.38 0.39 0.27]);
m_proj('robinson','lon',[-180 180]);
m_pcolor(lon,lat,T2);shading flat
m_coast('patch',[1 1 1],'edgecolor','k');
m_grid('XTick',[],'YTick',[]);
caxis([-7 7])
load('roma.mat')
colormap(roma);
ylabel('Kalman Smoother','Fontsize',10)

T2 = NaN(size(L1,1),1);
T2 = (TD(idx,1))*1000;
T2(sdis) = NaN;
T2 = reshape(T2,size(lon));
[x1,y1] = m_ll2xy(L1(:,1),L1(:,2));

axes3 = axes('Parent',figure1,...
    'Position',[0.1 0.12 0.39 0.27]);
m_proj('robinson','lon',[-180 180]);
m_pcolor(lon,lat,T2);shading flat
m_coast('patch',[1 1 1],'edgecolor','k');
m_grid('XTick',[],'YTick',[]);
caxis([-7 7])
load('roma.mat')
colormap(roma);
cb = colorbar('YTick',[-6:2:6],'Fontsize',8,'Location','Southoutside')
ylabel('Differences','Fontsize',10)

T2 = NaN(size(L1,1),1);
T2 = TACC(idx,1)*1000;
T2(sdis) = NaN;
T2 = reshape(T2,size(lon));

axes4 = axes('Parent',figure1,...
    'Position',[0.5 0.64 0.39 0.27]);
m_proj('robinson','lon',[-180 180]);
m_pcolor(lon,lat,T2);shading flat
m_coast('patch',[1 1 1],'edgecolor','k');
m_grid('XTick',[],'YTick',[]);
caxis([-.7 .7])
load('roma.mat')
colormap(roma);
title('Acceleration [mmyr^{-2}]','Fontsize',10)

T2 = NaN(size(L1,1),1);
T2 = TACC(idx,2)*1000;
T2(sdis) = NaN;
T2 = reshape(T2,size(lon));

axes5 = axes('Parent',figure1,...
    'Position',[0.5 0.38 0.39 0.27]);
m_proj('robinson','lon',[-180 180]);
m_pcolor(lon,lat,T2);shading flat
m_coast('patch',[1 1 1],'edgecolor','k');
m_grid('XTick',[],'YTick',[]);
caxis([-.7 .7])
load('roma.mat')
colormap(roma);

T2 = NaN(size(L1,1),1);
T2 = (TACCD(idx,1))*1000;
T2(sdis) = NaN;
T2 = reshape(T2,size(lon));

axes6 = axes('Parent',figure1,...
    'Position',[0.5 0.12 0.39 0.27]);
m_proj('robinson','lon',[-180 180]);
m_pcolor(lon,lat,T2);shading flat
hold on
m_coast('patch',[1 1 1],'edgecolor','k');
m_grid('XTick',[],'YTick',[]);
caxis([-.7 .7])
load('roma.mat')
colormap(roma);
cb = colorbar('YTick',[-.6:.2:.6],'Fontsize',8,'Location','Southoutside')




end