Permalink
Browse files

添加图片绘制代码

  • Loading branch information...
1 parent 54086c2 commit 8b800ee2e51637801d04f98e7b86511610367397 jingxu committed Jan 9, 2017
@@ -0,0 +1,4 @@
+#!/bin/bash
+gmt pscoast -R50/160/-15/15 -JM5.3i -Gburlywood -Sazure -A500 -K -P > GMT_contour-anno1.ps
+gmt grdcontour geoid.nc -J -O -B20f10 -BWSne -C10 -A20+f8p -Gd1.5i -S10 -T+lLH >> GMT_contour-anno1.ps
+rm gmt.*
@@ -0,0 +1,4 @@
+#!/bin/bash
+gmt pscoast -R50/160/-15/15 -JM5.3i -Gburlywood -Sazure -A500 -K -P > GMT_contour-anno2.ps
+gmt grdcontour geoid.nc -J -O -B20f10 -BWSne -C10 -A20+f8p -Gn1/1i -S10 -T+lLH >> GMT_contour-anno2.ps
+rm gmt.*
@@ -0,0 +1,10 @@
+#!/bin/bash
+cat << EOF > fix.txt
+80 -8.5
+55 -7.5
+102 0
+130 10.5
+EOF
+gmt pscoast -R50/160/-15/15 -JM5.3i -Gburlywood -Sazure -A500 -K -P > GMT_contour-anno3.ps
+gmt grdcontour geoid.nc -J -O -B20f10 -BWSne -C10 -A20+d+f8p -Gffix.txt/0.1i -S10 -T+lLH >> GMT_contour-anno3.ps
+rm gmt.*
@@ -0,0 +1,4 @@
+#!/bin/bash
+gmt pscoast -R50/160/-15/15 -JM5.3i -Gburlywood -Sazure -A500 -K -P > GMT_contour-anno4.ps
+gmt grdcontour geoid.nc -J -O -B20f10 -BWSne -C10 -A20+d+f8p -GLZ-/Z+ -S10 -T+lLH >> GMT_contour-anno4.ps
+rm gmt.*
@@ -0,0 +1,4 @@
+#!/bin/bash
+gmt pscoast -R50/160/-15/15 -JM5.3i -Gburlywood -Sazure -A500 -K -P > GMT_contour-anno5.ps
+gmt grdcontour geoid.nc -J -O -B20f10 -BWSne -C10 -A20+d+f8p -GXcross.txt -S10 -T+lLH >> GMT_contour-anno5.ps
+rm gmt.*
@@ -0,0 +1,5 @@
+#!/bin/bash
+gmt pscoast -R50/160/-15/15 -JM5.3i -Gburlywood -Sazure -A500 -K -P > GMT_contour-anno6.ps
+gmt grdcontour geoid.nc -J -O -K -B20f10 -BWSne -C10 -A20+d+f8p -Gl50/10S/160/10S -S10 -T+l >> GMT_contour-anno6.ps
+gmt psxy -R -J -O -SqD1000k:+g+LD+an+p -Wthick transect.txt >> GMT_contour-anno6.ps
+rm gmt.*
@@ -0,0 +1,5 @@
+#!/bin/bash
+gmt pscoast -R50/160/-15/15 -JM5.3i -Gburlywood -Sazure -A500 -K -P > GMT_contour-anno7.ps
+gmt grdcontour geoid.nc -J -O -K -B20f10 -BWSne -C10 -A20+d+u" m"+f8p -Gl50/10S/160/10S -S10 -T+l >> GMT_contour-anno7.ps
+gmt psxy -R -J -O -SqD15d:+gblack+fwhite+LD+o+u\\260 -Wthick transect.txt >> GMT_contour-anno7.ps
+rm gmt.*
@@ -0,0 +1,6 @@
+#!/bin/bash
+gmt convert -i0,1,4 -Em150 transect.txt | gawk '{print $1,$2,int($3)}' > fix2.txt
+gmt pscoast -R50/160/-15/15 -JM5.3i -Gburlywood -Sazure -A500 -K -P > GMT_contour-anno8.ps
+gmt grdcontour geoid.nc -J -O -K -B20f10 -BWSne -C10 -A20+d+u" m"+f8p -Gl50/10S/160/10S -S10 -T+l >> GMT_contour-anno8.ps
+gmt psxy -R -J -O -Sqffix2.txt:+g+an+p+Lf+u" m"+f8p -Wthick transect.txt >> GMT_contour-anno8.ps
+rm gmt.*
@@ -0,0 +1,27 @@
+#!/bin/bash
+R=-R-85/5/10/55
+gmt grdgradient topo5.nc -Nt1 -A45 -Gtopo5_int.nc
+gmt set FORMAT_GEO_MAP ddd:mm:ssF FONT_ANNOT_PRIMARY +9p FONT_TITLE 22p
+gmt project -E-74/41 -C-17/28 -G10 -Q > great_NY_Canaries.txt
+gmt project -E-74/41 -C2.33/48.87 -G100 -Q > great_NY_Paris.txt
+km=`echo -17 28 | gmt mapproject -G-74/41/k -fg --FORMAT_FLOAT_OUT=%.0f -o2`
+cat << EOF > ttt.cpt
+0 lightred 3 lightred
+3 lightyellow 6 lightyellow
+6 lightgreen 100 lightgreen
+EOF
+gmt grdimage ttt_atl.nc -Itopo5_int.nc -Cttt.cpt $R -JM5.3i -P -K -nc+t1 > GMT_contour-anno9.ps
+gmt grdcontour ttt_atl.nc -R -J -O -K -C0.5 -A1+u" hour"+v+f8p,Bookman-Demi -GL80W/31N/17W/26N,17W/28N/17W/50N -S2 >> GMT_contour-anno9.ps
+gmt psxy -R -J -Wfatter,white great_NY_Canaries.txt -O -K >> GMT_contour-anno9.ps
+gmt pscoast -R -J -B20f5 -BWSne+t"Tsunami travel times from the Canaries" -N1/thick -O -K -Glightgray -Wfaint -A500 >> GMT_contour-anno9.ps
+gmt convert great_NY_*.txt -E | gmt psxy -R -J -O -K -Sa0.15i -Gred -Wthin >> GMT_contour-anno9.ps
+gmt psxy -R -J -Wthick great_NY_Canaries.txt -O -K -Sqn1:+f8p,Times-Italic+l"Distance Canaries to New York = $km km"+ap+v >> GMT_contour-anno9.ps
+
+gmt psxy -R -J -Wthinner great_NY_Paris.txt -SqD1000k:+an+o+gblue+LDk+f7p,Helvetica-Bold,white -O -K >> GMT_contour-anno9.ps
+cat << EOF | gmt pstext -R -J -O -K -Gwhite -Wthin -Dj0.1i/0.1i -F+f8p,Bookman-Demi+j >> GMT_contour-anno9.ps
+74W 41N RT New York
+2.33E 48.87N CT Paris
+17W 28N CT Canaries
+EOF
+gmt psxy -R -J -O -T >> GMT_contour-anno9.ps
+rm gmt.*
View
@@ -0,0 +1,35 @@
+>
+58 -11
+63 -6
+
+>
+63 -6
+68 -2
+
+>
+68 -2
+72 -1
+
+>
+72 -1
+78 3
+
+>
+150 3
+158 13
+
+>
+95 -8
+102 -7
+
+>
+102 -7
+110 -8
+
+>
+110 -8
+120 -2
+
+>
+120 -2
+123 0
View
@@ -0,0 +1,4 @@
+80 -8.5
+55 -7.5
+102 0
+130 10.5
View
@@ -0,0 +1,6 @@
+77 4 -3595
+90.3631916207 2.02193063653 -2836
+103.693238887 -0.0652707067064 -6
+117.024382171 -2.14895473069 -9
+130.390638631 -4.11662676531 -4527
+143.818818989 -5.86098872838 2574
View
Binary file not shown.
Oops, something went wrong.

0 comments on commit 8b800ee

Please sign in to comment.