# Marmousi velocity model
# ~~~~~~~~~~~~~~~~~~~~~~~
# History file to convert original binary velocities 'velocity.h@'
# gridded with intervals 4m*4m into ascii velocities 'mar-vel.dat'
# gridded with intervals 40m*80m.  The slowness is averaged over grid
# cells 40m*80m.
# Since both original binary files 'velocity.h@' (nearly 7MB for each
# little and big endian) and the corresponding ascii file 'velocity.out'
# (nearly 14MB) are very large, only file 'mar-vel.dat' is stored in
# subdirectory 'mar' of package DATA for further processing.
# Binary file 'velocity.h@' of the Marmousi model and data set is
# required to execute this history file.

  n1=751      n2=2301     n3=1
  o1=0.0      o2=0.0      o3=0.0
  d1=4.0      d2=4.0      d3=1.0
  esize=4
  data_format="xdr_float"
  in="velocity.h@"

# Conversion of gridded velocity into an ascii formatted file
  GRD='velocity.out'  NDIG=0
  binasc: 'mar-vel.h'  /

# Checking the velocity conversion (disabled by #)
  GRD='velocity.out'  IN='velocity.tmp'
 #ascbin: 'mar-vel.h'  /
# Binary files 'velocity.h@' and 'velocity.tmp' should be identical.

# Calculating slowness on the dense grid
  CAL='inv.cal' GRD1='velocity.out' GRD2='slowness.out'
  grdcal: 'mar-vel.h' /

# Plotting velocity and slowness on the dense grid
  HSIZE=7.51  VSIZE=23.01
  VCIRC=-4800  VREF=1500  CREF=0.666667
# 1500 blue, 2300 cyan, 3100 green, 3900 yellow, 4700 red, 5500 magenta
  GRD='velocity.out' PS='velocity.ps'
  grdps: 'mar-vel.h' /
  VCIRC=.000581818  VREF=.000666667  CREF=0.666667  # 1/1500 blue,
# 1/1755 cyan, 1/2115 green, 1/2661 yellow, 1/3587 red, 1/5500 magenta
  GRD='slowness.out' PS='slowness.ps'
  grdps: 'mar-vel.h' /

# Average slowness on a coarse grid and the corresponding velocity
  N1NEW=76  N2NEW=116  D1NEW=40  D2NEW=80  O1NEW=0  O2NEW=0  GNORM=1
  GRD='slowness.out' GRDNEW='mar-slo.out'
  grdnorm: 'mar-vel.h' /
  N1=76  N2=116  D1=40  D2=80  O1=0  O2=0
  CAL='inv.cal' GRD1='mar-slo.out' GRD2='mar-vel.dat'
  grdcal: 'mar-vel.h' /