Matlab海洋绘图：海面热通量分布——以南海为例

7 K; W0 j2 n7 W, K4 O* J

图片赏析

& W+ J% p" { s) t2 \7 F& T: |

, o) O: y' z4 J% L9 ]4 g( ~9 s/ K

前言

7 w0 \/ E8 N) d/ s+ U/ D

南海的经纬度：南海范围北纬2°至北纬25° 东经104°至东经124°

数据来源：

% `$ |8 R: L& ^. n- E+ `$ Y! `7 c) ^+ H

National Climatic Data Center​www.ncdc.noaa.gov/cdr/atmospheric/ocean-heat-fluxes

8 W( }0 u3 P6 ?7 [; Q

数据链接提供者：广东海洋大学海洋与气象学院海洋科学系纪同学

% Z% D; c5 C" \8 ?1 V

数据来源机构：NOAA美国国家海洋与大气管理局-NCDC环境信息国家中心

产品名称：Ocean Heat Fluxes海洋热通量

, Z# r+ L% c. _* o

时间精度：每日3-hourly记录

空间精度：0.25°

. o; E: \* [$ F

注：本程序需调用UBC英属哥伦比亚大学Rich Pawlowicz 开发的程序包m_map；

亦需调用gshhs高精度海岸线数据。

源代码

%南海1998年12月31日某时的海面热通量数据及分布图 %% 3 L/ j0 x- o# s- H %第一步：数据的展示 tic;Path="E:\OneDrive - Inc\Share_science\NC文件\";%文件夹地址 ' ^# f$ A9 { u# y+ j filename="SEAFLUX-OSB-CDR_V02R00_FLUX_D19981231_C20160824.nc";%文件名 8 T g5 K/ ~# R4 {; K2 r7 P) k abs_path=fullfile(Path,filename);%合并，绝对路径文件 ' p6 U5 c) X0 f. K5 x clearvars -except abs_path;%清除除了abs_path以外所有变量 ; v; {' L, R. A ncdisp(abs_path);%读取所选文件 * m: {6 f7 {0 a* k; ~- S/ h) ` whos;toc; ) J" Z* q! n. A( t, t6 V4 ^+ Z %% 8 s, {9 a8 T8 H! ? %第二步：数据读取 8 k* z& ?* V6 x* i# i, Q tic;Lat_ini=ncread(abs_path,lat); & F6 p; ~% | {" M ?- T, s# B Lon_ini=ncread(abs_path,lon); $ q: t& _/ h3 K N Time_ini=ncread(abs_path,time); , O, Q2 @8 @+ M! f9 O6 h3 _; [2 U Lflx_ini=ncread(abs_path,surface_upward_latent_heat_flux); 2 h! Q( b* I5 B" e, v+ `, c Sflx_ini=ncread(abs_path,surface_upward_sensible_heat_flux); + i, A5 I Z+ G% u/ ^0 w: c toc;clear abs_path; + V4 q( ]! X. u+ [4 Z' { %% %第三步：确定边界/截取数据 ; ~, t2 x3 P* z* o1 n$ d %空间边界：南海范围北纬2°至北纬25° 东经104°至东经124° tic;Boundary=[104,124;2,25;]%修改范围[经度;纬度],注意顺序 Spatial_Resolution=abs(Lat_ini(1)-Lat_ini(2));Sp_Re=Spatial_Resolution;%自动确定精度 Y, o$ J* Y; M {1 P Lon_start_No=find(Lon_ini-(Boundary(1))<=Sp_Re&Lon_ini-(Boundary(1))>=0);%找到离104°最近格点的编号（行号） Lon_end_No=find(Lon_ini-(Boundary(3))<=Sp_Re&Lon_ini-(Boundary(3))>=0);%找到离124°最近格点的编号（行号） 5 d1 {5 k& K q0 ^- Q: b2 @ Lat_start_No=find(Lat_ini-(Boundary(2))<=Sp_Re&Lat_ini-(Boundary(2))>=0);%找到离2°最近格点的编号（行号） ; [# m! S) E& T! d Lat_end_No=find(Lat_ini-(Boundary(4))<=Sp_Re&Lat_ini-(Boundary(4))>=0);%找到离25°最近格点的编号（行号） Lon=Lon_ini(Lon_start_No:Lon_end_No); Lat=Lat_ini(Lat_start_No:Lat_end_No);%北半球可能要颠倒一下编号，因为是从北极往赤道计数的 clear Lat_ini;clear Lon_ini;%回收内存 1 ]$ E+ G% B6 w2 q Lat=double(Lat);Lon=double(Lon);%双精度化 ) N8 I: }2 E0 Y$ q$ {- s latlim=[min(Lat),max(Lat)];lonlim=[min(Lon),max(Lon)]; ; `" f. ~! \9 u+ f+ T% O% d Time=9%输入时间，注意三小时的倍数 Time_No=Time/3+1; 9 C' @ J l; x) t. m %截取数据并降维度 Lflx=squeeze(Lflx_ini(Lon_start_No:Lon_end_No,Lat_start_No:Lat_end_No,Time_No));;%北半球有时可能要颠倒一下编号，因为是从北极往赤道计数的 Sflx=squeeze(Sflx_ini(Lon_start_No:Lon_end_No,Lat_start_No:Lat_end_No,Time_No));;%北半球有时可能要颠倒一下编号，因为是从北极往赤道计数的 8 |2 c: o' Q! |- U/ r clear Sp_Re;toc;%清除不再使用的变量 ) Z0 k9 U+ E2 C1 Y. ^2 T %% # G" H. O. ^& r3 a. i %第四步：绘图 tic;figure; . c, l7 }5 C0 j8 o m_proj(mercator,lat,latlim,lon,lonlim); [longrid,latgrid]=meshgrid(Lon,Lat); m_pcolor(longrid,latgrid,Sflx); %pcolor或contourf $ i, i! V+ W6 |1 O* Q, r tic;m_gshhs_f(patch,k);toc;%添加海岸线：c<l<h<f m_grid(box,fancy,tickdir,in);%网格化 %brighten(.5);%亮化 $ d. T8 A+ M! [8 W# N+ K colormap(hsv); h = colorbar(h);%色标 h.Label.String = W·m^{-2} 2 G/ g9 p) h' ^2 p h.Location = eastoutside;%色标位置 title([1998年12月31日,num2str(Time),:00,南海海面感热通量],... FontName,黑体,fontsize,12,FontWeight,"bold") %标题 ( X. f" C* q5 \( x/ }* m+ c, q xlabel(Longitude);ylabel(Latitude) ;set(gca,fontsize,9); toc; 3 k) F% w, M. m# ~( E* X- s. u3 K

+ _, p. J6 |" G

附上另外俩残次品绘图结果，数据来源：IOCADS

* Q: Q8 A' y3 h2 S% p% M& u, @

6 N- j+ r7 X, t( {7 H. O

7 b- P6 ?( z% t3 Y5 ]8 N 9 b; H# l0 ^0 V: ~' @" o9 o