%GBMAP georges bank map for vcr movie
% R. Limeburner Jan 95
% Color basemap of georges bank

clear;clg
load map

cg1b=[.97 .97 .97];
cg1a=[.94 .94 .94];

cgw=[.9 1 1];
cg1aa=[.6 1 1];
cg0=[.70 1 1];
cg1=[.1 .8 1];
cg2=[.1 .6 1];
cg3=[.21 .5 1];
cg4=[.21 .41 .95];
cg5=[.1 .1 .9];
cg6=[.05 .05 .75];
cgb=[0 0 0];
cgb2=[1 .92 .64];
%axes('position',[.16,.27,.70,.6]);
axes('position',[.11,.14,.83,.76]);
axis([-71 -65 40 43]);
%set(gcf,'position',[256 256 641 481]);
hold on

%a=fill(m4000(:,2),m4000(:,1),cg5);set(a,'edgecolor','none');
aaa=fill(m4000(:,2),m4000(:,1),cg6);set(aaa,'edgecolor','none');
a=fill(m4000nw(:,2),m4000nw(:,1),cg5);set(a,'edgecolor','none');
%a=fill(isobath4000m1b(:,4),isobath4000m1b(:,3),cg5);set(a,'edgecolor','none');
%a=fill(isobath4000m1a(:,4),isobath4000m1a(:,3),cg5);set(a,'edgecolor','none');
%a=fill(isobath4000m2(:,4),isobath4000m2(:,3),cg5);set(a,'edgecolor','none');
%clear isobath4000m1a isobath4000m1b isobath4000m2
a=fill(isobath2000m1(:,4),isobath2000m1(:,3),cg4);set(a,'edgecolor','none');
a=fill(isobath2000m2(:,4),isobath2000m2(:,3),cg4);set(a,'edgecolor','none');
a=fill(isobath2000m3(:,4),isobath2000m3(:,3),cg4);set(a,'edgecolor','none');
a=fill(isobath2000m4(:,4),isobath2000m4(:,3),cg4);set(a,'edgecolor','none');
a=fill(isobath2000m5(:,4),isobath2000m5(:,3),cg4);set(a,'edgecolor','none');
clear isobath2000m1 isobath2000m2 isobath2000m3 isobath2000m4 isobath2000m5
a=fill(isobath200m1(:,4),isobath200m1(:,3),cg3);set(a,'edgecolor','none');
a=fill(isobath200m2(:,4),isobath200m2(:,3),cg3);set(a,'edgecolor','none');
a=fill(isobath200m3(:,4),isobath200m3(:,3),cg3);set(a,'edgecolor','none');
a=fill(isobath200m4(:,4),isobath200m4(:,3),cg3);set(a,'edgecolor','none');
a=fill(isobath200m5(:,4),isobath200m5(:,3),cg3);set(a,'edgecolor','none');
a=fill(isobath200m6(:,4),isobath200m6(:,3),cg3);set(a,'edgecolor','none');
a=fill(isobath200m7(:,4),isobath200m7(:,3),cg3);set(a,'edgecolor','none');
a=fill(isobath200m8(:,4),isobath200m8(:,3),cg3);set(a,'edgecolor','none');
a=fill(isobath200m9(:,4),isobath200m9(:,3),cg3);set(a,'edgecolor','none');
a=fill(isobath200m10(:,4),isobath200m10(:,3),cg3);set(a,'edgecolor','none');
a=fill(isobath200m11(:,4),isobath200m11(:,3),cg3);set(a,'edgecolor','none');
a=fill(isobath200m12(:,4),isobath200m12(:,3),cg3);set(a,'edgecolor','none');
a=fill(isobath200m13(:,4),isobath200m13(:,3),cg3);set(a,'edgecolor','none');
a=fill(isobath200m14(:,4),isobath200m14(:,3),cg3);set(a,'edgecolor','none');
a=fill(isobath200m15(:,4),isobath200m15(:,3),cg3);set(a,'edgecolor','none');
a=fill(isobath200m16(:,4),isobath200m16(:,3),cg3);set(a,'edgecolor','none');
clear isobath200m1 isobath200m2 isobath200m3 isobath200m4 isobath200m5
clear isobath200m6 isobath200m7 isobath200m8 isobath200m9 isobath200m10
clear isobath200m11 isobath200m12 isobath200m13 isobath200m14 isobath200m15
clear isobath200m16
a=fill(isobath100m1(:,4),isobath100m1(:,3),cg2);set(a,'edgecolor','none');
a=fill(isobath100m2(:,4),isobath100m2(:,3),cg2);set(a,'edgecolor','none');
a=fill(isobath100m3(:,4),isobath100m3(:,3),cg2);set(a,'edgecolor','none');
a=fill(isobath100m4(:,4),isobath100m4(:,3),cg2);set(a,'edgecolor','none');
a=fill(isobath100m5(:,4),isobath100m5(:,3),cg2);set(a,'edgecolor','none');
a=fill(isobath100m6(:,4),isobath100m6(:,3),cg2);set(a,'edgecolor','none');
a=fill(isobath100m7(:,4),isobath100m7(:,3),cg2);set(a,'edgecolor','none');
a=fill(isobath100m8(:,4),isobath100m8(:,3),cg2);set(a,'edgecolor','none');
a=fill(isobath100m9(:,4),isobath100m9(:,3),cg2);set(a,'edgecolor','none');
a=fill(isobath100m10(:,4),isobath100m10(:,3),cg2);set(a,'edgecolor','none');
a=fill(isobath100m11(:,4),isobath100m11(:,3),cg2);set(a,'edgecolor','none');
a=fill(isobath100m12(:,4),isobath100m12(:,3),cg2);set(a,'edgecolor','none');
a=fill(isobath100m13(:,4),isobath100m13(:,3),cg2);set(a,'edgecolor','none');
a=fill(isobath100m14(:,4),isobath100m14(:,3),cg2);set(a,'edgecolor','none');
clear isobath100m1 isobath100m2 isobath100m3 isobath100m4 isobath100m5
clear isobath100m6 isobath100m7 isobath100m8 isobath100m9 isobath100m10
clear isobath100m11 isobath100m12 isobath100m13 isobath100m14 

a=fill(isobath60m1(:,4),isobath60m1(:,3),cg1);set(a,'edgecolor','none');
a=fill(isobath60m2(:,4),isobath60m2(:,3),cg1);set(a,'edgecolor','none');
a=fill(isobath60m3(:,4),isobath60m3(:,3),cg1);set(a,'edgecolor','none');
a=fill(isobath60m4(:,4),isobath60m4(:,3),cg1);set(a,'edgecolor','none');
a=fill(isobath60m5(:,4),isobath60m5(:,3),cg1);set(a,'edgecolor','none');
clear isobath60m1 isobath60m2 isobath60m3 isobath60m4 isobath60m5
a=fill(isobath40m1(:,4),isobath40m1(:,3),cg0);set(a,'edgecolor','none');
a=fill(isobath40m2(:,4),isobath40m2(:,3),cg0);set(a,'edgecolor','none');
a=fill(isobath40m3(:,4),isobath40m3(:,3),cg0);set(a,'edgecolor','none');
a=fill(isobath40m4(:,4),isobath40m4(:,3),cg0);set(a,'edgecolor','none');
a=fill(isobath40m5(:,4),isobath40m5(:,3),cg0);set(a,'edgecolor','none');
a=fill(isobath40m6(:,4),isobath40m6(:,3),cg0);set(a,'edgecolor','none');
a=fill(isobath40m7(:,4),isobath40m7(:,3),cg0);set(a,'edgecolor','none');
a=fill(isobath40m8(:,4),isobath40m8(:,3),cg0);set(a,'edgecolor','none');
a=fill(isobath40m9(:,4),isobath40m9(:,3),cg0);set(a,'edgecolor','none');
a=fill(isobath40m10(:,4),isobath40m10(:,3),cg0);set(a,'edgecolor','none');
clear isobath40m1 isobath40m2 isobath40m3 isobath40m4 isobath40m5 isobath40m6
clear isobath40m7
clear isobath40m8 isobath40m9 isobath40m10

a=fill(isobath20m1(:,4),isobath20m1(:,3),cgw);set(a,'edgecolor','none');
a=fill(isobath20m2(:,4),isobath20m2(:,3),cgw);set(a,'edgecolor','none');
a=fill(isobath20m3(:,4),isobath20m3(:,3),cgw);set(a,'edgecolor','none');
a=fill(isobath20m4(:,4),isobath20m4(:,3),cgw);set(a,'edgecolor','none');
a=fill(isobath20m5(:,4),isobath20m5(:,3),cgw);set(a,'edgecolor','none');
a=fill(isobath20m6(:,4),isobath20m6(:,3),cgw);set(a,'edgecolor','none');
a=fill(isobath20m7(:,4),isobath20m7(:,3),cgw);set(a,'edgecolor','none');
clear isobath20m1 isobath20m2 isobath20m3 isobath20m4 isobath20m5 isobath20m6
clear isobath20m7


%a=fill(newcoastline1(:,4),newcoastline1(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline2(:,4),newcoastline2(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline3(:,4),newcoastline3(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline4(:,4),newcoastline4(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline5(:,4),newcoastline5(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline6(:,4),newcoastline6(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline7(:,4),newcoastline7(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline8(:,4),newcoastline8(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline9(:,4),newcoastline9(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline10(:,4),newcoastline10(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline11(:,4),newcoastline11(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline12(:,4),newcoastline12(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline13(:,4),newcoastline13(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline14(:,4),newcoastline14(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline15(:,4),newcoastline15(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline16(:,4),newcoastline16(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline17(:,4),newcoastline17(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline18(:,4),newcoastline18(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline19(:,4),newcoastline19(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline20(:,4),newcoastline20(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline21(:,4),newcoastline21(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline22(:,4),newcoastline22(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline23(:,4),newcoastline23(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newcoastline24(:,4),newcoastline24(:,3),cgb2);set(a,'edgecolor','none');
a=fill(newfound(:,4),newfound(:,3),cgb2);set(a,'edgecolor','none');

%a=plot(newcoastline1(:,4),newcoastline1(:,3));set(a,'color',cgb);
a=plot(newcoastline2(:,4),newcoastline2(:,3));set(a,'color',cgb);
%a=plot(newcoastline3(:,4),newcoastline3(:,3));set(a,'color',cgb);
a=plot(newcoastline4(:,4),newcoastline4(:,3));set(a,'color',cgb);
a=plot(newcoastline5(:,4),newcoastline5(:,3));set(a,'color',cgb);
a=plot(newcoastline6(:,4),newcoastline6(:,3));set(a,'color',cgb);
a=plot(newcoastline7(:,4),newcoastline7(:,3));set(a,'color',cgb);
a=plot(newcoastline8(:,4),newcoastline8(:,3));set(a,'color',cgb);
a=plot(newcoastline9(:,4),newcoastline9(:,3));set(a,'color',cgb);
a=plot(newcoastline10(:,4),newcoastline10(:,3));set(a,'color',cgb);
a=plot(newcoastline11(:,4),newcoastline11(:,3));set(a,'color',cgb);
a=plot(newcoastline12(:,4),newcoastline12(:,3));set(a,'color',cgb);
a=plot(newcoastline13(:,4),newcoastline13(:,3));set(a,'color',cgb);
a=plot(newcoastline14(:,4),newcoastline14(:,3));set(a,'color',cgb);
a=plot(newcoastline15(:,4),newcoastline15(:,3));set(a,'color',cgb);
a=plot(newcoastline16(:,4),newcoastline16(:,3));set(a,'color',cgb);
a=plot(newcoastline17(:,4),newcoastline17(:,3));set(a,'color',cgb);
a=plot(newcoastline18(:,4),newcoastline18(:,3));set(a,'color',cgb);
a=plot(newcoastline19(:,4),newcoastline19(:,3));set(a,'color',cgb);
a=plot(newcoastline20(:,4),newcoastline20(:,3));set(a,'color',cgb);
a=plot(newcoastline21(:,4),newcoastline21(:,3));set(a,'color',cgb);
a=plot(newcoastline22(:,4),newcoastline22(:,3));set(a,'color',cgb);
a=plot(newcoastline23(:,4),newcoastline23(:,3));set(a,'color',cgb);
a=plot(newcoastline24(:,4),newcoastline24(:,3));set(a,'color',cgb);
a=plot(newfound(:,4),newfound(:,3));set(a,'color',cgb);

clear newcoastline1 newcoastline2 newcoastline3 newcoastline4 newcoastline5
clear newcoastline6 newcoastline7 newcoastline8 newcoastline9 newcoastline10
clear newcoastline11 newcoastline12 newcoastline13 newcoastline14 newcoastline15
clear newcoastline16 newcoastline17 newcoastline18 newcoastline19 newcoastline20
clear newcoastline21 newcoastline22 newcoastline23 

% Improve axes labels
set(gca,'TickDir','out');
%yt=[40 41 42 43 ];
%set(gca,'Ytick',yt);
%xt=[-71 -70 -69 -68 -67 -66 -65];
%set(gca,'Xtick',xt);
set(gca,'fontname','times');
set(gca,'fontsize',22);
set(gca,'fontw','bold');

% Add title
%t1=text(-71.4,39.5,'U.S. GLOBEC Georges Bank 1995 Surface Drifters');
%set(t1,'fontsize',25);
%set(t1,'fontw','bold');
%set(t1,'fontname','times');

% Add labels

%t1=text(-65.3,40.45,'4000');
t1=text(-65.3,40.45,'');
set(t1,'Rotation',35)
set(t1,'fontname','times')
set(t1,'fontsize',10)
set(t1,'fontw','bold')
set(t1,'col',cgb)

%add thick border line
box=[-71 43; -71 40;-65 40;-65 43;-71 43];
%h=plot(box(:,1),box(:,2));%set(h,'linewidth',3)
clear box h
 set(gca,'fontsize',18)
clear t1  cg0 cgb cg1 cg2 cg3 cg4 cg5 cg6
clear a fourcorners xt yt ans cg1a cg1aa
clear cg1b cgb2 cgw
end

