*define
model	Example-1	# model name

*info
        Ray tracing in model {model} from Rayinvr
######################################################################

*setup
noapp
units	km	km

### some symbolic definitions

*define
vcont	6.0			# velocity level contoured in the model
fname	Example-1.v.in.save	# file name to save the model in when "Save" is pressed


### Generate a simple "zero-offset time" table

*call	readtab	map	table-zoffset
x	float	1
twtt	float	

argrang	x	0	1.0	300

const	
twtt	8.0

*nocall	pritab	*
table-zoffset

### Rayinvr-like model:

*call	rayinvr	gui	ps	Example-1
#vin	vp			dbas:v.in	# dbas:accrete.v.in1
#column	vp	100
#contour	vp	6.0
#zoffset	x		twtt	table-zoffset	

prefix	dbas:./

# bodies with names, densities, and optional "lithology" 
# classification types

grbody	-0.4		Root
-40	25
40	25
0	50

#grbody	2.67		Mountain
#-40	0
#0	-7.49
#40	0

# observed gravity profile. It can also be loaded from 
# database tables, but we are using a simple job line description here

grobs
-70	0.000
-60	0.000
-50	0.000
-40	0.000
-30	0.000
-20	0.000
-10	0.000
0	0.000
10	0.000
20	0.000
30	0.000
40	0.000
50	0.000
60	0.000
70	0

summary		# print brief model summary

### Define various graphical elements

*call	graphic
backgr	=same=
foregr	=same=

#palette	pal-vp		buor			5.5	8
#palette	pal-vp		gmtcopper			1	2	# 4.5	8.2
palette	pal-vp		gmtsplit			1	2	# 4.5	8.2
palhue	pal-vs		color	cont	2	5.	4	0.05	.95
palette	pal-raycodes	gmtjet			0	3
palette	pal-arrivals	gmtjet			0	4
palette	pal-ampl	gmtseis			-1	1 

palette	pal-density	gmtseis			2.5	3.3
palette	pal-adensity	gmtjet			-0.1	0.1

#palproc	edit	loop			pal-vp	# add a reversed copy of itself to the palette
#palproc	edit	repeat	2		pal-vp	# repeat the palette

palproc	edit	scale			pal-vp	# scale the Vp palette by Vp in the model
vp	Example-1

#palproc	edit	hbal			pal-vp		# histogram-balance the Vp palette
#vp	accrete

line	solid	1	melon		line-grcalc	# line for calculated gravity anomaly 
line	solid	1	yellow		line-observed

fill	solid	green		velocity

axis	axis-tt		bar	left
_label		T [s]
_ticks

axis	axis-grav	bar	left
_label		mGal
_ticks	

axis	axis-dist	bar	top
_label		X [km]
_ticks	

axis	axis-depth	bar	left
_label		Z [km]
_ticks	

# buttons and the corresponding commands to the ray tracer

button	button2		Reload	
Example-1	load-s	{fname}

button	button3		Save		
Example-1	save-s	{fname}

button	button4		Print		
Example-1	ps

### Create cross-section image
 
*call	image	im-section		##### plot the model and rays
1	1.	6	6
project	xy	width

range-x	-70	70
range-y	-15	60	yes

velmod	Example-1		
rho			shade			pal-density
a-rho			shade			pal-adensity
rho			paint

object	axis-depth
object	axis-dist

cscale	vert	0.8	-0.7	0.1	0.2	Density		# colour bar
pal-density

cscale	vert	0.6	-0.7	0.1	0.2	Anomaly		# colour bar
pal-adensity

# buttons to control modelling

object	button2
object	button3
object	button4

### Create gravity image

*call	image	im-gravity
1	1.	6	6
project	xy	width

backgr	gray

range-x	-70	70
range-y	-150	20

gravmod	Example-1			
calgrav						line-grcalc
observed					line-observed
# difference					line:red

object	axis-grav

### Display the interactive images

*call	gui
im-gravity			0	10
im-section

# make axes X coordinated in all sections

layout	vert				x	zoom

dockers
otree
#otree	graph
prop