|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
clear. U, {7 A( g' y. a6 X
ncinfo('rhum.mon.mean.nc');
' }4 U6 b/ b- p6 y2 ^7 yy_len=60;
- u8 v" }& f R" st_len=854;
5 v* d% Z3 n* [; O) m2 Nz_index=7
0 ~- L! ~3 U% \' i8 T4 Vx_index=120
# _. F% X# u7 p4 u( Px = ncread('rhum.mon.mean.nc','lon');%x) i3 D0 f* N5 ^2 k3 J3 c
z = ncread('rhum.mon.mean.nc', 'level');%z
^2 T! s; H T# zy = ncread('rhum.mon.mean.nc', 'lat');%y
$ {- X4 c- j7 Z/ s. y9 Pt = ncread('rhum.mon.mean.nc', 'time');%t& [- O& r+ `" U- B/ f$ b8 L5 I
rhum = ncread('rhum.mon.mean.nc', 'rhum');
( l/ G4 u- z) I" M) O9 x+ i6 Wfigure; hold on;
. J. C. L2 K- J' y- Gdraw_y = 30:1:60;2 Z! B9 {7 B5 |+ S9 v7 {
draw_t = 1:1:854;
3 W9 K* [9 L1 y) p$ Uf = rhum(x_index,draw_y,z_index,draw_t)
) O, T; v% a- m; f! t9 Ifor ii = 1:1:numel(draw_y)& M: O: h% ]* W, {# c- K; ^5 P8 K
for jj = 1:1:numel(draw_t)
/ i4 K3 |' |+ z+ w% B; G# |5 ] Z(jj,z_index,ii,x_index) = f(x_index,draw_y(ii),z_index,draw_t(jj));. \# ~' f( l) T* b+ p$ M
end
2 z+ _) ]; D8 lend/ b K0 g+ L# o. Q
[draw_y,draw_t] = meshgrid(draw_y,draw_t);
8 U) e/ w. ]0 ?& ccontour (draw_y,draw_t,Z);' R" C: m$ E: f: d8 j2 w% ~4 c
/ `4 _# h, F, m5 c4 [) P' U错误显示 试图访问 f(120,30,7,1);由于 size(f)=[1,31,1,854],索引超出范围。) k0 X3 Y W- Q2 i
改了循环里面之后会显示 试图访问 f(1,32,1,1);由于 size(f)=[1,31,1,854],索引超出范围。3 z+ {1 ~+ t/ `
请教大神,非常感谢!: ^4 ^$ }* L7 \+ q0 N: o' A- ]; T2 }
|
|