% jarkovskeho_jev_T_x_t_regolit_rocni.m
% Jarkovskeho jev 1-D: graf teplota(hloubka, cas) v regolitu pri rocnim jevu.
% Miroslav Broz (miroslav.broz@email.cz), Jan 28th 2006

global T_eq T_1_0_ f phi_th_ chi;

jarkovskeho_jev_params;

K = K_regolit;
rho = rho_regolit;
chi = K/(rho*C);

P = 4*365.25*24*3600.;	% seasonal
f = 1./P;
T_eq = 150;
A = 0.1;
delta = sqrt(chi/(pi*f));
phi_th_ = phi_th(f,T_eq,K);
T_1_0_ = T_1_0(f,T_eq,K);

x = 0:1e-2:1.5;
x = 0:1e-2:3*delta;
t = 0:1e-2:2;		% in periods

[X,Y] = meshgrid(x,t.*P);
Z = T_x_t(X,Y);
h1 = surf(X,Y./P,Z);

ylabel([ 'time {\itt} / periods {\itP} = ' sprintf('%.1f', P/3600/24/365.25) ' yr' ]);

jarkovskeho_jev_cb;

view([-30 30]);
print('-r192', '-dpng', 'jarkovskeho_jev_T_x_t_regolit_rocni.png');

