Plotting a Map of Greenland with Locations of All the Data Marked

Define the latitude and longitude of the selected region.

clc;
clear;
lat_start = 56;
lat_end = 86;
lon_start = -90;
lon_end = 0;

Start to creat a map of Greenland (this requires the m_map toolbox).

f = figure;
set(f,'units','centimeters','position',[1,1,12*1.5,9*1.5],'color','w')
m_proj('lambert','lon',[lon_start lon_end],'lat',[lat_start lat_end]);
hold on
m_grid('box','on','linest','none','fontname','times','linewidth',0.1,...
    'tickdir','out','backcolor',[1 1 1],'Fontsize',20);
m_gshhs_l('patch',[0.6 0.6 0.6]);

Set six different color for each of the year from 2016 to 2021.

color = jet(6);

Read latitude and longitude (they are stored in a '.mat' file) of all the profiles and plot them on the map. Colors are assigned according to the year when the data were collected.

for year = 2016:2021
    load(['../data/lonlat_',num2str(year),'.mat'])    
    m_plot(-lon,lat,'.','color',color(year-2015,:),'markersize',10)
end

Colorbar settings

colormap(color)

colorbar('ytick',linspace(0.1,0.9,6),'yTicklabel',2016:2021,'Position',...

[0.79 0.56 0.02 0.32],'fontname','times','fontsize',15)