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

hpucq

# n0 A* Y. L( S% v8 Oclear
ncinfo('rhum.mon.mean.nc');
. y9 W: Y( I6 ?7 X$ Yy_len=60;
" h$ `* y. O; Qt_len=854;
2 r3 s3 f+ ~1 W5 }z_index=7
3 D" }1 p  r) o; s; f" f: Mx_index=120
5 Q- U; m1 m7 O9 T  g) j% Y% Kx = ncread('rhum.mon.mean.nc','lon');%x
: }8 ?8 f* F/ Y8 P2 ]z = ncread('rhum.mon.mean.nc', 'level');%z
y = ncread('rhum.mon.mean.nc', 'lat');%y
t = ncread('rhum.mon.mean.nc', 'time');%t
rhum = ncread('rhum.mon.mean.nc', 'rhum');
/ B" w6 C& {( x8 afigure; hold on;
for i = 30:y_len
  for j=10:t_len
7 ?9 e( O" \7 Z. M. b[ii,jj]=meshgrid(j,i)
* ?. L7 Q$ E- k, jarray(i,j)=rhum(x_index,i,z_index,j);
contour(i,j,z,[rhum(x_index,i,z_index,j)])
  end
+ P) a2 @0 S3 O9 I2 g' Qend

5 }$ V. e6 w8 a错误使用 contour (line 84)
Z 的大小必须为 2x2 或更大。

5 f& O% d5 R8 F& }8 N% C! F1 j
' D7 k6 S+ h7 X6 h5 v% n$ P' r这个contour函数里 i,j很好定义，主要这是这个rhum，也就是contour里的矩阵z，rhum现在是固定两维，两个自变量。麻烦大神指点一下如何定义这个rhum， 成为一个合适的矩阵z？非常感谢！

Colbie

for i = 30:y_len
  for j=10:t_len
[ii,jj]=meshgrid(j,i)
6 F# q2 I) N" Xarray(i,j)=rhum(x_index,i,z_index,j);
contour(i,j,z,[rhum(x_index,i,z_index,j)])
  end
% d+ h& Z# J5 m7 c$ rend

这里不对，不应在循环里面嵌套contour，因为对循环每一步而言，i和j都只有一个数值，meshgrid不能把它们展开成矩阵，而且contour(i,j,z)里面的这个z也没赋过值。
举个正确用法的例子，循环每一步计算Z的一个元素，全部计算完毕后，用meshgrid把横纵坐标向量展开成矩阵，最后在循环外面用contour绘图
draw_x = 1:1/4:7;
draw_y = 2:1/3:6;
+ z  R4 }1 W( kf = @(x,y) x.^(exp( -y.^2/50 ));
! q- G+ U' ~! |$ U; }6 vfor ii = 1:1:numel(draw_x)
    for jj = 1:1:numel(draw_y)
: {$ v' d5 M8 r' [' S5 y  i        Z(jj,ii) = f(draw_x(ii),draw_y(jj));
    end
end
[draw_x,draw_y] = meshgrid(draw_x,draw_y);
- J0 ~  M2 Y5 ]contour (draw_x,draw_y,Z);