clear all;format compact

% coast
load coastlines;
lon2=coastlon;
lat2=coastlat;

lon1=coastlon;
lat1=coastlat;
for i=1:length(lon1)
  if(lon1(i)<0)
    lon1(i)=lon1(i)+360;
  end
  if((lon1(i)<=1)||(lon1(i)>=359))
    lon1(i)=NaN;
  end
end
clear coastlon coastlat lon2 lat2
% coast


file1 = '../seeding_masks/mask_0375_025-200_SSP245.nc';

xx = {'7.5 PCT';'10 PCT';'12.5 PCT';'15 PCT'};

f0 = '../seeding_masks/b.e21.BSSP245cmip6.f09_g17.CMIP6-MCB-cntl.000.SWCF.nc';
lon = ncread(f0,'lon');
lat = ncread(f0,'lat');

tmp=ncread(file1,'mask');
ii = [1:5:36];
for i=1:length(ii);
    mask(i,:,:,:) = squeeze(tmp(:,:,:,ii(i)));
end
clear tmp


ctb1=colormap(jet(13));
ctb1(1,:) = [1 1 1];


ctb1(13,:) = [.7 0 0 ];
ctb1(12,:) = [1 0 1];
ctb1(11,:) = [1 .2 0];

tt={'(a) 2.5% ocean area';'(b) 5.0% ocean area';'(c) 7.5% ocean area';...
    '(d) 10.0% ocean area';'(e) 12.5% ocean area';'(f) 15.0% ocean area';...
    '(g) 17.5% ocean area';'(h) 20.0% ocean area'}; 

ppp = [.1 .4 .2 .25;.32 .4 .2 .25;.54 .4 .2 .25;.76 .4 .2 .25;...
       .1 .1 .2 .25;.32 .1 .2 .25;.54 .1 .2 .25;.76 .1 .2 .25];

v=[0:1:13];

figure;orient tall;orient landscape;colormap(ctb1)
for i=1:8
subplot('position',ppp(i,:))
[c,h]=contourf(lon,lat,squeeze(sum(mask(i,:,:,:),4))',v);
set(gca,'clim',[0 13]);
    set(h,'edgecolor','none')
    hold on;plot(lon1,lat1,'k-')
    axis([0 360 -90 90])
    set(gca,'xtick',[0:60:360],'xminortick','on','xticklabel',[])
    set(gca,'ytick',[-90:30:90],'yminortick','on','yticklabel',[])
    set(title(char(tt(i))),'fontsize',8)      
end
h=colorbar('horiz');
set(h,'position',[0.2  0.07 0.6  0.01])