|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
clear
6 h, `; F! K( x1 k; @( R6 [ncinfo('rhum.mon.mean.nc');# v. M+ O2 I5 j' u: k b. i
y_len=60;; I$ v* a! z, `! j( v$ S, z
t_len=854;
7 L6 }( T, ?$ Ez_index=7 Y3 c. @) D& h8 o
x_index=1205 r7 k$ V& P- _* v4 `
x = ncread('rhum.mon.mean.nc','lon');%x/ a9 S/ v- @( V
z = ncread('rhum.mon.mean.nc', 'level');%z9 H# k$ g8 n, T' B: d
y = ncread('rhum.mon.mean.nc', 'lat');%y% E0 `; H2 d! t* \6 p; L
t = ncread('rhum.mon.mean.nc', 'time');%t
! \7 N; s# v) l* H2 t0 r2 frhum = ncread('rhum.mon.mean.nc', 'rhum');
4 {# }: a* j3 j/ i7 J4 _figure; hold on;
, l/ t/ S5 n* V- E! `) Z9 Adraw_y = 30:1:60;
0 A" q" `1 _- s" \; ~. edraw_t = 1:1:854;: i7 B- ^( G" Y' i% E* ^
f = rhum(x_index,draw_y,z_index,draw_t)
6 Y3 k; }* |- i( a/ I5 R+ E# Rfor ii = 1:1:numel(draw_y)
, |2 F( {& Y' q; q% f0 d for jj = 1:1:numel(draw_t)
9 ?$ `1 O7 ]0 X" J4 V6 u' ] Z(jj,z_index,ii,x_index) = f(x_index,draw_y(ii),z_index,draw_t(jj));
/ _9 [: u2 T" v* i) R) \ end/ J" \) k! U, T/ d- z- i* u
end8 J9 } b7 ^( B
[draw_y,draw_t] = meshgrid(draw_y,draw_t);) C' X4 {2 ^2 K* E3 \& n5 R8 `
contour (draw_y,draw_t,Z);
: A" c- ^. t6 E+ _. g5 _; m; \0 }. q W( w
错误显示 试图访问 f(120,30,7,1);由于 size(f)=[1,31,1,854],索引超出范围。- H5 M9 O* f, U1 P7 n
改了循环里面之后会显示 试图访问 f(1,32,1,1);由于 size(f)=[1,31,1,854],索引超出范围。4 [0 b% O' `+ h4 f
请教大神,非常感谢!
) Z- E6 E. G5 P' X- F$ [' R- F |
|