Tmax below 100m (Part 2)

clc;

clear;

% Loading the previously saved data

load('WHW_10.dat')

lat = WHW_10(:,1);

lon = WHW_10(:,2);

year = WHW_10(:,3);

Tmax = WHW_10(:,4);

depth_Tmax =WHW_10(:,5);

S_Tmax = WHW_10(:,6);

% median in each year

nyear = length(2016:2021);

Tmax_me = zeros(1,nyear);

for n = 1:nyear

Tmax_me(n) = median(Tmax(year==2015+n));

end

% Temporal variation of Tmax and Tmax_avg

f = figure;

set(f,'units','centimeters','position',[1,1,13*1.5,7*1.5],'color','w')

FontSize = 15;

plot(year,Tmax,'.','MarkerSize',15)

hold on

plot(2016:2021,Tmax_me,'o-r','LineWidth',2,'MarkerFaceColor','r')

xlim([2015.5 2021.5])

ylim([-1.2 1.2])

xlabel('Year','FontWeight','bold')

ylabel('Tmax','FontWeight','bold')

legend({'Stations','Median'},'Location','northeast')

set(gca,'FontSize',FontSize,'tickdir','out')