Detailed plasma part 3

process/data_structure/physics_variables.py

@@ -1308,6 +1308,12 @@
 len_plasma_debye_electron_profile: list[float] = None
 """Profile of electron Debye length in plasma (m)"""
 
+radius_plasma_deuteron_toroidal_larmor_isotropic_profile: list[float] = None
+"""Profile of deuteron toroidal Larmor radius in plasma, assuming equal speeds in all directions (m)"""
+
+radius_plasma_triton_toroidal_larmor_isotropic_profile: list[float] = None
+"""Profile of triton toroidal Larmor radius in plasma, assuming equal speeds in all directions (m)"""
+
 len_plasma_debye_electron_vol_avg: float = None
 """Volume averaged electron Debye length in plasma (m)"""
 
@@ -1320,6 +1326,12 @@
 vel_plasma_triton_profile: list[float] = None
 """Profile of triton thermal velocity in plasma (m/s)"""
 
+vel_plasma_alpha_thermal_profile: list[float] = None
+"""Profile of thermal alpha particle velocity in plasma (m/s)"""
+
+vel_plasma_alpha_birth: float = None
+"""Birth velocity of alpha particles in plasma (m/s)"""
+
 plasma_coulomb_log_electron_electron_profile: list[float] = None
 """Profile of electron-electron Coulomb logarithm in plasma"""
 
@@ -1329,9 +1341,18 @@
 plasma_coulomb_log_electron_triton_profile: list[float] = None
 """Profile of electron-triton Coulomb logarithm in plasma"""
 
+plasma_coulomb_log_deuteron_triton_profile: list[float] = None
+"""Profile of deuteron-triton Coulomb logarithm in plasma"""
+
+plasma_coulomb_log_electron_alpha_thermal_profile: list[float] = None
+"""Profile of electron-alpha Coulomb logarithm in plasma"""
+
 freq_plasma_electron_profile: list[float] = None
 """Electron plasma frequency profile (Hz)"""
 
+freq_plasma_deuteron_profile: list[float] = None
+"""Deuteron plasma frequency profile (Hz)"""
+
 freq_plasma_larmor_toroidal_electron_profile: list[float] = None
 """Profile of electron Larmor frequency in plasma due to toroidal magnetic field (Hz)"""
 
@@ -1341,6 +1362,9 @@
 freq_plasma_larmor_toroidal_triton_profile: list[float] = None
 """Profile of triton Larmor frequency in plasma due to toroidal magnetic field (Hz)"""
 
+freq_plasma_upper_hybrid_profile: list[float] = None
+"""Profile of upper hybrid frequency in plasma (Hz)"""
+
 
 def init_physics_module():
     """Initialise the physics module"""
@@ -1659,17 +1683,25 @@ def init_physics_variables():
         n_charge_plasma_effective_mass_weighted_vol_avg, \
         j_plasma_bootstrap_sauter_profile, \
         len_plasma_debye_electron_profile, \
+        radius_plasma_deuteron_toroidal_larmor_isotropic_profile, \
+        radius_plasma_triton_toroidal_larmor_isotropic_profile, \
         len_plasma_debye_electron_vol_avg, \
         vel_plasma_electron_profile, \
         vel_plasma_deuteron_profile, \
         vel_plasma_triton_profile, \
+        vel_plasma_alpha_thermal_profile, \
+        vel_plasma_alpha_birth, \
         plasma_coulomb_log_electron_electron_profile, \
         plasma_coulomb_log_electron_deuteron_profile, \
         plasma_coulomb_log_electron_triton_profile, \
+        plasma_coulomb_log_deuteron_triton_profile, \
+        plasma_coulomb_log_electron_alpha_thermal_profile, \
         freq_plasma_electron_profile, \
+        freq_plasma_deuteron_profile, \
         freq_plasma_larmor_toroidal_electron_profile, \
         freq_plasma_larmor_toroidal_deuteron_profile, \
-        freq_plasma_larmor_toroidal_triton_profile
+        freq_plasma_larmor_toroidal_triton_profile, \
+        freq_plasma_upper_hybrid_profile
 
     m_beam_amu = 0.0
     m_fuel_amu = 0.0
@@ -1928,14 +1960,22 @@ def init_physics_variables():
     n_charge_plasma_effective_profile = []
     n_charge_plasma_effective_mass_weighted_vol_avg = 0.0
     len_plasma_debye_electron_profile = []
+    radius_plasma_deuteron_toroidal_larmor_isotropic_profile = []
+    radius_plasma_triton_toroidal_larmor_isotropic_profile = []
     len_plasma_debye_electron_vol_avg = 0.0
     vel_plasma_electron_profile = []
     vel_plasma_deuteron_profile = []
     vel_plasma_triton_profile = []
+    vel_plasma_alpha_thermal_profile = []
+    vel_plasma_alpha_birth = 0.0
     plasma_coulomb_log_electron_electron_profile = []
     plasma_coulomb_log_electron_deuteron_profile = []
     plasma_coulomb_log_electron_triton_profile = []
+    plasma_coulomb_log_deuteron_triton_profile = []
+    plasma_coulomb_log_electron_alpha_thermal_profile = []
     freq_plasma_electron_profile = []
+    freq_plasma_deuteron_profile = []
     freq_plasma_larmor_toroidal_electron_profile = []
     freq_plasma_larmor_toroidal_deuteron_profile = []
     freq_plasma_larmor_toroidal_triton_profile = []
+    freq_plasma_upper_hybrid_profile = []