|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
clear$ N) j' I$ r7 l' @; t4 U
ncinfo('rhum.mon.mean.nc');! e1 Z8 U o9 ]
y_len=60;% ]( ]0 Q c0 h6 D/ V1 u7 r
t_len=854;( x% `! ^- x- [: ]& d$ }7 a
z_index=74 {- s3 q! o7 C3 E6 {' s! e0 q" B5 o: K
x_index=120
, S% F. b, i( p+ Vx = ncread('rhum.mon.mean.nc','lon');%x
; N# g* z! C$ iz = ncread('rhum.mon.mean.nc', 'level');%z) X( q6 o+ ^/ E. A) |1 @' w
y = ncread('rhum.mon.mean.nc', 'lat');%y: c+ n) a( a2 L9 Z$ s& K
t = ncread('rhum.mon.mean.nc', 'time');%t
7 G& K! `8 }3 w- Q d) Qrhum = ncread('rhum.mon.mean.nc', 'rhum');
6 _% C* o$ d" e% Q% @4 P0 y" @figure; hold on;$ i, R9 L( }; a9 k: ], T
draw_y = 30:1:60;( F3 `# i( Z$ ]: J, d
draw_t = 1:1:854;8 N/ l& I- t4 ?; M. E$ E
f = rhum(x_index,draw_y,z_index,draw_t): M$ t9 \+ m- M r" `6 H l
for ii = 1:1:numel(draw_y)( b; `% p# z+ N
for jj = 1:1:numel(draw_t)
X* p2 v! |: Z/ L Z(jj,z_index,ii,x_index) = f(x_index,draw_y(ii),z_index,draw_t(jj));0 ?# F) g% a* J1 v$ Z3 A
end1 v, h b# {* i
end+ v# W2 b) z* e
[draw_y,draw_t] = meshgrid(draw_y,draw_t);6 U4 K7 h! {, E/ E- \
contour (draw_y,draw_t,Z);
3 f1 u, d' w* Z5 |. @. u* J) c) U3 ]4 ~* O, ?
错误显示 试图访问 f(120,30,7,1);由于 size(f)=[1,31,1,854],索引超出范围。
6 f- b+ R4 O) g4 j改了循环里面之后会显示 试图访问 f(1,32,1,1);由于 size(f)=[1,31,1,854],索引超出范围。
* I/ ^+ H1 P/ y* G# \请教大神,非常感谢!
/ W; L, O1 s7 u; w |
|