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

hpucq

. a) k( ~8 l( j$ {8 E8 f( Rclear
7 z/ Z+ N# Q  x( Sncinfo('rhum.mon.mean.nc');
y_len=60;
& L, d% b$ x! H! ~0 xt_len=854;
2 L# F/ C) N. _4 Wz_index=7
# G" m0 J: {3 E# j0 T2 o4 j* jx_index=120
! F; l( M. E% ~  M/ rx = ncread('rhum.mon.mean.nc','lon');%x
8 u7 D3 H( B) J' o- az = ncread('rhum.mon.mean.nc', 'level');%z
* }! [, \& A, q; |- }y = ncread('rhum.mon.mean.nc', 'lat');%y
) z& H/ q" m8 gt = ncread('rhum.mon.mean.nc', 'time');%t
+ s4 q; B! T6 nrhum = ncread('rhum.mon.mean.nc', 'rhum');
figure; hold on;
for i = 30:y_len
  for j=10:t_len
0 g5 ^, ]) T/ W6 f2 E[ii,jj]=meshgrid(j,i)
4 O0 M3 l1 u& P, K- u2 E1 harray(i,j)=rhum(x_index,i,z_index,j);
contour(i,j,z,[rhum(x_index,i,z_index,j)])
  end
% ]6 _" {$ C/ kend

6 F+ s+ y. m0 K  G* ^错误使用 contour (line 84)
3 J0 M- f7 y3 H* D6 gZ 的大小必须为 2x2 或更大。

& P# b% v7 D$ P5 M
这个contour函数里 i,j很好定义，主要这是这个rhum，也就是contour里的矩阵z，rhum现在是固定两维，两个自变量。麻烦大神指点一下如何定义这个rhum， 成为一个合适的矩阵z？非常感谢！

Colbie

for i = 30:y_len
: Z5 d5 N% I2 a& a  for j=10:t_len
[ii,jj]=meshgrid(j,i)
array(i,j)=rhum(x_index,i,z_index,j);
contour(i,j,z,[rhum(x_index,i,z_index,j)])
  end
end
( {$ a  m9 ]& Y9 r* T' f! E8 K; C
这里不对，不应在循环里面嵌套contour，因为对循环每一步而言，i和j都只有一个数值，meshgrid不能把它们展开成矩阵，而且contour(i,j,z)里面的这个z也没赋过值。
. m+ |5 G+ ?7 V4 s0 v4 U% t举个正确用法的例子，循环每一步计算Z的一个元素，全部计算完毕后，用meshgrid把横纵坐标向量展开成矩阵，最后在循环外面用contour绘图
draw_x = 1:1/4:7;
9 P+ |# A8 g$ s: R& s! rdraw_y = 2:1/3:6;
f = @(x,y) x.^(exp( -y.^2/50 ));
for ii = 1:1:numel(draw_x)
/ z8 l' [8 I1 p1 g, B6 O0 p    for jj = 1:1:numel(draw_y)
1 U4 n2 ~# S6 V( K' E; y        Z(jj,ii) = f(draw_x(ii),draw_y(jj));
    end
end
$ |& L8 D5 s( \/ p3 N1 R/ K[draw_x,draw_y] = meshgrid(draw_x,draw_y);
  B' L5 ~; f5 r* ]contour (draw_x,draw_y,Z);