how to get the data of temperature pressure humidity at 20 km - program

clear all;fclose all;clc;close all;clear a=fopen('files.dat','r'); gg=0; while ~feof(a) fgets(a); gg=gg+1; end frewind(a) for i=1:gg b=fscanf(a,'%s',1); c=load(b); high(:,i)=c(:,1); pres(:,i)=c(:,2); temp(:,i)=c(:,3); rh(:,i)=c(:,4); u(:,i)=c(:,5); v(:,i)=c(:,6); ws(:,i)=c(:,7); wd(:,i)=c(:,8); end

 for i=1:400
t(i,:)=nanmean(tempe(i,:));
pres(i,:)=nanmean(pres(i,:));
rhs(i,:)=nanmean(rh(i,:));
uw(i,:)=nanmean(u(i,:));
vw(i,:)=nanmean(v(i,:));
wsp(i,:)=nanmean(ws(i,:));
wdr(i,:)=nanmean(wd(i,:));

 end
subplot(241);plot(t,hight);title('temperature');xlabel('temp');ylabel('high')
subplot(242);plot(pre,hight);title('pressure');
subplot(243);plot(rhs,hight);title('relative humidity');
subplot(244);plot(uw,hight);title('u-zonal');
subplot(245);plot(vw,hight);title('v-meriodanal');
subplot(246);plot(wsp,hight);title('wind speed');
subplot(247);plot(wdr,hight);title('wind direction');
for i=1:400
tdif(i,:)=t(i,:)-temp(i,:);
presdif(i,:)=pres(i,:)-pres(i,:);
rhsdif(i,:)=rhs(i,:)-rh(i,:);
uwdif(i,:)=uw(i,:)-u(i,:);
vwdif(i,:)=vw(i,:)-v(i,:);
wspdif(i,:)=wsp(i,:)-ws(i,:);
wdrdif(i,:)=wdr(i,:)-wd(i,:);

 end