|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
clear
. e, X8 P3 W* [3 w8 zncinfo('rhum.mon.mean.nc');$ ^' D( Z" m+ U$ H% u S2 U! }5 U
y_len=60;8 _5 e( R$ ]9 j: f' J
t_len=854;
2 M) t" d$ ^0 w& z0 F0 s! ~2 Wz_index=7+ r4 Q5 R1 @: u, P, b# m- ^! p
x_index=1205 i4 E+ s6 q2 y% G
x = ncread('rhum.mon.mean.nc','lon');%x$ Y& @) N% f) F) s. ~
z = ncread('rhum.mon.mean.nc', 'level');%z @- j/ v( F7 C- F; L
y = ncread('rhum.mon.mean.nc', 'lat');%y) Q4 B4 a i$ Y9 \' B7 }
t = ncread('rhum.mon.mean.nc', 'time');%t% n) [: |, ^4 Y! _2 y1 ^
rhum = ncread('rhum.mon.mean.nc', 'rhum');3 j4 e8 e' F: ^( h- x$ _
figure; hold on;5 P3 n+ Y' B+ n: t R6 [
draw_y = 30:1:60;5 @9 D X B. w& t# o3 k Q1 b
draw_t = 1:1:854;
6 f- O2 s' `! V9 y9 Of = rhum(x_index,draw_y,z_index,draw_t); O* U/ I. N& q9 N
for ii = 1:1:numel(draw_y)# w. N- t. R. Q: r g
for jj = 1:1:numel(draw_t) ?, x4 j4 b# ^/ N
Z(jj,z_index,ii,x_index) = f(x_index,draw_y(ii),z_index,draw_t(jj));
: {! X# ^7 N9 v$ B2 j end4 O4 n9 i) B/ H% D5 @/ j
end( L9 r1 s& V# ^+ F' C" L
[draw_y,draw_t] = meshgrid(draw_y,draw_t);
7 H3 c- B O9 e: n- j: v2 n+ ocontour (draw_y,draw_t,Z);& }' G1 e5 X. u& W
1 }# u' c( `& h& {0 |) i
错误显示 试图访问 f(120,30,7,1);由于 size(f)=[1,31,1,854],索引超出范围。4 r6 O/ h& w$ O# t- M
改了循环里面之后会显示 试图访问 f(1,32,1,1);由于 size(f)=[1,31,1,854],索引超出范围。: f0 Z: `) ^- Q: `, Z# G
请教大神,非常感谢!% C- @# l, C6 `
|
|