|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
clear
, ?8 Z# ]0 s/ f9 z ^/ U( V$ qncinfo('rhum.mon.mean.nc');
# n. C9 @/ O6 Zy_len=60;
$ d& d7 V3 d9 g v; j- Z) D3 ut_len=854;
5 g M" {3 |7 C$ @z_index=7$ Q+ ^4 M L d' J: p* _( R3 w7 _
x_index=120! O: ^$ O2 x: l) i3 m
x = ncread('rhum.mon.mean.nc','lon');%x- x3 s: w& c5 j7 r" g$ @" f
z = ncread('rhum.mon.mean.nc', 'level');%z
, N" ?: J+ j$ Z) n9 ly = ncread('rhum.mon.mean.nc', 'lat');%y% |/ X9 v2 G. e) L" |
t = ncread('rhum.mon.mean.nc', 'time');%t' U2 h! z! }- O
rhum = ncread('rhum.mon.mean.nc', 'rhum');
) s" m% D- a* gfigure; hold on;
( P) M6 c9 n( L" Sdraw_y = 30:1:60;7 \3 F2 c' t# U: Y0 y8 ^: {
draw_t = 1:1:854;
7 v6 l+ X$ o- x% D+ Z9 W9 B8 ]f = rhum(x_index,draw_y,z_index,draw_t)
% Y) A9 _1 b$ v- }$ V0 Vfor ii = 1:1:numel(draw_y); v7 e/ _. X: o
for jj = 1:1:numel(draw_t)
. ?1 P2 G# f* G2 [ Z(jj,z_index,ii,x_index) = f(x_index,draw_y(ii),z_index,draw_t(jj));
6 S+ z4 E1 b9 f5 f0 M end
6 w3 \& E* P$ _/ nend
! {% E8 O- l6 [% a5 U4 I& x[draw_y,draw_t] = meshgrid(draw_y,draw_t);$ X4 y9 o6 V$ U/ J. z
contour (draw_y,draw_t,Z);
$ y8 f# o2 o& r7 L3 G
; K) h, K" e3 k, i J错误显示 试图访问 f(120,30,7,1)；由于 size(f)=[1,31,1,854]，索引超出范围。8 c" f o" d6 Y+ P
改了循环里面之后会显示 试图访问 f(1,32,1,1)；由于 size(f)=[1,31,1,854]，索引超出范围。
& K) q3 k# d4 O. U. p' T请教大神，非常感谢！. b" R8 B! q$ L9 V+ ^: [
|
|
/1 
发表于 2020-6-2 16:07
收藏
淘帖
支持！
反对！