Reservoir model with thin layers and output of velocity and Q columns
 
######### common parameters are convenient to give by symbolic substitutions:

*define
model	majin1	# model name

qpval	30				# Qp (constant in this model)
qsval	20				# Qs (constant in this model)
minr	0.0001				# starting receiver distance
incr	0.02				# increment in receiver distance
maxr	1.01				# maximum receiver distance

#### "reflectivity" modeling. 
#### Note that in 'reflect', the depths and receiver coordinates (x) are in km

*call	reflect	displ	2	4096		x		{model}
modes	pp	# noconv	upgo
fwin	2	10	20	60	100	0.		# frequency window
phvwin	2	0.5	2.0	40.0	60.0	3000		# phase-velocity window
source	expl	0.0001	spike	0.	3.		# 5.
receiv	0.	minr	incr	maxr				# receiver line
reduce	2.5							# travel-time reduction (optional)
refldep	0							# depth of start of reflectivity

trials	1	run	0	1	no_q		# make three runs idenfified by header 'run': 
							# run = 1 in the model as described
							# run = 2 with fine layering
							# run = 3 with fine layering without Q	

layers	flat	0.02	0.04	alt	abs	-1	0.03		# apply alternating 40-m thick Vp variations by 3%

0		2.0		0.5	qpval		qsval			2.
0.6	2.0	2.438	0.5	1.025	qpval		qsval		2.	2.14	1
0.75	2.438	2.48	1.025	1.232	qpval		qsval		2.14	2.8	1
0.95	2.48	2.638	1.232	1.125	qpval		qsval		2.16	2.8	1
1.13	2.638	2.438	1.125	1.625	qpval		qsval		2.23	3.2	1
1.25	2.438	2.7	1.625	1.825	qpval		qsval		2.16	2.25	1
1.35	2.7	2.5	1.825	1.332	qpval		qsval		2.25	2.17	1
1.47	2.5	2.738	1.332	1.625	qpval		qsval		2.17	2.26	1
1.57	2.738	3.048	1.625	1.244	qpval		qsval		2.26	2.40	1	
1.69	3.048	2.838	1.244	1.650	qpval		qsval		2.40	2.28	1
1.79	2.838	3.148	1.650	1.750	qpval		qsval		2.28	2.43	1
1.91	3.148	3.148	1.750	1.750	qpval		qsval		2.43	2.43	1

	# output tables of model parameters vs depth and time:

column	column	twtt	stack	vel-stacking	
twtt	vp	vs	qp	qs	rho

column	column	twtt	int	vel-twtt
twtt	vp	vs	qp	qs	rho

column	column	depth	int	vel-depth
z	vp	vs	qp	qs	rho

######### info reporting the elapsed time

*info	eof	time