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

图片赏析

前言

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

数据来源：

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

7 O2 Y' {1 l5 [$ m

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

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

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

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

空间精度：0.25°

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

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

7 s- c1 v( h/ d$ c) M" M. \

源代码

, I* P5 X% O: M3 b7 \

%南海1998年12月31日某时的海面热通量数据及分布图 7 k) `# ]5 v8 D %% ) R0 B# ^& \# y4 x2 v* D$ C- w %第一步：数据的展示 tic;Path="E:\OneDrive - Inc\Share_science\NC文件\";%文件夹地址 " K: u! S4 @4 i" p- x* e! q filename="SEAFLUX-OSB-CDR_V02R00_FLUX_D19981231_C20160824.nc";%文件名 ' g/ z* e$ x+ h" [2 }1 ~& | abs_path=fullfile(Path,filename);%合并，绝对路径文件 7 M9 [- p- q8 \5 P clearvars -except abs_path;%清除除了abs_path以外所有变量 6 y7 ?( Z6 t5 }9 y ncdisp(abs_path);%读取所选文件 4 H0 H* s/ n. v: ?( M9 K whos;toc; %% 4 R$ U3 ~$ V7 X x %第二步：数据读取 1 a, `' m& P, J! |6 _ tic;Lat_ini=ncread(abs_path,lat); Lon_ini=ncread(abs_path,lon); Time_ini=ncread(abs_path,time); Lflx_ini=ncread(abs_path,surface_upward_latent_heat_flux); 8 ]+ R5 M* p3 g Sflx_ini=ncread(abs_path,surface_upward_sensible_heat_flux); E4 G2 J3 t }/ p8 Y toc;clear abs_path; / \) X6 {; c2 R5 W! W7 i H$ a %% %第三步：确定边界/截取数据 %空间边界：南海范围北纬2°至北纬25° 东经104°至东经124° tic;Boundary=[104,124;2,25;]%修改范围[经度;纬度],注意顺序 Spatial_Resolution=abs(Lat_ini(1)-Lat_ini(2));Sp_Re=Spatial_Resolution;%自动确定精度 5 V! G$ ]* N5 b1 P3 O0 }- n4 Q3 ] Lon_start_No=find(Lon_ini-(Boundary(1))<=Sp_Re&Lon_ini-(Boundary(1))>=0);%找到离104°最近格点的编号（行号） 5 A9 S5 R h/ ~* e# d Lon_end_No=find(Lon_ini-(Boundary(3))<=Sp_Re&Lon_ini-(Boundary(3))>=0);%找到离124°最近格点的编号（行号） Lat_start_No=find(Lat_ini-(Boundary(2))<=Sp_Re&Lat_ini-(Boundary(2))>=0);%找到离2°最近格点的编号（行号） 9 @ {3 l, b! d0 x4 K Lat_end_No=find(Lat_ini-(Boundary(4))<=Sp_Re&Lat_ini-(Boundary(4))>=0);%找到离25°最近格点的编号（行号） 0 \4 l- u; }+ k$ j$ G% @& w Lon=Lon_ini(Lon_start_No:Lon_end_No); : H1 N& M7 v: \ Lat=Lat_ini(Lat_start_No:Lat_end_No);%北半球可能要颠倒一下编号，因为是从北极往赤道计数的 clear Lat_ini;clear Lon_ini;%回收内存 # a Y" {/ I8 F+ I; Q Lat=double(Lat);Lon=double(Lon);%双精度化 latlim=[min(Lat),max(Lat)];lonlim=[min(Lon),max(Lon)]; * I6 W; I5 V" R6 U+ w1 o3 ~ Time=9%输入时间，注意三小时的倍数 , V8 G3 `3 M4 r; T Time_No=Time/3+1; n. Y/ Q3 H7 ^ %截取数据并降维度 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));;%北半球有时可能要颠倒一下编号，因为是从北极往赤道计数的 clear Sp_Re;toc;%清除不再使用的变量 %% %第四步：绘图 tic;figure; m_proj(mercator,lat,latlim,lon,lonlim); [longrid,latgrid]=meshgrid(Lon,Lat); 4 O+ K! N0 X/ v% Z2 u/ B) v" g m_pcolor(longrid,latgrid,Sflx); %pcolor或contourf ) |( T. t* a+ A, W9 [5 M0 Z1 S tic;m_gshhs_f(patch,k);toc;%添加海岸线：c<l<h<f : [# b+ x9 e, k% s: J2 k) e m_grid(box,fancy,tickdir,in);%网格化 %brighten(.5);%亮化 1 s: d& j' Y8 z+ l" w0 f$ q) L8 } colormap(hsv); h = colorbar(h);%色标 h.Label.String = W·m^{-2} h.Location = eastoutside;%色标位置 1 |. ?2 h9 \5 c title([1998年12月31日,num2str(Time),:00,南海海面感热通量],... 0 x w7 T3 n) g+ W& \+ H FontName,黑体,fontsize,12,FontWeight,"bold") %标题 xlabel(Longitude);ylabel(Latitude) ;set(gca,fontsize,9); toc;

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

3 F/ g9 [8 X" g& f5 h$ Z

& Y8 E+ G# M; i- g, ]: n ; J5 I F1 r( N& v