# A cookbook for using adiabatic conditions and a gravity 
# profile generated by a mineral physics toolkit and read in
# through an ascii data file. In this example, material properties
# are based on a Birch-Murnaghan equation of state. 

set Dimension                              = 2
set Use years in output instead of seconds = true
set End time                               = 2.63e8
set Output directory                       = output-burnman

subsection Material model
  set Model name = ascii reference profile

  subsection Ascii reference profile
    set Thermal viscosity exponent = 10.0
    set Viscosity prefactors = 1.0, 0.1, 1.0, 10.0

    subsection Ascii data model
      set Data directory = $ASPECT_SOURCE_DIR/data/adiabatic-conditions/ascii-data/
      set Data file name = example_isentrope.txt
    end
  end

  set Material averaging = harmonic average
end

# The geometry is a spherical shell with the inner and 
# outer radius of the Earth's mantle. 
subsection Geometry model
  set Model name = spherical shell

  subsection Spherical shell
    set Inner radius  = 3481000
    set Outer radius  = 6336000
  end
end

subsection Adiabatic conditions model
  set Model name = ascii data

  subsection Ascii data model
    set Data directory = $ASPECT_SOURCE_DIR/data/adiabatic-conditions/ascii-data/
    set Data file name = example_isentrope.txt
  end
end


# The gravity model reads an ascii data file that contains a gravity
# profile consistent with the material properties computed using the 
# Birch-Murnaghan equation of state. It automatically uses the same file 
# specified in the adiabatic conditions model.
subsection Gravity model
  set Model name = ascii data
  subsection Ascii data model
    set Data directory = $ASPECT_SOURCE_DIR/data/adiabatic-conditions/ascii-data/
    set Data file name = example_isentrope.txt
  end
end

# The model uses the anelastic liquid approximation.
subsection Formulation
  set Formulation = anelastic liquid approximation
end

# The present-day plate velocities are imposed on the upper model
# boundary.
subsection Boundary velocity model
  set Prescribed velocity boundary indicators = top:gplates
  set Tangential velocity boundary indicators = bottom

  subsection GPlates model
    set Data directory = $ASPECT_SOURCE_DIR/data/boundary-velocity/gplates/
    set Velocity file name = current_day.gpml
    set Data file time step = 1e6
    set Point one = 1.5708,4.87
    set Point two = 1.5708,5.24
  end
end


subsection Boundary temperature model
  set Fixed temperature boundary indicators = top, bottom
  set List of model names = spherical constant

  subsection Spherical constant
    set Inner temperature = 3500 
    set Outer temperature = 273
  end
end

# The initial temperature is an adiabatic profile, which is
# taken from the adabatic conditions model. In this case, this is
# is the ascii data model that loads a file containing profiles
# corresponding to the Birch-Murnaghan equation of state.
# In addition to the adiabatic profile, there are thermal boundary 
# layers at the sirface and the core-mantle boundary.
subsection Initial temperature model
  set Model name = adiabatic

  subsection Adiabatic
    set Age top boundary layer = 5e7
  end
end


subsection Mesh refinement
  set Refinement fraction                = 0.4
  set Coarsening fraction                = 0.05
  set Initial adaptive refinement        = 1
  set Initial global refinement          = 5
  set Strategy                           = temperature
  set Time steps between mesh refinement = 5
end


subsection Heating model
  set List of model names = adiabatic heating, shear heating

  subsection Adiabatic heating
    set Use simplified adiabatic heating = true
  end
end

subsection Postprocess
  set List of postprocessors = visualization, velocity statistics, temperature statistics, heat flux statistics, depth average

  subsection Visualization
    set Output format                 = vtu
    set Time between graphical output = 0
    set Number of grouped files       = 0
    set List of output variables      = material properties, gravity, adiabat, nonadiabatic temperature, nonadiabatic pressure
  end

  subsection Depth average
    set Time between graphical output = 1e6
  end
end

subsection Checkpointing
  set Time between checkpoint  = 1800
end

subsection Solver parameters
  subsection Stokes solver parameters
    set Linear solver tolerance = 1e-6
  end
end