|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
clear* s9 F; J( p* Y0 Y) ^& v7 K
ncinfo('rhum.mon.mean.nc');
% v+ c/ x& R$ Q* g9 U+ ~" n5 Ny_len=60;* o8 h7 C* ]8 ~8 z0 d$ |1 {
t_len=854;8 o, f; c- y1 @! {, `8 B4 v$ p
z_index=7
3 i* H" a; K$ |# t3 Qx_index=120
/ q, Q+ x$ ~/ jx = ncread('rhum.mon.mean.nc','lon');%x
5 `4 K' k: R2 ~' f, Az = ncread('rhum.mon.mean.nc', 'level');%z
) g+ O6 S; q' c+ t5 z4 b. f" F+ Zy = ncread('rhum.mon.mean.nc', 'lat');%y
& b' {6 E* n2 e V' Wt = ncread('rhum.mon.mean.nc', 'time');%t
& @& Z( K0 D/ `/ prhum = ncread('rhum.mon.mean.nc', 'rhum');( U$ {: m" k* \) k$ I5 i% U2 H
figure; hold on;
7 a7 S5 S1 y) n( |/ wdraw_y = 30:1:60;
6 m& U5 A L0 H9 `0 i5 }, Qdraw_t = 1:1:854;
& k' }$ |" n/ O) ^f = rhum(x_index,draw_y,z_index,draw_t)! R9 P* |; {! K# Y8 h
for ii = 1:1:numel(draw_y)
! a) G1 t& i; p, a8 t& S/ s5 m" [0 O1 H+ f for jj = 1:1:numel(draw_t)
9 r- K! v" l3 B6 S$ c3 z Z(jj,z_index,ii,x_index) = f(x_index,draw_y(ii),z_index,draw_t(jj));, A- @2 M7 Y0 Q
end3 f \2 n& h7 C" n2 i; R
end
6 g7 c+ ~- D6 ?. M6 N[draw_y,draw_t] = meshgrid(draw_y,draw_t);
- e& h, g2 B, u" ], dcontour (draw_y,draw_t,Z);
2 y& V! f; z' b- z
, |6 a% S' w `* c' x( \错误显示 试图访问 f(120,30,7,1)；由于 size(f)=[1,31,1,854]，索引超出范围。+ n3 c# A6 O- w& t
改了循环里面之后会显示 试图访问 f(1,32,1,1)；由于 size(f)=[1,31,1,854]，索引超出范围。
5 V4 r; u/ l1 ^( @0 y+ ], U请教大神，非常感谢！
1 f# m7 [* `5 h* ^: o" l |
|
/1 
发表于 2020-6-2 16:07
收藏
淘帖
支持！
反对！