|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
clear. u: |# w0 g. x$ Q) q C) U/ X
ncinfo('rhum.mon.mean.nc');# T" i- l( x" l9 W( F( M
y_len=60;2 h& T+ j9 y( `2 p" l% i0 M% x- G' r, \
t_len=854;7 O& E5 e9 _7 C' R% L
z_index=7
3 E6 X8 F' a3 D T8 g+ E# a1 Vx_index=120
, H+ i; R; g% ^5 P2 bx = ncread('rhum.mon.mean.nc','lon');%x2 o" W9 N- z# S7 W0 C, B
z = ncread('rhum.mon.mean.nc', 'level');%z1 h! T6 e$ a3 g
y = ncread('rhum.mon.mean.nc', 'lat');%y9 G9 b. D" I8 O0 j9 |, y2 D
t = ncread('rhum.mon.mean.nc', 'time');%t- ~2 d F$ r! {' P
rhum = ncread('rhum.mon.mean.nc', 'rhum');# k; t& g4 H; \( }
figure; hold on;: I; U7 I5 i9 y _% l! F
draw_y = 30:1:60;7 ?" ^/ ] u9 m2 J& o
draw_t = 1:1:854; g, Z" F6 `7 \) L
f = rhum(x_index,draw_y,z_index,draw_t)9 x) A& z Y' f# t$ S! t6 w
for ii = 1:1:numel(draw_y)1 A7 M. K5 o; K
for jj = 1:1:numel(draw_t)7 C& \& p" `3 z- p
Z(jj,z_index,ii,x_index) = f(x_index,draw_y(ii),z_index,draw_t(jj));
% ?, Z* i, U7 T- ~$ c7 b end) [! A. I# J8 T! p) b! q5 `, E
end
* s. ?( d J! p& Z[draw_y,draw_t] = meshgrid(draw_y,draw_t);
, e& F4 B0 d d) _! I0 Dcontour (draw_y,draw_t,Z);
# X/ u! X# f1 `; Q; N" Q
8 d) ]; f. o; j* Y% B5 ^; b/ s错误显示 试图访问 f(120,30,7,1);由于 size(f)=[1,31,1,854],索引超出范围。
' m( D) h$ ?& S' G改了循环里面之后会显示 试图访问 f(1,32,1,1);由于 size(f)=[1,31,1,854],索引超出范围。' D$ w3 ?( y/ A% }/ Z$ c
请教大神,非常感谢!" T5 Z& u/ N- F3 |4 e! \0 I
|
|