% jarkovskeho_jev_cb.m
% Jarkovskeho jev 1-D: spolecna nastaveni grafu.
% Miroslav Broz (miroslav.broz@email.cz), Jan 28th 2006

set(h1, 'EdgeColor', 'none');
cameratoolbar;
axis tight;

xlabel('depth {\itx} / m');
zlabel('temperature {\itT} / K');

%F_eq = 4*sigma*T_eq^4/(1-A);
%F = F_eq + F_1 * cos(2*pi*f*t);	% this is a real flux value
F = T_eq + T_1_0_ * cos(2*pi*f* (t*P)); % but plot the same amplitude as T
hold on
plot3(0.*t, t, F, 'k-');
h3 = text(x(1), t(length(t)), T_eq+T_1_0_, '  \propto{\itF}({\itt})', 'HorizontalAlignment', 'left');

%delta = sqrt(chi/(pi*f));
tmp=3*delta;
plot3([tmp tmp],[t(1) t(length(t))],[T_eq T_eq],'k-');
text(tmp, t(length(t)), T_eq, '  3\delta');

ha = gca;
cb = colorbar;
axes(cb);
set(cb,'Position',[0.875 0.11 0.03 0.815]);
ylabel('temperature {\it T} / K');
axes(ha);

