contour函数报错怎么办？

panphp

clear
1 n5 C$ S$ y$ X+ Y$ Kncinfo('rhum.mon.mean.nc');
$ o, i6 g' z% w# zy_len=60;
t_len=854;
2 S! J7 q2 C, d9 }7 P; Az_index=7
4 {% c8 E6 n- X& t+ t' V% ax_index=120
  [$ M: U; L( px = ncread('rhum.mon.mean.nc','lon');%x
z = ncread('rhum.mon.mean.nc', 'level');%z
' ~$ m! C: L/ i1 E2 ry = ncread('rhum.mon.mean.nc', 'lat');%y
( F) D0 o+ l: }t = ncread('rhum.mon.mean.nc', 'time');%t
& @. ~& x; p3 i/ V: prhum = ncread('rhum.mon.mean.nc', 'rhum');
9 [, {- T: G2 ^" W3 h. N4 U% afigure; hold on;
draw_y = 30:1:60;
draw_t = 1:1:854;
8 Q2 N5 o+ K1 w  x  H) ^f = rhum(x_index,draw_y,z_index,draw_t)
9 q2 s( s/ a  Qfor ii = 1:1:numel(draw_y)
- ?6 j& G. v  V' a1 a  R% O    for jj = 1:1:numel(draw_t)
, |: z* s+ y) W2 O        Z(jj,z_index,ii,x_index) = f(x_index,draw_y(ii),z_index,draw_t(jj));
5 R- P& J9 _( ?  `4 M    end
end
[draw_y,draw_t] = meshgrid(draw_y,draw_t);
contour (draw_y,draw_t,Z);

错误显示 试图访问 f(120,30,7,1)；由于 size(f)=[1,31,1,854]，索引超出范围。
8 F5 A% m* A0 u) J1 X改了循环里面之后会显示 试图访问 f(1,32,1,1)；由于 size(f)=[1,31,1,854]，索引超出范围。
$ i+ o1 X+ X1 y( I( L7 Y请教大神，非常感谢！
' \/ C, B# ]! s

Sky_wm

可能参数超出了范围