%load NASAs OMG data 2020
load output.mat % [NASA (2021)]
lon = output(:,1); % Longitude
lat = output(:,2); % Latitude
ind_disko = find(lat >= 68 & lat <= 72 & lon <= 90 & lon >= 50);
lon_disko = lon(ind_disko);
lat_disko = lat(ind_disko);
depth_disko = depth(ind_disko);
temp_disko = temp(ind_disko);
sal_disko = sal(ind_disko);
%% MAP Greenland with 2020 OMG stations (red) plus highlighted Disko Bay stations (blue)
clf, axis([-80 -5 55 90])
topoplot(axis,-5:1/8:0),hold on % [Lilly (2021)]
topoplot continents, latratio(55) % [Lilly (2021)]
a = plot(lon.*(-1),lat,'wo','markerfacecolor','r');
b = plot(lon_disko.*(-1),lat_disko,'wo','markerfacecolor','b');
hc=colorbar;hc.Label.String='Bathymetry [km]';
title('NASAs OMG Data 2020 (stations)');
legend([a b],'Stations','Disko Bay','location','southeast')
set(gca,'fontsize',15,'linewidth',2)
clf, axis([-65 -50 66 73])
topoplot(axis,-5:1/8:0) % [Lilly (2021)]
topoplot continents, latratio(55) % [Lilly (2021)]
b = plot(lon.*(-1),lat,'wo','markerfacecolor','r')
b =
Line with properties:
Color: [1 1 1]
LineStyle: 'none'
LineWidth: 0.5000
Marker: 'o'
MarkerSize: 6
MarkerFaceColor: [1 0 0]
XData: [1×1488442 double]
YData: [1×1488442 double]
ZData: [1×0 double]
Show all properties
c = plot(lon_disko.*(-1),lat_disko,'wo','markerfacecolor','b');
hc=colorbar;hc.Label.String='Bathymetry [km]';
title('Disko Bay - NASAs OMG Data 2020 (stations)');
set(gca,'fontsize',20,'linewidth',2);