%process ciras data clear din = 'C:\Users\cd833\Desktop\docs\newpapers\girdle\dataz\photo\C8_Licor\'; %din = 'Z:\Desktop\newpapers\girdle\dataz\\'; pp = 1; w=1; r=1; A = dir(din); datasxfall = zeros(1,36); for i = 3:5530 cfname = A(i).name; subdin = [din cfname ]; %[sta] = xlsread(subdin, '%s', 'headerlines', 0, 'delimiter', ' '); datasxf = xlsread(subdin); xx = size(datasxf); if xx(1)*xx(2)>100 An(pp) = nanmean(datasxf(20:end,5)); Astd(pp) = nanstd(datasxf(20:end,5)); gs(pp) = nanmean(datasxf(20:end,6)); gsstd(pp) = nanstd(datasxf(20:end,6)); par(pp) = nanmean(datasxf(20:end,25)); parstd(pp) = nanstd(datasxf(20:end,25)); temp(pp) = nanmean(datasxf(20:end,16)); tempstd(pp) = nanstd(datasxf(20:end,16)); CO2(pp) = nanmean(datasxf(20:end,18)); %to do Ci(pp) = nanmean(datasxf(20:end,7)); Cistd(pp) = nanstd(datasxf(20:end,7)); Trmmol(pp) = nanmean(datasxf(20:end,8)); Trmmolstd(pp) = nanstd(datasxf(20:end,8)); VPD(pp) = nanmean(datasxf(20:end,9)); name(pp) = {cfname}; sxcv = size(datasxf); if par(pp)>800 if An(pp)>0 Anf(w) = An(pp); Astdf(w) = Astd(pp); namef(w) = name(pp); parf(w) = par(pp); tempf(w) = temp(pp); gsf(w) = gs(pp); Cif(w) = Ci(pp); Trmmolf(w) = Trmmol(pp); VPDf(w) = VPD(pp); else Anf(w) = NaN; Astdf(w) = NaN; namef(w) = name(pp); parf(w) = NaN; tempf(w) = NaN; gsf(w) = NaN; Cif(w) = NaN; Trmmolf(w) = NaN; VPDf(w) = NaN; end w=w+1; else if An(pp)<-3.5 Res(r) = NaN; Resstdf(r) = Astd(pp); namer(r) = name(pp); parr(r) = NaN; tempr(r) = NaN; gsr(r) = NaN; Cir(r) = NaN; Trmmolr(r) = NaN; VPDr(r) = NaN; else Res(r) = An(pp); Resstdf(r) = Astd(pp); namer(r) = name(pp); parr(r) = par(pp); tempr(r) = temp(pp); gsr(r) = gs(pp); Cir(r) = Ci(pp); Trmmolr(r) = Trmmol(pp); VPDr(r) = VPD(pp); end r=r+1; end pp=pp+1; end end zzz=[Anf; Astdf; parf; tempf; gsf; Cif; Trmmolf; VPDf ]; zzzr=[Res; Resstdf; parr; tempr; gsr; Cir; Trmmolr; VPDr];