matlab ncep,Matlab全球风场图—以ERA5、NCEP和ECMWF再分析数据为例

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: