parameter_example.yml 5.74 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
  nr_queues:        0        # (Optional) The number of task queues to use. Use 0  to let the system decide.
12
13
  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).
14
  cell_split_size:  400      # (Optional) Maximal number of particles per cell (this is the default value).
Matthieu Schaller's avatar
Matthieu Schaller committed
15

16
# Parameters governing the time integration
Matthieu Schaller's avatar
Matthieu Schaller committed
17
TimeIntegration:
18
19
20
21
22
  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).

23
24
# Parameters governing the snapshots
Snapshots:
25
26
27
  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)
28
  compression: 0          # (Optional) Set the level of compression of the HDF5 datasets [0-9]. 0 does no compression.
29
30
31
32
33
  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)
34

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

41
# Parameters for the hydrodynamics scheme
Matthieu Schaller's avatar
Matthieu Schaller committed
42
SPH:
43
44
  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.
45
  max_ghost_iterations:  30       # (Optional) Maximal number of iterations allowed to converge towards the smoothing length.
46
47
  max_smoothing_length:  0.1      # Maximal smoothing length allowed (in internal units).
  CFL_condition:         0.1      # Courant-Friedrich-Levy condition for time integration.
48
  max_volume_change:     2.       # (Optional) Maximal allowed change of kernel volume over one time-step
Matthieu Schaller's avatar
Matthieu Schaller committed
49

50
# Parameters related to the initial conditions
Matthieu Schaller's avatar
Matthieu Schaller committed
51
InitialConditions:
52
  file_name:  SedovBlast/sedov.hdf5 # The file to read
53
54
  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
55
56
57
  shift_y:    0.
  shift_z:    0.

58
# Parameters governing domain decomposition
Matthieu Schaller's avatar
Matthieu Schaller committed
59
DomainDecomposition:
60
61
  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.
62
63
  initial_grid_y:    10
  initial_grid_z:    10
64
  repartition_type:   b     # (Optional) The re-decomposition strategy ("n", "b", "v", "e" or "x").
65
 
66
# Parameters related to external potentials --------------------------------------------
67
  
68
# Point mass external potentials
69
PointMassPotential:
70
  position_x:      50.      # location of external point mass (internal units)
71
  position_y:      50.
72
  position_z:      50.
73
74
  mass:            1e10     # mass of external point mass (internal units)
  timestep_mult:   0.03     # Dimensionless pre-factor for the time-step condition
75

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

84
# Disk-patch potential parameters
85
86
87
DiscPatchPotential:
  surface_density: 10.      # Surface density of the disc (internal units)
  scale_height:    100.     # Scale height of the disc (internal units)
88
  z_disc:          200.     # Position of the disc along the z-axis (internal units)
89
  timestep_mult:   0.03     # Dimensionless pre-factor for the time-step condition
90
  growth_time:     5.       # (Optional) Time for the disc to grow to its final size (multiple of the dynamical time)
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107

# 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