WikiStart: scaling5h_intrepid.par

File scaling5h_intrepid.par, 7.0 KB (added by daley, 10 years ago)

flash.par used on BG/P for 6144 processor (VN mode) run in weak scaling curve

Line 
1# Number of processors: 1536*4
2#####################################################################
3#
4# FLASH file for a Type Ia supernova in 3-d cartesian geometry
5#
6# off-set ignition whole-star model without tracer particles
7#
8# This file is for the weak convergence tests of July 2008
9#
10# This file is for a 0.2 km resolution run on Seaborg, with a centered
11# bubble (with varying radius)
12#
13#   study evaluates convergence with the sharpened KPP ADR flame model,
14#   also first large set of runs with quiet, realistic energy release
15#   and neutronization
16#
17# configuration:
18#       block of 16
19#       "size 16" model
20#       domain [-2457.6,2457.6]^3 km
21#       refine_min=1
22#       at min refinement, base of 3 blocks in each dimension (27 total)
23#
24#
25#  setup with:
26#  ./setup WD_def -3d +cube16 -maxblocks=150 +noio -auto
27#
28#-----------------------------
29# refine_max resolution  km
30#-----------------------------
31#       1        48   102.4
32#       2        96    51.2
33#       3       192    25.6
34#       4       384    12.8
35#       5       768     6.4
36#       6      1536     3.2
37#       7      3072     1.6
38#       8      6144     0.8
39#       9     12288     0.4
40#       10    24576     0.2
41#       11    49152     0.1
42#       12    98304     0.05
43#-----------------------------
44#
45#####################################################################
46
47##################
48# model parameters
49##################
50
51# amplitude of random velocity perturbations
52# not used
53v_pert                  = 0.e5
54n_pert                  = 0
55
56# central velocity dipole size and speed
57# not used
58r_conv                  = 0.e5
59v_conv                  = 0.e5
60
61# Euler angles of the dipole in degrees
62# not used
63phi_conv                = 0.e0
64theta_conv              = 0.e0
65
66# flame capturing and nuclear network
67useFlame                = .true.
68useBurn                 = .true.
69
70# ignition specific
71# position and size of spherical region set as burned material
72# at time zero
73ignite                  = .true.
74x_match                 =  0.e5
75y_match                 =  0.e5
76z_match                 =  0.e5
77r_match                 = 90.5000e5
78
79# damping
80# not used
81damp_method             = 4
82damp_const1             = 1.e-3
83damp_const2             = 1.e-3
84damp_const3             = 1.8e8
85damp_time               = 0.5e0
86
87#################
88# flame capturing
89#################
90
91#eos_coulombMult = 0.0
92
93
94# front tracking stuff
95# saves flame speed
96fl_debug                = .true.
97fl_epsilon              = 1.e-12
98fl_epsilon_0            = 1.e-3
99fl_epsilon_1            = 1.e-3
100fl_kpp_fact             = 1.309e0
101fl_b                    = 3.2e0
102
103
104# use an approximate atwood number for calculating
105# turbulent flame speed
106approx_atwood           = .true.
107
108# turbulent flame speed parameters
109# turbulent flame speed is
110# 0.5*sqrt(A g dx* fl_subgrid_s1)
111fl_subgrid_s1           = 0.06e0
112fl_subgrid_s2           = 0.e0
113flame_speed_mult        = 10.e0
114pbFqdotMult     = 1.e0
115
116# flame quenching at low densities
117# TODO do we want to change this to 1e6-1e5?
118use_quenching           = .true.
119quenching_dens0         = 5.e6
120quenching_dens1         = 1.e7
121
122#################
123# nuclear network
124#################
125
126# not used (iburn=0)
127# except xrenorm_burning is checked
128
129dnucmin                 = 1.e6
130dnucmax                 = 1.e14
131tnucmin                 = 1.e8
132tnucmax                 = 1.e12
133#shock_burning          = .false.
134#xrenorm_burning                = .true.
135#ode_steper             = 1
136
137#################
138# mesh definition
139#################
140
141# max and min levels of refinement
142lrefine_min             = 1
143lrefine_max             = 10
144lrefine_del             = 0
145
146# use custom refinement criteria
147refine_var_1            = "none"
148refine_var_2            = "none"
149refine_var_3            = "none"
150refine_var_4            = "none"
151
152refine_ddens            = 0.5e0
153refine_dtvel            = 0.2e0
154refine_xtvel            = 100.e5
155refine_dphi1            = 0.1e0
156refine_xphi1            = 1.e-15
157refine_sphi1            = 1.e-20
158refine_xenuc            = 1.e18
159
160# minimum density at which refinement criteria other than
161# flame front are checked
162refine_dens_min         = 5.e5
163# max density for weighting refinement
164# weighted by log(dens/mindens)/log(maxdens/mindens)
165refine_dens_max         = 1.e10
166# above uni dens, don't let resolution go coarser than uni_dx
167refine_uni_dens         = 3.e16
168refine_uni_dx           = 0.1e5
169# within uni_radius, don't let resolution go coarser than uni_dx
170refine_uni_radius       = 84.1000e5
171# outside max_radius force resolution to be coarser than 2*uni_dx
172refine_max_radius       = 2500.e5
173
174# enhancement of refinement in inner region
175# ( help decrease number of particles per processor)
176refine_inner_dens_min   = 2.5e10
177refine_inner_dens_dx    = 10.e5
178
179refine_ignition_time    = 10.e0
180#The following is unused!? - KW
181refine_ignition_radius  = 30.e5
182lrefine_ignition_del    = 0
183
184#gcd refinement values
185 
186gcd_focus_time          = 1.5
187gcd_focus_angle         = 30.0
188gcd_focus_min_radius    = 1500.e5
189gcd_focus_max_radius    = 3000.e5
190
191
192# physical domain
193geometry                = "cartesian"
194
195xmin                    =  -2457.6e5
196xmax                    =  2457.6e5
197ymin                    =  -2457.6e5
198ymax                    =  2457.6e5
199zmin                    =  -2457.6e5
200zmax                    =  2457.6e5
201
202quadrant                = .false.
203octant                  = .false.
204
205nblockx                 = 3
206nblocky                 = 3
207nblockz                 = 3
208
209# boundary conditions for the whole star
210xl_boundary_type        = "diode"
211xr_boundary_type        = "diode"
212
213yl_boundary_type        = "diode"
214yr_boundary_type        = "diode"
215
216zl_boundary_type        = "diode"
217zr_boundary_type        = "diode"
218
219#bc_velocity_type       = "diode"
220
221grav_boundary_type      = "isolated"
222
223#########
224# physics
225#########
226
227# self-gravity
228##igrav                 = 1
229mpole_lmax              = 4
230grav_3daxisymmetric     = .true.
231grav_dumpMoments        = .false.
232grav_subsamplepotential = 1
233
234# eos
235eos_coulombMult         = 1.e0
236
237#####
238# I/O
239#####
240
241# names of files
242basenm          = "conv_bgp_S5h_"
243
244# for starting a new run
245# remember to set enuc_tint below!!!
246restart         = .false.
247
248# checkpointFileNumber    = 1   
249# plotFileNumber                = 1
250
251# dump checkpoint files every trstrt seconds
252checkpointFileIntervalTime                      = 1.0E-01
253checkpointFileIntervalStep                      = 0
254
255# dump plot files every tplot seconds
256plotFileIntervalTime                    = 5.0e-3
257
258# go for nend steps or tmax seconds, whichever comes first
259nend                    = 10
260tmax                    = 0.1e0
261wall_clock_checkpoint   = 21600.e0 # 6 hrs
262#wall_clock_checkpoint   = 85700.e0 # 23.81 hrs
263
264# variables for plotting
265# TODO do these work without the spaces in them?
266plot_var_1              = "dens"
267plot_var_2              = "temp"
268plot_var_10             = "rpv1"
269
270corners                 = .false.
271
272##########
273# timestep
274##########
275
276# CFL limit
277cfl                     = 0.80e0
278cfl_ini                 = 0.1e0
279
280# initial and minimum
281dtinit                  = 1.59e-5
282dtmin                   = 1.e-20
283dtmax                   = 8.e-4
284tstep_change_factor     = 1.2e0
285#enuc_factor            = 1.e0
286
287# density, temperature and composition of the fluff
288dens_fluff              = 1.e-3
289temp_fluff              = 3.e7
290xc12_fluff              = 0.5e0
291xo16_fluff              = 0.5e0
292
293##################
294# tracer particles
295##################
296
297# ParticleFileNumber       = 1
298
299##ParticleFileIntervalTime = 5.0e-3
300##ParticleFileIntervalStep = 0
301
302useParticles = .false.
303pt_numParticlesWanted = 0
304pt_maxPerProc =         0
305
306#############
307# miscellaneous
308#############
309
310# floor values
311
312small                   = 1.e-100
313smalle                  = 1.e-100
314smallt                  = 1.e7
315smallu                  = 1.e-10
316# set to 1/5 of dens_fluff
317smlrho                  = 2.e-4
318smallp                  = 1.e15
319smallx                  = 1.e-100
320
321# misc (many are same as defaults)
322
323eintSwitch              = 0.e0
324cvisc                   = 0.1e0
325iplm                    = 0
326use_steepening          = .false.
327ppm_modifystates        = .false.
328leveque                 = .false.
329hybrid_riemann          = .true.
330convertToConsvdForMeshCalls     = .false.
331converttoconsvdinmeshinterp     = .true.
332use_cma_advection       = .false.
333dp_sh_md                = 1.e0
334dp_sh                   = 1.e0
335
336enuc_tint               = 0.e48  # this is a integral quantities thing.
337
338#checkpoint_timers       = .false.
339
340memory_stat_freq         = 1
341run_number = 2