探空仪测得数据按照分层报表要求进行处理（程序）

地面接收到的是GPS电子探空仪测得的是每个采样点的数据，而大气研究一般是按照高度对大气进行分层，研究每个高度区间大气的温度、湿度、气压、风场等参数的分布情况。提供给应用部门的气象产品往往也需要使用分层垂直剖面大气数据。因此，需要对探空仪测得的数据按高度区间进行处理，得到分层数据报表。

获得的GPS探空仪探测数据，需要对数据按照分层数据报表要求处理。由于低空大气活动比较强烈，高度2000m以下按照50m步长进行分层；高空大气活动相对平缓，因此2000m～10000m按照200m步长进行分层，10000m～30000m按照500m步长进行分层。

clear  
clc  
load D:\MATLAB7\work\qiu\snd1.tpu %存放有时间、温度、湿度、气压、高度  
load D:\MATLAB7\work\qiu\snd2.wnd %存放有风向、风速  

temp=snd1(:,2);%温度  
shi=snd1(:,3);%湿度  
pres=snd1(:,4);%压强  
h=snd1(:,6);%高度  
fx=snd2(:,3)*pi/180;%风向  
fs=snd2(:,2);%风速  

%高度0~2000m,40组数据，间隔50m  
nn=1;  
for i=1:40  
xuhao=find(h<=(50*i) & h>(50*(i-1)));  
n=length(xuhao);