EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
2 M6 B* g4 Y4 D* \! X8 r& L一、matlab读取NCEP再分析数据并绘制风场
/ U- h7 c) Q% g. [4 Z9 g
) k P( w- N) ?, O" l%该程序用于求水汽通量散度
# z: T5 d6 Y) x3 T9 Y%注意!我们这里读到的u是四维矩阵,分别是lon*lat*level*time,1 e, |( l# I* B5 `: A- p" N6 g
clc;clear;close all
( Q; w5 a5 Y. {# P& S& k7 xf_hgt = 'ps_level_20170121_0130.nc';! k& g! b8 M: ~2 w4 b9 c2 I K, @
% ncdisp(f_hgt);* X8 }% C4 z( U, p& L
time=ncread(f_hgt,'time');/ m" p4 L9 m1 h2 u5 r8 h
level=ncread(f_hgt,'level');+ j0 E$ W D* ]& ~: X- z
lon=ncread(f_hgt,'longitude');8 @0 u2 e( _' d0 |
lat=ncread(f_hgt,'latitude');
1 ]" ]2 d/ z7 u5 s7 e%%%%%%时间转换
4 |' [9 T. j+ T3 Ltime = double(time);
' M4 e1 \ C% z4 } N- H! eformat = 'mm dd, yyyy HH:MM:SS.FFF AM';%转换格式% y# x6 g# Q$ i! h7 ]
dstr = datestr((datenum('1900-01-01') + time./24),format);%转换后时间字符串存储
5 E/ A6 ]5 z8 B+ ATM = datevec(dstr);%将时间字符数组转化为数值数组
6 E- }9 ^+ Z% ?5 Xtidx=find(TM(:,2)==1 & TM(:,3)==28 & TM(:,4)==00);%筛选7月25日08时(世界时加8)/ ~% m. U( r* K( O- V5 G! Y! w
ps_lev=find(level ==850);%%删选出850hPa高度
6 P. H! B: R( U- O! Y( V; wstart=[1,1,ps_lev,tidx];%所指定变量的每一维的开始读取的位置
a+ i* \) G: