# A simple setup for convection in a 3d spherical shell.

set Dimension                              = 2
set Use years in output instead of seconds = true
set Start time                             = 0
set End time                               = 1.0e10 #0.236196
set CFL number                             = 1.0
set Timing output frequency                = 10
set Resume computation                     = false
set Use conduction timestep                = false
set Adiabatic surface temperature          = 1000.0 #0.5
set Pressure normalization                 = surface
set Surface pressure                       = 0.0
#set Maximum time step                      = 3.99e7
set Output directory                       = output_dimensional

subsection Solver parameters
  set Temperature solver tolerance = 1e-12
  subsection Advection solver parameters
    set GMRES solver restart length = 50
  end
  subsection Stokes solver parameters
    set Linear solver tolerance     = 1e-3
    set Number of cheap Stokes solver steps = 200
    set Maximum number of expensive Stokes solver steps = 50000
    set GMRES solver restart length = 50
  end
end

subsection Checkpointing
  set Steps between checkpoint = 50
end

subsection Formulation
  set Formulation          = Boussinesq approximation
#  set Mass conservation    = incompressible
#  set Temperature equation = reference density profile
end

subsection Adiabatic conditions model
  set Model name = function
  subsection Function
    set Function constants = 
    set Variable names = depth
    set Function expression = 298.15;9.81*depth;35000.0 #0.5;2770245845.206189*depth;1.0 
  end 
end


subsection Material model
  set Model name = depth dependent
  subsection Depth dependent model
    set Base model = simple
    set Depth dependence method = File
    set Data directory = /home/egrant93/newriver/Hackathon_2019/aspect/data/material-model/depth-dependent-visc/
    set Viscosity depth file = venus_like_dim.txt  
  end
  subsection Simple model
    #set Composition viscosity prefactor     = 1.0
    #set Density differential for compositional field 1 = 0
    set Maximum thermal prefactor           = 1.0e6
    set Minimum thermal prefactor           = 1.0e-6
    set Reference density                   = 3300.0 #1.0   
    set Reference specific heat             = 1250.0 #1.0
    set Reference temperature               = 1600.0 #0.5
    set Thermal conductivity                = 1.0    #1.0
    set Thermal expansion coefficient       = 2.0e-5 #1.0
    #set Thermal viscosity exponent          = 1.0    #5.756465 
    set Viscosity                           = 1.0 #1.0
  end
end

subsection Discretization
  set Stokes velocity polynomial degree = 2
  set Temperature polynomial degree     = 2

  subsection Stabilization parameters
#    set Stabilization method = SUPG
    set alpha = 2
    set beta  = 0.0
    set cR    = 0.0
  end
end

subsection Geometry model
  set Model name = spherical shell
  subsection Spherical shell
    set Inner radius  = 3450000.0 #0.514
    set Outer radius  = 6371000.0 #1.0
    set Opening angle = 360
    set Cells along circumference = 12
  end
end

subsection Boundary velocity model
  set Tangential velocity boundary indicators = top
  set Prescribed velocity boundary indicators =
  set Zero velocity boundary indicators       = bottom
end

subsection Boundary temperature model
  set Fixed temperature boundary indicators   = top, bottom
  set Model name = spherical constant
  subsection Spherical constant
    set Inner temperature = 2000.0 #1.0
    set Outer temperature = 290.0
  end
end

subsection Heating model
  set List of model names = constant heating
  subsection Constant heating
    set Radiogenic heating rate = 406.442 #4e-23 #96 #406.4421074  # 96 * (R/d)^2
  end
end

subsection Melt settings
  set Include melt transport                  = false
end

subsection Nullspace removal
  set Remove nullspace = net rotation
end

subsection Initial temperature model
  set Model name = harmonic perturbation
  subsection Harmonic perturbation
    set Vertical wave number = 1
    set Lateral wave number one = 8
    set Lateral wave number two = 4
    set Magnitude = 0.02
    #set Reference temperature = 2000.0 #0.5
  end
end

subsection Gravity model
  set Model name = radial constant
  subsection Radial constant
    set Magnitude = 9.81 #2770245845.206189
   # set Magnitude = 3.18e8
  end
end

subsection Termination criteria
  set Termination criteria = end time
end

subsection Mesh refinement
  set Initial global refinement          = 5
  set Initial adaptive refinement        = 0
  set Time steps between mesh refinement = 0
end


subsection Postprocess
  set List of postprocessors = visualization, velocity statistics, temperature statistics, heat flux statistics, heat flux densities, depth average, basic statistics, entropy viscosity statistics, material statistics

  subsection Visualization
    set List of output variables            = viscosity, density
    set Output format                       = vtu
    #set Time between graphical output       = 1000.0 #0.0118098
    set Time steps between graphical output = 25
    set Number of grouped files             = 1
  end

  subsection Depth average
    set Number of zones               = 96
    set Time between graphical output = 1000.0 #9.44784e-03
    set Output format                 = txt
  end
end