%% Placera .m-filen tillsammans med data-filerna från "Datastudio". Se till att inga viktiga filer med filnamnet "korrdata1.txt, korrdata2.txt, korrdata3.txt" befinner sig i samma mapp.

%Input : [t, W, TH, TC] = varmepump('TH.txt', 'TC.txt', 'Geigerräknare.txt')


function [t, W, TH, TC] = varmepump(varargin)
folder=fileparts(which('varmepump.m')); % Se till att MATLAB arbetar i mappen med .m-filen och de exporterade .txt-filerna.
cd(folder)
for i =1:3
% Anpassa värden till MATLABs decimalform.
temps=strcat('korrdata',int2str(i),'.txt');
fin = fopen(varargin{i},'rt');
fout = fopen(temps,'w');
while ~feof(fin)
   s = fgetl(fin);
   s = strrep(s, ',', '.');
   fprintf(fout,'%s\n',s);
   %disp(s)
end
fin  = fclose(fin);
fout = fclose(fout);
end


% Importera värden från korrigerade .txt-filer. 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
data1 = dlmread('korrdata1.txt', '\t', 2,0);
    t = data1(:,1); % Behöver bara extrahera samplingstiderna från en .txt-fil.
    TempA = data1(:,2);
data2 = dlmread('korrdata2.txt', '\t', 2,0);
    TempB = data2(:,2);
data3 = dlmread('korrdata3.txt', '\t', 2,0);
    Wtemp = data3(:,2); 
    W = 4500*vertcat(0,cumsum(Wtemp)); % Den tillförda energin är 4500 Ws för varje puls från geigerräknaren. Inför ett initialt värde vid t=0.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Släng de korrigerade filerna. 
delete('korrdata1.txt', 'korrdata2.txt','korrdata3.txt')
    
    
%Skriv om temperaturerna till Kelvin
TH = TempA + 273.15; 
TC = TempB + 273.15; 

end

