4 ^5 N% @7 ~7 u5 q% B. k
图片赏析
* }" c# W+ X' O; G( d7 I$ n/ ]+ q$ y ! s5 z d" w% p: s
7 i/ L# B" v$ e: p
0 W0 x. y: z2 R! t: \. o! A0 W
前言 - j# d. j" L# S i5 C
, s$ Q! {# J; n9 O( l
南海的经纬度:南海范围北纬2°至北纬25° 东经104°至东经124° 9 |7 N9 O# C3 N0 o; y" r/ N: \
数据来源:
& o! h% \5 E" A O) v/ |5 q8 N; e Y0 }( ?8 s& F+ v& T
数据链接提供者:广东海洋大学海洋与气象学院海洋科学系纪同学 5 [: e6 ?1 V7 r$ x
数据来源机构:NOAA美国国家海洋与大气管理局-NCDC环境信息国家中心 ; B0 i: W9 k1 ?
产品名称:Ocean Heat Fluxes海洋热通量
: u, X, c' |5 S5 p0 Z) I8 D, P 时间精度:每日3-hourly记录
! Y& _% j* e2 q 空间精度:0.25° : x4 z( b, [+ k7 I3 G
注:本程序需调用UBC英属哥伦比亚大学Rich Pawlowicz 开发的程序包m_map;
0 x ^, f$ Y1 H# `4 |# } 亦需调用gshhs高精度海岸线数据。
6 d& W6 S! D f# m) h 源代码
" u' r- i0 `3 l. \4 g, V %南海1998年12月31日某时的海面热通量数据及分布图! d9 Y8 e5 y" |. U" m, @8 v( G
%%
0 \+ N& d2 K+ i0 i3 y) w: ~ %第一步:数据的展示' p& O8 F! l9 Z2 J
tic;Path="E:\OneDrive - Inc\Share_science\NC文件\";%文件夹地址
4 y6 y" D7 m$ f- u8 k4 ? filename="SEAFLUX-OSB-CDR_V02R00_FLUX_D19981231_C20160824.nc";%文件名
# g- ~0 ^7 G' W4 ^ abs_path=fullfile(Path,filename);%合并,绝对路径文件5 O( {0 s8 X) n1 p/ W. v
clearvars -except abs_path;%清除除了abs_path以外所有变量
% ^+ a: a3 I# H1 S( Z4 k ncdisp(abs_path);%读取所选文件
. f% u( b% f0 k; H+ w" H whos;toc;, K5 R, v, q, [ V
%% }; I; ~5 r9 T0 Q4 k) ?
%第二步:数据读取
; g$ ` u& y! O1 k# v$ |& x& D2 q tic;Lat_ini=ncread(abs_path,lat);
% z$ I5 X! J' C0 d Lon_ini=ncread(abs_path,lon);' G: A1 Z2 k& c8 ?) L7 |
Time_ini=ncread(abs_path,time);- r/ R7 v& b( U% x8 X/ [
Lflx_ini=ncread(abs_path,surface_upward_latent_heat_flux);
& a; n1 ?7 `2 W Sflx_ini=ncread(abs_path,surface_upward_sensible_heat_flux);8 l4 H& r; W; }1 J
toc;clear abs_path;" K4 A" E7 S$ k- f! f: G6 ?2 u
%%6 d7 F4 U9 J0 @' |
%第三步:确定边界/截取数据/ j2 J/ o: E5 I9 W
%空间边界:南海范围北纬2°至北纬25° 东经104°至东经124°& C4 p; m' E9 R& _) s
tic;Boundary=[104,124;2,25;]%修改范围[经度;纬度],注意顺序; \$ d; f: x' m' ~" S% G1 n
Spatial_Resolution=abs(Lat_ini(1)-Lat_ini(2));Sp_Re=Spatial_Resolution;%自动确定精度9 j1 z1 P5 G/ X, V( q* D j
Lon_start_No=find(Lon_ini-(Boundary(1))<=Sp_Re&Lon_ini-(Boundary(1))>=0);%找到离104°最近格点的编号(行号)- V, M2 t4 e4 O2 l, z7 h0 M. O
Lon_end_No=find(Lon_ini-(Boundary(3))<=Sp_Re&Lon_ini-(Boundary(3))>=0);%找到离124°最近格点的编号(行号)0 i/ F6 |- s0 L% m6 j
Lat_start_No=find(Lat_ini-(Boundary(2))<=Sp_Re&Lat_ini-(Boundary(2))>=0);%找到离2°最近格点的编号(行号)0 s' Y# d. Y1 t" ?* ^. M9 d) j
Lat_end_No=find(Lat_ini-(Boundary(4))<=Sp_Re&Lat_ini-(Boundary(4))>=0);%找到离25°最近格点的编号(行号)
& [) o- g/ A) v Lon=Lon_ini(Lon_start_No:Lon_end_No);: k9 E4 z4 ~( ]7 R$ e( \
Lat=Lat_ini(Lat_start_No:Lat_end_No);%北半球可能要颠倒一下编号,因为是从北极往赤道计数的- W5 {5 s( q$ J2 A; |; c' d
clear Lat_ini;clear Lon_ini;%回收内存. c0 |* V2 A& @3 S. O, |( l
Lat=double(Lat);Lon=double(Lon);%双精度化 _# o8 w* u4 d( ^5 F
latlim=[min(Lat),max(Lat)];lonlim=[min(Lon),max(Lon)];. F. X" n( V$ C4 g9 ?
Time=9%输入时间,注意三小时的倍数& I0 R9 q2 {4 l1 Q5 i
Time_No=Time/3+1;
: r/ u1 b* `% q6 n %截取数据并降维度; U8 o+ F# X4 E. I
Lflx=squeeze(Lflx_ini(Lon_start_No:Lon_end_No,Lat_start_No:Lat_end_No,Time_No));;%北半球有时可能要颠倒一下编号,因为是从北极往赤道计数的
! [# W1 a N9 f7 {+ l Sflx=squeeze(Sflx_ini(Lon_start_No:Lon_end_No,Lat_start_No:Lat_end_No,Time_No));;%北半球有时可能要颠倒一下编号,因为是从北极往赤道计数的1 F& M- i; q0 G/ d0 P7 \. t
clear Sp_Re;toc;%清除不再使用的变量0 @* X; L- U/ A; L/ O) m
%%1 Z0 C+ M) q4 N, i6 a
%第四步:绘图! \) A/ B* P0 {8 W* K
tic;figure;. P" o' h- p6 d
m_proj(mercator,lat,latlim,lon,lonlim);
, c Z. |" z# f: e [longrid,latgrid]=meshgrid(Lon,Lat);
. Z" {4 Z) B0 N6 f m_pcolor(longrid,latgrid,Sflx); %pcolor或contourf5 m2 C) ]2 B- O
tic;m_gshhs_f(patch,k);toc;%添加海岸线:c<l<h<f$ g- }. G6 ? v4 x
m_grid(box,fancy,tickdir,in);%网格化
}6 g" F/ L1 C; P y8 ` %brighten(.5);%亮化8 \7 a0 Y1 f9 o% W: m
colormap(hsv);3 H/ X8 [5 \% u) R" X
h = colorbar(h);%色标2 Q G, C0 L, f$ q8 x: q- g
h.Label.String = W·m^{-2}5 F) i% c) |7 ]
h.Location = eastoutside;%色标位置' m5 G3 u z) r$ h, L7 r& h. d
title([1998年12月31日,num2str(Time),:00,南海海面感热通量],...
* ^# e1 s! I( d5 ]* Q7 s5 c FontName,黑体,fontsize,12,FontWeight,"bold") %标题( U! N- w, @& `1 r+ D# W( F
xlabel(Longitude);ylabel(Latitude) ;set(gca,fontsize,9);) g4 u& D/ c7 `' Q/ e
toc;
$ v5 p$ i3 @ j
l5 E+ m9 u1 {, {& F& c$ _ 附上另外俩残次品绘图结果,数据来源:IOCADS 3 K/ x: }( L& y% ^/ }
, @- p( _6 e- w _3 _ : [9 w- f2 P: K: y8 a
d- O y! R* u( l/ }5 Z3 p! t) ~, C
) ^) ]# ?, Y |+ b/ u+ |2 a3 P+ h7 |
7 @0 I9 o1 d v7 }+ d0 P" T$ f( R6 x0 S/ X* c! b. C9 d
|