/ t" p7 ~2 q4 t3 M' s 图片赏析
( N8 X) S: o3 w; x
* R8 Q9 j i* r! W8 U: D# q 2 k D* }% O" G! k- @" Q; k `: ]
2 J5 p; n- R! ?6 z1 M# {! w: K: V 前言 , v( c% o$ v+ ~7 T
7 Z7 K1 f5 {+ g. X+ q' g6 C 南海的经纬度：南海范围北纬2°至北纬25° 东经104°至东经124° & E# Z C) P( N* j7 X, A) T
数据来源：
& ? r! a. j1 E0 ] ' i2 K9 |" f' _% x, |* i5 E5 ]
数据链接提供者：广东海洋大学海洋与气象学院海洋科学系纪同学 " u- g# T' o3 z- P
数据来源机构：NOAA美国国家海洋与大气管理局-NCDC环境信息国家中心
8 S* p+ R* ^$ l# C0 l 产品名称：Ocean Heat Fluxes海洋热通量 : L" `" M L$ m/ N4 E
时间精度：每日3-hourly记录 : r! q, C( G% {( J; `" X
空间精度：0.25° & |- M7 A; h( d7 ~( P& t8 \& Z% C
注：本程序需调用UBC英属哥伦比亚大学Rich Pawlowicz 开发的程序包m_map； 5 Q0 w" c3 ]) _& o
亦需调用gshhs高精度海岸线数据。
( v4 o6 E% i& q; Z9 d& ?/ Q 源代码
0 C! `, j8 W; {. T9 q %南海1998年12月31日某时的海面热通量数据及分布图
7 o ~0 E$ w$ R0 Y %% |; ~% a* o1 p9 M1 Z! {* ~* ~
%第一步：数据的展示/ r7 x r* U" W) M2 F
tic;Path="E:\OneDrive - Inc\Share_science\NC文件\";%文件夹地址* p3 _) O2 ~2 }
filename="SEAFLUX-OSB-CDR_V02R00_FLUX_D19981231_C20160824.nc";%文件名
5 q: T7 E& C# d: ]8 y- T" g/ `/ h abs_path=fullfile(Path,filename);%合并，绝对路径文件
4 `4 q5 `6 Q6 {3 p3 f clearvars -except abs_path;%清除除了abs_path以外所有变量
" |' m' |$ |9 o6 \, ^# m9 ?0 M! ^! S1 m ncdisp(abs_path);%读取所选文件
8 O, M$ ~+ `, w whos;toc;
. N% M) k0 T- h' N* D' S' D %%
" _9 \& ~' I( j0 C %第二步：数据读取
1 y' T* ~) p# ]+ N8 U7 ^ tic;Lat_ini=ncread(abs_path,lat);
3 o. @$ z; X) Q3 p+ s Lon_ini=ncread(abs_path,lon);7 a! ]7 [5 U Q) ]5 g4 r
Time_ini=ncread(abs_path,time);$ z# g( y( [1 ?9 _6 q
Lflx_ini=ncread(abs_path,surface_upward_latent_heat_flux);; G, t% K9 r$ ` R
Sflx_ini=ncread(abs_path,surface_upward_sensible_heat_flux);
, M8 X- e6 h! U( @ toc;clear abs_path;, q$ C- ]( q- K# d; ?. c
%%1 V' }* h" c# g s
%第三步：确定边界/截取数据* J: I, G' d4 A4 ?
%空间边界：南海范围北纬2°至北纬25° 东经104°至东经124°1 }& \! s& m$ |5 r" j) ^- J- y$ H- r
tic;Boundary=[104,124;2,25;]%修改范围[经度;纬度],注意顺序
H% j+ m" Q' ? a Spatial_Resolution=abs(Lat_ini(1)-Lat_ini(2));Sp_Re=Spatial_Resolution;%自动确定精度
9 y6 M+ @& U2 @" u8 S& d Lon_start_No=find(Lon_ini-(Boundary(1))<=Sp_Re&Lon_ini-(Boundary(1))>=0);%找到离104°最近格点的编号（行号）
' t8 Q/ s5 Y5 C7 a7 J6 r Lon_end_No=find(Lon_ini-(Boundary(3))<=Sp_Re&Lon_ini-(Boundary(3))>=0);%找到离124°最近格点的编号（行号）2 M* Z; p3 N, D/ f- J, S2 M
Lat_start_No=find(Lat_ini-(Boundary(2))<=Sp_Re&Lat_ini-(Boundary(2))>=0);%找到离2°最近格点的编号（行号）- Z1 f1 G1 Z R+ ?" U9 P n5 p
Lat_end_No=find(Lat_ini-(Boundary(4))<=Sp_Re&Lat_ini-(Boundary(4))>=0);%找到离25°最近格点的编号（行号）
) S k. Y( R3 b) d0 Q+ x! ~1 i Lon=Lon_ini(Lon_start_No:Lon_end_No);* z" d* [* c, {& E( J
Lat=Lat_ini(Lat_start_No:Lat_end_No);%北半球可能要颠倒一下编号，因为是从北极往赤道计数的2 r1 O' H4 S4 ^2 k( v; K
clear Lat_ini;clear Lon_ini;%回收内存
* b3 c6 x+ ` Z' W# A5 ] Lat=double(Lat);Lon=double(Lon);%双精度化/ ]- r, a; z& e* ?+ I2 l
latlim=[min(Lat),max(Lat)];lonlim=[min(Lon),max(Lon)];$ _0 g4 k8 C" z% s, g/ C
Time=9%输入时间，注意三小时的倍数
9 m4 h2 `' d9 v4 Q3 Y1 o- \ Time_No=Time/3+1;
4 ]5 c3 t7 k8 W& L0 T; n %截取数据并降维度
8 o3 d/ m6 U$ q# ?) E Lflx=squeeze(Lflx_ini(Lon_start_No:Lon_end_No,Lat_start_No:Lat_end_No,Time_No));;%北半球有时可能要颠倒一下编号，因为是从北极往赤道计数的3 w8 y* }- w) n! F% W
Sflx=squeeze(Sflx_ini(Lon_start_No:Lon_end_No,Lat_start_No:Lat_end_No,Time_No));;%北半球有时可能要颠倒一下编号，因为是从北极往赤道计数的/ @7 [& O6 Z. M o( S" Z3 m
clear Sp_Re;toc;%清除不再使用的变量
]( p9 X) H! M. O7 M+ p %%: t- a. q, v3 F/ C7 a" `6 n' K2 m
%第四步：绘图- E2 s8 N) \5 f: _' }2 h. } j
tic;figure;& [8 z' G: D6 b" w! |* t
m_proj(mercator,lat,latlim,lon,lonlim);; g1 O3 E' l1 U( Q1 j
[longrid,latgrid]=meshgrid(Lon,Lat);; q8 D+ b5 e+ w4 e2 z
m_pcolor(longrid,latgrid,Sflx); %pcolor或contourf
2 i6 d0 g+ Y5 B+ V" i tic;m_gshhs_f(patch,k);toc;%添加海岸线：c<l<h<f
^: o- I1 L6 |: g m_grid(box,fancy,tickdir,in);%网格化 G$ O7 p# h) m
%brighten(.5);%亮化
1 q9 [0 E) Y# O5 I+ d% p& z6 x3 M colormap(hsv);* h/ ]5 b; d, z* i( ]3 o- h8 y
h = colorbar(h);%色标 B s! |# G: i% k* J2 b
h.Label.String = W·m^{-2}) ~$ S" R) r6 R' \$ W& Y, m0 K, x
h.Location = eastoutside;%色标位置- h, R6 b0 A) y/ V% l
title([1998年12月31日,num2str(Time),:00,南海海面感热通量],...1 v4 V" ]; M, U- N2 r* ]
FontName,黑体,fontsize,12,FontWeight,"bold") %标题. W8 L( m1 U6 x
xlabel(Longitude);ylabel(Latitude) ;set(gca,fontsize,9);4 v# Q# B* c. j
toc;: W$ E1 s) r o* j% \
: ]3 o0 r! Q3 ?/ \3 F ] 附上另外俩残次品绘图结果，数据来源：IOCADS 1 M: t* l- `* z. r3 E$ D" N' X+ i, K( C
- W0 c3 M% A* x z8 v
, k- i$ G- I" a9 E; C
. V' d% Q& m |- @5 H" l- L- A+ h9 o9 e, y$ g5 }2 g
# H/ ?4 g' E: _. Q
7 ?' M9 ~; v! M
|
