|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
clear% \$ a4 C5 o; M ?/ c. o
ncinfo('rhum.mon.mean.nc');0 C% O$ D2 |: z# C; s6 L/ V; y' J
y_len=60;
2 M/ } f/ }) Q+ {t_len=854;
. N/ s i+ j# K1 P! Wz_index=7$ X1 n4 r T' F" G9 T
x_index=120
. Z j. y. _# v( x7 wx = ncread('rhum.mon.mean.nc','lon');%x
0 C4 f4 z& `3 d- R0 E0 Fz = ncread('rhum.mon.mean.nc', 'level');%z
& q- ^2 L, J' d+ p% U& y+ l0 ^y = ncread('rhum.mon.mean.nc', 'lat');%y
9 [1 Z$ V" |+ ^t = ncread('rhum.mon.mean.nc', 'time');%t0 q! X) X4 J% j( @3 h
rhum = ncread('rhum.mon.mean.nc', 'rhum');
9 w# c. B+ Z! k: W9 b" G; k3 `figure; hold on;6 P% f q. W( P9 }
draw_y = 30:1:60;
; h$ R# [6 |3 a" m: Bdraw_t = 1:1:854;9 T4 U* u I" F4 o# R
f = rhum(x_index,draw_y,z_index,draw_t)( k. K; b7 U$ |7 C# U
for ii = 1:1:numel(draw_y)) Y. ~1 \$ }5 N% S' V
for jj = 1:1:numel(draw_t)- |' h5 [& @; @/ s- V+ k, s( g, `
Z(jj,z_index,ii,x_index) = f(x_index,draw_y(ii),z_index,draw_t(jj));7 D: l4 B4 O- Q, ^& i
end; O/ \ p* V [7 @) U9 a8 x1 h# H
end! g0 B6 U$ v% f9 p% P5 e# J5 {6 s
[draw_y,draw_t] = meshgrid(draw_y,draw_t);; [/ P& y& p4 P7 I# `/ x
contour (draw_y,draw_t,Z);: O, p6 `( N* F) \% R
5 v/ t1 }" L& N+ w8 H5 E
错误显示 试图访问 f(120,30,7,1);由于 size(f)=[1,31,1,854],索引超出范围。
" n' n+ q# m0 h# p& V改了循环里面之后会显示 试图访问 f(1,32,1,1);由于 size(f)=[1,31,1,854],索引超出范围。* l! B; e. v' h% H6 N
请教大神,非常感谢!
8 ^3 R- v( C( }% f4 C3 v |
|