海洋潮汐由诸多分潮组成,同潮图的绘制是潮汐研究的常见操作。这里介绍GMT
绘制同潮图的主要步骤,并导出振幅和相位的等值线,转为kml
格式,以便在google earth
查看。
数据下载
使用AVISO网站的FES2014潮汐模型,各分潮为单独的nc文件。
内容包含:
1 2 3 4
| $ 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
,代码如下。
1 2 3 4 5 6 7 8 9
| #!/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
|
![](/images/sd.png)
M2
![](/images/sd_s2.png)
S2
数据导出
首先使用grdcontour
的-D
输出等值线文件,再使用gmt2kml
制作kml文件,可进一步导入Google Earth。
1 2 3 4 5 6 7
| # 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
|
![](/images/ge.png)
导入Google Earth