@@ -23,6 +23,7 @@ def __init__(self,
2323 pedistal_radius = 0.8 ,
2424 ion_density_peaking_factor = 1 ,
2525 ion_temperature_peaking_factor = 8.06 ,
26+ ion_temperature_beta = 6.0 ,
2627 shafranov_shift = 0.0 ,
2728 number_of_bins = 100 ,
2829 plasma_type = 1 ,
@@ -44,6 +45,7 @@ def __init__(self,
4445 self .pedistal_radius = pedistal_radius # pedistal major rad
4546 self .ion_density_peaking_factor = ion_density_peaking_factor
4647 self .ion_temperature_peaking_factor = ion_temperature_peaking_factor
48+ self .ion_temperature_beta = ion_temperature_beta
4749 self .shafranov_shift = shafranov_shift
4850 self .number_of_bins = number_of_bins
4951 self .plasma_type = plasma_type # 0 = L mode anything else H/A mode
@@ -109,6 +111,17 @@ def ion_temperature_peaking_factor(self, ion_temperature_peaking_factor):
109111 else :
110112 self ._ion_temperature_peaking_factor = ion_temperature_peaking_factor
111113
114+ @property
115+ def ion_temperature_beta (self ):
116+ return self ._ion_temperature_beta
117+
118+ @ion_temperature_beta .setter
119+ def ion_temperature_beta (self , ion_temperature_beta ):
120+ if ion_temperature_beta < 0 :
121+ raise ValueError ('ion_temperature_beta is out of range' )
122+ else :
123+ self ._ion_temperature_beta = ion_temperature_beta
124+
112125 @property
113126 def ion_density_peaking_factor (self ):
114127 return self ._ion_density_peaking_factor
@@ -329,6 +342,7 @@ def export_plasma_source(self, output_filename):
329342 "pedistal_radius" : self .pedistal_radius ,
330343 "ion_density_peaking_factor" : self .ion_density_peaking_factor ,
331344 "ion_temperature_peaking_factor" : self .ion_temperature_peaking_factor ,
345+ "ion_temperature_beta" : self .ion_temperature_beta ,
332346 "minor_radius" : self .minor_radius ,
333347 "major_radius" : self .major_radius ,
334348 "elongation" : self .elongation ,
0 commit comments