( w" x1 P2 F( T9 Y) C
图片赏析
, X( N- S2 l0 k/ s/ N- U 4 X- N. R2 m* S8 e; P2 \- h, q: y0 m
' t& I8 F' o5 w7 @3 |$ M$ t
4 h6 |" j3 E/ b$ `1 \# F 前言 $ O3 q4 [* C7 K- V
& G g" ?0 R' P 南海的经纬度:南海范围北纬2°至北纬25° 东经104°至东经124°
; @8 P8 ~2 \1 Z' [/ G 数据来源: ) l& a- g9 }6 S3 j
3 e; {. P$ U) K+ w 数据链接提供者:广东海洋大学海洋与气象学院海洋科学系纪同学
% }" q; W! |! ?! M/ K# v6 { 数据来源机构:NOAA美国国家海洋与大气管理局-NCDC环境信息国家中心 0 ]( J' z) O+ h7 X4 J( w) C$ F
产品名称:Ocean Heat Fluxes海洋热通量 3 }! d. M8 _& W5 J8 ~2 `
时间精度:每日3-hourly记录
3 T% Y2 ^- q' a2 h& u/ U 空间精度:0.25° " p# g2 f( [ I5 U6 G
注:本程序需调用UBC英属哥伦比亚大学Rich Pawlowicz 开发的程序包m_map;
K' H! {1 ^8 d- n+ ]3 i1 e 亦需调用gshhs高精度海岸线数据。
" w; w% N. O& w3 ?1 V J& u+ z 源代码% q. Y5 p' ]/ n- x
%南海1998年12月31日某时的海面热通量数据及分布图. i5 s. n6 F" ~$ O
%%% e! O: ?& F9 I5 |% e
%第一步:数据的展示
- d P0 r N7 t' s tic;Path="E:\OneDrive - Inc\Share_science\NC文件\";%文件夹地址1 T! X+ Q9 @! h+ j& } V
filename="SEAFLUX-OSB-CDR_V02R00_FLUX_D19981231_C20160824.nc";%文件名
! |/ B' M" J e$ p/ Y5 x abs_path=fullfile(Path,filename);%合并,绝对路径文件) ?8 t2 M) F2 I1 T$ c9 b
clearvars -except abs_path;%清除除了abs_path以外所有变量( S" g$ `1 y5 W9 { g) ^& N4 L9 f
ncdisp(abs_path);%读取所选文件' Y; L( m+ o0 {! _
whos;toc;
, v k& r: e+ O; ~* ]. D %%
. c* D! j3 D0 ~* ]/ ^ %第二步:数据读取 d" I1 l1 ?2 P) `' g
tic;Lat_ini=ncread(abs_path,lat);' L$ ^4 c5 T) C& O' w* F, c# H
Lon_ini=ncread(abs_path,lon);# V" F, K$ r( W1 Y9 i
Time_ini=ncread(abs_path,time);
2 q' H0 [6 a+ V* U& b Lflx_ini=ncread(abs_path,surface_upward_latent_heat_flux);
+ J @ q7 n z; r# n2 o l Sflx_ini=ncread(abs_path,surface_upward_sensible_heat_flux);
) W: v+ ^3 P1 g8 b' Z toc;clear abs_path;( |/ n4 Q e. u5 r
%%
, \7 a& q! G- V8 I %第三步:确定边界/截取数据
' D7 `; `6 p: f9 Z %空间边界:南海范围北纬2°至北纬25° 东经104°至东经124°6 y& ~4 I7 p! J# e6 }$ J
tic;Boundary=[104,124;2,25;]%修改范围[经度;纬度],注意顺序- c3 u4 ]' p; S+ L0 o! O; k" y, h- g
Spatial_Resolution=abs(Lat_ini(1)-Lat_ini(2));Sp_Re=Spatial_Resolution;%自动确定精度6 y7 F* q7 E/ D
Lon_start_No=find(Lon_ini-(Boundary(1))<=Sp_Re&Lon_ini-(Boundary(1))>=0);%找到离104°最近格点的编号(行号)
) R$ k0 P/ _' T Lon_end_No=find(Lon_ini-(Boundary(3))<=Sp_Re&Lon_ini-(Boundary(3))>=0);%找到离124°最近格点的编号(行号)
' U- f3 Y) e8 M- q' {% _ Lat_start_No=find(Lat_ini-(Boundary(2))<=Sp_Re&Lat_ini-(Boundary(2))>=0);%找到离2°最近格点的编号(行号)& _0 ~( o3 N/ D) ?! ]- h3 Y
Lat_end_No=find(Lat_ini-(Boundary(4))<=Sp_Re&Lat_ini-(Boundary(4))>=0);%找到离25°最近格点的编号(行号)
$ N& w, ]! I) x- Z3 u Lon=Lon_ini(Lon_start_No:Lon_end_No);
2 H1 K6 t0 z% {+ [2 o Lat=Lat_ini(Lat_start_No:Lat_end_No);%北半球可能要颠倒一下编号,因为是从北极往赤道计数的
/ O3 n8 i+ T) Z- J' o4 A: H clear Lat_ini;clear Lon_ini;%回收内存
3 Y5 q0 J2 q9 ]6 |1 u/ _ Lat=double(Lat);Lon=double(Lon);%双精度化9 a7 w# j9 [) c$ q Z
latlim=[min(Lat),max(Lat)];lonlim=[min(Lon),max(Lon)];! D; k) D) J. G& c% Z/ I/ J4 l3 T. A
Time=9%输入时间,注意三小时的倍数3 Y# }- U B ]
Time_No=Time/3+1;& A7 p8 ?) Y. ]' {& A2 [" v6 l
%截取数据并降维度! P0 l6 E) l7 ?+ S
Lflx=squeeze(Lflx_ini(Lon_start_No:Lon_end_No,Lat_start_No:Lat_end_No,Time_No));;%北半球有时可能要颠倒一下编号,因为是从北极往赤道计数的
9 k- P0 Q9 y( @$ q, z3 A) H Sflx=squeeze(Sflx_ini(Lon_start_No:Lon_end_No,Lat_start_No:Lat_end_No,Time_No));;%北半球有时可能要颠倒一下编号,因为是从北极往赤道计数的! D4 e5 N; j2 I
clear Sp_Re;toc;%清除不再使用的变量: K' a$ ]' j/ ^
%%
' M9 W5 g& |" u. ?) d %第四步:绘图; \- V) O) |( d |4 o( d* f0 q
tic;figure;2 N, g0 U9 Q; O( K8 a: |
m_proj(mercator,lat,latlim,lon,lonlim);
+ _0 o# x6 ^. P a [longrid,latgrid]=meshgrid(Lon,Lat);" e6 N$ C; {% h
m_pcolor(longrid,latgrid,Sflx); %pcolor或contourf
|- m k( ^ X$ k' A tic;m_gshhs_f(patch,k);toc;%添加海岸线:c<l<h<f+ P0 Q' e! ]3 E; z: b
m_grid(box,fancy,tickdir,in);%网格化' v$ V$ x' j, S
%brighten(.5);%亮化
* Z/ R* w" V [) _ colormap(hsv);
& U3 H& A) v/ ]. J h = colorbar(h);%色标5 v$ | Z0 F+ c' ^# c- N
h.Label.String = W·m^{-2}( i# j) e7 V9 h
h.Location = eastoutside;%色标位置+ z' `7 m8 D. T. u0 Y
title([1998年12月31日,num2str(Time),:00,南海海面感热通量],.... V( @# \) e, W. _4 a
FontName,黑体,fontsize,12,FontWeight,"bold") %标题' y* a0 o( v) z' }$ X& F
xlabel(Longitude);ylabel(Latitude) ;set(gca,fontsize,9);
+ J/ ], ~0 I, K" x9 @2 ^ toc;
! k9 ~* R" O% T0 w% d ; C& _( S' G0 p, |7 W' J& D% W, W
附上另外俩残次品绘图结果,数据来源:IOCADS % B/ u+ a3 r5 F2 q' X( r& C
: K6 |) Z8 J. Y- j & W/ s/ T3 t1 S
* L6 c1 D; i4 D' N8 k! Y* u1 d3 L) `" p0 u9 j t) Y
5 V( l* u5 S$ C! @7 q8 b1 n! o8 q, @3 k$ s/ U' `
|