*params	
shotx	0.0		# shot position in the model
colx	25.0		# column position
ofact	-1		# factor to multiply offset (1 for forward, -1 for reverse)

*info
        Ray tracing in model Example 1 from Rayinvr
######################################################################

## set the distances and depths measured in km

*setup
noapp
units	km	km

### some symbolic definitions

*define
model	Example-1	# model name

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

*nocall	pritab	*
table-zoffset

### Rayinvr-like model:

*call	rayinvr	gui	ps	{model}
vin	vp			dbas:v.in	# dbas:accrete.v.in1
column	vp	colx
prefix	dbas:./
summary		# print brief model summary

#### reflectivity modeling

*call	reflect	displ	24	24000		offs		{model}-synth
modes	full	# noconv	upgo
fwin	2	2	4	8	12	0.
phvwin	2	2.	3.	80.	120.	800	50.e3
source	expl	0.0001	spike	0.	3.		# 5.
receiv	0.	0.01	0.5	80
refldep	0
reduce	27.0

layersm	sphere		{model}

*call	hdrmath
set	dummy	integer	0

*call	break	trace	dummy	break
*call	sortcut	break

####### an info reporting elapsed time

*info	eof	time

####### test printout of the resulting traces

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

#### multiply trace header 'offs' by 1000 
### to create "soffset" (in meters). Also create profile distance:
### x = shotx + x

*call	hdrmath
set	soffset	double
1000			offs
set	x	float	shotx
ofact			offs
set	xcmp	float
0.5			shotx	0.5			x

###### Write the output in SIA format:

*call	write		{model}-synth-{shotx}.sia
Reflectivity synthetics in model {model}

###### List trace headers in file 'synth-{model}.table'
###### Note how the name of 'model' is substituted

*call	table	*
x	soffset	xcmp	timstrt	lasttr

###### select only the traces with component number = 1 for plotting

*call	convert
extract	1

######## Travel-time reduction and scaling for the plot

*call	trcscal	0	20000

######## Plot the trace section in PostScript

*call	plot	WGVA	down	median	1.0	100
gain	.3
#setamp	rms		500	1200

ps	300	PS/{model}-synth-{shotx}.ps
show

window	6	6	.5	.5	0	8.5	11

offset	offs	4e-2
trange	0	20000

font	12				black
	
line	dash	0.1
1	3
timlab	1000	%1.1f
timtitl	left	TWTT [s]

ticks		0		0		0
margins	1.2	.8	.6

line	solid	0.1
scale	20	5		%1.0f
text	1010	0	0.1	0
\Offset [km]

font	10				black
title	Reflectivity synhetic