想请问一下这个怎么改 ？

panphp

/ k. X. R; N* f% I, R4 C* ]/ kclear;
x_len=144;
! d( u$ U* j, q7 K" n- |; py_len=60;
1 o4 H" W  D% `+ Xz_len=17;
3 {7 S0 B7 r' t3 L5 ct_len=854;
z_index=7;
x_index=120;

z = ncread('rhum.mon.mean.nc', 'level');%z
y = ncread('rhum.mon.mean.nc', 'lat');%y
1 M4 ?( D& f( H! B, u9 ?. ]' tx = ncread('rhum.mon.mean.nc', 'lon');%x
6 j0 o2 _$ s; C8 Mt = ncread('rhum.mon.mean.nc', 'time');%t
# X$ e# H4 D; T3 Yrhum = ncread('rhum.mon.mean.nc', 'rhum');

array = zeros(x_len,y_len);
for i = 30:y_len
    for j = 1:t_len
        array(i,j)=rhum(x_index,i,z_index,j);
3 Y! S- J8 p/ @    end
3 X. \1 a" z2 q1 p" r) h' {' Uend
: G9 m+ G/ R/ i" Omesh(t,y,array)
- `" U0 O: }7 e* W( B错误使用 mesh (line 79)
数据维度必须一致。
. L9 I* W; ^2 `9 R, I7 l想请问一下这个怎么改 ？

House

需要先用meshgrid根据横、纵坐标重新生成对应的矩阵
; |$ v/ b. S" q9 R/ z2 e* a另外注意循环嵌套的顺序，是Z(jj,ii)而不是Z(ii,jj)
draw_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)
    for jj = 1:1:numel(draw_y)
        Z(jj,ii) = f(draw_x(ii),draw_y(jj));
6 C& v# s5 v% W3 ?0 [. N; S    end
' n1 j0 {( L; A9 L# c6 h0 t1 P% _( Zend
[draw_x,draw_y] = meshgrid(draw_x,draw_y);
4 ?* X7 l3 @! r3 n+ Umesh(draw_x,draw_y,Z);