clear all
close all
% LEER DATOS
resultados= importdata('solucion.txt');
% VARIABLES ALMACENADAS EN VECTORES
NN=resultados.data(:,1);
x=resultados.data(:,2);
y=resultados.data(:,3);
z=resultados.data(:,4);
presion=resultados.data(:,5);
densidad=resultados.data(:,6);
velocidad=resultados.data(:,7);
x_vel=resultados.data(:,8);
y_vel=resultados.data(:,9);
z_vel=resultados.data(:,10);
temperatura=resultados.data(:,11);
y_star=resultados.data(:,12);
y_plus=resultados.data(:,13);
flujo_calor=resultados.data(:,14);
cell_volume=resultados.data(:,15);
% SISTEMA DE COORDENADAS CILÍNDRICO
[theta,rho,z]=cart2pol(x,y,z);
vector_operado=abs(unique(z)-5);
z_unicos=unique(z);
SC=z_unicos(find(min(vector_operado)==vector_operado));
a_z=find(z==SC);
x_corte=x(a_z);
y_corte=y(a_z);
z_corte=z(a_z);
presion_corte=presion(a_z);
densidad_corte=densidad(a_z);
velocidad_corte=velocidad(a_z);
x_vel_corte=x_vel(a_z);
y_vel_corte=y_vel(a_z);
z_vel_corte=z_vel(a_z);
temperatura_corte=temperatura(a_z);
y_star_corte=y_star(a_z);
y_plus_corte=y_plus(a_z);
flujo_calor_corte=flujo_calor(a_z);
cell_volume_corte=cell_volume(a_z);
theta_corte=theta(a_z);
rho_corte=rho(a_z);
% DIFERENCIAR ZONAS --> hay tubo, aire y sal
% DENSIDAD:
% sal--> densidad>1000
% tubo--> densidad>0.01 (es < que 0.01 porque aunque la densidad
% del haynes230 es 8970, la exporta como 0!!)
% aire--> 0.01<densidad<1000
% SAL:
x_sal=x_corte(find(densidad_corte>1000));
y_sal=y_corte(find(densidad_corte>1000));
z_sal=z_corte(find(densidad_corte>1000));
presion_sal=presion_corte(find(densidad_corte>1000));
densidad_sal=densidad_corte(find(densidad_corte>1000));
velocidad_sal=velocidad_corte(find(densidad_corte>1000));
x_vel_sal=x_vel_corte(find(densidad_corte>1000));
y_vel_sal=y_vel_corte(find(densidad_corte>1000));
z_vel_sal=z_vel_corte(find(densidad_corte>1000));
temperatura_sal=temperatura_corte(find(densidad_corte>1000));
y_star_sal=y_star_corte(find(densidad_corte>1000));
y_plus_sal=y_plus_corte(find(densidad_corte>1000));
flujo_calor_sal=flujo_calor_corte(find(densidad_corte>1000));
cell_volume_sal=cell_volume_corte(find(densidad_corte>1000));
theta_sal=theta_corte(find(densidad_corte>1000));
rho_sal=rho_corte(find(densidad_corte>1000));
% TUBO:
x_tubo=x_corte(find(densidad_corte<0.01));
y_tubo=y_corte(find(densidad_corte<0.01));
z_tubo=z_corte(find(densidad_corte<0.01));
presion_tubo=presion_corte(find(densidad_corte<0.01));
densidad_tubo=densidad_corte(find(densidad_corte<0.01));
velocidad_tubo=velocidad_corte(find(densidad_corte<0.01));
x_vel_tubo=x_vel_corte(find(densidad_corte<0.01));
y_vel_tubo=y_vel_corte(find(densidad_corte<0.01));
z_vel_tubo=z_vel_corte(find(densidad_corte<0.01));
temperatura_tubo=temperatura_corte(find(densidad_corte<0.01));
y_star_tubo=y_star_corte(find(densidad_corte<0.01));
y_plus_tubo=y_plus_corte(find(densidad_corte<0.01));
flujo_calor_tubo=flujo_calor_corte(find(densidad_corte<0.01));
cell_volume_tubo=cell_volume_corte(find(densidad_corte<0.01));
theta_tubo=theta_corte(find(densidad_corte<0.01));
rho_tubo=rho_corte(find(densidad_corte<0.01));
%AIRE:
x_aire=x_corte(find(densidad_corte>0.01 & densidad_corte<1000));
y_aire=y_corte(find(densidad_corte>0.01 & densidad_corte<1000));
z_aire=z_corte(find(densidad_corte>0.01 & densidad_corte<1000));
presion_aire=presion_corte(find(densidad_corte>0.01 & densidad_corte<1000));
densidad_aire=densidad_corte(find(densidad_corte>0.01 & densidad_corte<1000));
velocidad_aire=velocidad_corte(find(densidad_corte>0.01 & densidad_corte<1000));
x_vel_aire=x_vel_corte(find(densidad_corte>0.01 & densidad_corte<1000));
y_vel_aire=y_vel_corte(find(densidad_corte>0.01 & densidad_corte<1000));
z_vel_aire=z_vel_corte(find(densidad_corte>0.01 & densidad_corte<1000));
temperatura_aire=temperatura_corte(find(densidad_corte>0.01 & densidad_corte<1000));
y_star_aire=y_star_corte(find(densidad_corte>0.01 & densidad_corte<1000));
y_plus_aire=y_plus_corte(find(densidad_corte>0.01 & densidad_corte<1000));
flujo_calor_aire=flujo_calor_corte(find(densidad_corte>0.01 & densidad_corte<1000));
cell_volume_aire=cell_volume_corte(find(densidad_corte>0.01 & densidad_corte<1000));
theta_aire=theta_corte(find(densidad_corte>0.01 & densidad_corte<1000));
rho_aire=rho_corte(find(densidad_corte>0.01 & densidad_corte<1000));
% LIMITAR SUPERFICIES
% interior: rho_tubo=<0.01485
% exterior: rho_tubo=>0.0165
% Obtengo los puntos en los bordes del tubo: TE (tubo exterior)
x_TE=x_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
y_TE=y_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
z_TE=z_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
temperatura_TE=temperatura_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
y_plus_TE=y_plus_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
flujo_calor_TE=flujo_calor_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
cell_volume_TE=cell_volume_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
theta_TE=theta_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
rho_TE=rho_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo>=0.0165));
% Obtengo los puntos en los bordes del tubo: TI (tubo interior)
x_TI=x_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
y_TI=y_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
z_TI=z_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
temperatura_TI=temperatura_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
y_plus_TI=y_plus_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
flujo_calor_TI=flujo_calor_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
cell_volume_TI=cell_volume_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
theta_TI=theta_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
rho_TI=rho_tubo(find(flujo_calor_tubo~=0 & y_plus_tubo==0 & rho_tubo<=0.01485));
plot(temperatura_TE,theta_TE,'s');
posicion_max_temperatura_TE=find(temperatura_TE==max(max(temperatura_TE)));
posicion_0_theta_TE=find(theta_TE==min(abs(theta_TE)));