input.txt 7.54 KB
Newer Older
Lena Noack's avatar
Lena Noack committed
1
2
3
! Geometry	0
geometry	=	0 		! 0-box, 1-sph. annulus, 2-cylinder
dim		=	2 		! 2-2D, 3-3D
Oliver Henke-Seemann's avatar
Oliver Henke-Seemann committed
4
nl		=	40 	  	! number of grid points in lateral direction (2 ghost points are added when reflective boundary is used)
Lena Noack's avatar
Lena Noack committed
5
ny		=	0
Oliver Henke-Seemann's avatar
Oliver Henke-Seemann committed
6
nr		=	40 		! number of grid points in radial direction (+2 ghost points)
Lena Noack's avatar
Lena Noack committed
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
lmin		=	0.0		! only used when geometry=0
lmax		=	1.0		! only used when geometry=0
ymin		=	0.0
ymax		=	0.0
rmin		=	1.0 
rmax		=	2.0		! rmax=rmin+1.0
phi_factor	=	1		! only used when geometry>0, 1.0: full annulus, 0.5: halph annulus, ... 
CorrRot		=	0    


! Adaptivity    0
ad_factor       =       0		! 0-uniform box grid, 1-refine area (ad_min/ad_max values needed), 2-smallest/largest resolution (res_min/res_max) is set at (ad_l=0.5*lmax, ad_r=rmax)
res_ratio_l     =       0.1 
res_ratio_y	=	0.1
res_ratio_r     =       0.1 
Boissinot Pierre's avatar
Boissinot Pierre committed
22
23


Lena Noack's avatar
Lena Noack committed
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155

! Bnd_cond	0
Ttop		=	0.0		! temperature at surface, until now: has to be 0 or dimensional value (e.g. 288 for Earth)
Tbottom		=	1.0		! temperature at CMB, until now: has to be 1 or dimensional value (e.g. 3900 for Earth)
noslip_t	=	0		! 0: free-slip, 1: no-slip at surface boundary
noslip_b	=	0		! 0: free-slip, 1: no-slip at CMB boundary
inner_bound	=	0		! 1: isolated
outer_bound	=	0		! 1: isolated
open_b_bound	=	0		! 1: open boundary at bottom
periodic	=	0  		! 0: reflective boundary, 1: periodic boundary, is currently tested, normally: 0 for box and 1 for sphere


! Restart_TS 	0
use_snap	=	0		! insert number of TS from the snapshot, from which the simulation shall re-start
last_snap	=	0		! set to 1, then the last snapshot is automatically used (use_snap is ignored)
Cdt		=	10 		! Courant factor when positive value (e.g. 0.1), Delta criterium factor when negative (e.g. -10.0), dt_ini if 0
dt_ini		=	1.0e-6		! Initial time step or fixed time step if CourantCoeff=0
dt_max		=	1.0e-2		! Maximal time step
dt_min		=	1.0e-10		! Minimal time step
nbiter		=	0		! max number of time steps, not used if tmax>0
iter_out	=	0 		! output every ... time steps
tmax		=	3 		! maximal time, stop the simulation when reached; if tmax=0 then either nbiter or conv is used to stop the simulation
t_output	=	0.1		! output time
print_time	=	0


! Initial 	0
ampl		=	0.05 		! amplitude of initially applied sphercial harmonics
sph_l		=	-1 		! mode of spherical harmonics, chose +/-1.0 times the aspect ratio for box geometry and number of preferred plumes for polar geometry
sph_y		=	0
linear		=	3 		! initial temperature profile: 0: T=Tini everywhere; 1: linear, 2: conductive profile, 3: convective profile with TBLs
Tini		=	0.5		! initial mantle temperature for profile 3
s_bot		=	0.2		! lower thermal boundary layer (TBL) for profile 3
s_top		=	0.2		! upper thermal boundary layer (TBL) for profile 3
isotherm        =       0.7


! Viscosity     0
n_fac           =       0.0		! Arrhenius AND FK: non-Newtonian factor (1.0 - Newtonian, >1.0 non-Newtonian)
iniStrRate      =       0.0		! 10^-15 1/s ! Initial Strain rate
StartNN         =       0               ! Number of timesteps after which non-Newtonian law is applied
Tref            =       1.0		! reference temperature for Rayleigh number
zref            =       0.0		! reference depth for Rayleigh number
gam             =       0.0             ! FK Viscosity: gamma_T
gamP            =       0.0             ! FK Viscosity: gamma_p
e_gamT		=	1.0
e_gamP		=	1.0
E               =       0.0		! Arrhenius: actvation energy
V               =       0.0		! Arrhenius: activation volume
E_dis           =       0.0		! Arrhenius, non-Newtonian: activation energy
V_dis           =       0.0             ! Arrhenius, non-Newtonian: activation volume
A_dif           =       0.0             ! Arrhenius, mixed Newtonian: prefactor Newtonian
A_dis           =       0.0		! Arrhenius, mixed Newtonian: prefactor non-Newtonian: A' = A * eta_ref^n * (kappa/D^2)^(n-1)
T0              =       0.2		! Arrhenius: surface temperature (Tsurf/DeltaT)
AdVisc          =       0               ! adapt viscosity at boundaries, such that viscosity at ri/ro is correct
ViscMin		=	1.0e-10
ViscMax		=	1.0e+7


! Plasticity    0
plasticity      =       0               ! 0-Bingham, 1-nonsmooth, 2-angle-dependent
YieldStr0       =       0.0             ! plasticity = 0/1: surface yield stress (YS_0); plasticity = 2: surface cohesion
YieldStrZ       =       0               ! plasticity = 0/1: increase of yield stress with depth (YS_z*z); plasticity = 2: non-dim. factor for pressure-dependent part
friction        =       0.0		! non-dim factor: rho*g*D^3 / kappa / eta
mantle_fr_angl	=	0		! plasticity = 2: friction angle in degrees
mantle_fr_coh	=	0.0		! plasticity = 2: surface cohesion, 350 MPa


! Continents	0
cont_nr		=	0


! MantleRefVal	0
rho             =       1
rho_c		=	1
Tm_ref          =       1
k               =       1
Cp		=	1
Cp_c		=	1
g               =       1   	
alpha           =       1		
kappa           =       1		! actually not needed since rho, Cp and k are set and kappa=k/rho/Cp, but can be used to check if the other values are realistic or not...


! Ra_H0		0
Ra		=	1.0e+4		! Rayleigh number at Tref/zref
H0		=	0.0		! in W/m^3; initial heat source density, set to H0>0 for internal heating!  
lambda		=	0.0		! in 1/yr; decay of internal heat sources with time: H=H0*exp(-t*lambda)
CoreCooling	=	0


! Chemical	0
B		=	0.0		! Buoyancy ratio for chemical diffusion, needed in momentum equation
Le		=	0		! Lewis number for chemical diffusion, if 0, then no double-diffusive solver is used
LeV		=	0
comp_init	=	0 		! +-1 = linear, +-2 = rho value; negative: add cylinder, positive: add random perturbation
comp_ampl	=	0		! Amplitude of composition perturbation
C_ref		=	0 		! reference composition (should correspond to reference density)


! Particles 	0
use_part	=	0		! use particles instead of compositional field approach
part_n		=	0 		! total number of particles, ignored if part_n_cell>0
part_n_cell	=	10		! number of particles per cell
part_depth	=	1
restart_p	=	1


! Numerics 	0
debug		=	0		! set between 1 and 5 to obtain detailed information
max_outer       =       20		! number of max outer iterations (coupled energy-momentum solver), used if e_solver>0
max_inner       =       3               ! number of max inner iterations (coupled pressure-velocity solver), used if m_solver>0
noBouss		=	0		! if set to one then no Boussinesq is used in momentum solver
e_solver	=	3		! -1: no e-solver, 0: direct first-order explicit, 1: second-order iterative explicit, 2: first-oder iterative implicit, 3: second-order iterative implicit
c_solver	=	0
t_method	=	2		! 0=FDM, 1=FVM1 (+FDM for 2nd deriv.), 2=FVM2
s_method	=	1 		! 0=FDM, 1=FVM1 (+FDM for 2nd deriv.), 2=FVM2
no_m_solver	=	0
convT		=	1e-10		! criterion for convergence of outer iteration, used if e_solver>0 and number of outer iterations < max_outer
convV		=	1e-4 		! criterion for convergence of inner iteration, used if m_solver>0 and number of inner iterations < max_inner
conv		=	0		! stops the simulation if convergence of time steps is reached, not used if tmax>0 or nbiter is reached
test		=	0		! output of pressure, when used IDL keyword /test has to be used, as well
beta		=	1.0		! for flux limiter; 1: minmod, 2: superbee
average		=	1		! averaging scheme to be used for viscosity: 0=old, 1=harmonic, 2=arithmetic, 3=geometric
p_average	=	1		! averaging scheme to be used for particles: 1=harmonic, 2=arithmetic, 3=geometric
p_average2	=	2		! for density averaging; if 0 then set to p_average
p_penalty	=	1.0e-7		! this value shoud be 1.0e-7 and is divided by the maximal viscosity in the stokes equation
nmesures	=	0		! currently not used... should be included again for the cases when quasi-steady-state is reached to gain average values for the last X timesteps


EOF		=	1