REM Ehsan Momeni

​

gmtset proj_length_unit inch
set range=122/154/20/50
set filename=Ehsan_Momeni_Ex10.ps
set grdfile=ETOPO1_Bed_g_gmt4.grd
set Japan_grid=Japan_grid.grd
set volcano_data=volcanos_of_the_world.txt
set volcanos_name=Volacanos_name.txt
set coord_sys=M6i
set Japan_volcanos=Japan_volcanos.txt
set Japan_Volacanos_name=Japan_Volacanos_name.txt
set cpt=rainbow.cpt

gmt grdcut %grdfile% -G%Japan_grid% -R%range% -V
gmt grd2cpt %Japan_grid% -Z -C%cpt% -V > NewCpt.cpt
gmt psbasemap -J%coord_sys% -R%range% -B10f5 -Y3i -P -K -V -BWESN+t"Japan Volcanos" > %filename%
gmt grdimage %Japan_grid% -J%coord_sys% -R%range% -CNewCpt.cpt -P -O -K -V >> %filename%
gmt pscoast -R%range% -J%coord_sys% -Di -W0.5p -Na -K -O -P -V  -Lf130/48/40/1000+u>>%filename%

extract_volcanoes %volcano_data% -R%range%  -Z -V > volcanos.txt

gawk "{print $1, $2}" volcanos.txt > %Japan_volcanos% 
gmt psxy %Japan_volcanos% -Skvolcano/0.2 -GSIENNA3 -R -J -W0.1p -O -K -V >>%filename%
gmt pstext %Japan_Volacanos_name% -R -J -F+f14p,4,blue+a0+jTC -D1p/12p -O -K -V >>%filename%

gmt psxy ridge.gmt -J -R%range% -W3p,red -O -P -K -V >> %filename%
echo 125 27 CB ridge | pstext -J  -R -F+f14,1,red+a30+j -N -O -K -V >> %filename%

gmt psxy transform.gmt -J -R%range% -W3p,blue -O -P -K -V >> %filename%
echo 139 33 CB transform | pstext -J  -R -F+f14,1,blue+a30+j -N -O -K -V >> %filename%

gmt psxy trench.gmt -J -R%range% -W3p,black -O -P -K -V >> %filename%
echo 148 41 CB trench | pstext -J  -R -F+f14,1,black+a30+j -N -O -K -V >> %filename%

gmt psscale -CNewCpt.cpt -D3i/-0.5i/6i/0.12ih -Bx2000f1000+l"Elevation" -By+l"meter" -N -P -O -V >>%filename%

del gmt.conf gmt.history NewCpt.cpt Japan_grid.grd Japan_volcanos.txt volcanos.txt
%filename%