%% r        - polomer vaku v [cm]
%% h        - vyšška vaku v [cm]
%% V_ml     - objem krvi v [ml] vo vaku
%% Ro_krv   - hustota krvi [kg/m3], od 1048-1066

%% Typ plastu:
%%              1) Polyvinylchlorid (PVC)
%%              2) Polypropylén (PP)
%%              3) Polyetylén (PE)
%%              4) Etylenvinyl acetát (EVA)

%% T_krv_p - pociatocna teplota krvi v [°C]
%% T_krv_k - teplota na ktorú chcete krv ochladit v [°C]
%% T_chlad - teplota chladnicky v [°C]


function krv()
  prompt = {'T_krv_p', 'T_krv_k', 'T_chlad', 'r [cm]', 'h [cm]', 'V_ml [ml]', 'Ro_krv 1048-1066', 'Typ plastu'};
  defaults = {'36', '22.0', '20.0', '10', '18', '450', '1060', '1'};
  rowscols = [1,10; 1,10; 1,10; 1,10; 1,10; 1,10; 1,10; 1,10;];
  dims = inputdlg (prompt, 'Enter experiment params', rowscols, defaults);

  T_krv_p = str2num(dims{1});
  T_krv_k = str2num(dims{2});
  T_chlad = str2num(dims{3});
  r = str2num(dims{4});
  h = str2num(dims{5});
  V_ml = str2num(dims{6});
  Ro_krv = str2num(dims{7});
  typ_plast = str2num(dims{8});

  chladenie(T_krv_p, T_krv_k, T_chlad, r, h, V_ml, Ro_krv, typ_plast);
end

function chladenie( T_krv_p, T_krv_k, T_chlad, r, h, V_ml, Ro_krv, typ_plast)
   clc;
   lam_PVC=0.20;
   lam_PP=0.22;
   lam_PE=0.35;
   lam_EVA=0.34;
   c_krv = 3617; %% sú?inite? tepelnej vodivosti krvi
   x=0.00025; %% hrúbka plastu

   %% volanie funkcie
   lam_plast = get_lam_plastu(typ_plast);

   %% premeny jednotiek na zakladne
   r_vak = r*0.01;
   h_vak = h*0.01;
   T_krv_p_k=T_krv_p+273.15;
   T_krv_k_k=T_krv_k+273.15;
   T_chlad_k=T_chlad+273.15;
   V=V_ml*0.000001;

   C1=Ro_krv*V*c_krv; %% hmotnostna tepelna kapacita krvi
   S=4*((r_vak*h_vak)+(r_vak*x)+(h_vak*x)); %% plocha vaku
   T_roz1=(T_krv_p_k-T_chlad_k);
   %T_roz2=(T_krv_p_k-T_krv_k_k)

   Q1=C1*(T_krv_p_k-T_chlad_k); %% teplo, ktore moze krv odovzdat

   D1=(Q1*x)/(lam_plast*S);

   D_PVC=(Q1*x)/(lam_PVC*S);

   D_PP=(Q1*x)/(lam_PP*S);

   D_PE=(Q1*x)/(lam_PE*S);

   D_EVA=(Q1*x)/(lam_EVA*S);

   X1 = [];
   Y1 = [];
   X2 = [];
   Y2 = [];
   X3 = [];
   Y3 = [];
   X4 = [];
   Y4 = [];

   X5 = [];
   Y5 = [];
   X6 = [];
   Y6 = [];

   T_zmena1=T_krv_p_k;
   T_zmena2=T_krv_p_k;
       for i =1:T_roz1
            t1=D_PVC*(1/(T_zmena1-T_chlad_k)); %% vypocet casu
            t1=t1-(D_PVC*(1/(T_roz1)));

            t2=D_PP*(1/(T_zmena1-T_chlad_k)); %% vypocet casu
            t2=t2-(D_PP*(1/(T_roz1)));

            t3=D_PE*(1/(T_zmena1-T_chlad_k)); %% vypocet casu
            t3=t3-(D_PE*(1/(T_roz1)));

            t4=D_EVA*(1/(T_zmena1-T_chlad_k)); %% vypocet casu
            t4=t4-(D_EVA*(1/(T_roz1)));

            T_vys1=T_zmena1-273.15;
            T_zmena1=T_zmena1-1;

            Y1=[Y1;T_vys1]; %%

            %X1=[X1;datenum(0,0,0,0,0,t1)]; %% vyskresluje  vsetko do jedneho grafu
            X1=[X1;t1]; %% vyskresluje  vsetko do jedneho grafu

            %X2=[X2;datenum(0,0,0,0,0,t2)];%%
            X2=[X2;t2];%%

            %X3=[X3;datenum(0,0,0,0,0,t3)];%%
            X3=[X3;t3];%%

            %X4=[X4;datenum(0,0,0,0,0,t4)];%%
            X4=[X4;t4];%%
       end

       for i =1:T_roz1

            t2=D1*(1/(T_zmena2-T_chlad_k)); %% vypocet casu
            t2=t2-(D1*(1/(T_roz1)));
            Q_krv=(C1*T_zmena2)/1000000;
            T_vys2=T_zmena2-273.15;

            T_zmena2=T_zmena2-1;

            %X5=[X5;datenum(0,0,0,0,0,t2)];
            X5=[X5;t2];

            Y5=[Y5;T_vys2];

            %X6=[X6;datenum(0,0,0,0,0,t2)];
            X6=[X6;t2];

            Y6=[Y6;Q_krv];

       end


   figure;
   subplot(2,2,1);
   plot(X1, Y1, X2, Y1, X3, Y1, X4, Y1);
   grid on;
   title('Vplyv materialu na rychlost  chladenia');
   xlabel('t [s]');
   ylabel('T [°C]');
   %datetick('x','HH:MM:SS');
   xlim([0 inf]);
   legend('PVC', 'PP', 'PE', 'EVA');


   subplot(2,2,[3,4]);
   plot(X5, Y5);
   grid on;
   title('Teplota za cas');
   xlabel('t [s]');
   ylabel('T [°C]');
   %%datetick('x','HH:MM:SS');
   xlim([0 inf]);

   subplot(2,2,2);
   plot(X6, Y6);
   grid on;
   title('Teplo za cas');
   xlabel('t [s]');
   ylabel('Q [MJ]');
   %%datetick('x','HH:MM:SS');
   xlim([0 inf]);


end

%%volana funkcia, sucinitel prestupu tepla
function lam_plastu = get_lam_plastu(typ_plast)
   switch logical(true)
       case (typ_plast==1 ),  lam_plastu = 0.20;
       case (typ_plast==2 ),  lam_plastu = 0.22;
       case (typ_plast==3 ),  lam_plastu = 0.355;
       case (typ_plast==4 ),  lam_plastu = 0.34;
       otherwise,  lam_plastu = 0.20;
   end;
end