% jarkovskeho_jev_params.m
% Jarkovskeho jev 1-D: spolecne tepelne parametry.
% Miroslav Broz (miroslav.broz@email.cz), Jan 28th 2006

global sigma C rho eps F1 A chi;

sigma = 5.67e-8;        % W/m^2/K^4; Stefanova-Boltzmannova konstanta

K_bazalt = 1.;		% W/m/K; tepelna vodivost
K_regolit = 0.01;
rho_bazalt = 3500.;	% kg/m^3; hustota
rho_regolit = 1500.;
C = 680.;		% J/kg/K; merna tepelna kapacita
eps = 0.9;		% IR emisivita
A = 0.1;		% Bondovo albedo

K = K_bazalt;
rho = rho_bazalt;
chi = K/(rho*C);

F1 = 1360./(2.5^2)/2;	% W/m^2; amplituda toku slunecniho zareni
A = 0.1;		% Bondovo albedo

Pmin = 1e-2;
Pmax = 365.25*24;
fmin = 1./(Pmax.*3600.);
fmax = 1./(Pmin.*3600.);
logf = log10(fmin):(log10(fmax)-log10(fmin))/100:log10(fmax);
f = 10.^logf;
P = (1./f)./3600.;
T_eq = 100:2:200;

[X,Y] = meshgrid(logf,T_eq);

%set(gcf,'Renderer','opengl');
%colormap autumn(256);
%colormap jet(256);
colormap hot(256);
hold off;


