% This Programme is written for contouring the given data about the
% SNT 250-600 pump blade from the pre-design model output.
% Hamdi Nadir Tural
% 702061013
% ITU Informatics Institute
% Computational Science&Engineering
%

clear,clc;
x9=load('x9.txt');
y9=load('y9.txt');
z9=load('z9.txt');
z9_9=load('z9_9.txt');

figure %1. color contours
plot(x9(:,1),y9(:,1),'-k',x9(:,9),y9(:,9),'-k')
hold on
[c9 h9]=contour(x9,y9,z9,z9(:,1)); 
clabel(c9,h9,'manual');
hold off
axis equal
axis([-10 140 0 220])
colormap(jet(14))
xlabel('Z [mm]')
ylabel('R [mm]')
title('Iso-Power Contours on Blade Meridional Section')

figure %2. gray contours
plot(x9(:,1),y9(:,1),'-k',x9(:,9),y9(:,9),'-k')
hold on
[c9 h9]=contour(x9,y9,z9,z9(:,1)); 
clabel(c9,h9,'manual');
axis equal
axis([-10 140 0 220])
colormap gray(14)
colorbar
hold off
xlabel('Z [mm]')
ylabel('R [mm]')
title('Iso-Power Contours on Blade Meridional Section')

figure %3. color filled
h9=contourf(x9,y9,z9,z9(:,1)); 
axis equal
axis([-10 140 0 220])
colormap jet
colorbar
shading flat
xlabel('Z [mm]')
ylabel('R [mm]')
title('Iso-Power Contours on Blade Meridional Section')

figure %4. gray filled
h9=contourf(x9,y9,z9,z9(:,1)); 
hold on
plot(x9(:,9),y9(:,9),'-k',x9(:,1),y9(:,1),'-k')
hold off
axis equal
axis([-10 140 0 220])
colormap(gray(14))
shading flat
colorbar
xlabel('Z [mm]')
ylabel('R [mm]')
title('Iso-Power Contours on Blade Meridional Section')

figure %5. color matrix
h9=pcolor(x9,y9,z9_9);
hold on
plot(x9(:,9),y9(:,9),'-k',x9(:,1),y9(:,1),'-k')
hold off
axis equal
axis([-10 140 0 220])
colormap(jet(14))
colorbar
xlabel('Z [mm]')
ylabel('R [mm]')
title('Iso-Velocity Contours on Blade Meridional Section')

figure %6. gray matrix
h9=pcolor(x9,y9,z9_9);
hold on
plot(x9(:,9),y9(:,9),'-k',x9(:,1),y9(:,1),'-k')
hold off
axis equal
axis([-10 140 0 220])
colormap(gray(14))
colorbar
xlabel('Z [mm]')
ylabel('R [mm]')
title('Iso-Velocity Contours on Blade Meridional Section')

figure %7. color filled with gamma=0.4
h9=contourf(x9,y9,z9,z9(:,1)); 
axis equal
axis([-10 140 0 220])
colormap(jet(14))
colorbar
brighten(0.6) %gamma=1-beta
shading flat
xlabel('Z [mm]')
ylabel('R [mm]')
title('Iso-Power Contours on Blade Meridional Section (\gamma =0.4)')

figure %8. gray filled with gamma=0.5
h9=contourf(x9,y9,z9,z9(:,1)); 
hold on
plot(x9(:,9),y9(:,9),'-k',x9(:,1),y9(:,1),'-k')
hold off
axis equal
axis([-10 140 0 220])
colormap(gray(14))
brighten(0.5)
colorbar
shading flat
xlabel('Z [mm]')
ylabel('R [mm]')
title('Iso-Power Contours on Blade Meridional Section (\gamma =0.5)')