clear all; close all; clc


load('OSCAHR_TSG_corr.mat')
load('OSCAHR_STATIONS.mat')
load('OSCAHR_LAT_LON_MVP.mat')

% All
North = 44.5;
South = 42.5;
East = 9.5;
West = 6;

% %Leg1
North1 = 44;
South1 = 43;
East1 = 8;
West1 = 6;

%Leg2
North2 = 44.25;
South2 = 43.25;
East2 = 9.25;
West2 = 7.75;

figure(1)
set(gca,'FontSize',15,'FontName','Arial');hold on %Set font type and size.
m_proj('Equidistant Cylindrical','lon',[West  East],'lat',[South North]);%set up our projection and bourndaries.
hold on;
m_gshhs_h('patch',[.5 .5 .5],'edgecolor','k');%Set up our coastlines and continents.
m_grid('box','fancy','tickdir','in','fontsize',20);%Set up our fancy box.
m_scatter(LON,LAT,50,TRANSECT,'fill');
cbar1 = colorbar('location','EastOutside','FontSize',15);%position and fontsize of our colour bar.
colormap (hsv(19))
caxis([1 20])
m_plot(LON_STA,LAT_STA,'s','MarkerSize',12,'MarkerEdgeColor','k','MarkerFaceColor','k');
m_plot(LON_MVP,LAT_MVP,'+k','MarkerSize',5);
hold off;
print(gcf,'-dpng','OSCAHR_map_station_transect');

figure(2)
set(gca,'FontSize',15,'FontName','Arial');hold on %Set font type and size.
m_proj('Equidistant Cylindrical','lon',[West1  East1],'lat',[South1 North1]);%set up our projection and bourndaries.
hold on;
m_gshhs_h('patch',[.5 .5 .5],'edgecolor','k');%Set up our coastlines and continents.
m_grid('box','fancy','tickdir','in','fontsize',20);%Set up our fancy box.
m_scatter(LON,LAT,50,TRANSECT,'fill');
cbar2 = colorbar('location','EastOutside','FontSize',15);%position and fontsize of our colour bar.
colormap (hsv(19))
caxis([1 20])
m_plot(LON_STA,LAT_STA,'s','MarkerSize',12,'MarkerEdgeColor','k','MarkerFaceColor','k');
m_plot(LON_MVP,LAT_MVP,'+k','MarkerSize',5);
hold off;
print(gcf,'-dpng','OSCAHR_map_station_transect_LEG1');

figure(3)
set(gca,'FontSize',15,'FontName','Arial');hold on %Set font type and size.
m_proj('Equidistant Cylindrical','lon',[West2  East2],'lat',[South2 North2]);%set up our projection and bourndaries.
hold on;
m_gshhs_h('patch',[.5 .5 .5],'edgecolor','k');%Set up our coastlines and continents.
m_grid('box','fancy','tickdir','in','fontsize',20);%Set up our fancy box.
m_scatter(LON,LAT,50,TRANSECT,'fill');
cbar3 = colorbar('location','EastOutside','FontSize',15);%position and fontsize of our colour bar.
colormap (hsv(19))
caxis([1 20])
m_plot(LON_STA,LAT_STA,'s','MarkerSize',12,'MarkerEdgeColor','k','MarkerFaceColor','k');
m_plot(LON_MVP,LAT_MVP,'+k','MarkerSize',5);
hold off;
print(gcf,'-dpng','OSCAHR_map_station_transect_LEG2');
