function [figure1] = Plot_DataAvailability_Map(t,M,L,regionTG,LALT,regionSL)

regionTG(regionTG==7) = 3;
Avail = NaN(size(M,1),max(regionTG));
for i = 1:max(regionTG);
    s = find(regionTG==i);
    for j = 1:size(M,1);
        isn = find(~isnan(M(j,s)));
        if length(isn)>=1;
            Avail(j,i) = length(isn);
        end
    end
end

figure1 = figure;
axes1 = axes('Parent',figure1,...
    'Position',[0.1 0.7 0.68 0.239]);
m_proj('robinson','lon',[-180 180],'lat',[-90 90]);
[x1,y1] = m_ll2xy(LALT(:,1),LALT(:,2));
sreg = find(regionSL<7);
scatter(x1(sreg),y1(sreg),2,regionSL(sreg),'filled')
hold on
[x1,y1] = m_ll2xy(L(:,1),L(:,2));
load('roma.mat');
cm = colormap(roma(linspace(1,size(roma,1),6),:));
m_coast('patch',[1 1 1],'edgecolor',[0 0 0]);
m_grid('XTick',[],'YTick',[],'Fontsize',8);
Mar = {'o','square','diamond','^','v','<'}
for i = 1:max(regionTG);
    s = find(regionTG==i);
    scatter(x1(s),y1(s),12,regionTG(s),'filled','Marker',Mar{i},'MarkeredgeColor','k');
end
title('# of Tide Gauge: 516','Fontsize',10,'Fontweight','bold')
axes1a = axes('Parent',figure1,...
    'Position',[0.75 0.7 0.1 0.239],'visible','off','xlim',[0 13],'ylim',[0 13]);
hold all
Nregion = {'South Atlantic','Indo-Pacific','Northwest Pacific','Subpolar North Atlantic',...
    'Subtropical North Atlantic','East Pacific'};
for i = 1:max(regionTG);
    scatter(1,13-i,12,i,'filled','Marker',Mar{i},'MarkeredgeColor','k');
    hold on 
    text(2,13-i,Nregion(i),'Fontsize',8);
end

axes2 = axes('Parent',figure1,'Position',[0.2 0.55 0.55 0.1],...
    'xlim',[1900 2021],'Ylim',[0 500],'YTick',[100:100:400],'YTicklabel',[100:100:400],...
    'XTick',[1920:20:2000],'XTicklabel',[1920:20:2000],'Layer','top',...
    'Fontsize',8,'box','off','Yaxislocation','right');
ylabel('# Tide Gauges','Fontsize',10);
xlabel('Year','Fontsize',10);
title('Temporal Availability','Fontsize',10,'Fontweight','bold');
hold all
a = area(t,Avail)
for i = 1:6;
    a(i).FaceColor = cm(i,:);
end


end