Param
DTz = DTcrit_conv;
DTh = DTcrit_circ;

Ts1 = [200:.1:350];
Ts2 = Ts1-DTh;
t = (Ts1-Ts1(1))/10;
EPSA1= 1-exp(-t);
t = (Ts2-Ts2(1))/10;
EPSA2= 1-exp(-t);
OLR1 = EPSA1.*(Ts1-DTz).^4+(1-EPSA1).*Ts1.^4; 
OLR2 = EPSA2.*(Ts2-DTz).^4+(1-EPSA2).*Ts2.^4;
OLR = SIGMA*(OLR1+OLR2)/2;
ASR = SIGMA*(TEMI1^4+TEMI2^4)/2;

figure(1), clf,
plot(Ts1,OLR,'b',Ts1,ASR*ones(size(Ts1)),'g');

figure(2), clf,
plot(Ts1,EPSA1,'r',Ts2,EPSA2,'b');



function emiss