# History file '98-mod.h' to convert model '98' into VRML or GOCAD
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Input files required
  chk.pl: "data/98/" "98-mod.dat"
  chk.pl: "data/98/" "98-cam1.dat"
  chk.pl: "data/98/" "98-cam2.dat"
  chk.pl: "data/98/" "98-dli1.dat"
  chk.pl: "data/98/" "98-dli2.dat"
  chk.pl: "forms/"   "hsv.dat"

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

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

# Temporary files (several times overwritten)
  VRTX='vrtx.tmp'  TRGL='trgl.tmp'  PLGN='plgn.tmp'

# Setting of cameras (viewpoints) and lights, only VRML
  UP1=0  UP2=0  UP3=-1
  CAMERA='98-cam1.dat'
  DLIGHT='98-dli1.dat'  AMBIENT=0.00
  R=0.00  G=0.00  B=0.00  # Black background

# General colour properties
  COLORS='hsv.dat'
  TRANSP=0.00  SPECULAR=0.00  SHININESS=0.10

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

# Colours correspond to surface indices
  COLUMN07='ISRF'

# Triangulating and plotting interfaces
  N1=26  N2=26  N3=21  D1=0.2  D2=0.2  D3=0.2  O1=0.  O2=0.  O3=0.
  modsrf:
  trgl:
  KOLUMN=7
  VALUE=1  VRTXN='98-pts1.out'  TRGLN='98-trg1.out'  trglsort:
  VALUE=2  VRTXN='98-pts2.out'  TRGLN='98-trg2.out'  trglsort:
  VALUE=3  VRTXN='98-pts3.out'  TRGLN='98-trg3.out'  trglsort:
  KOLSRF=7  VPER=6  VREF=1  CREF=0.
  PROPERTIES='N1 N2 N3 ISRF'                       # just for GOCAD
  NAME='Surface1'                                  # just for GOCAD
  VRTX='98-pts1.out'  TRGL='98-trg1.out'  srfwrl:
  NAME='Surface2'                                  # just for GOCAD
  VRTX='98-pts2.out'  TRGL='98-trg2.out'  srfwrl:
  NAME='Surface3'                                  # just for GOCAD
  VRTX='98-pts3.out'  TRGL='98-trg3.out'  srfwrl:

# Adding model box wireframe
  LIN='98-bnd.out'
  KOLLIN=0  R=1.00  G=1.00  B=1.00
  NAME='ModelBox'  PROPERTIES=' '  # just for GOCAD
  bndlin:
  linwrl:

# Solid P-wave velocity model
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Setting of cameras and lights
  CAMERA='98-cam2.dat'  DLIGHT='98-dli2.dat'
  R=0.00  G=0.00  B=0.00  # Black background

# Name of output file
  WRL='98-vp.wrl'  # VRML
 #WRL='98-vp.mx'   # GOCAD
  iniwrl:
#
# Colours correspond to P-wave velocities
  COLUMN07='ISRF'  COLUMN08='+VP'
  KOLSRF=08  VPER=-6  VREF=1.5  CREF=0.666667
  PROPERTIES='N1 N2 N3 ISRF VP'               # just for GOCAD

# Triangulating individual velocity sections
# Section perpendicular to X3, minimum X3 coordinate, part1
  N1=2   N2=2   N3=1   D1=3.   D2=5.   D3=0.2  O1=0.   O2=0.   O3=0.
  NAME='Section3Min1'  # just for GOCAD
  modsrf:
  trgl:
  srfwrl:
# Section perpendicular to X3, minimum X3 coordinate, part2
                       D1=2.   D2=2.5          O1=3.
  NAME='Section3Min2'  # just for GOCAD
  modsrf:
  trgl:
  srfwrl:
# Section perpendicular to X3, maximum X3 coordinate, part1
                       D1=3.   D2=5.           O1=0.           O3=4.
  NAME='Section3Max1'  # just for GOCAD
  modsrf:
  trgl:
  srfwrl:
# Section perpendicular to X3, maximum X3 coordinate, part2
                       D1=2.   D2=2.5          O1=3.
  NAME='Section3Max2'  # just for GOCAD
  modsrf:
  trgl:
  srfwrl:
# Section perpendicular to X1, minimum X1 coordinate
  N1=1          N3=2           D2=5.   D3=-4.  O1=0.
  NAME='Section1Min'  # just for GOCAD
  modsrf:
  trgl:
  srfwrl:
# Section perpendicular to X1, maximum X1 coordinate
                               D2=2.5          O1=5.
  NAME='Section1Max'  # just for GOCAD
  modsrf:
  trgl:
  srfwrl:
# Section perpendicular to X1, inside model box
         N2=11  N3=21          D2=0.25 D3=-0.2 O1=3.   O2=2.5
  NAME='Section1Mid'  # just for GOCAD
  modsrf:
  trgl:
  srfwrl:
# Section perpendicular to X2, minimum X2 coordinate
  N1=26  N2=1   N3=2   D1=0.2          D3=-4.  O1=0.   O2=0.
  NAME='Section2Min'  # just for GOCAD
  modsrf:
  trgl:
  srfwrl:
# Section perpendicular to X2, maximum X2 coordinate
  N1=11                D1=0.3                          O2=5.
  NAME='Section2Max'  # just for GOCAD
  modsrf:
  trgl:
  srfwrl:
# Section perpendicular to X2, inside model box
                N3=21  D1=0.2          D3=-0.2 O1=3.   O2=2.5
  NAME='Section2Mid'  # just for GOCAD
  modsrf:
  trgl:
  srfwrl:

# 3D P-wave velocity cube (GOCAD Voxet), only GOCAD
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Name of output file
  WRL='98-vpv.mx'  # GOCAD
 #iniwrl:

# Data for gridding
  N1=126  N2=126  N3=101  D1=0.04  D2=0.04  D3=0.04  # sparser grid
# N1=251  N2=251  N3=201  D1=0.02  D2=0.02  D3=0.02  # minimum MRAM=13000000
  O1=0.0  O2=0.0  O3=0.0

  VEL='98-vpv.out'
 #grid:  # Gridding P velocities
  GRD='98-vpv.out' IN='98-vpv.bin'
 #ascbin:
 #swap:  # Gridded data for GOCAD must be big-endian (workstations)
         # Comment the above line on big-endian computers (workstations)
  NAME='VelocityCube'              # GOCAD object (and property) name
  PROPERTY='Velocity'              # GOCAD property name equal to NAME
  VPER=-6  VREF=1.5  CREF=0.666667 # Colour scale
 #grdwrl: