vel_relaxations=10 |
maximum number of velocity loops (default=10) |
piterations=100 |
maximum Uzawa iteration loops (default=100) |
viterations=20 |
number of velocity iterations before checking convergence (default=20) |
Solver=multigrid |
multigrid or (later, conjugate gradient) |
mg_cycle=1 |
style of multigrid cycle |
|
1 = V cycle, 2 = W cycle, ... (default=1) |
accuracy= |
desired accuracy of Uzawa algorithm (default=1.e-4) |
delta_accuracy_factor=1.0 |
change in accuracy level->level (> 1 for nonN, < 1 for Newt) |
|
0.001 < factor < 10.0 (default=0.2) |
gs_under_relax=1.0 |
Gauss-Seidel under-relaxation (default=1.0) |
|
# Grid Mesh and Levels |
mgunitx=5 |
(Integer) Coarsest (base) multigrid cell dimensions (elements) along the X axis |
mgunitz=3 |
(Integer) Multigrid cells along the Z axis. The default/minimum is 2. |
levels=4 |
(Integer) Number of multigrid levels. The default is 1. |
# Tracer setup |
Tracers=on |
initialise tracers (default=on) |
Tracer_appetite=0.5 |
( size(tracer1)+size(tracer2) ) x Tracer_appetite (default=0.5) |
Tracer_voids=off |
allow tracers to disappear (default=off) |
Tracer_rect=1 |
number of rectangular regions of different tracer densities (default=0) |
Tracer_rect_density=4 |
tracer density (N x N per finest element) (<=12) |
Tracer_rect_x1=0 |
the coordinate extent of region should be same size as the model box |
Tracer_rect_x2=dimenx |
|
Tracer_rect_z1=0 |
Tracer_rect_z2=dimenz |
|
|
Material_1_density= |
density (default=1.0) |
Material_1_porosity= |
initial porosity (default=0.0) |
|
NB: initial porosity = 0 ensures that Bulk_visc is constant |
Material_1_Bulk_visc= |
bulk visc ratio at initial porosity ( >1.0(2D), >2/3(3D) ) |
|
bulk visc = ratio*visc (default=-1.0=infinite) |
|
div(v) + p/(bulk_visc-2/3*visc) = 0 |
Material_1_Bulk_modulus= |
B in dp = B*div(v)*dt (slightly compressible formulation) (default=0.0) where dp is on tracers (isotropic stress) |
Material_1_reproduction=on |
allow tracer reproduction (default=on) |
Material_1_phases=1 |
number of unique phases (first phase is phase 0) (default=1) then visc = [ sum(1/visc_n) ] (-1) |
Material_1_T_block= |
(blocking) T above which phase change can occur (default=-1.e32) |
Material_1_rheol_cpts= |
number of rheological components (at least one per phase) (default=1) then visc = [ sum(1/visc_n) ] (-1) |
Material_1_Trange_min=-1.e8 |
temperature range to which rheology applies |
Material_1_Trange_max=1.e8 |
(default=-1.e8,1.e8) |
Material_1_rheol_phase=0 |
phase to which each rheology applies (start from 0) (default=0) |
|
# Rheological model |
Material_2_rheol_T_type=1 |
rheological temperature-dependence model (default=2) |
|
(1) visc=N0*exp(-T1*T) (Frank-Kamenetski) |
|
(2) visc=N0*exp (E+Z*z)/(T1*(T+T0)) (Arrhenius) where z=depth |
Material_2_viscN0= |
N0 in viscosity models (default=1.0) |
Material_2_viscT1= |
T1 in viscosity models (default=1.0) |
Material_2_viscT0= |
T0 in Arhennius viscosity model (default=0.0) |
Material_2_viscZ= |
Z in Arrhenius viscosity model (default=0.0) |
Material_2_viscE= |
E in Arrhenius viscosity model (default=0.0) |
Material_2_viscTmax= |
maximum and minimum T to use in calculating viscosity |
Material_2_viscTmin= |
(default=1.e32,0.0) |
Material_2_sdepv_expt= |
exponent "s" in stress dependance of viscosity |
|
|
# Stress-strain relationship |
Material_2_yield_stress_minimum= |
minimum yield stress for plastic deformation (default=1.e-32) |
Material_2_yield_stress_maximum= |
maximum yield stress for semi-brittle effect (default=1.e32) |
Material_2_yield_stress_B0= |
"cohesion" B0 in above eqn (default=1.e32) |
Material_2_yield_stress_Bp= |
"friction coefficient" Bp in above eqn (default=0.0) |
Material_2_yield_stress_Ea= |
ratio Ea = f(0,0)/f(E0,0) (default=1.0,range=[0,1])0 |
Material_2_yield_stress_E0= |
strain weakening E0 (default=1.e32) |
Material_2_yield_stress_En= |
exponent En in f(e), e<E0 (default=0.0) |
Material_2_yield_stress_Bz= |
"friction coefficient" Bz in above eqn (default=0.0) |
Material_2_yield_stress_Bc= |
tension cutoff Bc in above law (default=1.e32) |
Material_2_yield_stress_Edota= |
ratio Edota = f(0,0)/f(0,Edot0) (default=1.0,range=[0,1]) |
Material_2_yield_stress_Edot0= |
strain rate weakening Edot0 (default=0.0) |
Material_2_yield_stress_Edotn= |
exponent Edotn in f(e), edot<Edot0 (default=0.0) |
Material_2_yield_stress_ET= |
T above which strain weakening is reset (default=1.e32) |
Material_2_yield_stress_E0dt= |
time rate of strain reduction (healing) |
|
previous_temperature_file="" |
initial temperature configuration to use |
particle_input="" |
initial particle configuration (XDR binary) |
previous_particle_data="" |
initial particle properties to use |
|
possible keywords for binary_data and particle_data: |
|
Temp = temperature |
|
Pres, PorP = solid pressure, pore pressure |
|
Visc = viscosity |
|
Poro = porosity |
|
Edot = strain rate |
|
StrP, StrT = integrated plastic/total strain |
|
Grsz = grain size |
|
|
|
datatypes="Temp,Pres,Pstn" |
desired nodal output variables (ascii file) |
|
possible keywords for datatypes: |
|
Velx, Vely, Velz = x, y, z velocity |
|
PoVx, PoVz = x, z pore liquid velocity |
|
Pres, PorP = (nodal) solid pressure, pore pressure |
|
Temp = temperature, Comp = compression |
|
Strf = stream function (2D cartesian coords only) |
|
Pbdy = phase boundary |
|
|
|
observables="" |
surface observables (slice at z=0/zmax) (ascii) |
|
possible keywords for observables: |
|
Shfl, Bhfl = surface and basal heat fluxes |
|
Vxsf, Vysf, Vzsf = surface velocities |
|
Tpgx, Tpbx = surface and basal topography |
|
Tpgk, Tpbk = surface and basal topography wavenumbers |
|
Grvx, Grbx, Grtx = surface, basal, and topographical gravity signals |
|
Grvk, Grbk, Grtk = wavenumbers of above |
|
Geox, Gebx, Getx = surface, basal, and topographical geoid |
|
Geok, Gebk, Getk = wavenumbers of above |
|
PPM_files= |
number of PPM files at each output step (default=1) & first PPM file is *.ppm0, etc. |
PPM_height= |
vertical size of output PPM file (default=256) |
PPM_coloring= |
variable upon which to base colouring |
|
possible choices for PPM_coloring: (default=1) |
|
1=temperature, 2=viscosity, 3=stress=visc*edot, |
|
4=solid pressure, 5=grainsize, 6=compression, |
|
7=strain rate, 8=accumulated strain |
|
9=pore pressure, 10=permeability, 13=melt production |
|
|
PPM_coloring_autorange= |
automatically scale colour (default=1) |
PPM_coloring_min= |
min value for color scale (default=0.0) |
PPM_coloring_max= |
max value for color scale (default=1.0) |
PPM_show_strain= |
colour according to actual strength change (default=0.0) |
|
|
|
Sampling_tracers= |
number of sampling tracers (default=0) |
Sampling_lagrangian= |
(fixed) Eulerian=0 (Stay fixed with mesh) Lagrangian=1 (Move within mesh) (default=0) |
Sampling_plot_num= |
PPM file in which profile is stored (default=0) |
Sampling_x= |
initial x location (default=0.0) |
Sampling_z= |
initial z location (default=0.0) |
Sampling_y= |
initial z location (default=0.0) |
Sampling_field= |
field to sample (default=0) |
|
1=temperature, 2=x velocity, 3=z velocity |
|
4=nodal pressure, 5=strain rate, 6=stress=visc*edot 14=depl |
Sampling_dirn= |
profile direction (1=x, 2=z, 3=y) (default=0) |
Sampling_normalize= |
0=unnormalized, 1=normalized (default=0) |
Sampling_plot_min= |
If not autoranging, then need a max/min for the scale |
Sampling_plot_max= |
(default=1e4) |
Sampling_R= |
RGB colour mixes for profiles scaled between 0 and 1 (255 real RGB terms) (default=0.0) |
Sampling_G= |
Sampling_B= |