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...

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

Detailed Description

Contains all runtime parameters as well as phys and math constants.

Contains the following definitions:

all runtime simulation parameters;
mathematical and physical constants.

Author: Darko Mocelj

Date: 22.01.03

Editors:

01.12.11: Christian Winteler
22.07.13: Marcella Ugliano
11.01.14: Oleg Korobkin
11.01.21: Moritz Reichert

Function/Subroutine Documentation

◆ check_param()

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.

Author: : M. Reichert

Edited:

21.05.19, MR

Definition at line 1215 of file parameter_class.f90.

Here is the call graph for this function:

◆ levenshteindistance()

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

Author: M. Reichert

Date: 02.11.22

Parameters

[in]	par1	Input string 1
[in]	par2	Input string 2

Returns: Measure of the similarity

Definition at line 1339 of file parameter_class.f90.

Here is the call graph for this function:

◆ output_binary_parameter_data()

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.

Author: M. Reichert

Date: 21.07.23

Parameters

[in] path Path to the output directory

Definition at line 1654 of file parameter_class.f90.

Here is the call graph for this function:

◆ output_param()

subroutine parameter_class::output_param

Output parameters to a file.

Parameters are sorted in alphabetical order

Definition at line 1049 of file parameter_class.f90.

Here is the call graph for this function:

◆ output_param_prepared_network()

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

Parameters

[in] path Path to the output directory

Definition at line 1523 of file parameter_class.f90.

Here is the call graph for this function:

◆ read_float_param()

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.

Author: Moritz Reichert

Date: 22.01.21

Parameters

[in]	input_string	String from param file
[in]	param_name	Parameter name

Returns: Converted float value from input string

Definition at line 851 of file parameter_class.f90.

Here is the call graph for this function:

◆ read_integer_param()

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.

Author: Moritz Reichert

Date: 22.01.21

Parameters

[in]	input_string	String from param file
[in]	param_name	Parameter name

Returns: Converted integer value from input string

Definition at line 822 of file parameter_class.f90.

Here is the call graph for this function:

◆ read_param()

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:

01.02.14

Definition at line 267 of file parameter_class.f90.

Here is the call graph for this function:

◆ read_sanity_check_prepared_network()

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.

Note: This routine is only called if the parameter "use_prepared_network" is set to ".true."

Author: M. Reichert

Date: 21.07.23

Definition at line 1420 of file parameter_class.f90.

Here is the call graph for this function:

◆ search_close_parameter()

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.

Author: M. Reichert

Date: 02.11.22

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 1284 of file parameter_class.f90.

Here is the call graph for this function:

◆ set_default_param()

subroutine parameter_class::set_default_param

Sets default parameters Parameters are sorted in alphabetical order.

Definition at line 880 of file parameter_class.f90.

◆ set_param()

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:

M.R. 02.11.22, Implemented more meaningfull error msg, made use of *_params

Definition at line 329 of file parameter_class.f90.

Here is the call graph for this function:

◆ unit_define()

subroutine, public parameter_class::unit_define

Declares values for the elements in unit_type.

Definition at line 232 of file parameter_class.f90.

Variable Documentation

◆ adapt_stepsize_maxcount

integer parameter_class::adapt_stepsize_maxcount

max. iterations in adapting the stepsize

Definition at line 204 of file parameter_class.f90.

◆ alpha_decay_file

character(max_fname_len) parameter_class::alpha_decay_file

File with additional alpha decays.

Definition at line 182 of file parameter_class.f90.

◆ alpha_decay_ignore_all

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.

◆ alpha_decay_src_ignore

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.

◆ alpha_decay_zmax

integer parameter_class::alpha_decay_zmax

Maximum Z for additional alpha decay rates.

Definition at line 78 of file parameter_class.f90.

◆ alpha_decay_zmin

integer parameter_class::alpha_decay_zmin

Minimum Z for additional alpha decay rates.

Definition at line 77 of file parameter_class.f90.

◆ beta_decay_file

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.

◆ beta_decay_src_ignore

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.

◆ bfission_file

character(max_fname_len) parameter_class::bfission_file

Fission table for beta-delayed fission.

Definition at line 178 of file parameter_class.f90.

◆ calc_nsep_energy

logical parameter_class::calc_nsep_energy

calculate neutron separation energy?

Definition at line 121 of file parameter_class.f90.

◆ chem_pot_file

character(max_fname_len) parameter_class::chem_pot_file

tabulated chemical potential of electron gas

Definition at line 172 of file parameter_class.f90.

◆ custom_snapshots

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.

◆ dcc

real(r_kind), dimension(:), allocatable parameter_class::dcc

Definition at line 216 of file parameter_class.f90.

◆ detailed_balance_src_ignore

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.

◆ detailed_balance_src_q_reac

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.

◆ detailed_balance_src_q_winvn

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.

◆ engen

real(r_kind) parameter_class::engen

total energy generation [MeV/s]

Definition at line 109 of file parameter_class.f90.

◆ engen_alpha

real(r_kind) parameter_class::engen_alpha

energy generation from alpha-decays [MeV/s]

Definition at line 110 of file parameter_class.f90.

◆ engen_beta

real(r_kind) parameter_class::engen_beta

energy generation from beta-decays [MeV/s]

Definition at line 111 of file parameter_class.f90.

◆ engen_every

integer parameter_class::engen_every

Energy generation output frequency.

Definition at line 68 of file parameter_class.f90.

◆ engen_fiss

real(r_kind) parameter_class::engen_fiss

energy generation from fission [MeV/s]

Definition at line 112 of file parameter_class.f90.

◆ enue

character(max_fname_len) parameter_class::enue

average electron-neutrino energies [MeV]

Definition at line 192 of file parameter_class.f90.

◆ enuebar

character(max_fname_len) parameter_class::enuebar

average electron-antineutrino energies [MeV]

Definition at line 193 of file parameter_class.f90.

◆ enux

character(max_fname_len) parameter_class::enux

average Muon and Tauon neutrino energies [MeV]

Definition at line 196 of file parameter_class.f90.

◆ enuxbar

character(max_fname_len) parameter_class::enuxbar

average Muon and Tauon antineutrino energies [MeV]

Definition at line 197 of file parameter_class.f90.

◆ expansiontype

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.

◆ extrapolation_width

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.

◆ final_dens

real(r_kind) parameter_class::final_dens

termination density [g/cm3]

Definition at line 131 of file parameter_class.f90.

◆ final_temp

real(r_kind) parameter_class::final_temp

termination temperature [GK]

Definition at line 130 of file parameter_class.f90.

◆ final_time

real(r_kind) parameter_class::final_time

termination time in seconds

Definition at line 129 of file parameter_class.f90.

◆ fissflag

integer parameter_class::fissflag

defines type of fission fragment distribution used

Definition at line 86 of file parameter_class.f90.

◆ fission_format_beta_delayed

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.

◆ fission_format_n_induced

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.

◆ fission_format_spontaneous

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.

◆ fission_frag_beta_delayed

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.

◆ fission_frag_missing

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.

◆ fission_frag_n_induced

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.

◆ fission_frag_spontaneous

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.

◆ fission_rates_beta_delayed

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.

◆ fission_rates_n_induced

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.

◆ fission_rates_spontaneous

character(max_fname_len) parameter_class::fission_rates_spontaneous

reaction library for spontaneous fission rates

Definition at line 159 of file parameter_class.f90.

◆ flow_every

integer parameter_class::flow_every

flow output frequency

Definition at line 64 of file parameter_class.f90.

◆ freeze_rate_temp

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.

◆ gear_cfactor

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.

◆ gear_eps

real(r_kind) parameter_class::gear_eps

Abundance accuracy for gear solver.

Definition at line 136 of file parameter_class.f90.

◆ gear_escale

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.

◆ gear_ignore_adapt_stepsize

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.

◆ gear_nr_eps

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.

◆ gear_nr_maxcount

integer parameter_class::gear_nr_maxcount

Maximum newton-raphson iterations for gear solver.

Definition at line 141 of file parameter_class.f90.

◆ gear_nr_mincount

integer parameter_class::gear_nr_mincount

Minimum newton-raphson iterations for gear solver.

Definition at line 142 of file parameter_class.f90.

◆ gear_timestep_max

real(r_kind) parameter_class::gear_timestep_max

For gear solver.

Definition at line 140 of file parameter_class.f90.

◆ h_custom_snapshots

logical parameter_class::h_custom_snapshots

Same, but in hdf5 format.

Definition at line 108 of file parameter_class.f90.

◆ h_engen_detailed

logical parameter_class::h_engen_detailed

Output the energy per parent nucleus and reaction type.

Definition at line 122 of file parameter_class.f90.

◆ h_engen_every

integer parameter_class::h_engen_every

Energy generation output frequency in hdf5 format.

Definition at line 69 of file parameter_class.f90.

◆ h_finab

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.

◆ h_flow_every

integer parameter_class::h_flow_every

flow output frequency in hdf5 format

Definition at line 65 of file parameter_class.f90.

◆ h_mainout_every

integer parameter_class::h_mainout_every

HDF5 output frequency of the mainout.

Definition at line 72 of file parameter_class.f90.

◆ h_nu_loss_every

integer parameter_class::h_nu_loss_every

Output neutrino loss and gain in hdf5 format.

Definition at line 101 of file parameter_class.f90.

◆ h_snapshot_every

integer parameter_class::h_snapshot_every

snapshot output frequency in hdf5 format

Definition at line 63 of file parameter_class.f90.

◆ h_timescales_every

integer parameter_class::h_timescales_every

timescales output frequency in hdf5 format

Definition at line 67 of file parameter_class.f90.

◆ h_track_nuclei_every

integer parameter_class::h_track_nuclei_every

frequency of track nuclei output in hdf5 format

Definition at line 74 of file parameter_class.f90.

◆ heating_density

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.

◆ heating_frac

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.

◆ heating_mode

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.

◆ heating_t9_tol

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.

◆ htpf_file

character(max_fname_len) parameter_class::htpf_file

high-temperature partition functions (htpf.dat)

Definition at line 157 of file parameter_class.f90.

◆ initemp_cold

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.

◆ initemp_hot

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.

◆ initial_stepsize

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.

◆ interp_mode

integer parameter_class::interp_mode

Mode for interpolation of temperature and density.

Definition at line 151 of file parameter_class.f90.

◆ isotopes_file

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.

◆ iwformat

integer parameter_class::iwformat

defines format of the weak rates (0 = tabulated, 1 = log<ft>)

Definition at line 81 of file parameter_class.f90.

◆ iwinterp

integer parameter_class::iwinterp

defines the interpolation for the weak rates (0 = bilinear, 1 = bicubic)

Definition at line 82 of file parameter_class.f90.

◆ le

character(max_fname_len) parameter_class::le

electron-neutrino luminosities [erg/s]

Definition at line 190 of file parameter_class.f90.

◆ lebar

character(max_fname_len) parameter_class::lebar

electron-antineutrino luminosities [erg/s]

Definition at line 191 of file parameter_class.f90.

◆ lx

character(max_fname_len) parameter_class::lx

Muon and Tauon neutrino luminosities [erg/s].

Definition at line 194 of file parameter_class.f90.

◆ lxbar

character(max_fname_len) parameter_class::lxbar

Muon and Tauon antineutrino luminosities [erg/s].

Definition at line 195 of file parameter_class.f90.

◆ mainout_every

integer parameter_class::mainout_every

frequency of mainout output

Definition at line 71 of file parameter_class.f90.

◆ matcc

real(r_kind), dimension(:), allocatable parameter_class::matcc

matrices for Clenshaw-Curtis integration

Definition at line 216 of file parameter_class.f90.

◆ max_fname_len

integer, parameter, public parameter_class::max_fname_len = 400

maximum length of filenames

Definition at line 58 of file parameter_class.f90.

◆ mcc

real(r_kind), dimension(:,:), allocatable parameter_class::mcc

Definition at line 217 of file parameter_class.f90.

◆ ncc

integer parameter_class::ncc = 256

nr of points for Clenshaw-Curtis integration

Definition at line 215 of file parameter_class.f90.

◆ net_source

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.

◆ neutrino_loss_file

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.

◆ neutrino_mode

character(max_fname_len) parameter_class::neutrino_mode

Similar to trajectory mode.

Definition at line 185 of file parameter_class.f90.

◆ nfission_file

character(max_fname_len) parameter_class::nfission_file

Fission table for neutron-induced fission.

Definition at line 179 of file parameter_class.f90.

◆ nr_maxcount

integer parameter_class::nr_maxcount

– Newton-Raphson iterative loop parameters

no more that this many iterations in NR

Definition at line 201 of file parameter_class.f90.

◆ nr_mincount

integer parameter_class::nr_mincount

Minimum iterations in NR.

Definition at line 202 of file parameter_class.f90.

◆ nr_tol

real(r_kind) parameter_class::nr_tol

exit NR if tolerance less than this value

Definition at line 203 of file parameter_class.f90.

◆ nrdiag_every

integer parameter_class::nrdiag_every

frequency of NR loop diagnostic messages (wrt iteration counter)

Definition at line 70 of file parameter_class.f90.

◆ nse_calc_every

integer parameter_class::nse_calc_every

Compute NSE abundances every x step.

Definition at line 89 of file parameter_class.f90.

◆ nse_delt_t9min

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.

◆ nse_descend_t9start

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.

◆ nse_max_it

integer parameter_class::nse_max_it

Maximum amount of NSE iterations.

Definition at line 117 of file parameter_class.f90.

◆ nse_nr_tol

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.

◆ nse_solver

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.

◆ nsep_energies_file

character(max_fname_len) parameter_class::nsep_energies_file

neutron separation energies

Definition at line 173 of file parameter_class.f90.

◆ nsetemp_cold

real(r_kind) parameter_class::nsetemp_cold

T [GK] for the nse->network switch.

Definition at line 119 of file parameter_class.f90.

◆ nsetemp_hot

real(r_kind) parameter_class::nsetemp_hot

T [GK] for the nse<-network switch.

Definition at line 120 of file parameter_class.f90.

◆ nu_loss_every

integer parameter_class::nu_loss_every

Output neutrino loss and gain.

Definition at line 100 of file parameter_class.f90.

◆ nuchannel_file

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.

◆ nuflag

integer parameter_class::nuflag

defines type of neutrino reactions used

Definition at line 85 of file parameter_class.f90.

◆ nunucleo_rates_file

character(max_fname_len) parameter_class::nunucleo_rates_file

neutrino reaction rates on nucleons

Definition at line 174 of file parameter_class.f90.

◆ nurates_file

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.

◆ out_every

integer parameter_class::out_every

– runtime parameters set when reading parameter file

main output frequency

Definition at line 61 of file parameter_class.f90.

◆ par_binary_name

character(len=*), parameter, private parameter_class::par_binary_name ='parameter.windat'

private

other static variables

Definition at line 207 of file parameter_class.f90.

◆ param_name_len

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.

◆ prepared_network_path

character(max_fname_len) parameter_class::prepared_network_path

Prepared network folder.

Definition at line 198 of file parameter_class.f90.

◆ reaclib_file

character(max_fname_len) parameter_class::reaclib_file

reaction rate library (reaclib)

Definition at line 158 of file parameter_class.f90.

◆ read_initial_composition

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.

◆ rho_analytic

character(max_fname_len) parameter_class::rho_analytic

analytic density [g/cm3]

Definition at line 187 of file parameter_class.f90.

◆ rkm_analytic

character(max_fname_len) parameter_class::rkm_analytic

analytic radial scale [km]

Definition at line 188 of file parameter_class.f90.

◆ screening_mode

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.

◆ seed_file

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.

◆ seed_format

character(max_fname_len) parameter_class::seed_format

Seed format.

Definition at line 154 of file parameter_class.f90.

◆ sfission_file

character(max_fname_len) parameter_class::sfission_file

Fission table for spontaneous fission.

Definition at line 180 of file parameter_class.f90.

◆ snapshot_every

integer parameter_class::snapshot_every

snapshot output frequency

Definition at line 62 of file parameter_class.f90.

◆ snapshot_file

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.

◆ solver

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.

◆ t9_analytic

character(max_fname_len) parameter_class::t9_analytic

analytic temperature [T9]

Definition at line 186 of file parameter_class.f90.

◆ t_analytic

real(r_kind) parameter_class::t_analytic

for parameteric trajectories: initial time

Definition at line 126 of file parameter_class.f90.

◆ tabulated_rates_file

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.

◆ tabulated_temperature_file

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.

◆ temp_reload_exp_weak_rates

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.

◆ termination_criterion

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.

◆ timescales_every

integer parameter_class::timescales_every

timescales output frequency

Definition at line 66 of file parameter_class.f90.

◆ timestep_factor

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.

◆ timestep_hydro_factor

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.

◆ timestep_max

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.

◆ timestep_traj_limit

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.

◆ timestep_ymin

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.

◆ top_engen_every

integer parameter_class::top_engen_every

frequency of energy generators toplist

Definition at line 75 of file parameter_class.f90.

◆ track_nuclei_every

integer parameter_class::track_nuclei_every

frequency of track nuclei output

Definition at line 73 of file parameter_class.f90.

◆ track_nuclei_file

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.

◆ trajectory_file

character(max_fname_len) parameter_class::trajectory_file

name of trajectory data file

Definition at line 152 of file parameter_class.f90.

◆ trajectory_format

character(max_fname_len) parameter_class::trajectory_format

Format to read the trajectory.

Definition at line 184 of file parameter_class.f90.

◆ trajectory_mode

character*20 parameter_class::trajectory_mode

determines how trajectory is calculated (from_file|analytic|expand)

Definition at line 90 of file parameter_class.f90.

◆ unit

type(unit_type), public parameter_class::unit

constants and unit conversion factors

Definition at line 54 of file parameter_class.f90.

◆ use_alpha_decay_file

logical parameter_class::use_alpha_decay_file

Switch for using additional alpha decay rates.

Definition at line 76 of file parameter_class.f90.

◆ use_beta_decay_file

logical parameter_class::use_beta_decay_file

switch for using different format for beta decays

Definition at line 93 of file parameter_class.f90.

◆ use_detailed_balance

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.

◆ use_detailed_balance_q_reac

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.

◆ use_htpf

logical parameter_class::use_htpf

Use high temperature partition functions or not.

Definition at line 145 of file parameter_class.f90.

◆ use_neutrino_loss_file

logical parameter_class::use_neutrino_loss_file

Use a file with Qnu values?

Definition at line 105 of file parameter_class.f90.

◆ use_prepared_network

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.

◆ use_tabulated_rates

logical parameter_class::use_tabulated_rates

switch for using tabulated rates (e.g. talysNGrates.dat)

Definition at line 92 of file parameter_class.f90.

◆ use_thermal_nu_loss

logical parameter_class::use_thermal_nu_loss

Whether to include thermal neutrino loss or not.

Definition at line 99 of file parameter_class.f90.

◆ use_timmes_mue

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.

◆ weak_rates_file

character(max_fname_len) parameter_class::weak_rates_file

weak rates library (twr.dat)

Definition at line 169 of file parameter_class.f90.

◆ weights

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 211 of file parameter_class.f90.

◆ weightscalculated

logical parameter_class::weightscalculated = .false.

switch to calculated weights and nodes for [1,infinity]

Definition at line 213 of file parameter_class.f90.

◆ xnodes

real(r_kind), dimension(:), allocatable parameter_class::xnodes

corresponding nodes "

Definition at line 212 of file parameter_class.f90.

◆ ye_analytic

character(max_fname_len) parameter_class::ye_analytic

analytic electron fraction

Definition at line 189 of file parameter_class.f90.

Data Types

Functions/Subroutines

Variables

Detailed Description

Function/Subroutine Documentation

◆ check_param()

◆ levenshteindistance()

◆ output_binary_parameter_data()

◆ output_param()

◆ output_param_prepared_network()

◆ read_float_param()

◆ read_integer_param()

◆ read_param()

◆ read_sanity_check_prepared_network()

◆ search_close_parameter()

◆ set_default_param()

◆ set_param()

◆ unit_define()

Variable Documentation

◆ adapt_stepsize_maxcount

◆ alpha_decay_file

◆ alpha_decay_ignore_all

◆ alpha_decay_src_ignore

◆ alpha_decay_zmax

◆ alpha_decay_zmin

◆ beta_decay_file

◆ beta_decay_src_ignore

◆ bfission_file

◆ calc_nsep_energy

◆ chem_pot_file

◆ custom_snapshots

◆ dcc

◆ detailed_balance_src_ignore

◆ detailed_balance_src_q_reac

◆ detailed_balance_src_q_winvn

◆ engen

◆ engen_alpha

◆ engen_beta

◆ engen_every

◆ engen_fiss

◆ enue

◆ enuebar

◆ enux

◆ enuxbar

◆ expansiontype

◆ extrapolation_width

◆ final_dens

◆ final_temp

◆ final_time

◆ fissflag

◆ fission_format_beta_delayed

◆ fission_format_n_induced

◆ fission_format_spontaneous

◆ fission_frag_beta_delayed

◆ fission_frag_missing

◆ fission_frag_n_induced

◆ fission_frag_spontaneous

◆ fission_rates_beta_delayed

◆ fission_rates_n_induced

◆ fission_rates_spontaneous

◆ flow_every

◆ freeze_rate_temp

◆ gear_cfactor

◆ gear_eps

◆ gear_escale

◆ gear_ignore_adapt_stepsize

◆ gear_nr_eps

◆ gear_nr_maxcount

◆ gear_nr_mincount

◆ gear_timestep_max

◆ h_custom_snapshots

◆ h_engen_detailed

◆ h_engen_every

◆ h_finab

◆ h_flow_every

◆ h_mainout_every

◆ h_nu_loss_every

◆ h_snapshot_every

◆ h_timescales_every

◆ h_track_nuclei_every