T=linspace(0,45,100); P=101.325; cento=P_ws(T); omega_cento=0.62185.*cento./(P-cento); novanta=0.9.*P_ws(T); omega_novanta=0.62185.*novanta./(P-novanta); ottanta=0.8.*P_ws(T); omega_ottanta=0.62185.*ottanta./(P-ottanta); settanta=0.7.*P_ws(T); omega_settanta=0.62185.*settanta./(P-settanta); sessanta=0.6.*P_ws(T); omega_sessanta=0.62185.*sessanta./(P-sessanta); cinquanta=0.5.*P_ws(T); omega_cinquanta=0.62185.*cinquanta./(P-cinquanta); quaranta=0.4.*P_ws(T); omega_quaranta=0.62185.*quaranta./(P-quaranta); trenta=0.3.*P_ws(T); omega_trenta=0.62185.*trenta./(P-trenta); venti=0.2.*P_ws(T); omega_venti=0.62185.*venti./(P-venti); dieci=0.1.*P_ws(T); omega_dieci=0.62185.*dieci./(P-dieci); T_building=23; RH_building=0.7; omega_building=omega(T_building,RH_building,P); h_building=EnthalpyMoistAir(T_building,omega_building); a=omega_building*ones(size(T)); TOA=5; RHOA=0.8; omegaOA=omega(TOA,RHOA,P); Tiso=(h_building-omegaOA.*2501)./(1.006+1.805.*omegaOA); Tisenthalpic=linspace(19,45); omegaisenthalpic=(h_building-1.006.*Tisenthalpic)./(2501+1.805.*Tisenthalpic); h=plot(T_building,omega_building,'kp',T,omega_cento,'k',T,omega_novanta,'k:',T,omega_ottanta,'k:',T,omega_settanta,'k:',... T,omega_sessanta,'k:',T,omega_cinquanta,'k:',T,omega_quaranta,'k:',T,omega_trenta,'k:',... T,omega_venti,'k:',T,omega_dieci,'k:',TOA,omegaOA,'r*',T_building,omega_building,Tisenthalpic,omegaisenthalpic,'k--'); set(h,'MarkerSize',10) axis([0,45,0.0001,0.02]) TPB=[20 21 22 23 24 25 26 26 25 24 23 22 21 20]; RHPB=[.65 .65 .65 .65 .65 .65 .65 .75 .75 .75 .75 .75 .75 .75]; omegaPB=omega(TPB,RHPB,P); patch(TPB,omegaPB,'green','FaceAlpha',.2) TPB=[22 23 24 24 23 22]; RHPB=[.68 .68 .68 .72 .72 .72]; omegaPB=omega(TPB,RHPB,P); patch(TPB,omegaPB,'green','FaceAlpha',.4) %P1=[TOA omegaOA]; %P2=[Tiso omegaOA]; %dp=P2-P1; %hold on %quiver(P1(1),P1(2),dp(1),dp(2),0,'MaxHeadSize',0.5) xlabel('Temperature [\circC]','FontSize',20) ylabel('Moisture content [kg_H_2_O/kg_d_a]','FontSize',20) set(gca,'FontSize',18)