Define variables associated to the turbulence models.

Variables
logical	is_rans_fully_explicit = .false.

Variables associated to the RANS models
integer	rans_model = rans_model_k_omega
	RANS model choice.

logical	rans_is_k_omega_sst = .false.
	RANS SST option of k omega (false by default)

double precision	rans_time_step = 0d0
	Current time step.

double precision	rans_time_step_n = 0d0
	Previous time step.

double precision	rans_time_step_nm1 = 0d0
	Ante previous time step.

integer	rans_time_order_discretization = time_order_1
	Time order discretization.

integer	rans_explicit_time_order_discretization = time_order_nssp32
	Explicit time order discretization.

logical	rans_use_explicit_terms_accumulator = .false.
	A switch for when to use explicit accumulator.

Boundary condition RANS variables
type(t_boundary_condition)	turbulent_kinetic_energy_boundary_condition
	Boundary conditions.

type(t_boundary_condition)	specific_dissipation_rate_boundary_condition

type(t_boundary_condition)	turbulence_dissipation_boundary_condition

type(t_boundary_condition)	mean_squared_velocity_boundary_condition

type(t_boundary_condition)	f_22_boundary_condition

Initializer
class(t_scalar_initializer), pointer	turbulent_kinetic_energy_initializer => null()
	Initial condition initializer.

class(t_scalar_initializer), pointer	specific_dissipation_rate_initializer => null()

class(t_scalar_initializer), pointer	turbulence_dissipation_initializer => null()

class(t_scalar_initializer), pointer	mean_squared_velocity_initializer => null()

class(t_scalar_initializer), pointer	f_22_initializer => null()

class(t_cell_bc_initializer), pointer	turbulent_kinetic_energy_bc_initializer => null()
	Boundary conditions initializer.

class(t_cell_bc_initializer), pointer	specific_dissipation_rate_bc_initializer => null()

class(t_cell_bc_initializer), pointer	turbulence_dissipation_bc_initializer => null()

class(t_cell_bc_initializer), pointer	mean_squared_velocity_bc_initializer => null()

class(t_cell_bc_initializer), pointer	f_22_bc_initializer => null()

class(t_rans_ibc_initializer_item), dimension(:), allocatable	turbulent_kinetic_energy_ibc_initializer
	Immersed boundary conditions initializers.

class(t_rans_ibc_initializer_item), dimension(:), allocatable	specific_dissipation_rate_ibc_initializer

class(t_rans_ibc_initializer_item), dimension(:), allocatable	turbulence_dissipation_ibc_initializer

class(t_rans_ibc_initializer_item), dimension(:), allocatable	mean_squared_velocity_ibc_initializer

class(t_rans_ibc_initializer_item), dimension(:), allocatable	f_22_ibc_initializer

class(t_scalar_initializer), pointer	turbulent_kinetic_energy_source_term_initializer => null()
	Source term initializer.

class(t_scalar_initializer), pointer	specific_dissipation_rate_source_term_initializer => null()

class(t_scalar_initializer), pointer	turbulence_dissipation_source_term_initializer => null()

class(t_scalar_initializer), pointer	mean_squared_velocity_source_term_initializer => null()

class(t_scalar_initializer), pointer	f_22_source_term_initializer => null()

class(t_scalar_initializer), pointer	turbulent_kinetic_energy_linear_term_initializer => null()
	Linear term initializer.

class(t_scalar_initializer), pointer	specific_dissipation_rate_linear_term_initializer => null()

class(t_scalar_initializer), pointer	turbulence_dissipation_linear_term_initializer => null()

class(t_scalar_initializer), pointer	mean_squared_velocity_linear_term_initializer => null()

class(t_scalar_initializer), pointer	f_22_linear_term_initializer => null()

Logical list to activate terms of the RANS equation
logical	turbulent_kinetic_energy_has_source_term = .false.
	Has source term?

logical	specific_dissipation_rate_has_source_term = .false.

logical	turbulence_dissipation_has_source_term = .true.

logical	mean_squared_velocity_has_source_term = .false.

logical	f_22_has_source_term = .false.

logical	rans_has_advection_term = .true.
	Has advection term? Unique variable for (k,omega)

logical	k_omega_sst_has_advection_term = .true.

logical	v2_f_has_advection_term = .true.

logical	rans_has_diffusion_term = .true.
	Has diffusion term? Unique variable for (k,omega)

logical	k_omega_sst_has_diffusion_term = .true.

logical	v2_f_has_diffusion_term = .true.

logical	k_omega_sst_has_div_u_advection_term = .false.
	Has \( div(\mathbf{u}) \) advection term? Unique variable for (k,omega)

logical	v2_f_has_div_u_advection_term = .true.

logical	k_omega_sst_has_temporal_term = .true.
	Has temporal term? Unique variable for (k,omega)

logical	v2_f_has_temporal_term = .true.

logical	turbulent_kinetic_energy_has_linear_term = .false.
	Has linear term?

logical	specific_dissipation_rate_has_linear_term = .false.

logical	turbulence_dissipation_has_linear_term = .true.

logical	mean_squared_velocity_has_linear_term = .false.

logical	f_22_has_linear_term = .false.

Variables related to immersed boundaries
logical	rans_ib_has_one_sided_inner_discretization
	One sided discretization close to immersed boundaries.

integer, dimension(:), allocatable	rans_ib_inner_discretization_order

Variables associated to the discretization of the advection term
integer	rans_advection_term_discretization_type = implicit_discretization
	Advection term discretization type (implicit or explicit)

integer	rans_advection_term_scheme = cell_centered_o2_advection_scheme
	Advection term scheme.

integer	rans_advection_term_splitting_scheme = lie_trotter_splitting
	Advection term: Splitting Method for Lax-Wendroff TVD Superbee.

double precision, dimension(3)	rans_splitting_time_coeff =1d0

logical	rans_advection_term_dir_split = .false.
	Advection term: directional splitting option.

logical	is_rans_specify_advection_scheme_stability_factor = .false.
	Specify a CFL factor for explicit sub iterations.

double precision	rans_advection_scheme_stability_factor = -1d0
	CFL used for explicit sub iterations (if `is_rans_specify_advection_scheme_stability_factor`, otherwize, using default values)

type(t_fv_flux), parameter	default_rans_advection_explicit_flux = default_flux_godunov
	Fluxes type The default explicit cell advection flux.

type(t_fv_flux)	rans_advection_explicit_flux = default_rans_advection_explicit_flux
	Fluxes type The effective explicit cell advection flux.


integer	rans_diffusion_term_discretization_type = implicit_discretization
	Diffusion term discretization type (implicit or explicit)

integer	rans_diffusion_term_scheme = cell_centered_o2_diffusion_scheme
	Diffusion term scheme.

integer	rans_diffusion_explicit_time_order_discretization

Miscellaneous variables
logical	rans_has_ghost_boundary_cells = .false.
	Ghost boudary cell choice.

type(t_cell_stencil)	rans_stencil
	Stencil of the energy equation.

integer	rans_enforced_stencil = -1
	Enforced stencil.

Linear system solver variables
type(t_ls_map)	rans_ls_map
	Mapping variables.

class(t_linear_solver), allocatable	rans_solver
	Linear system solver.


double precision, parameter	min_turbulent_kinetic_energy = 1d-12
	Numerical thresholds.

double precision, parameter	min_specific_dissipation_rate = 1d-12

k-omega SST constants
double precision	sst_sigma_k1 = 0.85d0

double precision	sst_sigma_omega1 = 0.5d0

double precision	sst_beta1 = 0.075d0

double precision	sst_gamma1 = 0.553166667d0

double precision	sst_sigma_k2 = 1d0

double precision	sst_sigma_omega2 = 0.856d0

double precision	sst_beta2 = 0.0828d0

double precision	sst_gamma2 = 0.440354667d0

double precision	sst_beta_star = 0.09d0

double precision	sst_kappa = 0.41d0

double precision	sst_a1 = 0.31d0

k-omega constants
double precision	beta_star = 0.09d0

double precision	clim = 0.875d0

double precision	sigma_star = 0.6d0

double precision	sigma = 0.5d0

double precision	alpha = 0.52d0

double precision	beta_0 = 0.0708d0

double precision	sigma_d0 = 0.125d0

v2_f model constants
double precision	v2_f_sigma_k = 1d0

double precision	v2_f_sigma_e = 1.3d0

double precision	v2_f_c_mu = 0.19d0

double precision	v2_f_c_e_2 = 1.9d0

double precision	v2_f_c_1 = 1.4d0

double precision	v2_f_c_2 = 0.3d0

double precision	v2_f_c_l = 0.3d0

double precision	v2_f_c_eta = 70d0

Variables