matlab计算的柱相对湿度有很多大于1的值，找不到原因。。

蘑菇

以下是我的程序，经纬度是240*121，分为27层。1-21为中层，21-27为低层，此为上半年数据所以是181天，求得的柱相对湿度有很多大于1的值，不应该是大于0.85的就不多了吗？搞不明白问题出在哪里了，跪求指导啊{:soso_e196:}{:soso_e196:}



%由温度求饱和水汽压
[es]=GroffEw(temp);

%求饱和比湿
sshum=0.622*es/1013.2;

%积分得850h以上中层w*和850h以下低层w*
for i=1:240;
   for j=1:121;
      for t=1:181;
          aa=sshum(i,j,1:21,t);
          aa=squeeze(aa);
          wwm(i,j,t)=trapz(level(1:21),aa)/9.8;
          bb=sshum(i,j,21:27,t);
          bb=squeeze(bb);
          wwl(i,j,t)=trapz(level(21:27),bb)/9.8;
      end;
   end;  
end;

%比湿：积分得850h以上中层w和850h以下低层w
for i=1:240;
   for j=1:121;
      for t=1:181;
          aa=shummean(i,j,1:21,t);
          aa=squeeze(aa);
          wm(i,j,t)=trapz(level(1:21),aa)/9.8;
          bb=shummean(i,j,21:27,t);
          bb=squeeze(bb);
          wl(i,j,t)=trapz(level(21:27),bb)/9.8;
      end;
   end;  
end;

%求中层和低层的CRH
CRHm2010=wm./wwm;
CRHl2010=wl./wwl;