# Created by cfggui.py 20170412 at 14:56:38
# Description of this simulation
description = fast19
# Number of Tracers along X (width)
NumberOfTracersAlongX = 2000
# Number of Tracers along Y (depth)
NumberOfTracersAlongY = 800
# Viscosity limits for rocks, Pa
etamin = 1e17 pascal * second # Lower limit, Pa
# Viscosity maximum limit for rocks
etamax = 1.0e25 pascal * second # Upper limit, Pa
# Lower stress limit for power law, Pa
stressmin = 10.0e3 pascal
# Viscoelastic timestep, years
timemax = 1.0e4 year
# timemax = 30.0e3 * year
# Maximal tracers displacement step, number of gridsteps
tracerMaxMove = 0.3
# Moving Tracers:
tracerMoveAlg = runge-kutta
# Maximal temperature change, allowed for one timestep, K
tempmax = 20.0 kelvin
# Number of timesteps
stepmax = 4001
# Ouput file
output_file = fast19.hdf5
# Log file
log_file = fast19.log
num_threads = 1
# Visco-elasto-plastic iteration
vep_iteration = local
# Maximum number of Visco-elasto-plastic iterations
vep_max_iter = 10
# Tracers interpolation to image
tracers2image = False
# Temperature solver options
# temperatureSolver = '''
# ksp_type=preonly
# pc_type=lu
# pc_factor_mat_solver_package=superlu_dist
# mat_superlu_dist_colperm=PARMETIS
# mat_superlu_dist_parsymbfact=1
# '''
temperatureSolver = '''
ksp_type=preonly
pc_type=lu
pc_factor_mat_solver_package=mumps
mat_mumps_icntl_14=70
'''
#temperatureSolver = '''
#ksp_type=preonly
#pc_type=lu
#pc_factor_mat_solver_package=mkl_pardiso
#'''
# Stokes solver options
# stokesSolver = '''
# ksp_type=preonly
# pc_type=lu
# pc_factor_mat_solver_package=superlu_dist
# mat_superlu_dist_colperm=PARMETIS
# mat_superlu_dist_parsymbfact=1
# '''
stokesSolver = '''
ksp_type=preonly
pc_type=lu
pc_factor_mat_solver_package=mumps
mat_mumps_icntl_14=70
'''
#stokesSolver = '''
#ksp_type=preonly
#pc_type=lu
#pc_factor_mat_solver_package=mkl_pardiso
#'''
# Gas Constant
RGAS = 8.314 joule / kelvin / mole
# Universal gravity constant
UG = 66.72e-12 meter ** 2 * newton / kilogram ** 2
# Acceleration in x
gx = 0.0 meter / second ** 2
# Acceleration in y
gy = 9.81 meter / second ** 2
# Save every `save_every` evolution steps
save_every = 50
# Starting step
# step0 = 1331
# Use numerical Stress subgrid diffusion
stress_subgrid = True
# Numerical Stress subgrid value
stress_subgrid_value = 1.0 dimensionless
# Use numerical Temperature subgrid diffusion
temperature_subgrid = True
# Numerical Temperature subgrid value
temperature_subgrid_value = 1.0 dimensionless
# Apply shear heating
shear_heating = True
# Apply adiabatic heating
adiabatic_heating = True
# Pressure boundary conditions:
# onecell: pressure in one cell definition
# topbottom: pressure at the top and in the bottom of the channel
pressure_bc = onecell
# Cell boundary condition pressure
pressure_value = 100.0e3 pascal
# Maximum number of iterations for temperature
temp_iter_max = 100
# External file with Lagriangian Tracers
tracersFile = ""
# Starting step
# step0 = 501
# External configuration files blank separated
configFiles = ""
# Phase change switch
phase_change = False
[Melting]
# Enable melting
enabled = False
model = gerya
# # Polyon defining melting region
# polygon = ""
# # Melting fraction
# fraction = 0.25
# # Base rock lithology
# lithology = NoLithology
[Mesh]
[[X]]
width = 3000.0 kilometer
nodes = 201
[[[Distribution]]]
[[[[Variable]]]]
end = 800.0 kilometer
step = -8.0 kilometer
nodes = 25
# Starting point abscissa
start = -1.0 meter
[[[[Constant]]]]
end = 2000.0 kilometer
nodes = 150
# Starting point abscissa
start = -1.0 meter
# Step width/depth (0 for constant distribution, <0 for right to left)
step = 0.0 meter
[[[[Variable1]]]]
step = 8.0 kilometer
nodes = 25
# Starting point abscissa
start = -1.0 meter
# Ending point abscissa
end = -1.0 meter
[[Y]]
depth = 400.0 kilometer
nodes = 81
[[[Distribution]]]
[[[[Constant]]]]
end = 200.0 kilometer
nodes = 50
# Starting point abscissa
start = -1.0 meter
# Step width/depth (0 for constant distribution, <0 for right to left)
step = 0.0 meter
[[[[Variable]]]]
step = 4.0 kilometer
nodes = 30
# Starting point abscissa
start = -1.0 meter
# Ending point abscissa
end = -1.0 meter
[[Reference Point]]
# Abscissa
x = 0.0 meter
# Depth
y = 0.0 meter
# Lithologies at the reference point
lithologies = lithoalias1, lithoalias2
[Lithologies]
[[air]]
# AD
AD = 0.0 1 / pascal / second
# a
a = 0.0 watt / meter
# Layer alias
alias = air
# Cohesion 0
cohesion_0 = 0.0 pascal
# Cohesion 1
cohesion_1 = 0.0 pascal
# Compressibility
compressibility = 10.0e-12 1 / pascal
# Standard density
density = 1000.0 kilogram / meter ** 3
# Ea
Ea = 0.0 joule / mole
# Epsilon 0
epsilon_0 = 0.0 dimensionless
# Epsilon 1
epsilon_1 = 1.0 dimensionless
# Heat Capacity
heat_capacity = 3000.0 joule / kilogram
# Heterogeneity params: mean, std, min, max
heterogeneity = 0, 0, 0, 0
# Radiogenic heat production
hr = 0.0 watt / meter ** 3
# k0
k0 = 300.0 watt / kelvin / meter
# Melt density
melt_density = 0.0 kilogram / meter ** 3
# Melt Viscosity
melt_viscosity = 0.0 pascal * second
# n
n = 0.0 dimensionless
# Apply Peierls
peierls = False
# Phi 0
phi_0 = 0.0 degree
# Phi 1
phi_1 = 0.0 degree
# Use power law
power_law = False
# Rock type
rock_type = none
# Scaling Factor
scaling_factor = 0.0 dimensionless
# Scaling Factor Equation
scaling_factor_eq = none
# Shear modulus
shear_modulus = 100.0e18 pascal
# This equation controls melting
SolidusEquation = 0 dimensionless
# Thermal expansion
thermal_expansion = 30.0e-6 1 / kelvin
# Multiplyer factor for units in polygons
units_multiplyer = 1.0 kilometer
# Va
Va = 0.0 meter ** 3
# Viscosity
viscosity = 1.0e18 pascal * second
# Polygons defined by their vertices, one per line,
# separated by nl
polygons = '''
0.0 0.0
0.0 7.0
1450.0 7.0
1500.0 7.0
3000.0 7.0
3000.0 0.0'''
# Phase change file
phase_change = ""
# Alias rock of wet status
wet = ""
# Alias rock of dry status
dry = ""
# Pore fluid pressure factor
lambda = 1.0 dimensionless
# Fluid model for the layer
fluid = none
[[sediments]]
# Phi 0
phi_0 = 1.7 degree
# Phi 1
phi_1 = 1.7 degree
# Epsilon 1
epsilon_1 = 1.0 dimensionless
# Standard density
density = 2600.0 kilogram / meter ** 3
# AD
AD = 320.0e-6 1 / pascal / second
# a
a = 807.0 watt / meter
# Layer alias
alias = sediments
# Cohesion 0
cohesion_0 = 1.0e6 pascal
# Cohesion 1
cohesion_1 = 1.0e6 pascal
# Compressibility
compressibility = 10.0e-12 1 / pascal
# Ea
Ea = 154.0e3 joule / mole
# Epsilon 0
epsilon_0 = 0.0 dimensionless
# Heat Capacity
heat_capacity = 1000.0 joule / kilogram
# Heterogeneity params: mean, std, min, max
heterogeneity = 0, 0, 0, 0
# Radiogenic heat production
hr = 2.0e-6 watt / meter ** 3
# k0
k0 = 640.0e-3 watt / kelvin / meter
# Melt density
melt_density = 2400.0 kilogram / meter ** 3
# Melt Viscosity
melt_viscosity = 0.0 pascal * second
# n
n = 2.3 dimensionless
phase_change = sed300.dat
fluid = wet
# Apply Peierls
peierls = False
# Use power law
power_law = True
# Rock type
rock_type = none
# Scaling Factor
scaling_factor = 0.0 dimensionless
# Scaling Factor Equation
scaling_factor_eq = none
# Shear modulus
shear_modulus = 10.0e9 pascal
# This equation controls melting
# Thermal expansion
thermal_expansion = 30.0e-6 1 / kelvin
# Multiplyer factor for units in polygons
units_multiplyer = 1.0 kilometer
# Va
Va = 0.0 meter ** 3
# Polygons defined by their vertices, one per line,
# separated by nl
polygons = '''
1450.0 7.0
3000.0 7.0
3000.0 11.0
1540.0 11.0
1530.0 13.0
1520.0 18.0
1500.0 18.0
1500.0 10.0'''
SolidusEquation = 1 dimensionless
# Viscosity
viscosity = 0.0 pascal * second
# Alias rock of wet status
wet = ""
# Alias rock of dry status
dry = ""
# Pore fluid pressure factor
lambda = 1.0 dimensionless
[[Oceanic Upper crust]]
# AD
AD = 320.0e-6 1 / pascal / second
# a
a = 474.0 watt / meter
# Layer alias
alias = oceanic upper crust
# Cohesion 0
cohesion_0 = 1.0e6 pascal
# Cohesion 1
cohesion_1 = 1.0e6 pascal
# Compressibility
compressibility = 10.0e-12 1 / pascal
# Standard density
density = 3000.0 kilogram / meter ** 3
# Ea
Ea = 154.0e3 joule / mole
# Epsilon 0
epsilon_0 = 0.0 dimensionless
# Epsilon 1
epsilon_1 = 1.0 dimensionless
# Heat Capacity
heat_capacity = 1000.0 joule / kilogram
# Heterogeneity params: mean, std, min, max
heterogeneity = 0, 0, 0, 0
# Radiogenic heat production
hr = 250.0e-9 watt / meter ** 3
# k0
k0 = 1.18 watt / kelvin / meter
# Melt density
melt_density = 2400.0 kilogram / meter ** 3
# Melt Viscosity
melt_viscosity = 0.0 pascal * second
# n
n = 2.3 dimensionless
# Apply Peierls
peierls = False
phase_change = morb300.dat
fluid = wet
# Phi 0
phi_0 = 1.7 degree
# Phi 1
phi_1 = 1.7 degree
# Use power law
power_law = True
# Rock type
rock_type = none
# Scaling Factor
scaling_factor = 0.0 dimensionless
# Scaling Factor Equation
scaling_factor_eq = none
# Shear modulus
shear_modulus = 25.0e9 pascal
# This equation controls melting
SolidusEquation = 2 dimensionless
# Thermal expansion
thermal_expansion = 30.0e-6 1 / kelvin
# Multiplyer factor for units in polygons
units_multiplyer = 1.0 kilometer
# Va
Va = 0.0 meter ** 3
# Viscosity
viscosity = 0.0 pascal * second
# Polygons defined by their vertices, one per line,
# separated by nl
polygons = '''
1530.0 13.0
1540.0 11.0
3000.0 11.0
3000.0 13.0'''
# Alias rock of wet status
wet = ""
# Alias rock of dry status
dry = ""
# Pore fluid pressure factor
lambda = 1.0 dimensionless
[[Oceanic Lower crust]]
# AD
AD = 330.0e-6 1 / pascal / second
# a
a = 474.0 watt / meter
# Layer alias
alias = oceanic lower crust
# Cohesion 0
cohesion_0 = 1.0e6 pascal
# Cohesion 1
cohesion_1 = 1.0e6 pascal
# Compressibility
compressibility = 10.0e-12 1 / pascal
# Standard density
density = 3200.0 kilogram / meter ** 3
# Ea
Ea = 238.0e3 joule / mole
# Epsilon 0
epsilon_0 = 0.0 dimensionless
# Epsilon 1
epsilon_1 = 1.0 dimensionless
# Heat Capacity
heat_capacity = 1000.0 joule / kilogram
# Heterogeneity params: mean, std, min, max
heterogeneity = 0, 0, 0, 0
# Radiogenic heat production
hr = 250.0e-9 watt / meter ** 3
# k0
k0 = 1.18 watt / kelvin / meter
# Melt density
melt_density = 2700.0 kilogram / meter ** 3
# Melt Viscosity
melt_viscosity = 0.0 pascal * second
# n
n = 3.2 dimensionless
# Apply Peierls
peierls = False
fluid = wet
# Phi 0
phi_0 = 11.5 degree
# Phi 1
phi_1 = 11.5 degree
# Use power law
power_law = True
# Rock type
rock_type = none
# Scaling Factor
scaling_factor = 0.0 dimensionless
# Scaling Factor Equation
scaling_factor_eq = none
# Shear modulus
shear_modulus = 25.0e9 pascal
# This equation controls melting
SolidusEquation = 2 dimensionless
# Thermal expansion
thermal_expansion = 30.0e-6 1 / kelvin
# Multiplyer factor for units in polygons
units_multiplyer = 1.0 kilometer
# Va
Va = 0.0 meter ** 3
# Viscosity
viscosity = 0.0 pascal * second
# Polygons defined by their vertices, one per line,
# separated by nl
polygons = '''
1520.0 18.0
1530.0 13.0
3000.0 13.0
3000.0 18.0'''
# Phase change file
phase_change = gabbro300.dat
# Alias rock of wet status
wet = ""
# Alias rock of dry status
dry = ""
# Pore fluid pressure factor
lambda = 1.0 dimensionless
[[Lid]]
# AD
AD = 25.0e3 1 / pascal / second
# a
a = 1293.0 watt / meter
# Layer alias
alias = lid
# Cohesion 0
cohesion_0 = 1.0e6 pascal
# Cohesion 1
cohesion_1 = 1.0e6 pascal
# Compressibility
compressibility = 10.0e-12 1 / pascal
# Standard density
density = 3300.0 kilogram / meter ** 3
# Ea
Ea = 532.0e3 joule / mole
# Epsilon 0
epsilon_0 = 0.0 dimensionless
# Epsilon 1
epsilon_1 = 1.0 dimensionless
# Heat Capacity
heat_capacity = 1000.0 joule / kilogram
# Heterogeneity params: mean, std, min, max
heterogeneity = 0, 0, 0, 0
# Radiogenic heat production
hr = 22.0e-9 watt / meter ** 3
# k0
k0 = 730.0e-3 watt / kelvin / meter
# Melt density
melt_density = 2700.0 kilogram / meter ** 3
# Melt Viscosity
melt_viscosity = 0.0 pascal * second
# n
n = 3.5 dimensionless
# Apply Peierls
peierls = False
# Phi 0
phi_0 = 37.0 degree
# Phi 1
phi_1 = 37.0 degree
# Use power law
power_law = True
# Rock type
rock_type = none
# Scaling Factor
scaling_factor = 0.0 dimensionless
# Scaling Factor Equation
scaling_factor_eq = none
# Shear modulus
shear_modulus = 67.0e9 pascal
# This equation controls melting
SolidusEquation = 3 dimensionless
# Thermal expansion
thermal_expansion = 30.0e-6 1 / kelvin
# Multiplyer factor for units in polygons
units_multiplyer = 1.0 kilometer
# Va
Va = 10.0 centimeter ** 3
# Viscosity
viscosity = 0.0 pascal * second
# Polygons defined by their vertices, one per line,
# separated by nl
polygons = '''
0.0 42.0
1462.0 42.0
1350.0 137.0
0.0 137.0
1550.0 18.0
3000.0 18.0
3000.0 111.0
1520.0 111.0
1450.0 137.0
1365.0 137.0
1520.0 18.0'''
wet = hmantle
fluid = dry/wet
# Phase change file
phase_change = peridotite300.dat
# Alias rock of dry status
dry = ""
# Pore fluid pressure factor
lambda = 1.0 dimensionless
[[Asthenosphere]]
# AD
AD = 25.0e3 1 / pascal / second
# a
a = 1293.0 watt / meter
# Layer alias
alias = asthenosphere
# Cohesion 0
cohesion_0 = 1.0e6 pascal
# Cohesion 1
cohesion_1 = 1.0e6 pascal
# Compressibility
compressibility = 10.0e-12 1 / pascal
# Standard density
density = 3300.0 kilogram / meter ** 3
# Ea
Ea = 532.0e3 joule / mole
# Epsilon 0
epsilon_0 = 0.0 dimensionless
# Epsilon 1
epsilon_1 = 1.0 dimensionless
# Heat Capacity
heat_capacity = 1000.0 joule / kilogram
# Heterogeneity params: mean, std, min, max
heterogeneity = 0, 0, 0, 0
# Radiogenic heat production
hr = 22.0e-9 watt / meter ** 3
# k0
k0 = 730.0e-3 watt / kelvin / meter
# Melt density
melt_density = 2700.0 kilogram / meter ** 3
# Melt Viscosity
melt_viscosity = 0.0 pascal * second
# n
n = 3.5 dimensionless
# Apply Peierls
peierls = False
fluid = dry/wet
# Phi 0
phi_0 = 37.0 degree
# Phi 1
phi_1 = 37.0 degree
# Use power law
power_law = True
# Rock type
rock_type = none
# Scaling Factor
scaling_factor = 0.0 dimensionless
# Scaling Factor Equation
scaling_factor_eq = none
# Shear modulus
shear_modulus = 67.0e9 pascal
# This equation controls melting
SolidusEquation = 3 dimensionless
# Thermal expansion
thermal_expansion = 30.0e-6 1 / kelvin
# Multiplyer factor for units in polygons
units_multiplyer = 1.0 kilometer
# Va
Va = 10.0 centimeter ** 3
# Viscosity
viscosity = 0.0 pascal * second
# Polygons defined by their vertices, one per line,
# separated by nl
polygons = '''
0.0 137.0
1350.0 137.0
1450.0 137.0
1520.0 111.0
3000.0 111.0
3000.0 400.0
0.0 400.0'''
# Phase change file
phase_change = peridotite300.dat
# Alias rock of wet status
wet = hmantle
# Alias rock of dry status
dry = ""
# Pore fluid pressure factor
lambda = 1.0 dimensionless
[[Weak Zone]]
# AD
AD = 2000.0 1 / pascal / second
# a
a = 1293.0 watt / meter
# Layer alias
alias = weak
# Cohesion 0
cohesion_0 = 1.0e6 pascal
# Cohesion 1
cohesion_1 = 1.0e6 pascal
# Compressibility
compressibility = 10.0e-12 1 / pascal
# Standard density
density = 3200.0 kilogram / meter ** 3
# Ea
Ea = 471.0e3 joule / mole
# Epsilon 0
# Epsilon 0
epsilon_0 = 0.0 dimensionless
# Epsilon 1
epsilon_1 = 1.0 dimensionless
# Heat Capacity
heat_capacity = 1000.0 joule / kilogram
# Heterogeneity params: mean, std, min, max
heterogeneity = 0, 0, 0, 0
# Radiogenic heat production
hr = 22.0e-9 watt / meter ** 3
# k0
k0 = 730.0e-3 watt / kelvin / meter
# Melt density
melt_density = 2700.0 kilogram / meter ** 3
# Melt Viscosity
melt_viscosity = 0.0 pascal * second
# n
n = 4.0 dimensionless
# Apply Peierls
peierls = False
# Phi 0
phi_0 = 0.0 degree
# Phi 1
phi_1 = 0.0 degree
# Use power law
power_law = True
# Rock type
rock_type = none
# Scaling Factor
scaling_factor = 0.0 dimensionless
# Scaling Factor Equation
scaling_factor_eq = none
# Shear modulus
shear_modulus = 67.0e9 pascal
# This equation controls melting
SolidusEquation = 0 dimensionless
# Thermal expansion
thermal_expansion = 30.0e-6 1 / kelvin
# Multiplyer factor for units in polygons
units_multiplyer = 1.0 kilometer
# Va
Va = 0.0 meter ** 3
# Viscosity
viscosity = 0.0 pascal * second
# Polygons defined by their vertices, one per line,
# separated by nl
polygons = '''
1520.0 18.0
1500.0 18.0
1462.0 42.0
1350.0 137.0
1365.0 137.0'''
# Phase change file
phase_change = ""
# Alias rock of wet status
wet = ""
# Alias rock of dry status
dry = ""
# Pore fluid pressure factor
lambda = 1.0 dimensionless
# Fluid model for the layer
fluid = none
[[Continental Upper crust]]
# AD
AD = 320.0e-6 1 / pascal / second
# a
a = 807.0 watt / meter
# Layer alias
alias = continental upper crust
# Cohesion 0
cohesion_0 = 1.0e6 pascal
# Cohesion 1
cohesion_1 = 1.0e6 pascal
# Compressibility
compressibility = 10.0e-12 1 / pascal
# Standard density
density = 2700.0 kilogram / meter ** 3
# Ea
Ea = 154.0e3 joule / mole
# Epsilon 0
epsilon_0 = 0.0 dimensionless
# Epsilon 1
epsilon_1 = 1.0 dimensionless
# Heat Capacity
heat_capacity = 1000.0 joule / kilogram
# Heterogeneity params: mean, std, min, max
heterogeneity = 0, 0, 0, 0
# Radiogenic heat production
hr = 1.0e-6 watt / meter ** 3
# k0
k0 = 640.0e-3 watt / kelvin / meter
# Melt density
melt_density = 2400.0 kilogram / meter ** 3
# Melt Viscosity
melt_viscosity = 0.0 pascal * second
# n
n = 2.3 dimensionless
# Apply Peierls
peierls = False
# Phi 0
phi_0 = 11.5 degree
# Phi 1
phi_1 = 11.5 degree
# Use power law
power_law = True
# Rock type
rock_type = none
# Scaling Factor
scaling_factor = 0.0 dimensionless
# Scaling Factor Equation
scaling_factor_eq = none
# Shear modulus
shear_modulus = 10.0e9 pascal
# This equation controls melting
SolidusEquation = 1 dimensionless
# Thermal expansion
thermal_expansion = 30.0e-6 1 / kelvin
# Multiplyer factor for units in polygons
units_multiplyer = 1.0 kilometer
# Va
Va = 0.0 meter ** 3
# Viscosity
viscosity = 0.0 pascal * second
# Polygons defined by their vertices, one per line,
# separated by nl
polygons = '''
0.0 7.0
0.0 27.0
1450.0 27.0
1500.0 13.0
1500.0 10.0
1450.0 7.0'''
# Phase change file
phase_change = ""
# Alias rock of wet status
wet = ""
# Alias rock of dry status
dry = ""
# Pore fluid pressure factor
lambda = 1.0 dimensionless
# Fluid model for the layer
fluid = none
[[Continental Lower crust]]
# AD
AD = 330.0e-6 1 / pascal / second
# a
a = 474.0 watt / meter
# Layer alias
alias = continental lower crust
# Cohesion 0
cohesion_0 = 1.0e6 pascal
# Cohesion 1
cohesion_1 = 1.0e6 pascal
# Compressibility
compressibility = 10.0e-12 1 / pascal
# Standard density
density = 2900.0 kilogram / meter ** 3
# Ea
Ea = 238.0e3 joule / mole
# Epsilon 0
epsilon_0 = 0.0 dimensionless
# Epsilon 1
epsilon_1 = 1.0 dimensionless
# Heat Capacity
heat_capacity = 1000.0 joule / kilogram
# Heterogeneity params: mean, std, min, max
heterogeneity = 0, 0, 0, 0
# Radiogenic heat production
hr = 500.0e-9 watt / meter ** 3
# k0
k0 = 1.18 watt / kelvin / meter
# Melt density
melt_density = 2700.0 kilogram / meter ** 3
# Melt Viscosity
melt_viscosity = 0.0 pascal * second
# n
n = 3.2 dimensionless
# Apply Peierls
peierls = False
# Phi 0
phi_0 = 11.5 degree
# Phi 1
phi_1 = 11.5 degree
# Use power law
power_law = True
# Rock type
rock_type = none
# Scaling Factor
scaling_factor = 0.0 dimensionless
# Scaling Factor Equation
scaling_factor_eq = none
# Shear modulus
shear_modulus = 25.0e9 pascal
# This equation controls melting
SolidusEquation = 2 dimensionless
# Thermal expansion
thermal_expansion = 30.0e-6 1 / kelvin
# Multiplyer factor for units in polygons
units_multiplyer = 1.0 kilometer
# Va
Va = 0.0 meter ** 3
# Viscosity
viscosity = 0.0 pascal * second
# Polygons defined by their vertices, one per line,
# separated by nl
polygons = '''
0.0 27.0
0.0 42.0
1462.0 42.0
1500.0 18.0
1500.0 13.0
1450.0 27.0'''
# Phase change file
phase_change = ""
# Alias rock of wet status
wet = ""
# Alias rock of dry status
dry = ""
# Pore fluid pressure factor
lambda = 1.0 dimensionless
# Fluid model for the layer
fluid = none
[[Hydrated mantle]]
# AD
AD = 2000.0 1 / pascal / second
# a
a = 1293.0 watt / meter
# Layer alias
alias = hmantle
# Cohesion 0
cohesion_0 = 1.0e6 pascal
# Cohesion 1
cohesion_1 = 1.0e6 pascal
# Compressibility
compressibility = 10.0e-12 1 / pascal
# Standard density
density = 3200.0 kilogram / meter ** 3
# Ea
Ea = 471.0e3 joule / mole
# Epsilon 0
# Epsilon 0
epsilon_0 = 0.0 dimensionless
# Epsilon 1
epsilon_1 = 1.0 dimensionless
# Heat Capacity
heat_capacity = 1000.0 joule / kilogram
# Heterogeneity params: mean, std, min, max
heterogeneity = 0, 0, 0, 0
# Radiogenic heat production
hr = 22.0e-9 watt / meter ** 3
# k0
k0 = 730.0e-3 watt / kelvin / meter
# Melt density
melt_density = 2700.0 kilogram / meter ** 3
# Melt Viscosity
melt_viscosity = 0.0 pascal * second
# n
n = 4.0 dimensionless
# Apply Peierls
peierls = False
# Phi 0
phi_0 = 1.7 degree
# Phi 1
phi_1 = 1.7 degree
# Use power law
power_law = True
# Rock type
rock_type = none
# Scaling Factor
scaling_factor = 0.0 dimensionless
# Scaling Factor Equation
scaling_factor_eq = none
# Shear modulus
shear_modulus = 67.0e9 pascal
# This equation controls melting
SolidusEquation = 4 dimensionless
# Thermal expansion
thermal_expansion = 30.0e-6 1 / kelvin
# Multiplyer factor for units in polygons
units_multiplyer = 1.0 kilometer
# Va
Va = 0.0 meter ** 3
# Viscosity
viscosity = 0.0 pascal * second
# Phase change file
phase_change = peridotite300.dat
# Alias rock of wet status
wet = ""
# Alias rock of dry status
dry = ""
# Pore fluid pressure factor
lambda = 1.0 dimensionless
# Fluid model for the layer
fluid = wet
# Polygons defined by their vertices, one per line,
# separated by nl
polygons = ""
[Geothermic Model]
# Bottom temperature
T1 = 1465.0 degC
# Surface temperature
T0 = 0.0 degC
[[geotherm0]]
# Geothermic Model type
model = constant
# Constant layer temperature
value = 0.0 degC
# Polygon defined by their vertices, one per line, separated by nl
polygon = '''
0. 0.
0. 7.
1450. 7.
1520. 11.
3000. 11.
3000. 0.
'''
# Layer age
age = 0.0 year
# Thermal gradient
dtz = 0.0 kelvin / kilometer
# Thermal disturbance
dt = 0.0 kelvin
# Thermal diffusivity
kappa = 0.0 meter ** 2 / second
# Radius of anomaly
radius = 0.0 kilometer
# Direction of anomaly is right to left
rtol = False
# Surface temperature of layer
T0 = 0.0 degC
# Bottom temperature of layer
T1 = 0.0 degC
# List / tuple of temperature values with units enclosed in quotes
Ti = 0.0 degC, 100.0 degC
# Center position of thermal anomaly: x
x0 = 0.0 kilometer
# Center position of thermal anomaly: y
y0 = 0.0 kilometer
# Thermal Plate thickness
yL = 0.0 kilometer
# Radiogenic production
H = 0.0 meter ** 2 * watt
# Length scale of radiogenic production
hr = 0.0 kilometer
# Thermal conductivity
kd = 0.0 watt / kelvin / meter
# Surface heat flow
Qs = 0.0 watt / meter ** 2
# Min oceanic depth
ocemin = 0.0 meter
[[geotherm1]]
# Geothermic Model type
model = half-space
# Layer age
age = 60.0e6 year
# Thermal Plate thickness
yL = 111.0 kilometer
# Thermal diffusivity
kappa = 1.0e-6 meter ** 2 / second
# Surface temperature of layer
T0 = 0.0 degC
# Polygon defined by their vertices, one per line, separated by nl
polygon = '''
1520. 11.
1520. 111.
3000. 111.
3000. 11.
'''
# Thermal gradient
dtz = 0.0 kelvin / kilometer
# Thermal disturbance
dt = 0.0 kelvin
# Radius of anomaly
radius = 0.0 kilometer
# Direction of anomaly is right to left
rtol = False
# Bottom temperature of layer
T1 = 0.0 degC
# List / tuple of temperature values with units enclosed in quotes
Ti = 0.0 degC, 100.0 degC
# Constant layer temperature
value = 0.0 degC
# Center position of thermal anomaly: x
x0 = 0.0 kilometer
# Center position of thermal anomaly: y
y0 = 0.0 kilometer
# Radiogenic production
H = 0.0 meter ** 2 * watt
# Length scale of radiogenic production
hr = 0.0 kilometer
# Thermal conductivity
kd = 0.0 watt / kelvin / meter
# Surface heat flow
Qs = 0.0 watt / meter ** 2
# Min oceanic depth
ocemin = 0.0 meter
[[geotherm2]]
model = continent
polygon = '''
0. 7.
1450. 7.
1450. 42
0. 42.
'''
T0 = 0.0 degC
yL = 137.0 kilometer
# Layer age
age = 0.0 year
# Thermal gradient
dtz = 0.0 kelvin / kilometer
# Thermal disturbance
dt = 0.0 kelvin
# Thermal diffusivity
kappa = 0.0 meter ** 2 / second
# Radius of anomaly
radius = 0.0 kilometer
# Direction of anomaly is right to left
rtol = False
# Bottom temperature of layer
T1 = 500.0 degC
# List / tuple of temperature values with units enclosed in quotes
Ti = 0.0 degC, 100.0 degC
# Constant layer temperature
value = 0.0 degC
# Center position of thermal anomaly: x
x0 = 0.0 kilometer
# Center position of thermal anomaly: y
y0 = 0.0 kilometer
# Radiogenic production
H = 0.0 meter ** 2 * watt
# Length scale of radiogenic production
hr = 0.0 kilometer
# Thermal conductivity
kd = 0.0 watt / kelvin / meter
# Surface heat flow
Qs = 0.0 watt / meter ** 2
# Min oceanic depth
ocemin = 0.0 meter
[[geotherm3]]
model = continent
polygon = '''
0. 42.
1450. 42.
1450. 137
0. 137.
'''
T0 = 500.0 degC
yL = 137.0 kilometer
# Layer age
age = 0.0 year
# Thermal gradient
dtz = 0.0 kelvin / kilometer
# Thermal disturbance
dt = 0.0 kelvin
# Thermal diffusivity
kappa = 0.0 meter ** 2 / second
# Radius of anomaly
radius = 0.0 kilometer
# Direction of anomaly is right to left
rtol = False
# Bottom temperature of layer
T1 = 0.0 degC
# List / tuple of temperature values with units enclosed in quotes
Ti = 0.0 degC, 100.0 degC
# Constant layer temperature
value = 0.0 degC
# Center position of thermal anomaly: x
x0 = 0.0 kilometer
# Center position of thermal anomaly: y
y0 = 0.0 kilometer
# Radiogenic production
H = 0.0 meter ** 2 * watt
# Length scale of radiogenic production
hr = 0.0 kilometer
# Thermal conductivity
kd = 0.0 watt / kelvin / meter
# Surface heat flow
Qs = 0.0 watt / meter ** 2
# Min oceanic depth
ocemin = 0.0 meter
[[geotherm4]]
model = interpolated
Ti = 0.0 degC, 500.0 degC, 562.18 degC, 0.0 degC
polygon = '''
1450. 7.
1450. 42.
1520. 42.
1520. 11.
'''
[[geotherm5]]
model = interpolated
Ti = 500.0 degC, 1330.0 degC, 1330.0 degC, 1079.11 degC, 562.18 degC
polygon = '''
1450. 42.
1450. 137.
1520. 111.
1520. 80.0
1520. 42.0
'''
[Boundary Conditions]
[[Left]]
vx = 0.0 centimeter / year
vy = 0.0 centimeter / year
# All velocity components 0 on the boundary
noslip = False
# Constant Temperature
T = 0.0 degC
# Simmetric Temperature (no heat flow)
simmetricT = True
# Constant Pressure
P = 0.0 pascal
# Starting evolution step for this set of parameters
start = 0
[[Right]]
vx = -2.0 centimeter / year
vy = 0.0 centimeter / year
# All velocity components 0 on the boundary
noslip = False
# Constant Temperature
T = 0.0 degC
# Simmetric Temperature (no heat flow)
simmetricT = True
# Constant Pressure
P = 0.0 pascal
# Starting evolution step for this set of parameters
start = 0
[[Top]]
vx = 0.0 centimeter / year
vy = 0.0 centimeter / year
# All velocity components 0 on the boundary
noslip = False
# Constant Temperature
T = 0.0 degC
# Simmetric Temperature (no heat flow)
simmetricT = False
# Constant Pressure
P = 0.0 pascal
# Starting evolution step for this set of parameters
start = 0
[[Bottom]]
vx = 0.0 centimeter / year
vy = 8.45016e-11 meter / second
T = 1465.0 degC
# All velocity components 0 on the boundary
noslip = False
# Simmetric Temperature (no heat flow)
simmetricT = False
# Constant Pressure
P = 0.0 pascal
# Starting evolution step for this set of parameters
start = 0
[Topography]
# Enable topography
enabled = False
# Water level ???
waterlevel = 0.0 meter
# Topography resolution (number of nodes)
nx = 0
# Length of topography
size = 0.0 meter
[Peierls]
# Sigma ???
sigma = 9.1e9 pascal
# A ???
A = 63.0e-6 1 / pascal ** 2 / second
# m ???
m = 1.0 dimensionless
# n ???
n = 2.0 dimensionless
# Critical Pressure
Pr_crit = 200.0e6 pascal
# Critical Temperature
T_crit = 1100.0 degC
[Fluid]
enabled = True
threshold = 0.1
velocity_model = Darcy
rheology_model = DryWet
# minimum depth of fluid tracers
minDepth = 7.0e3 meter
[[Darcy]]
# Porosity
porosity = 0.01 dimensionless
# Permeability
permeability = 1e-18 meter ** 2
# Fluid phase density
fluid_density = 1000.0 kilogram / meter ** 3
# Fluid phase viscosity
fluid_viscosity = 1.0e-4 pascal * second