station='UNEC';

if ~isfolder([station '_trop_hourly_elev']), mkdir([station '_trop_hourly_elev']); end
lis1=ls([station '_trop\*.cmn']);
for i=1:size(lis1,1)
disp(['Processing,... ' num2str(i) ' of ' num2str(size(lis1,1))]);
    filename=lis1(i,:);
    fid2=fopen([station '_trop/' filename]);
    data1=textscan(fid2, '%f %f %f %f %f %f %f %f %f %f %f', 'headerlines', 1);
    hour1=data1{2}; el1=data1{4}; trop1=data1{11};
fido=fopen([station '_trop_hourly_elev/' filename], 'wt');
fprintf(fido, 'Hour\t Elevation\t Trop_delay\n');
forhh=0:23
        trop2=ones(1,9)*nan;
foreell=10:10:90
ff=find(floor(hour1)==hh& el1>eell-10 & el1<eell);
            trop2(round(eell/10))=nanmedian(trop1(ff));
end
fprintf(fido, '%2.0f\t %2.0f\t %3.4f\n', [ones(1,9)*hh; 10:10:90; trop2]);
end
fcloseall;
end