|
海洋潮汐由诸多分潮组成，同潮图的绘制是潮汐研究的常见操作。这里介绍GMT绘制同潮图的主要步骤，并导出振幅和相位的等值线，转为kml格式，以便在google earth查看。 数据下载使用AVISO网站的FES2014潮汐模型，各分潮为单独的nc文件。内容包含： [C] 纯文本查看 复制代码 $ ls
2n2.nc k2.nc m3.nc m8.nc mm.nc msf.nc nu2.nc q1.nc s2.nc ssa.nc
j1.nc l2.nc m4.nc mf.nc mn4.nc mu2.nc o1.nc r2.nc s4.nc t2.nc
k1.nc m2.nc m6.nc mks2.nc ms4.nc n2.nc p1.nc s1.nc sa.nc 以山东周边海域M2和S2分潮为例绘图。 主要用到了grdcontour绘制等值线，在相位一行代码中，注意加入-Z+p，代码如下。 [C] 纯文本查看 复制代码 #!/usr/bin/env bash
ps=sd.ps
gmt pscoast -R115/125/34/42 -JM6i -Baf -Gburlywood -Di -A10 -P -K --FORMAT_GEO_MAP=dddF > $ps
gmt grdimage "m2.nc?amplitude" -J -R -O -K -Q >> $ps
gmt pscoast -R -J -Gwhite -Df -K -O >> $ps
gmt grdcontour -R "m2.nc?phase" -J -Gl123/37/120/34,122/39/125/39 -C20 -A40+f7p+a0+u@. -Z+p -Wc0.25p,white -Wa0.75p,white -O -K >> $ps
gmt grdcontour -R "m2.nc?amplitude" -J -GlBL/TR -C20 -A40+f7p+a0 -Wc0.25p,black -Wa0.75p,black -O -K >> $ps M2 ( [/ I6 M+ Q! X% N1 I
S2 , ?5 X* G( b: m* l: D5 \
0 R. t; h7 p# ]" O; g) F
; r$ x1 w% }! Q3 A% N数据导出首先使用grdcontour的-D输出等值线文件，再使用gmt2kml制作kml文件，可进一步导入Google Earth。 * d1 ~) O" t' ^: y
[C] 纯文本查看 复制代码 # Output data
gmt grdcontour -R "m2.nc?phase" -J -C10 -A20+f7p+a0 -Wc0.25p,gray -Wa0.75p,gray -Dcontours_p.txt
gmt grdcontour -R "m2.nc?amplitude" -J -C10 -A20+f7p+a0 -Wc0.25p,red -Wa0.75p,red -Dcontours_a.txt
# Transform to kml
gmt 2kml contours_p.txt -Wthick,black -Fl > contours_p.kml
gmt 2kml contours_a.txt -Wthick,white -Fl > contours_a.kml 导入Google Earth + ?2 s1 U* d5 k) x) d2 ]
1 y* X C% J+ Q+ [ Y) K | 
