Matlab海洋绘图：海面热通量分布——以南海为例 - 海洋地质热流测量

/ j! d c( i6 c/ N5 {9 I& p* j9 f

图片赏析

& ~+ ^6 w9 b% s$ x) j

* i+ i; q* E$ g( X+ E

前言

5 O/ T; d2 h0 a' z* E& }

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

数据来源：

& f; x; U( P( }7 u$ h

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

( F' z1 P3 ~% g/ E. }! T

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

5 T# x4 b! @8 j& q7 |

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

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

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

空间精度：0.25°

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

/ D1 k3 \$ D8 J" V, C3 _$ W# L

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

0 V" u/ m) j& T. }

源代码

. j( ~9 G- o% C7 q' k/ F7 }, N

%南海1998年12月31日某时的海面热通量数据及分布图 ; W9 l3 X. f/ r5 C %% %第一步：数据的展示 tic;Path="E:\OneDrive - Inc\Share_science\NC文件\";%文件夹地址 filename="SEAFLUX-OSB-CDR_V02R00_FLUX_D19981231_C20160824.nc";%文件名 abs_path=fullfile(Path,filename);%合并，绝对路径文件 ( e- q! p* z& ]% F" d clearvars -except abs_path;%清除除了abs_path以外所有变量 ( N# z2 C0 d$ F" j) f( o ncdisp(abs_path);%读取所选文件 7 i" V* U. `1 \. E: ~6 E whos;toc; - X+ H w! }9 P! M4 ] %% * L H* V: [5 \1 |+ B4 s %第二步：数据读取 4 W; N$ Q& [, m" Y6 R3 h P8 V4 } tic;Lat_ini=ncread(abs_path,lat); Lon_ini=ncread(abs_path,lon); * d" P6 B0 p$ [ Time_ini=ncread(abs_path,time); 8 _4 O& u* ?! R3 e4 \0 V: a Lflx_ini=ncread(abs_path,surface_upward_latent_heat_flux); 6 s: {: ?) Y& i: d' {! I5 s Sflx_ini=ncread(abs_path,surface_upward_sensible_heat_flux); 3 _3 j3 {0 ]0 j) t, \ toc;clear abs_path; %% %第三步：确定边界/截取数据 $ U) @/ R* S4 a %空间边界：南海范围北纬2°至北纬25° 东经104°至东经124° 0 r N8 e2 D8 |5 B- ^" E tic;Boundary=[104,124;2,25;]%修改范围[经度;纬度],注意顺序 7 K$ s. ^1 W4 r+ P5 x. ^3 o h Spatial_Resolution=abs(Lat_ini(1)-Lat_ini(2));Sp_Re=Spatial_Resolution;%自动确定精度 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°最近格点的编号（行号） Lat_start_No=find(Lat_ini-(Boundary(2))<=Sp_Re&Lat_ini-(Boundary(2))>=0);%找到离2°最近格点的编号（行号） ; C" A1 w [/ d$ {2 w 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);%北半球可能要颠倒一下编号，因为是从北极往赤道计数的 7 ]) b2 J" g# E. { clear Lat_ini;clear Lon_ini;%回收内存 : {% I% n: _8 r1 I" e2 A Lat=double(Lat);Lon=double(Lon);%双精度化 latlim=[min(Lat),max(Lat)];lonlim=[min(Lon),max(Lon)]; - @: H. Y3 e* u a, [8 K X- n Time=9%输入时间，注意三小时的倍数 Time_No=Time/3+1; %截取数据并降维度 Lflx=squeeze(Lflx_ini(Lon_start_No:Lon_end_No,Lat_start_No:Lat_end_No,Time_No));;%北半球有时可能要颠倒一下编号，因为是从北极往赤道计数的 & n/ i& a# a4 Z Sflx=squeeze(Sflx_ini(Lon_start_No:Lon_end_No,Lat_start_No:Lat_end_No,Time_No));;%北半球有时可能要颠倒一下编号，因为是从北极往赤道计数的 clear Sp_Re;toc;%清除不再使用的变量 %% * l" Q1 ^8 D7 l6 U" a! W, e) v/ M %第四步：绘图 * |; L( f& a, w: j" Z0 }" N2 l6 ]6 R6 P tic;figure; % i3 P2 ]- e( O7 }* G m_proj(mercator,lat,latlim,lon,lonlim); [longrid,latgrid]=meshgrid(Lon,Lat); m_pcolor(longrid,latgrid,Sflx); %pcolor或contourf tic;m_gshhs_f(patch,k);toc;%添加海岸线：c<l<h<f m_grid(box,fancy,tickdir,in);%网格化 %brighten(.5);%亮化 . r0 U( l( t/ G0 x8 A& R, E colormap(hsv); h = colorbar(h);%色标 h.Label.String = W·m^{-2} 9 Z6 L! y5 I) ~- ]3 }# |2 v h.Location = eastoutside;%色标位置 3 @2 r2 ]. H: ]) O$ h title([1998年12月31日,num2str(Time),:00,南海海面感热通量],... & O H3 g2 X T; ?0 }: Z FontName,黑体,fontsize,12,FontWeight,"bold") %标题 . C# d H3 @+ C$ @ xlabel(Longitude);ylabel(Latitude) ;set(gca,fontsize,9); / J% r; O* L8 o8 e* H toc;

# r) D) _4 b3 K/ `- v, _" k ^6 l

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

) Z4 ~* }' v! k1 h/ }3 X2 ^; I+ T

) l; s+ f {8 }* z6 u 4 M+ J. _4 z+ e/ s. x8 K, ^ * ?' a9 ~" D: F0 P' Q ! I A- r4 t+ o3 X ' X- u1 ~. @; M2 o