|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
clear
+ m; T( T6 M& Yncinfo('rhum.mon.mean.nc');
2 w6 C' Y# ^. N/ v( h( Wy_len=60;
. ?: y( V% ?" k% H1 S+ Dt_len=854;
) d! m3 j5 e7 V, cz_index=7/ A. o. l4 |4 B
x_index=120
* O& }' W5 _! Q9 Yx = ncread('rhum.mon.mean.nc','lon');%x M3 L- _4 K2 U5 ^, n1 h | ?
z = ncread('rhum.mon.mean.nc', 'level');%z
6 m' i* k0 L1 }2 Iy = ncread('rhum.mon.mean.nc', 'lat');%y/ N+ w5 H: f- H1 E! x
t = ncread('rhum.mon.mean.nc', 'time');%t: ]1 }9 O' O0 J4 y
rhum = ncread('rhum.mon.mean.nc', 'rhum');
' K# M) Y2 E. tfigure; hold on;
6 f) e, o- g6 Z" R5 fdraw_y = 30:1:60;
8 e7 ^; E: G; k0 A9 ^draw_t = 1:1:854;& f+ Q" w, m% H8 a
f = rhum(x_index,draw_y,z_index,draw_t)+ K$ U4 L/ _. e! Z
for ii = 1:1:numel(draw_y)
}5 t& H& F/ U; } for jj = 1:1:numel(draw_t)
. l& |1 c7 ^% r& C2 ^0 i Z(jj,z_index,ii,x_index) = f(x_index,draw_y(ii),z_index,draw_t(jj));
, C) K, F1 M/ q2 h6 k3 ^6 X end7 n% T9 ` c4 y. v9 D/ M; H1 q
end
& B' D5 K# {, v V4 [[draw_y,draw_t] = meshgrid(draw_y,draw_t);' l: ~! f9 A C
contour (draw_y,draw_t,Z);
- {8 t. R/ J! t* ~; Q5 `' m3 V- g
0 l4 L" F2 l& Q% i错误显示 试图访问 f(120,30,7,1);由于 size(f)=[1,31,1,854],索引超出范围。" B7 ?( p% S4 G
改了循环里面之后会显示 试图访问 f(1,32,1,1);由于 size(f)=[1,31,1,854],索引超出范围。+ h$ G, A5 o0 M' @) c5 }' @2 D9 K }
请教大神,非常感谢!. Z: I) i1 c" c* ]
|
|