1- # settings file for BOUT++
2- #
3- # Blob simulation in a 2D slab
4- #
5- # This case has blob size
6- #
7- # delta = 0.3*256 ~ 10 * delta_*
8-
9-
10- # settings used by the core code
11-
12- nout = 3 # number of time-steps
13- timestep = 50 # time between outputs [1/wci]
141
152
163MXG = 2 # Number of X guard cells
@@ -20,7 +7,7 @@ MYG = 0 # No y derivatives, so no guard cells needed in y
207
218nx = 260 # Note: 4 guard cells
229ny = 1
23- nz = 256
10+ nz = 256
2411
2512dx = 0.3 # Grid spacing [rho_s]
2613dz = 0.3
@@ -30,7 +17,7 @@ dz = 0.3
3017
3118[mesh:ddx]
3219
33- first = C2
20+ first = C2
3421second = C2
3522upwind = W3
3623
@@ -51,9 +38,22 @@ upwind = W3
5138
5239[solver]
5340
54- atol = 1.0e -10 # absolute tolerance
55- rtol = 1.0e-5 # relative tolerance
41+ atol = 1e -10 # absolute tolerance
42+ rtol = 1e-05 # relative tolerance
5643mxstep = 10000 # Maximum internal steps per output
44+ output_step = 50 # time between outputs [1/wci]
45+ # settings file for BOUT++
46+ #
47+ # Blob simulation in a 2D slab
48+ #
49+ # This case has blob size
50+ #
51+ # delta = 0.3*256 ~ 10 * delta_*
52+
53+
54+ # settings used by the core code
55+
56+ nout = 3 # number of time-steps
5757
5858###################################################
5959# Electrostatic potential solver
@@ -66,19 +66,17 @@ pctype = user # Preconditioning type
6666fourth_order = true # 4th order or 2nd order
6767
6868flags = 0 # inversion flags for phi
69- # 0 = Zero value
70- # 10 = Zero gradient AC inner & outer
71- # 15 = Zero gradient AC and DC
72- # 768 = Zero laplace inner & outer
69+ # 0 = Zero value
70+ # 10 = Zero gradient AC inner & outer
71+ # 15 = Zero gradient AC and DC
72+ # 768 = Zero laplace inner & outer
7373
7474[phiSolver:precon] # Preconditioner (if pctype=user)
75- filter = 0. # Must not filter solution
76- flags = 49152 # set_rhs i.e. identity matrix in boundaries
75+ filter = 0.0 # Must not filter solution
76+ flags = 49152 # set_rhs i.e. identity matrix in boundaries
7777
7878###################################################
7979# Electrostatic potential solver (Boussinesq)
80-
81- [phiBoussinesq]
8280# By default type is tri (serial) or spt (parallel)
8381#flags = 0
8482
@@ -89,16 +87,16 @@ flags = 49152 # set_rhs i.e. identity matrix in boundaries
8987
9088Te0 = 5 # Electron Temperature (eV)
9189
92- n0 = 2e18 # Background plasma density (m^-3)
90+ n0 = 2e+18 # Background plasma density (m^-3)
9391
94- compressible = false # Compressibility?
92+ compressible = false # Compressibility?
9593
9694boussinesq = true # Boussinesq approximation (no perturbed n in vorticity)
9795
98- D_vort = 1e-6 # Viscosity
99- D_n = 1e-6 # Diffusion
100-
101- R_c = 1.5 # Radius of curvature (m)
96+ D_vort = 1e-06 # Viscosity
97+ D_n = 1e-06 # Diffusion
98+
99+ R_c = 1.5 # Radius of curvature (m)
102100
103101# settings for individual variables
104102# The section "All" defines default settings for all variables
@@ -117,5 +115,3 @@ height = 0.5
117115width = 0.05
118116
119117function = 1 + height * exp(-((x-0.25)/width)^2 - ((z/(2*pi) - 0.5)/width)^2)
120-
121-
0 commit comments