5 [( O8 p( Q/ l" k, P& `
图片赏析
. @+ N; U9 J* D% @# H# y5 b0 `5 `
3 H7 @/ S) V* s3 x$ X& S* k
/ i/ e, Z: A2 w; l
# u; \& \0 ]2 l- j! w1 V# o8 F 前言 ; h3 L3 {0 S: C
& N' }/ x3 U8 n9 S 南海的经纬度：南海范围北纬2°至北纬25° 东经104°至东经124° 5 R {; A2 z9 F9 D) ^
数据来源： 8 q3 B/ L8 L) h6 t7 g, N
' U& C6 K& Y3 w! t. |7 ^
数据链接提供者：广东海洋大学海洋与气象学院海洋科学系纪同学 + e7 r1 \/ m2 U- `( I; O
数据来源机构：NOAA美国国家海洋与大气管理局-NCDC环境信息国家中心 & u: p' a3 g7 l" t( S7 J
产品名称：Ocean Heat Fluxes海洋热通量 7 @1 v: U9 Z. B* {
时间精度：每日3-hourly记录 % j+ \+ [0 I$ G
空间精度：0.25° 2 b/ q# @* }/ ?% o( c# w
注：本程序需调用UBC英属哥伦比亚大学Rich Pawlowicz 开发的程序包m_map；
# h! h5 M) ?/ ]2 b0 O7 W9 ^ 亦需调用gshhs高精度海岸线数据。 ! R3 i- f O7 x2 }
源代码4 V X/ M2 ]2 T( `2 v
%南海1998年12月31日某时的海面热通量数据及分布图4 ~( e8 C7 z- [+ h
%%
+ z. m' m+ }# N$ e" n %第一步：数据的展示
$ n2 q4 d( a7 v4 _: j tic;Path="E:\OneDrive - Inc\Share_science\NC文件\";%文件夹地址$ H4 M" Y4 L( |, D% S; T6 K. G+ U. t
filename="SEAFLUX-OSB-CDR_V02R00_FLUX_D19981231_C20160824.nc";%文件名
U2 ^ E }9 B! \0 j& K abs_path=fullfile(Path,filename);%合并，绝对路径文件* r: B/ x5 o8 ~ y7 F+ [
clearvars -except abs_path;%清除除了abs_path以外所有变量
6 i" |' X# ^# |' P ncdisp(abs_path);%读取所选文件# h* y4 C8 A1 e9 ?" F
whos;toc;( t9 J: ~2 r6 v0 `& b z# h" f6 H
%%
8 K, U' o! x* A3 N% A0 ` %第二步：数据读取8 k9 P- y) X4 ~) l
tic;Lat_ini=ncread(abs_path,lat);
3 O \) ?: A; g# J' B Lon_ini=ncread(abs_path,lon);
y: |- {$ U; ~: \ Time_ini=ncread(abs_path,time);
5 Z5 Q1 O5 J+ l) [+ _8 i" ] Lflx_ini=ncread(abs_path,surface_upward_latent_heat_flux);
, ~; [. f; `1 B0 O Sflx_ini=ncread(abs_path,surface_upward_sensible_heat_flux);4 b4 |& @2 B* Y
toc;clear abs_path;( `' Z" D' P7 {/ I: X2 a" F' }
%%# k4 F) S2 Z V1 ^- ?
%第三步：确定边界/截取数据% O+ k% F4 ?! @: ?8 y' \
%空间边界：南海范围北纬2°至北纬25° 东经104°至东经124°
1 J2 {4 j: X! \. [ tic;Boundary=[104,124;2,25;]%修改范围[经度;纬度],注意顺序) J% ]4 f D' r O4 V& k
Spatial_Resolution=abs(Lat_ini(1)-Lat_ini(2));Sp_Re=Spatial_Resolution;%自动确定精度/ \: e/ f! r1 L4 k% }8 {
Lon_start_No=find(Lon_ini-(Boundary(1))<=Sp_Re&Lon_ini-(Boundary(1))>=0);%找到离104°最近格点的编号（行号）
4 k! ?( ]: g" J9 P Lon_end_No=find(Lon_ini-(Boundary(3))<=Sp_Re&Lon_ini-(Boundary(3))>=0);%找到离124°最近格点的编号（行号）
4 X0 q f& q1 Y9 O Lat_start_No=find(Lat_ini-(Boundary(2))<=Sp_Re&Lat_ini-(Boundary(2))>=0);%找到离2°最近格点的编号（行号）& G4 h& Y# K/ e0 g$ V
Lat_end_No=find(Lat_ini-(Boundary(4))<=Sp_Re&Lat_ini-(Boundary(4))>=0);%找到离25°最近格点的编号（行号）; b* ?5 l3 ^/ h) J; d8 l- h
Lon=Lon_ini(Lon_start_No:Lon_end_No);- C( y5 a3 P' Z6 E+ v( E
Lat=Lat_ini(Lat_start_No:Lat_end_No);%北半球可能要颠倒一下编号，因为是从北极往赤道计数的
" ^% I* Q" y2 { clear Lat_ini;clear Lon_ini;%回收内存3 A+ p6 z: X$ \! L+ k& M
Lat=double(Lat);Lon=double(Lon);%双精度化6 h, u2 n/ C X' F5 _8 w
latlim=[min(Lat),max(Lat)];lonlim=[min(Lon),max(Lon)];1 Z' b6 Y( p: m' t" f3 m) b
Time=9%输入时间，注意三小时的倍数
1 W/ C2 R$ U; o) G3 e& q Time_No=Time/3+1;
: J+ m% m/ e- h: U4 ` %截取数据并降维度; {8 n+ N3 B& b) C P' \ _
Lflx=squeeze(Lflx_ini(Lon_start_No:Lon_end_No,Lat_start_No:Lat_end_No,Time_No));;%北半球有时可能要颠倒一下编号，因为是从北极往赤道计数的6 j* H' k* |6 _6 n
Sflx=squeeze(Sflx_ini(Lon_start_No:Lon_end_No,Lat_start_No:Lat_end_No,Time_No));;%北半球有时可能要颠倒一下编号，因为是从北极往赤道计数的% m* w0 a5 T- s
clear Sp_Re;toc;%清除不再使用的变量: a9 U- V, W% _! Y+ @
%%
; v9 h& j- |8 @ %第四步：绘图7 u$ \/ `: Z7 A1 X
tic;figure; r: Q! f' W; v
m_proj(mercator,lat,latlim,lon,lonlim);# A) R# S; ]( j$ D
[longrid,latgrid]=meshgrid(Lon,Lat);
: ^7 F- `! \7 C. O m_pcolor(longrid,latgrid,Sflx); %pcolor或contourf
( c8 ~, K$ i2 n% O tic;m_gshhs_f(patch,k);toc;%添加海岸线：c<l<h<f
5 h( h3 ^6 L/ |% {9 c7 Z K m_grid(box,fancy,tickdir,in);%网格化
: a9 `8 G- e$ q# i5 W4 b %brighten(.5);%亮化4 e1 u2 F, \9 K$ V g P% t7 k
colormap(hsv);( k/ v/ r0 H% n; A% ]
h = colorbar(h);%色标
' |: q' w5 T, ]8 ~! M) D h.Label.String = W·m^{-2}# }: _6 ^( e% X' v9 M/ B: V/ s
h.Location = eastoutside;%色标位置
% r/ i# T2 j& J. h+ ]- W title([1998年12月31日,num2str(Time),:00,南海海面感热通量],...3 z$ S8 J. P) C" C2 T3 T8 F N
FontName,黑体,fontsize,12,FontWeight,"bold") %标题
5 ~7 j8 Q0 M: C, d xlabel(Longitude);ylabel(Latitude) ;set(gca,fontsize,9);2 Z( S( C8 Z) ]2 ^' @4 V
toc; }+ }9 q# u D+ W3 d
6 N5 \- Q3 z8 X% x 附上另外俩残次品绘图结果，数据来源：IOCADS
8 k) G. t9 ~1 ^% v1 I& n: v 1 r' `8 ^) k j2 u$ h( B
$ F1 F2 k: q" @5 U2 I
$ S# G8 Y5 _2 C& R5 }3 S
5 m8 |% r; F6 n4 G- x6 N9 v9 I
v% d8 }2 z. E+ Q
' U$ P: U1 n, e/ W | 
