🎨 Pressure profile plot_proc

process/data_structure/physics_variables.py

@@ -806,10 +806,23 @@
 pres_plasma_on_axis: float = None
 """central total plasma pressure (Pa)"""
 
+pres_plasma_total_profile: list[float] = None
+"""Profile of total pressure in plasma (Pa)"""
+
+pres_plasma_electron_profile: list[float] = None
+"""Profile of electron pressure in plasma (Pa)"""
+
+pres_plasma_ion_total_profile: list[float] = None
+"""Profile of ion pressure in plasma (Pa)"""
+
+pres_plasma_fuel_profile: list[float] = None
+"""Profile of fuel pressure in plasma (Pa)"""
 
 j_plasma_on_axis: float = None
 """Central plasma current density (A/m2)"""
 
+n_plasma_profile_elements: int = None
+"""Number of elements in plasma profile"""
 
 pres_plasma_vol_avg: float = None
 """Volume averaged plasma pressure (Pa)"""
@@ -1476,7 +1489,12 @@ def init_physics_variables():
     global nd_plasma_ions_on_axis
     global m_s_limit
     global pres_plasma_on_axis
+    global pres_plasma_total_profile
+    global pres_plasma_electron_profile
+    global pres_plasma_ion_total_profile
+    global pres_plasma_fuel_profile
     global j_plasma_on_axis
+    global n_plasma_profile_elements
     global f_dd_branching_trit
     global pden_plasma_alpha_mw
     global pden_alpha_total_mw
@@ -1723,7 +1741,12 @@ def init_physics_variables():
     nd_plasma_ions_on_axis = 0.0
     m_s_limit = 0.3
     pres_plasma_on_axis = 0.0
+    pres_plasma_total_profile = []
+    pres_plasma_electron_profile = []
+    pres_plasma_ion_total_profile = []
+    pres_plasma_fuel_profile = []
     j_plasma_on_axis = 0.0
+    n_plasma_profile_elements = 500
     f_dd_branching_trit = 0.0
     pden_plasma_alpha_mw = 0.0
     pden_alpha_total_mw = 0.0