maxrjones
/
gmt-fork
mirror of https://github.com/meghanrjones/gmt.git


  
1

	
2

	
3

	
4

	
5

	
6

	
7

	
8

	
9

	
10

	
11

	
12

	
13

	
14

	
15

	
16

	
17

	
            #!/usr/bin/env bash
#
#
# script to plot mendo wa1 combined data
#
# cdp = 40 * coordinate on line, 30km max depth
#

ps=segy_04.ps
area1=-R-35/6/0/30
proj1="-Jx0.15i/-0.15i"

gmt makecpt -T-5/5 -Cpolar > test.cpt
gmt segy2grd $area1 -I0.1/0.1 "${src:-.}"/wa1_mig13.segy -Gtest.nc
gmt grdimage $area1 $proj1 -K test.nc -Ctest.cpt -Xc -P > $ps
gmt psbasemap -R -J -Baf -O  >> $ps
rm -f test.nc test.cpt