thermalblock_output_function.m

function s=thermalblock_output_function(model_data)

% output computation: Integral over boundary segment


%shortcuts
p=model_data.grid.p;
e=model_data.grid.e;
t=model_data.grid.t;

zero_y_ind=find(p(2,:)==0);             %find indices with y=0

x_values=p(1,:);
x_value=x_values(zero_y_ind);   %corresponding x values with y=0

[x_value_sort,i]=sort(x_value); %sort vector

u_zero_y=u(zero_y_ind);         %u values at y=0
u_zero_y_sort=u_zero_y(i);              %corresponding u values sorted

%the following section tests, if the vector x_value contains 0,1, because
%we want to compute the integral in the intervall [0,1]. 
if min(x_value) ~= 0
        x_value=[0,x_value];
        u_zero=interp1(x_value_sort, u_zero_y_sort,0, 'cubic');
        u_zero_y_sort=[u_zero,u_zero_y_sort];
end

if max(x_value) ~= 1
        x_value=[x_value ,1];
        u_one = interp1(x_value_sort,u_zero_y_sort,1,'cubic');
        u_zero_y_sort=[u_zero_y_sort,u_one];
end

s = trapz(x_value_sort, u_zero_y_sort);

if model.plot_output == 1
        figure;
        plot(x_value_sort, u_zero_y_sort)
end