*define
model	Example-1	# model name

*info
        Ray tracing in model {model} from Rayinvr 
	with Receiver Functions
######################################################################

*setup
noapp
units	km	km

### some symbolic definitions

*define
minx	28			# end of X range displayed
maxx	46			# end of X range displayed
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 synthetic Receiver Functions for three locations

*call	genrate	50	10000	-1000
gauss	1.0	0.0	200

hlist	x	real
30
hlist	p	real
0.015

### generate synthetic RF

*call	genrf	rf-rays			Example-1
ray	1.0	p	r	0	1	x	p

*call	table	*
x	p	timstrt	_tnrec_	lasttr

*nocall	print	test.print
samples	0	1000
traces	0	100

### plot the synthetic Receiver Functions along S-wave rays

*call	plot	wgva+-	up	median	1.0	0.05
name	rf-section			0	0

gain	200.0
setamp	peak

offset	x	1.0

trange	-1000	5000.0

*end		# remove all traces

### Rayinvr-like model:

*call	rayinvr	gui	ps	Example-1
*incl	Example-1.igm

summary		# print brief model summary

### Define various graphical elements

*call	graphic
backgr	=same=
foregr	=same=

palette	pal-vp1
-0.3	red		0	blue
0	white		0.3	red

#palette	pal-vp		buor			5.5	8
palette	pal-vp		gmtcopper			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 

#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	red		l-layers
line	solid	1	blue		l-columns
line	solid	1	cyan		l-contour
line	solid	1	blue		cells-vert
line	solid	1	yellow		cells-diag
line	solid	1	emerald		l-rays1
line	solid	1	plum		l-rays2
line	solid	1	brown		l-tau
line	solid	1	melon		line-grcalc	# line for calculated gravity anomaly 

line	solid	1	blue		ttimes-forw
line	solid	1	red		ttimes-rev
line	solid	1	yellow		ttimes-observed

line	solid	1	yellow		wiggle-line	# line for the wiggles

fill	solid	green		velocity
fill	solid	blue		wet
fill	solid	blue		lakes

fill	solid	yellow		pos-fill
fill	solid	gray		neg-fill

marker	star	0.2						m-sources

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	button3		Reload	
Example-1	load-s	{fname}

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

button	button5		Print		
Example-1	ps

### Create travel-time image

*call	image	im-ttimes		##### travel times
1	1.	6	6
project	xy	width

action	timered			8.0	# travel-time reduction, initial V=8 km/s

range-x	minx	maxx
range-y	-2	6	yes

traces	rf-section
wiggle			line			wiggle-line
wiggle			fill-pos		pos-fill
wiggle			fill-neg		neg-fill

object	axis-tt

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

range-x	minx	maxx
range-y	-5	60	yes

velmod	Example-1		
layers						l-layers
#cells			vert			cells-vert
#cells			diag			cells-diag
#vs			shade			pal-vs
#vp			shade			pal-vp
vp			column			l-columns
vp			contour			l-contour
sources						m-sources
sources			label			Shot %d
sources			traces			source-%d	# interactive data input
rays			forward			l-rays1
#rays			reverse			l-rays2
#rays						pal-raycodes
tau			velocity		5.7
tau			line			l-tau	

object	rf-rays
object	axis-depth
object	axis-dist


# buttons to control modelling

object	button3
object	button4
object	button5
object	-spacer-

### Display the interactive images

*call	gui
im-ttimes
im-section

# add interactive editor for parameters of "genrate"

proprty	genrate.1

# make axes X coordinated in all sections

layout	vert				x	zoom

caption	Interactive modeling of Receiver Functions

dockers
otree
#otree	flow
prop