%%%% 9 model cloud Od and cloud Fraction % GC #1 K1 = ncread('Jstat_cloudODCF.nc','K1') ; P1 = ncread('Jstat_cloudODCF.nc','P1') ; dP1 = ncread('Jstat_cloudODCF.nc','dP1') ; CF1 = ncread('Jstat_cloudODCF.nc','CF1') ; OD1 = ncread('Jstat_cloudODCF.nc','OD1') ; ODp1 = ncread('Jstat_cloudODCF.nc','ODp1'); sum(OD1) % 1.8500 2.5500 % GFDL #2 K2 = ncread('Jstat_cloudODCF.nc','K2') ; P2 = ncread('Jstat_cloudODCF.nc','P2') ; dP2 = ncread('Jstat_cloudODCF.nc','dP2') ; CF2 = ncread('Jstat_cloudODCF.nc','CF2') ; OD2 = ncread('Jstat_cloudODCF.nc','OD2') ; ODp2 = ncread('Jstat_cloudODCF.nc','ODp2'); sum(OD2) % 10.7287 9.7873 % GISS #3 K3 = ncread('Jstat_cloudODCF.nc','K3') ; P3 = ncread('Jstat_cloudODCF.nc','P3') ; dP3 = ncread('Jstat_cloudODCF.nc','dP3') ; CF3 = ncread('Jstat_cloudODCF.nc','CF3') ; OD3 = ncread('Jstat_cloudODCF.nc','OD3') ; ODp3 = ncread('Jstat_cloudODCF.nc','ODp3'); sum(OD3) % 7.5419 10.8292 % GMI #4 K4 = ncread('Jstat_cloudODCF.nc','K4') ; P4 = ncread('Jstat_cloudODCF.nc','P4') ; dP4 = ncread('Jstat_cloudODCF.nc','dP4') ; CF4 = ncread('Jstat_cloudODCF.nc','CF4') ; OD4 = ncread('Jstat_cloudODCF.nc','OD4') ; ODp4 = ncread('Jstat_cloudODCF.nc','ODp4'); sum(OD4) % 1.8011 2.8352 % IFS-rev #5 K5 = ncread('Jstat_cloudODCF.nc','K5') ; P5 = ncread('Jstat_cloudODCF.nc','P5') ; dP5 = ncread('Jstat_cloudODCF.nc','dP5') ; CF5 = ncread('Jstat_cloudODCF.nc','CF5') ; OD5 = ncread('Jstat_cloudODCF.nc','OD5') ; ODp5 = ncread('Jstat_cloudODCF.nc','ODp5'); sum(OD5) % 0.7366 0.8643 % MOCA #6 K6 = ncread('Jstat_cloudODCF.nc','K6') ; P6 = ncread('Jstat_cloudODCF.nc','P6') ; dP6 = ncread('Jstat_cloudODCF.nc','dP6') ; CF6 = ncread('Jstat_cloudODCF.nc','CF6') ; OD6 = ncread('Jstat_cloudODCF.nc','OD6') ; ODp6 = ncread('Jstat_cloudODCF.nc','ODp6'); sum(OD6) % 0.4796 1.6332 % NCAR #7 K7 = ncread('Jstat_cloudODCF.nc','K7') ; P7 = ncread('Jstat_cloudODCF.nc','P7') ; dP7 = ncread('Jstat_cloudODCF.nc','dP7') ; CF7 = ncread('Jstat_cloudODCF.nc','CF7') ; OD7 = ncread('Jstat_cloudODCF.nc','OD7') ; ODp7 = ncread('Jstat_cloudODCF.nc','ODp7'); sum(OD7) % 4.4318 4.1618 % UCI #8 K8 = ncread('Jstat_cloudODCF.nc','K8') ; P8 = ncread('Jstat_cloudODCF.nc','P8') ; dP8 = ncread('Jstat_cloudODCF.nc','dP8') ; CF8 = ncread('Jstat_cloudODCF.nc','CF8') ; OD8 = ncread('Jstat_cloudODCF.nc','OD8') ; ODp8 = ncread('Jstat_cloudODCF.nc','ODp8'); sum(OD8) % 3.1733 5.2393 % UKCA #9 K9 = ncread('Jstat_cloudODCF.nc','K9') ; P9 = ncread('Jstat_cloudODCF.nc','P9') ; dP9 = ncread('Jstat_cloudODCF.nc','dP9') ; CF9 = ncread('Jstat_cloudODCF.nc','CF9') ; OD9 = ncread('Jstat_cloudODCF.nc','OD9') ; ODp9 = ncread('Jstat_cloudODCF.nc','ODp9'); sum(OD9) % 0.7384 4.3064 TOD(:,1) = sum(OD1); TOD(:,2) = sum(OD2); TOD(:,3) = sum(OD3); TOD(:,4) = sum(OD4); TOD(:,5) = sum(OD5); TOD(:,6) = sum(OD6); TOD(:,7) = sum(OD7); TOD(:,8) = sum(OD8); TOD(:,9) = sum(OD9); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %Mclr = linespecer(10,'sequential'); %Mclr = Mclr*0.9; %Mclr(5,:) = RGB('LightGray'); %Mclr(6,:) = RGB('Gold'); Mclr = ... [0.3318, 0.2788, 0.5718; 0.1805, 0.5034, 0.6642; 0.4102, 0.7108, 0.5812; 0.6855, 0.8152, 0.5669; 0.8242, 0.8242, 0.8242; 1.0000, 0.8398, 0; 0.8948, 0.6824, 0.3880; 0.8716, 0.4319, 0.2451; 0.7673, 0.2389, 0.2774; 0.5576, 0.0035, 0.2329]; Mclr(11,:) = [ 0 0 0]; Mclr(9,:) = Mclr(10,:); Mdls = 09; Mlb = ['GC '; 'GFDL'; 'GISS'; 'GMI '; 'IFS '; 'MOCA'; 'NCAR'; 'UCI '; 'UKCA']; Plb =['Tr.Pacific'; 'No.Pacific']; Xlb = ['Cloud Optical Depth (/100 hPa)'; 'Cloud Fraction ']; Xlm = [ 10, 0.4]; Xlm0= [ .01, 0.0]; Ylb = ['pressure (hPa)']; ODx = [0; 0]; Px = [1000; 999]; %%%% for plots, set min on ODp# so no line breaks ODp1 = max(ODp1,Xlm0(1)); % ODp2 = max(ODp2,Xlm0(1)); ODp3 = max(ODp3,Xlm0(1)); ODp4 = max(ODp4,Xlm0(1)); ODp5 = max(ODp5,Xlm0(1)); ODp6 = max(ODp6,Xlm0(1)); ODp7 = max(ODp7,Xlm0(1)); ODp8 = max(ODp8,Xlm0(1)); ODp9 = max(ODp9,Xlm0(1)); figure(225) hold on set(gcf, 'units','points','outerposition',[50 50 750 750]); figure(gcf) k=1; for b=1:2 subplot(2,2,b) hold on set(gca, 'YLim',[100, 1000], 'YTick', 100:100:1000, 'YDir', 'reverse') set(gca, 'XLim',[Xlm0(k), Xlm(k)], 'XScale', 'log') set(gca, 'FontSize',12, 'LineWidth',1) box on grid off plot(ODp1(:,b),P1,'Color',Mclr(1,:),'LineStyle','-','LineWidth',3) plot(ODp2(:,b),P2,'Color',Mclr(2,:),'LineStyle','-','LineWidth',3) plot(ODp3(:,b),P3,'Color',Mclr(3,:),'LineStyle','-','LineWidth',3) plot(ODp4(:,b),P4,'Color',Mclr(4,:),'LineStyle','-','LineWidth',3) plot(ODp5(:,b),P5,'Color',Mclr(5,:),'LineStyle','-','LineWidth',3) plot(ODp6(:,b),P6,'Color',Mclr(6,:),'LineStyle','-','LineWidth',3) plot(ODp7(:,b),P7,'Color',Mclr(7,:),'LineStyle','-','LineWidth',3) plot(ODp8(:,b),P8,'Color',Mclr(8,:),'LineStyle','-','LineWidth',3) plot(ODp9(:,b),P9,'Color',Mclr(9,:),'LineStyle','-','LineWidth',3) for m=1:9 text(+3.00,100+50*m, num2str(TOD(b,m),'% 10.2f'),... 'FontSize',13,'Color',Mclr(m,:)) end legend({Mlb(1:Mdls,:)},'FontSize',11,'Location','NorthWest') xlabel(Xlb(k,:)) ylabel(Ylb) text(0.2,150,Plb(b,:),'FontSize',12) end k=2; for b=1:2 subplot(2,2,b+2) hold on set(gca, 'YLim',[100, 1000], 'YTick', 100:100:1000, 'YDir', 'reverse') set(gca, 'XLim',[Xlm0(k), Xlm(k)]) set(gca, 'FontSize',12, 'LineWidth',1) box on grid off % for m=1:Mdls plot(CF1(:,b),P1,'Color',Mclr(1,:),'LineStyle','-','LineWidth',3) plot(CF2(:,b),P2,'Color',Mclr(2,:),'LineStyle','-','LineWidth',3) plot(CF3(:,b),P3,'Color',Mclr(3,:),'LineStyle','-','LineWidth',3) plot(CF4(:,b),P4,'Color',Mclr(4,:),'LineStyle','-','LineWidth',3) plot(CF5(:,b),P5,'Color',Mclr(5,:),'LineStyle','-','LineWidth',3) plot(CF6(:,b),P6,'Color',Mclr(6,:),'LineStyle','-','LineWidth',3) plot(CF7(:,b),P7,'Color',Mclr(7,:),'LineStyle','-','LineWidth',3) plot(CF8(:,b),P8,'Color',Mclr(8,:),'LineStyle','-','LineWidth',3) plot(CF9(:,b),P9,'Color',Mclr(9,:),'LineStyle','-','LineWidth',3) % end legend({Mlb(1:Mdls,:)},'FontSize',11,'Location','NorthWest') xlabel(Xlb(k,:)) ylabel(Ylb) text(0.21,150,Plb(b,:),'FontSize',12) end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%