Maps of CTD stations of NASA's OMG Project

Tianyu Zhou, 10/12/2021, at UDel

1. The Map of all the Stations around Greenland (2016-2021)

Define the latitude and longitude range of the mapped 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 (https://www.eoas.ubc.ca/~rich/map.html)

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

Read bathymetry data around Greenland which are provided by GEBCO (https://www.gebco.net/)

fbath = '../data/GEBCO/gebco_Greenland.nc';
lat_bath = ncread(fbath,'lat');
lon_bath = ncread(fbath,'lon');
elevation = ncread(fbath,'elevation');
lat_bath = lat_bath(1:6:end);
lon_bath = lon_bath(1:6:end);
elevation = elevation(1:6:end,1:6:end);
elevation(elevation>=0) = nan;
[lat_bath,lon_bath] = meshgrid(lat_bath,lon_bath);
m_pcolor(lon_bath,lat_bath,elevation)
shading interp
m_contour(lon_bath,lat_bath,elevation,[-3000 -3000],'k','linewidth',1)
m_contour(lon_bath,lat_bath,elevation,[-1000 -1000],'k','linewidth',1)
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]);

Read latitude and longitude (they are stored in a '.mat' file) of all the stations and plot them on the map

for year = 2016:2021
    load(['../data/lonlat_',num2str(year),'.mat'])    
    m_plot(-lon,lat,'o','markersize',3,'color',[0.8500 0.3250 0.0980],'markerfacecolor',[0.8500 0.3250 0.0980])
end

Colorbar settings

color = cmocean('deep');

color = flip(color,1);

colormap(color)

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

2. The Map of Stations in Nares Strait (2016-2021)

Define the latitude and longitude range of Nares Strait

lat_start = 77.5;
lat_end = 83.3;
lon_start = -85;
lon_end = -47;

Load data located within Nares Strait (pre-selected and stored in a '.mat' file)

load('../data/OMG_data_Nares')
profi_num = size(output,1);

Create a map of Nares Strait

f = figure;
set(f,'units','centimeters','position',[1,1,12*1.5,9*1.5],'color','w')
m_proj('mercator','lon',[lon_start lon_end],'lat',[lat_start lat_end]);
hold on
m_pcolor(lon_bath,lat_bath,elevation)
shading interp
m_contour(lon_bath,lat_bath,elevation,[-300 -300],'k','linewidth',1)
m_grid('box','on','linest','none','fontname','times','linewidth',0.1,'tickdir','out','backcolor',[1 1 1],'Fontsize',20);
m_gshhs_h('patch',[0.6 0.6 0.6]);
for n = 1:profi_num
    m_plot(-output{n,3},output{n,4},'.','color','r','markersize',20)  
end

Colorbar settings

color = cmocean('deep');

color = flip(color,1);

colormap(color)

colorbar

set(gca,'fontsize',16,'color','w');