|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
clear$ G' }' D7 g& R2 \$ t
ncinfo('rhum.mon.mean.nc');
* r' m8 b. Y% h+ fy_len=60;& V$ P% @! t7 G
t_len=854;
6 m3 w5 I; V6 t5 |& }z_index=7
- k }( G! e1 r+ u3 qx_index=120
+ {* ^, R3 N) D& b! w4 p7 H7 Wx = ncread('rhum.mon.mean.nc','lon');%x7 s2 S0 N) c+ \2 K5 }9 t2 X, T
z = ncread('rhum.mon.mean.nc', 'level');%z- J- H% w j+ h4 L* o
y = ncread('rhum.mon.mean.nc', 'lat');%y
9 @) M$ a( i, Q* v5 W' e, Rt = ncread('rhum.mon.mean.nc', 'time');%t' u- t% X Z. r- B
rhum = ncread('rhum.mon.mean.nc', 'rhum');9 ?/ \0 O' m, \) x- g
figure; hold on;4 g: y$ I0 S0 D
draw_y = 30:1:60;
9 g8 p- {- v8 udraw_t = 1:1:854;# n+ y, ]8 H" ^- o8 G3 E6 X
f = rhum(x_index,draw_y,z_index,draw_t)5 g0 B. s$ n+ V
for ii = 1:1:numel(draw_y)6 C: c) e( s" \' A0 q
for jj = 1:1:numel(draw_t)% E8 j( n3 O% x8 T# A" D$ l
Z(jj,z_index,ii,x_index) = f(x_index,draw_y(ii),z_index,draw_t(jj));
% Q% G% x) U' y* Y end7 m9 z+ n; H* G9 s; Q: x8 X
end1 {$ o5 a0 [4 X0 z: j
[draw_y,draw_t] = meshgrid(draw_y,draw_t);2 S( i( s8 d7 a9 R1 R4 k- }2 ~& ^
contour (draw_y,draw_t,Z);7 K8 [. b1 S6 U% G8 Z0 }' s. X
% v& j8 t6 J' F错误显示 试图访问 f(120,30,7,1)；由于 size(f)=[1,31,1,854]，索引超出范围。
$ C+ f' ?# R9 P3 E3 Z/ f% l改了循环里面之后会显示 试图访问 f(1,32,1,1)；由于 size(f)=[1,31,1,854]，索引超出范围。
6 W5 }7 G2 {# q3 u7 @" K$ b请教大神，非常感谢！
5 H) M! ^( s- S {' n' q |
|
/1 
发表于 2020-6-2 16:07
收藏
淘帖
支持！
反对！