麻烦大神指点一下如何定义这个rhum， 成为一个合适的矩阵z？非常感谢！

hpucq

clear
. H7 a/ }, c! Y" Zncinfo('rhum.mon.mean.nc');
y_len=60;
t_len=854;
z_index=7
3 q! w& Y0 `3 @x_index=120
9 Q! \. k# P0 Wx = ncread('rhum.mon.mean.nc','lon');%x
/ S4 j" I. m9 Qz = ncread('rhum.mon.mean.nc', 'level');%z
y = ncread('rhum.mon.mean.nc', 'lat');%y
t = ncread('rhum.mon.mean.nc', 'time');%t
2 N$ S/ i) _' v+ j7 yrhum = ncread('rhum.mon.mean.nc', 'rhum');
( \# k: _8 l! i6 h# tfigure; hold on;
) }) u! _% M5 w, p" Y( Nfor i = 30:y_len
  for j=10:t_len
[ii,jj]=meshgrid(j,i)
: b( v, X& ~2 Q9 v9 e. p  {5 darray(i,j)=rhum(x_index,i,z_index,j);
contour(i,j,z,[rhum(x_index,i,z_index,j)])
  end
end
3 b  T# S% C: x
) X$ F: i; X0 J- I% z0 h& r错误使用 contour (line 84)
Z 的大小必须为 2x2 或更大。

( I+ y5 k  u2 ?
9 l/ Z1 ?4 n% {这个contour函数里 i,j很好定义，主要这是这个rhum，也就是contour里的矩阵z，rhum现在是固定两维，两个自变量。麻烦大神指点一下如何定义这个rhum， 成为一个合适的矩阵z？非常感谢！

Colbie

for i = 30:y_len
2 |/ G$ [$ G3 j4 x) j$ ~  G  for j=10:t_len
[ii,jj]=meshgrid(j,i)
# @  U7 D' _) B; earray(i,j)=rhum(x_index,i,z_index,j);
* v, F2 n/ W+ c' H( ^% K+ Bcontour(i,j,z,[rhum(x_index,i,z_index,j)])
( \. m& E  i5 v  end
. R9 a* S' Q/ O' ~end

! z. m3 P. e) f, C- y' N5 h! k这里不对，不应在循环里面嵌套contour，因为对循环每一步而言，i和j都只有一个数值，meshgrid不能把它们展开成矩阵，而且contour(i,j,z)里面的这个z也没赋过值。
举个正确用法的例子，循环每一步计算Z的一个元素，全部计算完毕后，用meshgrid把横纵坐标向量展开成矩阵，最后在循环外面用contour绘图
2 d9 @7 p. e! Pdraw_x = 1:1/4:7;
draw_y = 2:1/3:6;
f = @(x,y) x.^(exp( -y.^2/50 ));
for ii = 1:1:numel(draw_x)
, b! T) c3 `6 Y( j    for jj = 1:1:numel(draw_y)
        Z(jj,ii) = f(draw_x(ii),draw_y(jj));
    end
end
[draw_x,draw_y] = meshgrid(draw_x,draw_y);
contour (draw_x,draw_y,Z);