contour函数报错怎么办？

panphp

clear
# d, P+ w3 l) K  y: vncinfo('rhum.mon.mean.nc');
0 B; S* {. k2 |) a: O* c6 Ky_len=60;
t_len=854;
z_index=7
x_index=120
8 _2 U1 p( j( p9 c7 n# Vx = ncread('rhum.mon.mean.nc','lon');%x
, v8 G9 p- |9 _5 mz = ncread('rhum.mon.mean.nc', 'level');%z
6 g3 K& g" }, D: |9 e1 A8 [y = ncread('rhum.mon.mean.nc', 'lat');%y
* D& `* h# x& e5 ]! e' J) c' i* ot = ncread('rhum.mon.mean.nc', 'time');%t
/ Z- {% A4 x0 U3 J! ?rhum = ncread('rhum.mon.mean.nc', 'rhum');
7 x0 F3 a7 Z5 `" l8 i: tfigure; hold on;
' i# H, c9 Q2 X6 r/ _4 t  gdraw_y = 30:1:60;
draw_t = 1:1:854;
1 m4 r! K8 G$ m% Q+ x# Wf = rhum(x_index,draw_y,z_index,draw_t)
for ii = 1:1:numel(draw_y)
    for jj = 1:1:numel(draw_t)
        Z(jj,z_index,ii,x_index) = f(x_index,draw_y(ii),z_index,draw_t(jj));
    end
$ ~% t7 F9 Y/ P; e0 z; F) Bend
[draw_y,draw_t] = meshgrid(draw_y,draw_t);
8 v6 o; R, @9 `% acontour (draw_y,draw_t,Z);
- K! v6 C4 R& U% N! R5 c4 F/ ]) h
4 b) o0 t+ t3 O  z  s% H7 G& g错误显示 试图访问 f(120,30,7,1)；由于 size(f)=[1,31,1,854]，索引超出范围。
改了循环里面之后会显示 试图访问 f(1,32,1,1)；由于 size(f)=[1,31,1,854]，索引超出范围。
请教大神，非常感谢！
! ^  C. H* m- e" r" H

Sky_wm

可能参数超出了范围