Prints out the default settings as a namelist
subroutine print_default_namelist
!+ Prints out the default settings as a namelist
write(*,'(a)') "!Default Atomic Table Settings"
write(*,'(a)') "&general_settings"
write(*,'(a)') "n_max = 12, !Number of initial atomic energy levels"
write(*,'(a)') "m_max = 12, !Number of final atomic energy levels"
write(*,'(a)') "tables_file = './atomic_tables.h5'"
write(*,'(a)') "/"
write(*,'(a)') "!Hydrogen-Hydrogen Cross Sections"
write(*,'(a)') "&H_H_cross"
write(*,'(a)') "calculate = T, !Calculate Table"
write(*,'(a)') "nenergy = 200, !Number of energy values"
write(*,'(a)') "emin = 1.0E-3, !Minimum energy/amu [keV/amu]"
write(*,'(a)') "emax = 8.0E2 !Maximum energy/amu [keV/amu]"
write(*,'(a)') "/"
write(*,'(a)') "!Hydrogen-Electron Cross Sections"
write(*,'(a)') "&H_e_cross"
write(*,'(a)') "calculate = T, !Calculate Table"
write(*,'(a)') "nenergy = 200, !Number of energy values"
write(*,'(a)') "emin = 1.0E-3, !Minimum energy/amu [keV/amu]"
write(*,'(a)') "emax = 8.0E2 !Maximum energy/amu [keV/amu]"
write(*,'(a)') "/"
write(*,'(a)') "!Hydrogen-Impurity Cross Sections. Up to 10 impurity charges"
write(*,'(a)') "&H_Aq_cross"
write(*,'(a)') "calculate = T, !Calculate Table"
write(*,'(a)') "q(1) = 5, !Impurity charge: Boron: 5, Carbon: 6, ..."
write(*,'(a)') "q(2) = 6, !Impurity charge: Boron: 5, Carbon: 6, ..."
write(*,'(a)') "nenergy = 200, !Number of energy values"
write(*,'(a)') "emin = 1.0E-3, !Minimum energy/amu [keV/amu]"
write(*,'(a)') "emax = 8.0E2 !Maximum energy/amu [keV/amu]"
write(*,'(a)') "/"
write(*,'(a)') "!Deuterium-Deuterium Nuclear Cross Sections"
write(*,'(a)') "&D_D_cross"
write(*,'(a)') "calculate = T, !Calculate Table"
write(*,'(a)') "nenergy = 200, !Number of energy values"
write(*,'(a)') "emin = 1.0E-3, !Minimum energy [keV]"
write(*,'(a)') "emax = 8.0E2 !Maximum energy [keV]"
write(*,'(a)') "/"
write(*,'(a)') "!Deuterium-Helium3 Nuclear Cross Sections"
write(*,'(a)') "&D_He3_cross"
write(*,'(a)') "calculate = T, !Calculate Table"
write(*,'(a)') "nenergy = 100, !Number of energy values"
write(*,'(a)') "emin = 1.0E-3, !Minimum energy [keV]"
write(*,'(a)') "emax = 4.0E2, !Maximum energy [keV]"
write(*,'(a)') "/"
write(*,'(a)') "!Hydrogen-Hydrogen Reaction Rates"
write(*,'(a)') "&H_H_rates"
write(*,'(a)') "calculate = T, !Calculate Table"
write(*,'(a)') "nenergy = 100, !Number of energy/amu values"
write(*,'(a)') "emin = 1.0E-3, !Minimum energy/amu [keV/amu]"
write(*,'(a)') "emax = 4.0E2, !Maximum energy/amu [keV/amu]"
write(*,'(a)') "ntemp = 100, !Number of temperature/amu values"
write(*,'(a)') "tmin = 1.0E-3, !Minimum ion temperature/amu [keV/amu]"
write(*,'(a)') "tmax = 2.0E1 !Maximum ion temperature/amu [keV/amu]"
write(*,'(a)') "/"
write(*,'(a)') "!Hydrogen-Electron Reaction Rates"
write(*,'(a)') "&H_e_rates"
write(*,'(a)') "calculate = T, !Calculate Table"
write(*,'(a)') "nenergy = 100, !Number of energy/amu values"
write(*,'(a)') "emin = 1.0E-3, !Minimum energy/amu [keV/amu]"
write(*,'(a)') "emax = 4.0E2, !Maximum energy/amu [keV/amu]"
write(*,'(a)') "ntemp = 100, !Number of temperature values"
write(*,'(a)') "tmin = 1.0E-3, !Minimum electron temperature [keV]"
write(*,'(a)') "tmax = 2.0E1 !Maximum electron temperature [keV]"
write(*,'(a)') "/"
write(*,'(a)') "!Hydrogen-Impurity Reaction Rates. Up to 10 impurity charges"
write(*,'(a)') "&H_Aq_rates"
write(*,'(a)') "calculate = T, !Calculate Table"
write(*,'(a)') "q(1) = 5, !Impurity charge: Boron: 5, Carbon: 6, ..."
write(*,'(a)') "q(2) = 6, !Impurity charge: Boron: 5, Carbon: 6, ..."
write(*,'(a)') "mass = 12.011, !Impurity mass [amu]"
write(*,'(a)') "nenergy = 100, !Number of energy/amu values"
write(*,'(a)') "emin = 1.0E-3, !Minimum energy/amu [keV/amu]"
write(*,'(a)') "emax = 4.0E2, !Maximum energy/amu [keV/amu]"
write(*,'(a)') "ntemp = 100, !Number of temperature values"
write(*,'(a)') "tmin = 1.0E-3, !Minimum ion temperature [keV]"
write(*,'(a)') "tmax = 2.0E1 !Maximum ion temperature [keV]"
write(*,'(a)') "/"
write(*,'(a)') "!Deuterium-Deuterium Nuclear Reaction Rates"
write(*,'(a)') "&D_D_rates"
write(*,'(a)') "calculate = T, !Calculate Table"
write(*,'(a)') "nenergy = 100, !Number of energy values"
write(*,'(a)') "emin = 1.0E-3, !Minimum energy [keV]"
write(*,'(a)') "emax = 4.0E2, !Maximum energy [keV]"
write(*,'(a)') "ntemp = 100, !Number of temperature values"
write(*,'(a)') "tmin = 1.0E-3, !Minimum deuterium temperature [keV]"
write(*,'(a)') "tmax = 2.0E1 !Maximum deuterium temperature [keV]"
write(*,'(a)') "/"
write(*,'(a)') "!Deuterium-Helium3 Nuclear Reaction Rates"
write(*,'(a)') "&D_He3_rates"
write(*,'(a)') "calculate = T, !Calculate Table"
write(*,'(a)') "nenergy = 100, !Number of energy values"
write(*,'(a)') "emin = 1.0E-3, !Minimum energy [keV]"
write(*,'(a)') "emax = 4.0E2, !Maximum energy [keV]"
write(*,'(a)') "ntemp = 100, !Number of temperature values"
write(*,'(a)') "tmin = 1.0E-3, !Minimum deuterium temperature [keV]"
write(*,'(a)') "tmax = 2.0E1 !Maximum deuterium temperature [keV]"
write(*,'(a)') "/"
end subroutine print_default_namelist