#############################################################
#        ACCREETE model in 3D
#############################################################

###### Job set-up. In particular, use kilometers for distance units
###################################################################

*setup
noapp
units	km	km

###### load the profile track line and midpoints from ASCII files
#################################################################

*call	readtab table   x-profile
x	real	1
lat	float
lon	float

file    all     pc-coords.txt

*call	readtab table   midpoints
lon	float
lat	float	1
z	float

file    all     ../MAP/midpoints-rev.table
	z	0.	# draw midpoints at z=0.0

### load the interpolated Moho depth grid
#####################################################

*call	readtab	map	moho
lon	float	1
lat	float	2
z	float

argrang	lon	-131.3	0.01	111
argrang	lat	54.6	0.01	61

file	all	moho_depth.xyz

###### create graphic elements
#####################################################

*call	graphic
backgr	black
foregr	black

# line and fill styles:

line	solid	1	gray		line-shore
line	solid	1	white		line-borders
line	solid	1	blue		line-rivers
line	dash	1	melon		line-grid

fill	none	blue		fill-wet
fill	none	blue		fill-lakes

# palette for Moho depth:

backgr	-same-
foregr	-same-		# foreground of the same color as top of the palette

#palette	pal-moho	gmtrainbow		25.0	31.0
#palette	pal-moho	gmtjet			25.0	31.0
palette	pal-moho	buor			25.0	31.0

# buttons for interactive operation:

button	button_trace	Trace rays
pc-section	exec	trace

button	button_save	Save model
pc-section	exec	save

button	button_summary	Print summary
pc-section	exec	print.summary

button	button_display	Display data
pc-section	print	This option only prints in this demo!

###### Define 3D projectors
#####################################################

# this will project 2D images onto profile cross-section:

*call	view3d	fence3D
fence	x	lon	lat	x-profile
geom	sphere				km

# this will project 2D images onto the surface:

*call	view3d	surface3D
geom	sphere				km

##### Create surface image
####################################################

*call	image	surface-maps
range-x	-131	-130		# longitudes
range-y	55	56		# latitudes

# map of the Moho

line	none

surface	lon	lat	z			moho
pal-moho		z		# fill colors based on z (depth)

# scatter plot of midpoints:

line	points	3	red		line-midpoints

3dtabs	lon	lat	z
midpoints

###### Create 2D Accrete Portland Canal section
####################################################

*incl	pc-image.inc

###### Assemble the 3D image
####################################################

*call	image	acc-3Dimage
rangec3	400			# 200-km viewing box

# list selectable views, the first is the default:
 
view	From above		-131.5	55	0.	-20	10.	1.0
view	Globe			0	0

# place 2-D cross-section in "fence" projection:

object	pc-section			fence3D

# place 2-D sirface maps in "surface" projection:

object	surface-maps			surface3D		

# draw in a GMT basemap at "high" resolution:

coast	-134	-129	52	58	high
surf	land			fill-dry
surf	sea			fill-wet
surf	lakes			fill-lakes
bndry	shore			line-shore
bndry	natnl			line-borders
bndry	state			line-borders
bndry	marine			line-borders
rivers	r			line-rivers
grid	line			line-grid	1.0	1.0

# Moho depth color bar:

cscale	vert	0.1	0.8	0.05	0.2	Moho depth (km)
pal-moho

# image lighting:

light	0	0	1000	0.9	Illumination
diff		white
amb		white
spec		cyan

# buttons for interactive operation:

object	button_trace
object	button_save
object	button_summary
object	button_display

# animation (autoplay):

a-rot	0.1	60	30	1.0

##### start the display server
####################################################

*call	gui
acc-3Dimage

# request docking windows:

dockers
otree		canvas		# object tree
otree		graph		# graphics DB tree
prop				# property editor