# History file 'prem-mod.h' transforms surfaces and velocity section
# in model PREM into VRML or GOCAD
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Input files required
  chk.pl: "data/prem/" "prem-mod.dat"
  chk.pl: "data/len/" "len-cam.dat"
  chk.pl: "data/len/" "len-dli1.dat"
  chk.pl: "data/len/" "len-dli2.dat"
  chk.pl: "forms/"    "hsv.dat"

# General settings
# ~~~~~~~~~~~~~~~~
# Form of output file
  VRML='VRML2'
 #VRML='GOCAD'

# Input data
  MODEL='prem-mod.dat'

# Colour properties
  COLORS='hsv.dat'
  TRANSP=0.00  SPECULAR=0.22  SHININESS=0.10

# Structural interfaces in the model
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Name of output file
  WRL='prem-mod.wrl'  # VRML
 #WRL='prem-mod.mx'   # GOCAD
  iniwrl:

# Colours correspond to surface indices
  COLUMN07='ISRF'
  KOLSRF=7  VPER=-4 VREF=1  CREF=0.6666667

# Triangulating and plotting interfaces
  O1=0.   O2=-1.55  O3=0.
  D1=0.1  D2=0.05   D3=10.
  N1=32   N2=63     N3=639
  VRTX='vrtx.tmp'  TRGL='trgl.tmp'  PLGN='plgn.tmp'
  modsrf:
  trgl:
  KOLUMN=7
  VALUE=1  VRTXN='prem-p01.out'  TRGLN='prem-t01.out'  trglsort:
  VALUE=2  VRTXN='prem-p02.out'  TRGLN='prem-t02.out'  trglsort:
  VALUE=3  VRTXN='prem-p03.out'  TRGLN='prem-t03.out'  trglsort:
  VALUE=4  VRTXN='prem-p04.out'  TRGLN='prem-t04.out'  trglsort:
  VALUE=5  VRTXN='prem-p05.out'  TRGLN='prem-t05.out'  trglsort:
  VALUE=6  VRTXN='prem-p06.out'  TRGLN='prem-t06.out'  trglsort:
  VALUE=7  VRTXN='prem-p07.out'  TRGLN='prem-t07.out'  trglsort:
  VALUE=8  VRTXN='prem-p08.out'  TRGLN='prem-t08.out'  trglsort:
  VALUE=9  VRTXN='prem-p09.out'  TRGLN='prem-t09.out'  trglsort:
  VALUE=10 VRTXN='prem-p10.out'  TRGLN='prem-t10.out'  trglsort:
  VALUE=11 VRTXN='prem-p11.out'  TRGLN='prem-t11.out'  trglsort:
  VALUE=12 VRTXN='prem-p12.out'  TRGLN='prem-t12.out'  trglsort:
  VALUE=13 VRTXN='prem-p13.out'  TRGLN='prem-t13.out'  trglsort:
  KOOR=2 TOCART=1
  PTS='prem-p01.out' PTSOUT='prem-v01.out'  coorchg:
  PTS='prem-p02.out' PTSOUT='prem-v02.out'  coorchg:
  PTS='prem-p03.out' PTSOUT='prem-v03.out'  coorchg:
  PTS='prem-p04.out' PTSOUT='prem-v04.out'  coorchg:
  PTS='prem-p05.out' PTSOUT='prem-v05.out'  coorchg:
  PTS='prem-p06.out' PTSOUT='prem-v06.out'  coorchg:
  PTS='prem-p07.out' PTSOUT='prem-v07.out'  coorchg:
  PTS='prem-p08.out' PTSOUT='prem-v08.out'  coorchg:
  PTS='prem-p09.out' PTSOUT='prem-v09.out'  coorchg:
  PTS='prem-p10.out' PTSOUT='prem-v10.out'  coorchg:
  PTS='prem-p11.out' PTSOUT='prem-v11.out'  coorchg:
  PTS='prem-p12.out' PTSOUT='prem-v12.out'  coorchg:
  PTS='prem-p13.out' PTSOUT='prem-v13.out'  coorchg:
  PTS=               PTSOUT=
  PROPERTIES='N1 N2 N3 ISRF'                         # just for GOCAD
  NAME='Surface01' VRTX='prem-v01.out'  TRGL='prem-t01.out'  srfwrl:
  NAME='Surface02' VRTX='prem-v02.out'  TRGL='prem-t02.out'  srfwrl:
  NAME='Surface03' VRTX='prem-v03.out'  TRGL='prem-t03.out'  srfwrl:
  NAME='Surface04' VRTX='prem-v04.out'  TRGL='prem-t04.out'  srfwrl:
  NAME='Surface05' VRTX='prem-v05.out'  TRGL='prem-t05.out'  srfwrl:
  NAME='Surface06' VRTX='prem-v06.out'  TRGL='prem-t06.out'  srfwrl:
  NAME='Surface07' VRTX='prem-v07.out'  TRGL='prem-t07.out'  srfwrl:
  NAME='Surface08' VRTX='prem-v08.out'  TRGL='prem-t08.out'  srfwrl:
  NAME='Surface09' VRTX='prem-v09.out'  TRGL='prem-t09.out'  srfwrl:
  NAME='Surface10' VRTX='prem-v10.out'  TRGL='prem-t10.out'  srfwrl:
  NAME='Surface11' VRTX='prem-v11.out'  TRGL='prem-t11.out'  srfwrl:
  NAME='Surface12' VRTX='prem-v12.out'  TRGL='prem-t12.out'  srfwrl:
  NAME='Surface13' VRTX='prem-v13.out'  TRGL='prem-t13.out'  srfwrl:


# P-wave velocity section
# ~~~~~~~~~~~~~~~~~~~~~~~

# Colours correspond to P-wave velocities
  COLUMN07='ISRF'  COLUMN08='+VP'
  KOLSRF=8  VPER=-15  VREF=1.5  CREF=0.6666667
  PROPERTIES='N1 N2 N3 ISRF VP'  # just for GOCAD

# Triangulating and plotting the velocity section
  O1=0.   O2=-1.55  O3=0.
  D1=0.1  D2=0.05   D3=10.
  N1=1    N2=63     N3=639
  NAME='VelocitySection'  # just for GOCAD
  VRTX='vrtx.tmp'  TRGL='trgl.tmp'  PLGN='plgn.tmp'
  modsrf:
  trgl:
  KOOR=2 TOCART=1
  PTS='vrtx.tmp' PTSOUT='vrtx.out'  coorchg:
  PTS=           PTSOUT=
  VRTX='vrtx.out'  TRGL='trgl.tmp'  srfwrl: