|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
clear9 H9 f3 S) h5 h' P0 \4 @. M8 q
ncinfo('rhum.mon.mean.nc');8 b, |' S/ |& b% }3 L7 Y
y_len=60;2 K. N' R: e% W7 b5 }5 b
t_len=854;) l( H, `. p6 s' K
z_index=7( \6 p+ F W5 w- e
x_index=120( X+ X0 b0 X8 Z& N2 S
x = ncread('rhum.mon.mean.nc','lon');%x
# {& `2 c N/ |5 w( a8 p2 ez = ncread('rhum.mon.mean.nc', 'level');%z
/ }9 n/ o! Y8 }" w1 w o8 oy = ncread('rhum.mon.mean.nc', 'lat');%y* @( i x" x/ L3 u$ r' n
t = ncread('rhum.mon.mean.nc', 'time');%t& @+ _; | A# P& U5 M, C
rhum = ncread('rhum.mon.mean.nc', 'rhum');* h5 s# w, A$ r0 ~0 N @
figure; hold on;& v4 b n" b' A/ T, h- S! t3 T
draw_y = 30:1:60;
# D$ t; ?" R+ xdraw_t = 1:1:854;
9 t. N. y/ r! P/ {f = rhum(x_index,draw_y,z_index,draw_t)
& y& C; C7 E4 J* N. Vfor ii = 1:1:numel(draw_y)4 q% i# y. E: x% Q* G* {1 ]
for jj = 1:1:numel(draw_t)! p; e) ?- W- Z4 y) `0 d
Z(jj,z_index,ii,x_index) = f(x_index,draw_y(ii),z_index,draw_t(jj));
( S2 Y4 k1 D" J/ n w8 c' n end. h2 p$ ~3 t2 P) j$ @' Z2 c
end; M" T2 w" O( w4 i8 y0 G
[draw_y,draw_t] = meshgrid(draw_y,draw_t);
1 y8 a4 g7 }5 H" [' s& wcontour (draw_y,draw_t,Z);
. P- l8 C- h: ^2 _+ Q& L8 S2 V' e: S2 }; \0 [) g
错误显示 试图访问 f(120,30,7,1);由于 size(f)=[1,31,1,854],索引超出范围。7 v3 G3 d/ a% A- q& Z
改了循环里面之后会显示 试图访问 f(1,32,1,1);由于 size(f)=[1,31,1,854],索引超出范围。
/ B2 H# Z: f3 O# W0 c请教大神,非常感谢!
: p& W, ]* s9 T" b$ e |
|