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

& G' C3 t( R2 W. _

图片赏析

2 _6 Z$ _0 k6 w# T

; \8 G+ K% p. i

0 I5 H* h' T/ X. T: s

前言

( G$ p; F$ Q' a8 `( P

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

数据来源：

+ t# U2 M. t8 z! L: {0 M

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

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

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

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

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

空间精度：0.25°

1 a; b0 a$ v( x' t

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

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

6 Z! r0 t5 H5 W& A8 A/ R( v6 c

源代码

! N, ]; G" w' s* q$ ?

%南海1998年12月31日某时的海面热通量数据及分布图 %% %第一步：数据的展示 9 e# f6 N. `. V& p! d3 ]8 p tic;Path="E:\OneDrive - Inc\Share_science\NC文件\";%文件夹地址 filename="SEAFLUX-OSB-CDR_V02R00_FLUX_D19981231_C20160824.nc";%文件名 1 h! s( g+ l' B3 M abs_path=fullfile(Path,filename);%合并，绝对路径文件 7 A9 Q' J; w+ D! o( d clearvars -except abs_path;%清除除了abs_path以外所有变量 1 q6 U9 l! s3 A/ K9 | ncdisp(abs_path);%读取所选文件 7 A1 M5 m' |) M7 ]5 \7 d7 t whos;toc; %% ! h, H2 G; E3 H$ u: W %第二步：数据读取 tic;Lat_ini=ncread(abs_path,lat); Lon_ini=ncread(abs_path,lon); 4 W2 m- j9 d( I: T$ M6 h9 k Time_ini=ncread(abs_path,time); Lflx_ini=ncread(abs_path,surface_upward_latent_heat_flux); Sflx_ini=ncread(abs_path,surface_upward_sensible_heat_flux); toc;clear abs_path; %% - X. u- B6 Q$ K0 e$ ~4 C %第三步：确定边界/截取数据 : k7 t5 q8 j {! b. H2 U3 p' N6 R %空间边界：南海范围北纬2°至北纬25° 东经104°至东经124° / V% ~ ]) x: j/ B! F tic;Boundary=[104,124;2,25;]%修改范围[经度;纬度],注意顺序 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°最近格点的编号（行号） T, P- X! ` k% W% _2 J# `, k Lat_start_No=find(Lat_ini-(Boundary(2))<=Sp_Re&Lat_ini-(Boundary(2))>=0);%找到离2°最近格点的编号（行号） $ g8 x0 h# a' i2 [% h3 K5 [' y Lat_end_No=find(Lat_ini-(Boundary(4))<=Sp_Re&Lat_ini-(Boundary(4))>=0);%找到离25°最近格点的编号（行号） ; ?3 B' h7 g- Q! } Lon=Lon_ini(Lon_start_No:Lon_end_No); , l5 v$ V: N8 G0 W Lat=Lat_ini(Lat_start_No:Lat_end_No);%北半球可能要颠倒一下编号，因为是从北极往赤道计数的 8 O* a: m( I3 m1 W; k clear Lat_ini;clear Lon_ini;%回收内存 / z# X1 B$ `0 B& K4 T+ V& q( H Lat=double(Lat);Lon=double(Lon);%双精度化 9 O- ]! _ `& t5 U3 f0 U: O# W latlim=[min(Lat),max(Lat)];lonlim=[min(Lon),max(Lon)]; 7 P6 e/ B8 g( K$ `6 W L. V Time=9%输入时间，注意三小时的倍数 # _$ I: i- ~0 ` Time_No=Time/3+1; %截取数据并降维度 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;%清除不再使用的变量 %% # U% b$ |) e# @# @" R& I- p %第四步：绘图 tic;figure; ; }( _# n) m6 B2 P6 C) D m_proj(mercator,lat,latlim,lon,lonlim); [longrid,latgrid]=meshgrid(Lon,Lat); 5 m& @; l' Z( c& W' U m_pcolor(longrid,latgrid,Sflx); %pcolor或contourf $ c( E% t' e6 p0 w8 J. J tic;m_gshhs_f(patch,k);toc;%添加海岸线：c<l<h<f / {; S0 q; ]6 | ~ m_grid(box,fancy,tickdir,in);%网格化 %brighten(.5);%亮化 , E! W' P+ V9 Q1 g1 T" w* @ colormap(hsv); 3 f$ H" X: o+ i7 @ h = colorbar(h);%色标 * O7 q" ~. Y% b5 G' f- m/ Z# A9 T+ ^ h.Label.String = W·m^{-2} 9 T2 r8 Q4 i9 [3 [9 P! S7 I r; f h.Location = eastoutside;%色标位置 title([1998年12月31日,num2str(Time),:00,南海海面感热通量],... FontName,黑体,fontsize,12,FontWeight,"bold") %标题 ; X8 A* p1 \; B) g' R2 m xlabel(Longitude);ylabel(Latitude) ;set(gca,fontsize,9); toc; 1 @) @5 P; N) C4 F# Z: @

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

, \ m7 o2 Z# O- }+ C, u+ M

' w4 P: s$ g6 h1 r6 c& _5 L2 j, r8 _$ l 6 v U# r) ^! \: E1 A1 g; a