parameter_example.yml 6.42 KB
Newer Older
1
# Define the system of units to use internally. 
2
InternalUnitSystem:
3
4
5
6
7
  UnitMass_in_cgs:     1   # Grams
  UnitLength_in_cgs:   1   # Centimeters
  UnitVelocity_in_cgs: 1   # Centimeters per second
  UnitCurrent_in_cgs:  1   # Amperes
  UnitTemp_in_cgs:     1   # Kelvin
Matthieu Schaller's avatar
Matthieu Schaller committed
8

9
# Parameters for the task scheduling
Matthieu Schaller's avatar
Matthieu Schaller committed
10
Scheduler:
11
12
13
14
15
16
  nr_queues:             0        # (Optional) The number of task queues to use. Use 0  to let the system decide.
  cell_max_size:         8000000  # (Optional) Maximal number of interactions per task if we force the split (this is the default value).
  cell_sub_size:         64000000 # (Optional) Maximal number of interactions per sub-task  (this is the default value).
  cell_split_size:       400      # (Optional) Maximal number of particles per cell (this is the default value).
  cell_max_count:        10000    # (Optional) Maximal number of particles per cell allowed before triggering a sanitizing (this is the default value).
  max_top_level_cells:   12       # (Optional) Maximal number of top-level cells per dimension. The number of top-level cells will be the cube of this.
Matthieu Schaller's avatar
Matthieu Schaller committed
17

18
# Parameters governing the time integration
Matthieu Schaller's avatar
Matthieu Schaller committed
19
TimeIntegration:
20
21
22
23
24
  time_begin: 0.    # The starting time of the simulation (in internal units).
  time_end:   1.    # The end time of the simulation (in internal units).
  dt_min:     1e-6  # The minimal time-step size of the simulation (in internal units).
  dt_max:     1e-2  # The maximal time-step size of the simulation (in internal units).

25
26
# Parameters governing the snapshots
Snapshots:
27
28
29
  basename:   output      # Common part of the name of output files
  time_first: 0.          # Time of the first output (in internal units)
  delta_time: 0.01        # Time difference between consecutive outputs (in internal units)
30
  compression: 0          # (Optional) Set the level of compression of the HDF5 datasets [0-9]. 0 does no compression.
31
32
33
34
35
  UnitMass_in_cgs:     1  # (Optional) Unit system for the outputs (Grams)
  UnitLength_in_cgs:   1  # (Optional) Unit system for the outputs (Centimeters)
  UnitVelocity_in_cgs: 1  # (Optional) Unit system for the outputs (Centimeters per second)
  UnitCurrent_in_cgs:  1  # (Optional) Unit system for the outputs (Amperes)
  UnitTemp_in_cgs:     1  # (Optional) Unit system for the outputs (Kelvin)
36

37
38
# Parameters governing the conserved quantities statistics
Statistics:
Matthieu Schaller's avatar
Matthieu Schaller committed
39
40
  delta_time:          1e-2      # Time between statistics output
  energy_file_name:    energy    # (Optional) File name for energy output
41
  timestep_file_name:  timesteps # (Optional) File name for timing information output. Note: No underscores "_" allowed in file name 
42

43
# Parameters for the hydrodynamics scheme
Matthieu Schaller's avatar
Matthieu Schaller committed
44
SPH:
45
46
47
  resolution_eta:        1.2348   # Target smoothing length in units of the mean inter-particle separation (1.2348 == 48Ngbs with the cubic spline kernel).
  delta_neighbours:      0.1      # The tolerance for the targetted number of neighbours.
  CFL_condition:         0.1      # Courant-Friedrich-Levy condition for time integration.
48
  max_ghost_iterations:  30       # (Optional) Maximal number of iterations allowed to converge towards the smoothing length.
49
  max_volume_change:     2.       # (Optional) Maximal allowed change of kernel volume over one time-step
Matthieu Schaller's avatar
Matthieu Schaller committed
50

51
# Parameters related to the initial conditions
Matthieu Schaller's avatar
Matthieu Schaller committed
52
InitialConditions:
53
  file_name:  SedovBlast/sedov.hdf5 # The file to read
54
55
  h_scaling:  1.                    # (Optional) A scaling factor to apply to all smoothing lengths in the ICs.
  shift_x:    0.                    # (Optional) A shift to apply to all particles read from the ICs (in internal units).
Matthieu Schaller's avatar
Matthieu Schaller committed
56
57
58
  shift_y:    0.
  shift_z:    0.

59
# Parameters governing domain decomposition
Matthieu Schaller's avatar
Matthieu Schaller committed
60
DomainDecomposition:
61
62
  initial_type:       m     # (Optional) The initial strategy ("g", "m", "w", or "v").
  initial_grid_x:    10     # (Optional) Grid size if the "g" strategy is chosen.
63
64
  initial_grid_y:    10
  initial_grid_z:    10
65
  repartition_type:   b     # (Optional) The re-decomposition strategy ("n", "b", "v", "e" or "x").
66
 
67
# Parameters related to external potentials --------------------------------------------
68
  
69
# Point mass external potentials
70
PointMassPotential:
71
  position_x:      50.      # location of external point mass (internal units)
72
  position_y:      50.
73
  position_z:      50.
74
75
  mass:            1e10     # mass of external point mass (internal units)
  timestep_mult:   0.03     # Dimensionless pre-factor for the time-step condition
76

77
# Isothermal potential parameters
78
IsothermalPotential:
79
  position_x:      100.     # Location of centre of isothermal potential with respect to centre of the box (internal units)
80
  position_y:      100.
81
  position_z:      100.
82
83
84
  vrot:            200.     # Rotation speed of isothermal potential (internal units)
  timestep_mult:   0.03     # Dimensionless pre-factor for the time-step condition

85
86
87
88
89
90
91
92
93
# External potential parameters
SoftenedIsothermalPotential:
  position_x:      0.       # Location of centre of isothermal potential with respect to centre of the box (internal units)
  position_y:      0.
  position_z:      0.	
  vrot:            200.     # rotation speed of isothermal potential (internal units)
  epsilon:         0.1      # Softening size (internal units)
  timestep_mult:   0.03     # controls time step
  
94
# Disk-patch potential parameters
95
96
97
DiscPatchPotential:
  surface_density: 10.      # Surface density of the disc (internal units)
  scale_height:    100.     # Scale height of the disc (internal units)
98
  z_disc:          200.     # Position of the disc along the z-axis (internal units)
99
  timestep_mult:   0.03     # Dimensionless pre-factor for the time-step condition
100
  growth_time:     5.       # (Optional) Time for the disc to grow to its final size (multiple of the dynamical time)
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117

# Parameters related to cooling function  ----------------------------------------------

# Constant du/dt cooling function
ConstCooling:
  cooling_rate: 1.          # Cooling rate (du/dt) (internal units)
  min_energy:   1.          # Minimal internal energy per unit mass (internal units)
  cooling_tstep_mult: 1.    # Dimensionless pre-factor for the time-step condition

# Constant lambda cooling function
LambdaCooling:
  lambda:                      2.0   # Cooling rate (in cgs units)
  minimum_temperature:         1.0e4 # Minimal temperature (Kelvin)
  mean_molecular_weight:       0.59  # Mean molecular weight
  hydrogen_mass_abundance:     0.75  # Hydrogen mass abundance (dimensionless)
  cooling_tstep_mult:          1.0   # Dimensionless pre-factor for the time-step condition