![]() |
Contains all runtime parameters as well as phys and math constants. More...
Data Types | |
type | unit_type |
constants and unit conversion factors More... | |
Functions/Subroutines | |
subroutine, public | unit_define () |
Declares values for the elements in unit_type. More... | |
subroutine | read_param (parfile) |
This function reads the parameter file. More... | |
subroutine | set_param (param_name, param_value) |
Sets a global parameter param_name to the value, given by its string representation param_value. More... | |
integer function | read_integer_param (input_string, param_name) |
Converts a string to an integer. More... | |
real(r_kind) function | read_float_param (input_string, param_name) |
Converts a string to an float. More... | |
subroutine | set_default_param |
Sets default parameters Parameters are sorted in alphabetical order. More... | |
subroutine | output_param |
Output parameters to a file. More... | |
subroutine | check_param |
Check for consistency of parameter. More... | |
subroutine | search_close_parameter (possible_pars, input_par, score, cl_par) |
Search for similar existing parameters. More... | |
real(r_kind) function | levenshteindistance (par1, par2) |
Calculates the Levenshtein distance between two strings. More... | |
subroutine | read_sanity_check_prepared_network |
Read important parameters from a binary file and check for consistency. More... | |
subroutine | output_param_prepared_network (path) |
Output relevant parameters to a file in the prepared network path. More... | |
subroutine | output_binary_parameter_data (path) |
Output the parameter data to a binary file. More... | |
Variables | |
type(unit_type), public | unit |
constants and unit conversion factors More... | |
integer, parameter, public | param_name_len = 200 |
– hardcoded parameters More... | |
integer, parameter, public | max_fname_len = 400 |
maximum length of filenames More... | |
integer | out_every |
– runtime parameters set when reading parameter file More... | |
integer | snapshot_every |
snapshot output frequency More... | |
integer | h_snapshot_every |
snapshot output frequency in hdf5 format More... | |
integer | flow_every |
flow output frequency More... | |
integer | h_flow_every |
flow output frequency in hdf5 format More... | |
integer | timescales_every |
timescales output frequency More... | |
integer | h_timescales_every |
timescales output frequency in hdf5 format More... | |
integer | engen_every |
Energy generation output frequency. More... | |
integer | h_engen_every |
Energy generation output frequency in hdf5 format. More... | |
integer | nrdiag_every |
frequency of NR loop diagnostic messages (wrt iteration counter) More... | |
integer | mainout_every |
frequency of mainout output More... | |
integer | h_mainout_every |
HDF5 output frequency of the mainout. More... | |
integer | track_nuclei_every |
frequency of track nuclei output More... | |
integer | h_track_nuclei_every |
frequency of track nuclei output in hdf5 format More... | |
integer | top_engen_every |
frequency of energy generators toplist More... | |
logical | use_alpha_decay_file |
Switch for using additional alpha decay rates. More... | |
integer | alpha_decay_zmin |
Minimum Z for additional alpha decay rates. More... | |
integer | alpha_decay_zmax |
Maximum Z for additional alpha decay rates. More... | |
character(max_fname_len) | alpha_decay_src_ignore |
Source flag(s) to ignore within the alpha decay rates. More... | |
logical | alpha_decay_ignore_all |
Flag whether rates should actually get replaced or only added. More... | |
integer | iwformat |
defines format of the weak rates (0 = tabulated, 1 = log<ft>) More... | |
integer | iwinterp |
defines the interpolation for the weak rates (0 = bilinear, 1 = bicubic) More... | |
real(r_kind) | temp_reload_exp_weak_rates |
temperature below which one should not use theoretical weak rates so they are replaced with exp. weak rates (min 1.d-2) More... | |
real(r_kind) | freeze_rate_temp |
Tmperature at which rates get frozen (for reacl. rates this should be 1d-2GK) More... | |
integer | nuflag |
defines type of neutrino reactions used More... | |
integer | fissflag |
defines type of fission fragment distribution used More... | |
integer | expansiontype |
defines prescription used for parametrized expansion after the last timestep of the hydro input More... | |
integer | extrapolation_width |
how many points from the end of trajectory to use when computing residual expansion More... | |
integer | nse_calc_every |
Compute NSE abundances every x step. More... | |
character *20 | trajectory_mode |
determines how trajectory is calculated (from_file|analytic|expand) More... | |
logical | read_initial_composition |
specify whether initial distribution of abundances should be read from file More... | |
logical | use_tabulated_rates |
switch for using tabulated rates (e.g. talysNGrates.dat) More... | |
logical | use_beta_decay_file |
switch for using different format for beta decays More... | |
logical | use_prepared_network |
Use a prepared folder with all necessary data in binary format. More... | |
character(max_fname_len) | beta_decay_src_ignore |
Source flag(s) to ignore within the beta decay file. More... | |
logical | use_timmes_mue |
Use electron chemical potentials from timmes EOS for theoretical weak rates. More... | |
logical | use_detailed_balance |
Calculate the inverse reactions via detailed balance rather than using them form file. More... | |
logical | use_detailed_balance_q_reac |
Use Q-value from reaclib for the calculation of detailed balance. More... | |
logical | use_thermal_nu_loss |
Whether to include thermal neutrino loss or not. More... | |
integer | nu_loss_every |
Output neutrino loss and gain. More... | |
integer | h_nu_loss_every |
Output neutrino loss and gain in hdf5 format. More... | |
character(max_fname_len) | detailed_balance_src_ignore |
Source flag(s) to ignore within calculated detailed balance. More... | |
character(max_fname_len) | detailed_balance_src_q_reac |
Source flag(s) to use q-value from rate file for inverse reaction. More... | |
character(max_fname_len) | detailed_balance_src_q_winvn |
Source flag(s) to use q-value from winvn file for inverse reaction. More... | |
logical | use_neutrino_loss_file |
Use a file with Qnu values? More... | |
character(max_fname_len) | neutrino_loss_file |
Path to a file containing Qnu values. More... | |
logical | custom_snapshots |
If true, a file must be provided with numbers in days. Snapshots will be created for these points in time. More... | |
logical | h_custom_snapshots |
Same, but in hdf5 format. More... | |
real(r_kind) | engen |
total energy generation [MeV/s] More... | |
real(r_kind) | engen_alpha |
energy generation from alpha-decays [MeV/s] More... | |
real(r_kind) | engen_beta |
energy generation from beta-decays [MeV/s] More... | |
real(r_kind) | engen_fiss |
energy generation from fission [MeV/s] More... | |
real(r_kind) | initemp_cold |
T [GK] lowest allowed temperature to start the calculation from. More... | |
real(r_kind) | initemp_hot |
T [GK] for the starting point of the trajectory: =0: from the beginning; >0: from the last T>initemp. More... | |
real(r_kind) | nse_delt_t9min |
Minimum temperature [GK] when descending to desired temperature in NSE. More... | |
real(r_kind) | nse_nr_tol |
Tolerance for the NR loop in the NSE calculation. More... | |
integer | nse_max_it |
Maximum amount of NSE iterations. More... | |
integer | nse_solver |
Solver for calculating NSE. 0: Newton-Raphson, 1: Powell's hybrid method. More... | |
real(r_kind) | nsetemp_cold |
T [GK] for the nse->network switch. More... | |
real(r_kind) | nsetemp_hot |
T [GK] for the nse<-network switch. More... | |
logical | calc_nsep_energy |
calculate neutron separation energy? More... | |
logical | h_engen_detailed |
Output the energy per parent nucleus and reaction type. More... | |
real(r_kind) | heating_frac |
use this fraction of nuclear-generated energy for heating More... | |
real(r_kind) | heating_density |
Density at which nuclear heating will be switched on (-1) to always include heating. More... | |
real(r_kind) | nse_descend_t9start |
high initial temperature in GK for winnse_descend subroutine More... | |
real(r_kind) | t_analytic |
for parameteric trajectories: initial time More... | |
integer | termination_criterion |
condition to terminate the simulation ([0]=trajectory_file, 1=final_time, 2=final_temp, 3=final_dens, 4=neutron freeze-out) More... | |
real(r_kind) | initial_stepsize |
this value is used as a stepsize at initial step More... | |
real(r_kind) | final_time |
termination time in seconds More... | |
real(r_kind) | final_temp |
termination temperature [GK] More... | |
real(r_kind) | final_dens |
termination density [g/cm3] More... | |
real(r_kind) | timestep_max |
Maximum factor for the change of the timestep. The new timestep is only allowed to be timestep_max * old_timestep. Default value is 2. More... | |
real(r_kind) | timestep_factor |
Factor for the change of the timestep (see nu in Winteler 2012 Eq. 2.49). Default value is 1.0d-1. More... | |
real(r_kind) | timestep_hydro_factor |
Factor for the maximum change of the hydrodynamic quantities (density and temperature) More... | |
real(r_kind) | timestep_ymin |
Lower limit of the abundance to contribute to the timestep calculation, default value is 1.0d-10. More... | |
real(r_kind) | gear_eps |
Abundance accuracy for gear solver. More... | |
real(r_kind) | gear_escale |
Normalization cutoff for gear solver, similar to timestep_Ymin for Euler. More... | |
real(r_kind) | gear_cfactor |
Conservative timestep factor for gear solver [0.1, ... , 0.4]. More... | |
real(r_kind) | gear_nr_eps |
Convergence criterion for the newton-raphson of the gear solver. More... | |
real(r_kind) | gear_timestep_max |
For gear solver. More... | |
integer | gear_nr_maxcount |
Maximum newton-raphson iterations for gear solver. More... | |
integer | gear_nr_mincount |
Minimum newton-raphson iterations for gear solver. More... | |
logical | gear_ignore_adapt_stepsize |
Flag whether gear should ignore the adapt stepsize loop. More... | |
logical | timestep_traj_limit |
Should the timestep be limited by the timestep of the trajectory. More... | |
logical | use_htpf |
Use high temperature partition functions or not. More... | |
logical | h_finab |
Store the finab in hdf5 format rather than in ascii format. More... | |
integer | solver |
solver flag (0 - implicit Euler, 1 - Gear's method, ...), is integer as it is faster than comparing strings every timestep More... | |
integer | heating_mode |
Mode for heating: 0 - no heating, 1 - heating using an entropy equation, 2 - heating from the energy generation and trajectory changes. More... | |
real(r_kind) | heating_t9_tol |
Convergence parameter for the temperature in the heating mode. More... | |
integer | screening_mode |
Mode for coulomb corrections: 0 - no screening, 1 - screening using the prescription of Kravchuk & Yakovlev 2014. More... | |
integer | interp_mode |
Mode for interpolation of temperature and density. More... | |
character(max_fname_len) | trajectory_file |
name of trajectory data file More... | |
character(max_fname_len) | seed_file |
name of file with initial seeds for trajectory run More... | |
character(max_fname_len) | seed_format |
Seed format. More... | |
character(max_fname_len) | net_source |
list of isotopes included in the network (sunet) More... | |
character(max_fname_len) | isotopes_file |
properties of all isotopes in the network: masses, partition functions etc. (winvn) More... | |
character(max_fname_len) | htpf_file |
high-temperature partition functions (htpf.dat) More... | |
character(max_fname_len) | reaclib_file |
reaction rate library (reaclib) More... | |
character(max_fname_len) | fission_rates_spontaneous |
reaction library for spontaneous fission rates More... | |
character(max_fname_len) | fission_rates_n_induced |
reaction library for neutron induced fission rates More... | |
character(max_fname_len) | fission_rates_beta_delayed |
reaction library for beta delayed fission rates More... | |
integer | fission_format_beta_delayed |
Format of beta-delayed fission rates (0: Off, 1: Reaclib, 2: Halflifes, 3: Probability) More... | |
integer | fission_format_n_induced |
Format of neutron-induced fission rates (0: Off, 1: Reaclib) More... | |
integer | fission_format_spontaneous |
Format of spontaneous fission rates (0: Off, 1: Reaclib, 2: Halflifes) More... | |
integer | fission_frag_spontaneous |
Fragment distribution of spontaneous fission rates in case of custom fragments (fissflag=4) More... | |
integer | fission_frag_n_induced |
Fragment distribution of n-induced fission rates in case of custom fragments (fissflag=4) More... | |
integer | fission_frag_beta_delayed |
Fragment distribution of beta-delayed fission rates in case of custom fragments (fissflag=4) More... | |
integer | fission_frag_missing |
Fragment distribution in case of missing fragments in case of custom fragments (fissflag=4) More... | |
character(max_fname_len) | weak_rates_file |
weak rates library (twr.dat) More... | |
character(max_fname_len) | tabulated_rates_file |
tabulated rates library (e.g. talysNGrates.dat) More... | |
character(max_fname_len) | tabulated_temperature_file |
file containing grid of tabulated temperature file ("default" for default grid) More... | |
character(max_fname_len) | chem_pot_file |
tabulated chemical potential of electron gas More... | |
character(max_fname_len) | nsep_energies_file |
neutron separation energies More... | |
character(max_fname_len) | nunucleo_rates_file |
neutrino reaction rates on nucleons More... | |
character(max_fname_len) | nuchannel_file |
Contains neutrino channel information as in Sieverding et al. 2018. More... | |
character(max_fname_len) | nurates_file |
Neutrino reactions on heavy nuclei as in Sieverding et al. 2018. More... | |
character(max_fname_len) | snapshot_file |
File that contains days, where a snapshot should be made. More... | |
character(max_fname_len) | bfission_file |
Fission table for beta-delayed fission. More... | |
character(max_fname_len) | nfission_file |
Fission table for neutron-induced fission. More... | |
character(max_fname_len) | sfission_file |
Fission table for spontaneous fission. More... | |
character(max_fname_len) | track_nuclei_file |
File of nuclei to track. Gives an output similar to mainout.dat. More... | |
character(max_fname_len) | alpha_decay_file |
File with additional alpha decays. More... | |
character(max_fname_len) | beta_decay_file |
File for reading in beta decays in different format. More... | |
character(max_fname_len) | trajectory_format |
Format to read the trajectory. More... | |
character(max_fname_len) | neutrino_mode |
Similar to trajectory mode. More... | |
character(max_fname_len) | t9_analytic |
analytic temperature [T9] More... | |
character(max_fname_len) | rho_analytic |
analytic density [g/cm3] More... | |
character(max_fname_len) | rkm_analytic |
analytic radial scale [km] More... | |
character(max_fname_len) | ye_analytic |
analytic electron fraction More... | |
character(max_fname_len) | le |
electron-neutrino luminosities [erg/s] More... | |
character(max_fname_len) | lebar |
electron-antineutrino luminosities [erg/s] More... | |
character(max_fname_len) | enue |
average electron-neutrino energies [MeV] More... | |
character(max_fname_len) | enuebar |
average electron-antineutrino energies [MeV] More... | |
character(max_fname_len) | lx |
Muon and Tauon neutrino luminosities [erg/s]. More... | |
character(max_fname_len) | lxbar |
Muon and Tauon antineutrino luminosities [erg/s]. More... | |
character(max_fname_len) | enux |
average Muon and Tauon neutrino energies [MeV] More... | |
character(max_fname_len) | enuxbar |
average Muon and Tauon antineutrino energies [MeV] More... | |
real(r_kind) | nu_max_time |
Maximum time for neutrino reactions to react. More... | |
real(r_kind) | nu_min_t |
Neutrino temperature cutoff for neutrino reactions [MeV]. More... | |
real(r_kind) | nu_min_l |
Neutrino luminosity cutoff for neutrino reactions [erg/g/s]. More... | |
character(max_fname_len) | prepared_network_path |
Prepared network folder. More... | |
integer | nr_maxcount |
– Newton-Raphson iterative loop parameters More... | |
integer | nr_mincount |
Minimum iterations in NR. More... | |
real(r_kind) | nr_tol |
exit NR if tolerance less than this value More... | |
integer | adapt_stepsize_maxcount |
max. iterations in adapting the stepsize More... | |
character(len= *), parameter, private | par_binary_name ='parameter.windat' |
other static variables More... | |
real(r_kind), dimension(:), allocatable | weights |
– parameters for efficient numerical integration of effphase in the interval [1,infinity] More... | |
real(r_kind), dimension(:), allocatable | xnodes |
corresponding nodes " More... | |
logical | weightscalculated = .false. |
switch to calculated weights and nodes for [1,infinity] More... | |
integer | ncc = 256 |
nr of points for Clenshaw-Curtis integration More... | |
real(r_kind), dimension(:), allocatable | dcc |
real(r_kind), dimension(:), allocatable | matcc |
matrices for Clenshaw-Curtis integration More... | |
real(r_kind), dimension(:,:), allocatable | mcc |
Contains all runtime parameters as well as phys and math constants.
Contains the following definitions:
Editors:
subroutine parameter_class::check_param |
Check for consistency of parameter.
Some parameters might not work together or might cause inconsistencies. This routine checks the user input and complains.
Edited:
Definition at line 1233 of file parameter_class.f90.
real(r_kind) function parameter_class::levenshteindistance | ( | character(*), intent(in) | par1, |
character(*), intent(in) | par2 | ||
) |
Calculates the Levenshtein distance between two strings.
The Levenshtein distance gives a measure of the similarity between two strings. The similarity is expressed in terms of a float afterwards. This routine follows the pseudocode here
[in] | par1 | Input string 1 |
[in] | par2 | Input string 2 |
Definition at line 1357 of file parameter_class.f90.
subroutine parameter_class::output_binary_parameter_data | ( | character(len=*), intent(in) | path | ) |
Output the parameter data to a binary file.
This routine outputs important parameters to a binary file. This is exclusively done to check later when reading the binary rate files.
[in] | path | Path to the output directory |
Definition at line 1672 of file parameter_class.f90.
subroutine parameter_class::output_param |
Output parameters to a file.
Parameters are sorted in alphabetical order
Definition at line 1064 of file parameter_class.f90.
subroutine parameter_class::output_param_prepared_network | ( | character(len=*), intent(in) | path | ) |
Output relevant parameters to a file in the prepared network path.
Parameters are sorted according to category
[in] | path | Path to the output directory |
Definition at line 1541 of file parameter_class.f90.
real(r_kind) function parameter_class::read_float_param | ( | character(len=*), intent(in) | input_string, |
character(len=*), intent(in) | param_name | ||
) |
Converts a string to an float.
If the string is not a valid float, an error message is raised.
[in] | input_string | String from param file |
[in] | param_name | Parameter name |
Definition at line 863 of file parameter_class.f90.
integer function parameter_class::read_integer_param | ( | character(len=*), intent(in) | input_string, |
character(len=*), intent(in) | param_name | ||
) |
Converts a string to an integer.
If the string is not a valid integer, an error message is raised.
[in] | input_string | String from param file |
[in] | param_name | Parameter name |
Definition at line 834 of file parameter_class.f90.
subroutine parameter_class::read_param | ( | character(*), intent(in) | parfile | ) |
This function reads the parameter file.
Blank lines as well as comment lines are skipped. The parameters are read into the proper variable that are declarated at the beginning of the file.
Edited:
Definition at line 270 of file parameter_class.f90.
subroutine parameter_class::read_sanity_check_prepared_network |
Read important parameters from a binary file and check for consistency.
This routine reads important parameters from a binary file and checks if they are consistent with the current run. This is done to check if the binary files are compatible with the current parameter file.
Definition at line 1438 of file parameter_class.f90.
subroutine parameter_class::search_close_parameter | ( | character(*), intent(in) | possible_pars, |
character(*), intent(in) | input_par, | ||
real(r_kind), intent(out) | score, | ||
character(*), intent(out) | cl_par | ||
) |
Search for similar existing parameters.
Takes a list of parameters that are separated by ":" and a comparison parameter and returns the most similar parameter from the list. Furthermore, a score is returned that gives a measure of the similarity of the strings. This function is only called if a user specified parameter was not included in the list of possible parameters.
[in] | possible_pars | All possible parameter, long string |
[in] | input_par | Parameter to compare |
[out] | cl_par | Most similar parameter |
[out] | score | Minimum distance between parameters |
Definition at line 1302 of file parameter_class.f90.
subroutine parameter_class::set_default_param |
Sets default parameters Parameters are sorted in alphabetical order.
Definition at line 892 of file parameter_class.f90.
subroutine parameter_class::set_param | ( | character(*), intent(in) | param_name, |
character(*), intent(in) | param_value | ||
) |
Sets a global parameter param_name to the value, given by its string representation param_value.
Edited:
Definition at line 332 of file parameter_class.f90.
subroutine, public parameter_class::unit_define |
Declares values for the elements in unit_type.
Definition at line 235 of file parameter_class.f90.
integer parameter_class::adapt_stepsize_maxcount |
max. iterations in adapting the stepsize
Definition at line 207 of file parameter_class.f90.
character(max_fname_len) parameter_class::alpha_decay_file |
File with additional alpha decays.
Definition at line 182 of file parameter_class.f90.
logical parameter_class::alpha_decay_ignore_all |
Flag whether rates should actually get replaced or only added.
Definition at line 80 of file parameter_class.f90.
character(max_fname_len) parameter_class::alpha_decay_src_ignore |
Source flag(s) to ignore within the alpha decay rates.
Definition at line 79 of file parameter_class.f90.
integer parameter_class::alpha_decay_zmax |
Maximum Z for additional alpha decay rates.
Definition at line 78 of file parameter_class.f90.
integer parameter_class::alpha_decay_zmin |
Minimum Z for additional alpha decay rates.
Definition at line 77 of file parameter_class.f90.
character(max_fname_len) parameter_class::beta_decay_file |
File for reading in beta decays in different format.
Definition at line 183 of file parameter_class.f90.
character(max_fname_len) parameter_class::beta_decay_src_ignore |
Source flag(s) to ignore within the beta decay file.
Definition at line 95 of file parameter_class.f90.
character(max_fname_len) parameter_class::bfission_file |
Fission table for beta-delayed fission.
Definition at line 178 of file parameter_class.f90.
logical parameter_class::calc_nsep_energy |
calculate neutron separation energy?
Definition at line 121 of file parameter_class.f90.
character(max_fname_len) parameter_class::chem_pot_file |
tabulated chemical potential of electron gas
Definition at line 172 of file parameter_class.f90.
logical parameter_class::custom_snapshots |
If true, a file must be provided with numbers in days. Snapshots will be created for these points in time.
Definition at line 107 of file parameter_class.f90.
real(r_kind), dimension(:), allocatable parameter_class::dcc |
Definition at line 219 of file parameter_class.f90.
character(max_fname_len) parameter_class::detailed_balance_src_ignore |
Source flag(s) to ignore within calculated detailed balance.
Definition at line 102 of file parameter_class.f90.
character(max_fname_len) parameter_class::detailed_balance_src_q_reac |
Source flag(s) to use q-value from rate file for inverse reaction.
Definition at line 103 of file parameter_class.f90.
character(max_fname_len) parameter_class::detailed_balance_src_q_winvn |
Source flag(s) to use q-value from winvn file for inverse reaction.
Definition at line 104 of file parameter_class.f90.
real(r_kind) parameter_class::engen |
total energy generation [MeV/s]
Definition at line 109 of file parameter_class.f90.
real(r_kind) parameter_class::engen_alpha |
energy generation from alpha-decays [MeV/s]
Definition at line 110 of file parameter_class.f90.
real(r_kind) parameter_class::engen_beta |
energy generation from beta-decays [MeV/s]
Definition at line 111 of file parameter_class.f90.
integer parameter_class::engen_every |
Energy generation output frequency.
Definition at line 68 of file parameter_class.f90.
real(r_kind) parameter_class::engen_fiss |
energy generation from fission [MeV/s]
Definition at line 112 of file parameter_class.f90.
character(max_fname_len) parameter_class::enue |
average electron-neutrino energies [MeV]
Definition at line 192 of file parameter_class.f90.
character(max_fname_len) parameter_class::enuebar |
average electron-antineutrino energies [MeV]
Definition at line 193 of file parameter_class.f90.
character(max_fname_len) parameter_class::enux |
average Muon and Tauon neutrino energies [MeV]
Definition at line 196 of file parameter_class.f90.
character(max_fname_len) parameter_class::enuxbar |
average Muon and Tauon antineutrino energies [MeV]
Definition at line 197 of file parameter_class.f90.
integer parameter_class::expansiontype |
defines prescription used for parametrized expansion after the last timestep of the hydro input
Definition at line 87 of file parameter_class.f90.
integer parameter_class::extrapolation_width |
how many points from the end of trajectory to use when computing residual expansion
Definition at line 88 of file parameter_class.f90.
real(r_kind) parameter_class::final_dens |
termination density [g/cm3]
Definition at line 131 of file parameter_class.f90.
real(r_kind) parameter_class::final_temp |
termination temperature [GK]
Definition at line 130 of file parameter_class.f90.
real(r_kind) parameter_class::final_time |
termination time in seconds
Definition at line 129 of file parameter_class.f90.
integer parameter_class::fissflag |
defines type of fission fragment distribution used
Definition at line 86 of file parameter_class.f90.
integer parameter_class::fission_format_beta_delayed |
Format of beta-delayed fission rates (0: Off, 1: Reaclib, 2: Halflifes, 3: Probability)
Definition at line 162 of file parameter_class.f90.
integer parameter_class::fission_format_n_induced |
Format of neutron-induced fission rates (0: Off, 1: Reaclib)
Definition at line 163 of file parameter_class.f90.
integer parameter_class::fission_format_spontaneous |
Format of spontaneous fission rates (0: Off, 1: Reaclib, 2: Halflifes)
Definition at line 164 of file parameter_class.f90.
integer parameter_class::fission_frag_beta_delayed |
Fragment distribution of beta-delayed fission rates in case of custom fragments (fissflag=4)
Definition at line 167 of file parameter_class.f90.
integer parameter_class::fission_frag_missing |
Fragment distribution in case of missing fragments in case of custom fragments (fissflag=4)
Definition at line 168 of file parameter_class.f90.
integer parameter_class::fission_frag_n_induced |
Fragment distribution of n-induced fission rates in case of custom fragments (fissflag=4)
Definition at line 166 of file parameter_class.f90.
integer parameter_class::fission_frag_spontaneous |
Fragment distribution of spontaneous fission rates in case of custom fragments (fissflag=4)
Definition at line 165 of file parameter_class.f90.
character(max_fname_len) parameter_class::fission_rates_beta_delayed |
reaction library for beta delayed fission rates
Definition at line 161 of file parameter_class.f90.
character(max_fname_len) parameter_class::fission_rates_n_induced |
reaction library for neutron induced fission rates
Definition at line 160 of file parameter_class.f90.
character(max_fname_len) parameter_class::fission_rates_spontaneous |
reaction library for spontaneous fission rates
Definition at line 159 of file parameter_class.f90.
integer parameter_class::flow_every |
flow output frequency
Definition at line 64 of file parameter_class.f90.
real(r_kind) parameter_class::freeze_rate_temp |
Tmperature at which rates get frozen (for reacl. rates this should be 1d-2GK)
Definition at line 84 of file parameter_class.f90.
real(r_kind) parameter_class::gear_cfactor |
Conservative timestep factor for gear solver [0.1, ... , 0.4].
Definition at line 138 of file parameter_class.f90.
real(r_kind) parameter_class::gear_eps |
Abundance accuracy for gear solver.
Definition at line 136 of file parameter_class.f90.
real(r_kind) parameter_class::gear_escale |
Normalization cutoff for gear solver, similar to timestep_Ymin for Euler.
Definition at line 137 of file parameter_class.f90.
logical parameter_class::gear_ignore_adapt_stepsize |
Flag whether gear should ignore the adapt stepsize loop.
Definition at line 143 of file parameter_class.f90.
real(r_kind) parameter_class::gear_nr_eps |
Convergence criterion for the newton-raphson of the gear solver.
Definition at line 139 of file parameter_class.f90.
integer parameter_class::gear_nr_maxcount |
Maximum newton-raphson iterations for gear solver.
Definition at line 141 of file parameter_class.f90.
integer parameter_class::gear_nr_mincount |
Minimum newton-raphson iterations for gear solver.
Definition at line 142 of file parameter_class.f90.
real(r_kind) parameter_class::gear_timestep_max |
For gear solver.
Definition at line 140 of file parameter_class.f90.
logical parameter_class::h_custom_snapshots |
Same, but in hdf5 format.
Definition at line 108 of file parameter_class.f90.
logical parameter_class::h_engen_detailed |
Output the energy per parent nucleus and reaction type.
Definition at line 122 of file parameter_class.f90.
integer parameter_class::h_engen_every |
Energy generation output frequency in hdf5 format.
Definition at line 69 of file parameter_class.f90.
logical parameter_class::h_finab |
Store the finab in hdf5 format rather than in ascii format.
Definition at line 146 of file parameter_class.f90.
integer parameter_class::h_flow_every |
flow output frequency in hdf5 format
Definition at line 65 of file parameter_class.f90.
integer parameter_class::h_mainout_every |
HDF5 output frequency of the mainout.
Definition at line 72 of file parameter_class.f90.
integer parameter_class::h_nu_loss_every |
Output neutrino loss and gain in hdf5 format.
Definition at line 101 of file parameter_class.f90.
integer parameter_class::h_snapshot_every |
snapshot output frequency in hdf5 format
Definition at line 63 of file parameter_class.f90.
integer parameter_class::h_timescales_every |
timescales output frequency in hdf5 format
Definition at line 67 of file parameter_class.f90.
integer parameter_class::h_track_nuclei_every |
frequency of track nuclei output in hdf5 format
Definition at line 74 of file parameter_class.f90.
real(r_kind) parameter_class::heating_density |
Density at which nuclear heating will be switched on (-1) to always include heating.
Definition at line 124 of file parameter_class.f90.
real(r_kind) parameter_class::heating_frac |
use this fraction of nuclear-generated energy for heating
Definition at line 123 of file parameter_class.f90.
integer parameter_class::heating_mode |
Mode for heating: 0 - no heating, 1 - heating using an entropy equation, 2 - heating from the energy generation and trajectory changes.
Definition at line 148 of file parameter_class.f90.
real(r_kind) parameter_class::heating_t9_tol |
Convergence parameter for the temperature in the heating mode.
Definition at line 149 of file parameter_class.f90.
character(max_fname_len) parameter_class::htpf_file |
high-temperature partition functions (htpf.dat)
Definition at line 157 of file parameter_class.f90.
real(r_kind) parameter_class::initemp_cold |
T [GK] lowest allowed temperature to start the calculation from.
Definition at line 113 of file parameter_class.f90.
real(r_kind) parameter_class::initemp_hot |
T [GK] for the starting point of the trajectory: =0: from the beginning; >0: from the last T>initemp.
Definition at line 114 of file parameter_class.f90.
real(r_kind) parameter_class::initial_stepsize |
this value is used as a stepsize at initial step
Definition at line 128 of file parameter_class.f90.
integer parameter_class::interp_mode |
Mode for interpolation of temperature and density.
Definition at line 151 of file parameter_class.f90.
character(max_fname_len) parameter_class::isotopes_file |
properties of all isotopes in the network: masses, partition functions etc. (winvn)
Definition at line 156 of file parameter_class.f90.
integer parameter_class::iwformat |
defines format of the weak rates (0 = tabulated, 1 = log<ft>)
Definition at line 81 of file parameter_class.f90.
integer parameter_class::iwinterp |
defines the interpolation for the weak rates (0 = bilinear, 1 = bicubic)
Definition at line 82 of file parameter_class.f90.
character(max_fname_len) parameter_class::le |
electron-neutrino luminosities [erg/s]
Definition at line 190 of file parameter_class.f90.
character(max_fname_len) parameter_class::lebar |
electron-antineutrino luminosities [erg/s]
Definition at line 191 of file parameter_class.f90.
character(max_fname_len) parameter_class::lx |
Muon and Tauon neutrino luminosities [erg/s].
Definition at line 194 of file parameter_class.f90.
character(max_fname_len) parameter_class::lxbar |
Muon and Tauon antineutrino luminosities [erg/s].
Definition at line 195 of file parameter_class.f90.
integer parameter_class::mainout_every |
frequency of mainout output
Definition at line 71 of file parameter_class.f90.
real(r_kind), dimension(:), allocatable parameter_class::matcc |
matrices for Clenshaw-Curtis integration
Definition at line 219 of file parameter_class.f90.
integer, parameter, public parameter_class::max_fname_len = 400 |
maximum length of filenames
Definition at line 58 of file parameter_class.f90.
real(r_kind), dimension(:,:), allocatable parameter_class::mcc |
Definition at line 220 of file parameter_class.f90.
integer parameter_class::ncc = 256 |
nr of points for Clenshaw-Curtis integration
Definition at line 218 of file parameter_class.f90.
character(max_fname_len) parameter_class::net_source |
list of isotopes included in the network (sunet)
Definition at line 155 of file parameter_class.f90.
character(max_fname_len) parameter_class::neutrino_loss_file |
Path to a file containing Qnu values.
Definition at line 106 of file parameter_class.f90.
character(max_fname_len) parameter_class::neutrino_mode |
Similar to trajectory mode.
Definition at line 185 of file parameter_class.f90.
character(max_fname_len) parameter_class::nfission_file |
Fission table for neutron-induced fission.
Definition at line 179 of file parameter_class.f90.
integer parameter_class::nr_maxcount |
– Newton-Raphson iterative loop parameters
no more that this many iterations in NR
Definition at line 204 of file parameter_class.f90.
integer parameter_class::nr_mincount |
Minimum iterations in NR.
Definition at line 205 of file parameter_class.f90.
real(r_kind) parameter_class::nr_tol |
exit NR if tolerance less than this value
Definition at line 206 of file parameter_class.f90.
integer parameter_class::nrdiag_every |
frequency of NR loop diagnostic messages (wrt iteration counter)
Definition at line 70 of file parameter_class.f90.
integer parameter_class::nse_calc_every |
Compute NSE abundances every x step.
Definition at line 89 of file parameter_class.f90.
real(r_kind) parameter_class::nse_delt_t9min |
Minimum temperature [GK] when descending to desired temperature in NSE.
Definition at line 115 of file parameter_class.f90.
real(r_kind) parameter_class::nse_descend_t9start |
high initial temperature in GK for winnse_descend subroutine
Definition at line 125 of file parameter_class.f90.
integer parameter_class::nse_max_it |
Maximum amount of NSE iterations.
Definition at line 117 of file parameter_class.f90.
real(r_kind) parameter_class::nse_nr_tol |
Tolerance for the NR loop in the NSE calculation.
Definition at line 116 of file parameter_class.f90.
integer parameter_class::nse_solver |
Solver for calculating NSE. 0: Newton-Raphson, 1: Powell's hybrid method.
Definition at line 118 of file parameter_class.f90.
character(max_fname_len) parameter_class::nsep_energies_file |
neutron separation energies
Definition at line 173 of file parameter_class.f90.
real(r_kind) parameter_class::nsetemp_cold |
T [GK] for the nse->network switch.
Definition at line 119 of file parameter_class.f90.
real(r_kind) parameter_class::nsetemp_hot |
T [GK] for the nse<-network switch.
Definition at line 120 of file parameter_class.f90.
integer parameter_class::nu_loss_every |
Output neutrino loss and gain.
Definition at line 100 of file parameter_class.f90.
real(r_kind) parameter_class::nu_max_time |
Maximum time for neutrino reactions to react.
Definition at line 198 of file parameter_class.f90.
real(r_kind) parameter_class::nu_min_l |
Neutrino luminosity cutoff for neutrino reactions [erg/g/s].
Definition at line 200 of file parameter_class.f90.
real(r_kind) parameter_class::nu_min_t |
Neutrino temperature cutoff for neutrino reactions [MeV].
Definition at line 199 of file parameter_class.f90.
character(max_fname_len) parameter_class::nuchannel_file |
Contains neutrino channel information as in Sieverding et al. 2018.
Definition at line 175 of file parameter_class.f90.
integer parameter_class::nuflag |
defines type of neutrino reactions used
Definition at line 85 of file parameter_class.f90.
character(max_fname_len) parameter_class::nunucleo_rates_file |
neutrino reaction rates on nucleons
Definition at line 174 of file parameter_class.f90.
character(max_fname_len) parameter_class::nurates_file |
Neutrino reactions on heavy nuclei as in Sieverding et al. 2018.
Definition at line 176 of file parameter_class.f90.
integer parameter_class::out_every |
– runtime parameters set when reading parameter file
main output frequency
Definition at line 61 of file parameter_class.f90.
|
private |
other static variables
Definition at line 210 of file parameter_class.f90.
integer, parameter, public parameter_class::param_name_len = 200 |
– hardcoded parameters
maximum length of a parameter name
Definition at line 57 of file parameter_class.f90.
character(max_fname_len) parameter_class::prepared_network_path |
Prepared network folder.
Definition at line 201 of file parameter_class.f90.
character(max_fname_len) parameter_class::reaclib_file |
reaction rate library (reaclib)
Definition at line 158 of file parameter_class.f90.
logical parameter_class::read_initial_composition |
specify whether initial distribution of abundances should be read from file
Definition at line 91 of file parameter_class.f90.
character(max_fname_len) parameter_class::rho_analytic |
analytic density [g/cm3]
Definition at line 187 of file parameter_class.f90.
character(max_fname_len) parameter_class::rkm_analytic |
analytic radial scale [km]
Definition at line 188 of file parameter_class.f90.
integer parameter_class::screening_mode |
Mode for coulomb corrections: 0 - no screening, 1 - screening using the prescription of Kravchuk & Yakovlev 2014.
Definition at line 150 of file parameter_class.f90.
character(max_fname_len) parameter_class::seed_file |
name of file with initial seeds for trajectory run
Definition at line 153 of file parameter_class.f90.
character(max_fname_len) parameter_class::seed_format |
Seed format.
Definition at line 154 of file parameter_class.f90.
character(max_fname_len) parameter_class::sfission_file |
Fission table for spontaneous fission.
Definition at line 180 of file parameter_class.f90.
integer parameter_class::snapshot_every |
snapshot output frequency
Definition at line 62 of file parameter_class.f90.
character(max_fname_len) parameter_class::snapshot_file |
File that contains days, where a snapshot should be made.
Definition at line 177 of file parameter_class.f90.
integer parameter_class::solver |
solver flag (0 - implicit Euler, 1 - Gear's method, ...), is integer as it is faster than comparing strings every timestep
Definition at line 147 of file parameter_class.f90.
character(max_fname_len) parameter_class::t9_analytic |
analytic temperature [T9]
Definition at line 186 of file parameter_class.f90.
real(r_kind) parameter_class::t_analytic |
for parameteric trajectories: initial time
Definition at line 126 of file parameter_class.f90.
character(max_fname_len) parameter_class::tabulated_rates_file |
tabulated rates library (e.g. talysNGrates.dat)
Definition at line 170 of file parameter_class.f90.
character(max_fname_len) parameter_class::tabulated_temperature_file |
file containing grid of tabulated temperature file ("default" for default grid)
Definition at line 171 of file parameter_class.f90.
real(r_kind) parameter_class::temp_reload_exp_weak_rates |
temperature below which one should not use theoretical weak rates so they are replaced with exp. weak rates (min 1.d-2)
Definition at line 83 of file parameter_class.f90.
integer parameter_class::termination_criterion |
condition to terminate the simulation ([0]=trajectory_file, 1=final_time, 2=final_temp, 3=final_dens, 4=neutron freeze-out)
Definition at line 127 of file parameter_class.f90.
integer parameter_class::timescales_every |
timescales output frequency
Definition at line 66 of file parameter_class.f90.
real(r_kind) parameter_class::timestep_factor |
Factor for the change of the timestep (see nu in Winteler 2012 Eq. 2.49). Default value is 1.0d-1.
Definition at line 133 of file parameter_class.f90.
real(r_kind) parameter_class::timestep_hydro_factor |
Factor for the maximum change of the hydrodynamic quantities (density and temperature)
Definition at line 134 of file parameter_class.f90.
real(r_kind) parameter_class::timestep_max |
Maximum factor for the change of the timestep. The new timestep is only allowed to be timestep_max * old_timestep. Default value is 2.
Definition at line 132 of file parameter_class.f90.
logical parameter_class::timestep_traj_limit |
Should the timestep be limited by the timestep of the trajectory.
Definition at line 144 of file parameter_class.f90.
real(r_kind) parameter_class::timestep_ymin |
Lower limit of the abundance to contribute to the timestep calculation, default value is 1.0d-10.
Definition at line 135 of file parameter_class.f90.
integer parameter_class::top_engen_every |
frequency of energy generators toplist
Definition at line 75 of file parameter_class.f90.
integer parameter_class::track_nuclei_every |
frequency of track nuclei output
Definition at line 73 of file parameter_class.f90.
character(max_fname_len) parameter_class::track_nuclei_file |
File of nuclei to track. Gives an output similar to mainout.dat.
Definition at line 181 of file parameter_class.f90.
character(max_fname_len) parameter_class::trajectory_file |
name of trajectory data file
Definition at line 152 of file parameter_class.f90.
character(max_fname_len) parameter_class::trajectory_format |
Format to read the trajectory.
Definition at line 184 of file parameter_class.f90.
character*20 parameter_class::trajectory_mode |
determines how trajectory is calculated (from_file|analytic|expand)
Definition at line 90 of file parameter_class.f90.
type(unit_type), public parameter_class::unit |
constants and unit conversion factors
Definition at line 54 of file parameter_class.f90.
logical parameter_class::use_alpha_decay_file |
Switch for using additional alpha decay rates.
Definition at line 76 of file parameter_class.f90.
logical parameter_class::use_beta_decay_file |
switch for using different format for beta decays
Definition at line 93 of file parameter_class.f90.
logical parameter_class::use_detailed_balance |
Calculate the inverse reactions via detailed balance rather than using them form file.
Definition at line 97 of file parameter_class.f90.
logical parameter_class::use_detailed_balance_q_reac |
Use Q-value from reaclib for the calculation of detailed balance.
Definition at line 98 of file parameter_class.f90.
logical parameter_class::use_htpf |
Use high temperature partition functions or not.
Definition at line 145 of file parameter_class.f90.
logical parameter_class::use_neutrino_loss_file |
Use a file with Qnu values?
Definition at line 105 of file parameter_class.f90.
logical parameter_class::use_prepared_network |
Use a prepared folder with all necessary data in binary format.
Definition at line 94 of file parameter_class.f90.
logical parameter_class::use_tabulated_rates |
switch for using tabulated rates (e.g. talysNGrates.dat)
Definition at line 92 of file parameter_class.f90.
logical parameter_class::use_thermal_nu_loss |
Whether to include thermal neutrino loss or not.
Definition at line 99 of file parameter_class.f90.
logical parameter_class::use_timmes_mue |
Use electron chemical potentials from timmes EOS for theoretical weak rates.
Definition at line 96 of file parameter_class.f90.
character(max_fname_len) parameter_class::weak_rates_file |
weak rates library (twr.dat)
Definition at line 169 of file parameter_class.f90.
real(r_kind), dimension(:), allocatable parameter_class::weights |
– parameters for efficient numerical integration of effphase in the interval [1,infinity]
weights for the numerical integration
Definition at line 214 of file parameter_class.f90.
logical parameter_class::weightscalculated = .false. |
switch to calculated weights and nodes for [1,infinity]
Definition at line 216 of file parameter_class.f90.
real(r_kind), dimension(:), allocatable parameter_class::xnodes |
corresponding nodes "
Definition at line 215 of file parameter_class.f90.
character(max_fname_len) parameter_class::ye_analytic |
analytic electron fraction
Definition at line 189 of file parameter_class.f90.