figure(1), clf, 
subplot(2,2,1)
plot(Time/YEAR,Ta1,'r',Time/MONTH,Ta2,'b-'); hold on
plot(Time/YEAR,To1,'m-',Time/MONTH,To2,'c-'); hold on
plot(Time/YEAR,Ts,'g');
plot(Time/YEAR,TEMI1*ones(1,Nt),'r--',Time/YEAR,TEMI2*ones(1,Nt),'b--');
%plot(Time(1:Nt-1)/YEAR,PSIa/SV,'k');
axis([0 Time(end)/YEAR 225 325]);
subplot(2,2,2)
plot(Time(1:Nt-1)/YEAR,(Fa*pi*RADIUS^2)/PW,'g'); hold on
plot(Time(1:Nt-1)/YEAR,(Fo*pi*RADIUS^2)/PW,'y'); hold on
%axis([0 160 -300 300]);
subplot(2,2,3)
plot(Time(1:Nt-1)/YEAR,PSIa/SV,'k');
subplot(2,2,4)
plot(Time(1:Nt-1)/YEAR,Ft,'g');


figure(2), clf,
plot(Time/MONTH,Ts1-Ts2-DTcrit_circ,'g');

return

figure(3), clf,
for n = 1:200:Nt-1
    subplot(2,3,1)
    [h]=plot(Ta1(n),To1(n),'r.'); hold on;set(h,'markersize',10); 
    [h]=plot(Ta2(n),To2(n),'b.'); set(h,'markersize',10);
    %[h]=plot((Ta1(n)+Ta2(n))/2,(To1(n)+To2(n))/2,'g.'); set(h,'markersize',10);
    plot([Ta1_end Ta1_end],[0 400],'r');plot([0 400],[To1_end To1_end],'r');
    plot([Ta2_end Ta2_end],[0 400],'b');plot([0 400],[To2_end To2_end],'b');
    set(gca,'Fontsize',10);
    axis([225 325 225 325]);
    xlabel('TA (K)','Fontsize',15);
    ylabel('TO (K)','Fontsize',15);
    pause(1/5);
    
    subplot(2,3,2)
    [h]=plot(Fa(n)*pi*(RADIUS^2)/PW,Fo(n)*pi*(RADIUS^2)/PW,'m.'); 
    hold on;set(h,'markersize',10); 
    plot([Fa_end Fa_end],[-400 400],'m');plot([-400 400],[Fo_end Fo_end],'m');
    H=5; axis([-H H -H H]);
    xlabel('HA (PW)','Fontsize',15);
    ylabel('HO (PW)','Fontsize',15);
    pause(1/5);
    
    subplot(2,3,3)
    [h]=plot(Feva1(n),Feva2(n),'c.'); 
    hold on;set(h,'markersize',10); 
    H=200; axis([0 H 0 H]);
    xlabel('LvE1 (W m-2)','Fontsize',15);
    ylabel('LvE2 (W m-2)','Fontsize',15);
    pause(1/5);
    
    subplot(2,3,4)
    [h]=plot(Ts1(n),Ts2(n),'g.'); 
    hold on;set(h,'markersize',10);
    plot([200:1:400],[200:1:400]-DTcrit_circ,'g--');
    plot([Ts1_end Ts1_end],[0 400],'g');plot([0 400],[Ts2_end Ts2_end],'g');
    axis([240 320 240 320]);
    xlabel('Ts1 (K)','Fontsize',15);
    ylabel('Ts2 (K)','Fontsize',15);
    pause(1/5);
    
    subplot(2,3,5)
    [h]=plot(Ts(n),Ft(n),'k.'); 
    hold on;set(h,'markersize',10);
    plot([Ts_end Ts_end],[-500 500],'k');plot([250 400],[Ft_end Ft_end],'k');
    axis([250 325 -200 200]);
    xlabel('Ts (K)','Fontsize',15);
    ylabel('Ft (Wm-2)','Fontsize',15);
    pause(1/5);
end