% A few plots to illustrate the structure of the data files


load RDM_global_05degree.mat


%% Plot global roughness map

figure
pcolor(lon_ave,lat_ave,roughness)
shading interp
set(gca,'dataaspectratio',[1 1 1])
hold on %overlay coastline
load coast
plot(long,lat,'k')
set(gca,'clim',[0 1000])
colorbar
fontsize(14,14,14,14)
xlabel('Longitude ','fontsize',14)
ylabel('Latitude ','fontsize',14)
title('Topographic Roughness [m] ', 'fontsize', 14)

%% Plot maps of diffusivity at various depths

i=50% choose depth level. depth goes from -7000 m (i=1) to -100 m (i=70)

figure
pcolor(lon_ave,lat_ave,log10(K_RDM(:,:,i)))
shading interp
set(gca,'dataaspectratio',[1 1 1])
hold on %overlay coastline
load coast
plot(long,lat,'k')
colorbar
fontsize(14,14,14,14)
xlabel('Longitude ','fontsize',14)
ylabel('Latitude ','fontsize',14)
title(['Diapycnal Diffusivity at ' num2str(depth(i)) ' m.   [log_{10} (K)]' ], 'fontsize', 14)

%% Plot Meridional sections

figure
lon_section=-150.25; %pick longitude. 
i=find(lon_ave==lon_section);
section=squeeze(K_RDM(:,i,:));
pcolor(lat_ave,depth,log10(section'))
shading interp
daspect([100 20000 1])
colorbar
fontsize(14,14,14,14)
xlabel('Longitude ','fontsize',14)
ylabel('Latitude ','fontsize',14)
title(['Meridional diffusivity section along ' num2str(lon_ave(i)) '   [log _{10} (K)]' ], 'fontsize', 14)


%% Plot Zonal sections

figure
lat=0.25; %pick latitude
i=find(lat_ave==lat);
section=squeeze(K_RDM(i,:,:));

pcolor(lon_ave,depth,log10(section'))
shading interp
daspect([100 5000 1])
colorbar
fontsize(14,14,14,14)
xlabel('Longitude ','fontsize',14)
ylabel('Latitude ','fontsize',14)
title(['Zonal diffusivity section along ' num2str(lat_ave(i)) '   [log _{10} (K)]' ], 'fontsize', 14)