matlabncep,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_';! 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/ A
6 ]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: [" h: H- Tcount=[41,31,1,1];%从start指定的开始位置算起，⼀共读取的每⼀维要素的数⽬1 |4 B. {! H2 I) }2 S
strip=[1,1,1,1];%从start开始，每⼀维读取的数⽬为count时，每⼀维的读取的步长
$ A6 C1 s( \) @. x8 ehgt=ncread(f_hgt,'z',start,count,strip);%读取温度值，单位K+ G5 Y/ G0 F) j3 K7 K; @; }
u=ncread(f_hgt,'u',start,count,strip);%读取温度值，单位K
/ ~2 v# V9 u4 n; Cv=ncread(f_hgt,'v',start,count,strip);%读取温度值，单位K
& _- k# b* d. v0 S+ u# E[X,Y]=meshgrid(lon,lat);
5 `' d) i2 Y% ]1 Y# Rfigure(1)
1 X; j9 e& b1 l! {1 Ym_proj('Mercator','lat',[25,35],'lon',[100,115]);, ]0 N: f9 d' h; M [
% m_grid('linestyle','none','tickdir','out','fontsize',12,'fontname','Times New Roman');$ i% n o' `, w* O9 j
m_grid('linestyle','none','box','fancy','fontsize',11,'tickdir','in','xtick',[100:3:115],'ytick',[25:2:35]);
. \& V; c" P/ P6 R2 U% e, Ahold on# ]6 m6 F5 `) a- c3 H5 p7 j$ l3 E
m_windbarb(X',Y',u,v,'color','k')' N7 y* g- |8 u' r' [: R) ~
%m_coast('patch',[.9 .9 .9],'edgecolor','none');3 h1 t6 l& n6 B/ a6 R6 ?% \( G
[C,h]=m_contour(X',Y',hgt/10 ,'color','k','LineWidth',1);%[5000:80:5900]%[1520:25:1680]) b* V) U* O, f& H/ T
%set(b,'ShowText','on','TextStep',get(b,'LevelStep'),'fontsize',12,'fontname','Times New Roman'); %
在等⾼线上叠加数值(⽂后详情)
. t0 L: O' ~, l K2 s8 A dclabel(C,h,'FontSize',13,'fontname','Times New Roman');% j; j9 k& {+ j# d
ma=shaperead('F:/RMeteoInfo/data/map/bou2_4l.shp');
( ?- M9 i, ~, h9 ?1 B5 I. T% m_line( [ma(：).X], [ma(：).Y],'color',[0.5,0.5,0.5]);%绘制范围内的地图
加热搅拌反应釜2 [& r. ~7 z
3 g# Z* o& U/ Mm_line([ma(：).X],[ma(：).Y],'color','r');%绘制范围内的地图% H/ }( R" f9 p! c9 \
m_plot(105.5,29.43,'marker','^','MarkerSize',7,'color','k','MarkeRFaceColor','k'). f- B3 P4 q7 q( j! g! }
5 A- Z; n$ W. o$ o
6 O' I) ]. g/ P1 k
3.png (579.26 KB, 下载次数: 0)
2021-3-3 09:46 上传
7 P$ t& m3 Q( j. y" e2 t% U
. ]3 ~; j0 k# I1 \' F* Q, R7 q g: K# ~
⼆、matlab读取ERA5并绘制全球风场图
5 X d. k$ @. K7 d! a
) z$ a( C7 o9 o! d; uclc;clear;close all
9 u( ?3 s( X0 cu0=ncread('','uas'); %读取其中⼀项# K0 L. z3 S- L) T
v0=ncread('','vas');- x* y# G* I1 g0 A3 c
time0=ncread('','time');新型高效台车式退火炉
7 u- ]- B7 b1 {2 f- Glat0=ncread('','lat');3 o9 A4 J1 r* I5 ^5 F
lon0=ncread('','lon');
新金瓶酶23 ~
防盗电子狗
4 y7 T/ \; W9 _) t& X2 y g8 Z[lat,lon]=meshgrid(lat0,lon0);
2 o8 G7 e/ g7 O' {
3 j6 zu=u0(:,:,2);% \. N9 T0 ~, n" T9 s j
v=v0(:,:,2);% h/ M I; t! U& j6 j" A
b=sqrt(u.^2+v.^2);: }& s7 r4 x1 Y \) M
figure(1);
: F5 y' F3 h! Lm_proj('Equidistant Cylindrical','long',[-180 180],'lat',[-90 90]);%矩形投影；取区域观察
% P7 n2 b2 K0 y. W[lon1,lat1]=meshgrid([-180:2.5:179.75],[-90:2.5:90]);
: ^# {. z# O% nu2=u(1:10:end,1:10:end)';
. t9 C* b I, ~( ]3 D$ a, Y9 H% u3=[u2,u2(:,end)];
! R' W& c2 z4 h: mv2=v(1:gap:end,1:gap:end)';. \; \0 ^; C5 x: ~, V" J0 l9 Z0 c
% v3=[v2,v2(:,end)];
2 X Z; N/ b小鼠步态分析实验
3 t7 I
4 j* q
5 X( fm_quiver(lon1,lat1,u2,v2,0);: N2 ~# L' @% x: `
hold off
* { A2 |( P' U0 M( bm_coast('patch',[1 .85 .7]);0 F; X/ p; Q7 K) c- ^
m_grid('box','on','tickdir','out');5 Z/ d O( s: f+ ^9 z
, m1 {4 @- {, H+ X$ P; G
4.png (1.48 MB, 下载次数: 0)
2021-3-3 09:49 上传
' r5 V# Q$ t8 X
4 Q7 R- X) g% X7 n8 W
! G) ?1 [6 d6 _0 _) L2 |
( E1 C# k& G, D. o z' ], |$ Mm_contourf(lon,lat,b);%在mmap基础上的画! e# v6 J* U, U7 S7 P shading interp;%使数据插值! w$ \0 I* H$ L5 g, _3 o
hold on;
3 @0 j9 u ~7 u1 Qm_quiver(lon1,lat1,u2,v2,0);: @9 g9 J! f7 [: M$ n. u* }9 W% n
hold off;: l: w* f; {* `! O
x: `2 d2 D4 }2 ~, J; C" E) x
人工鱼礁5.png (1000.81 KB, 下载次数: 0)
2021-3-3 09:46 上传
: |( J5 Y* s+ N3 j/ A- p# g* }/ m6 S8 c, F9 p. \
4 S9 l; t3 e. s" Q& T# K6 `
' j5 \) N; X* Zm_quiver(lon1,lat1,u2,v2,0);
/ j6 f7 g; ?. q, r/ f' `& z
$ m0 o; B& |* s7 S2 g* {) ^* [! C
6.png (929.8 KB, 下载次数: 0)
2021-3-3 09:46 上传
( s$ A7 o3 c- o) J. z) z5 B
+ C0 N$ j" o# M9 r5 Q5 E: h9 Q
" |) h x: d$ F1 c$ _# ^: }/ o
) `9 w# t' g7 ^9 d! g' @" I s
m_contourf(lon,lat,b,500,'linestyle','none') %在mmap基础上的画6 [* X$ [/ _# P
shading interp; %使数据插值
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议