Add PlasmaDensityLimit class and integrate into physics and stellarat…

documentation/source/physics-models/plasma_density.md

@@ -1,7 +1,7 @@
-# Density Limit
+# Density Limit | `PlasmaDensityLimit`
 
 Several density limit models are available in PROCESS. These are
-calculated in routine `calculate_density_limit()`, which is called by `physics`.
+calculated in the run routine `density_limit.run()`, which is called by `physics()`.
 
 This constraint can be activated by stating `icc = 5` in the input file.
 
@@ -13,7 +13,7 @@ For the models below $P_{\perp}$ is the mean heat flux density across the separa
 
 -----------------
 
-## ASDEX model
+## ASDEX model | `calculate_asdex_density_limit()`
 
 Switch value: `i_density_limit = 1`[^1][^2]
 
@@ -23,7 +23,7 @@ $$
 
 -----------------
 
-## Borrass model for ITER, I
+## Borrass model for ITER, I | `calculate_borrass_iter_i_density_limit()`
 
 Switch value: `i_density_limit = 2` [^1]
 
@@ -35,7 +35,7 @@ $C \approx$  1.8 for ITER-like conditions.
 
 -----------------
 
-## Borrass model for ITER, II 
+## Borrass model for ITER, II | `calculate_borrass_iter_ii_density_limit()`
 
 Switch value: `i_density_limit = 3` [^1]
 
@@ -45,7 +45,7 @@ $$
 
 -----------------
 
-## JET edge radiation model
+## JET edge radiation model | `calculate_jet_edge_radiation_density_limit()`
 
 Switch value: `i_density_limit = 4` [^1]
 
@@ -55,7 +55,7 @@ $$
 
 -----------------
 
-## JET simplified model
+## JET simplified model | `calculate_jet_simple_density_limit()`
 
 Switch value: `i_density_limit = 5` [^1]
 
@@ -73,7 +73,7 @@ where $\kappa \approx 1.5, \Delta \approx 0.1a$ has been taken from JET.
 
 -----------------
 
-## Hugill-Murakami model
+## Hugill-Murakami model | `calculate_hugill_murakami_density_limit()`
 
 Switch value: `i_density_limit = 6` [^2]
 
@@ -84,7 +84,7 @@ $$
 
 -----------------
 
-## Greenwald model
+## Greenwald model | `calculate_greenwald_density_limit()`
 
 Switch value: `i_density_limit = 7` [^3][^4]
 
@@ -96,7 +96,7 @@ For the Greenwald model the limit applies to the line-averaged electron density,
 
 ---------------------
 
-## ASDEX New model
+## ASDEX New model | `calculate_asdex_new_density_limit()`
 
 Switch value: `i_density_limit = 8` [^5][^6]
 

process/main.py

@@ -90,6 +90,7 @@
     LowerHybrid,
     NeutralBeam,
 )
+from process.models.physics.density_limit import PlasmaDensityLimit
 from process.models.physics.impurity_radiation import initialise_imprad
 from process.models.physics.physics import (
     DetailedPhysics,
@@ -698,11 +699,13 @@ def __init__(self):
         )
         self.plasma_beta = PlasmaBeta()
         self.plasma_inductance = PlasmaInductance()
+        self.plasma_density_limit = PlasmaDensityLimit()
         self.physics = Physics(
             plasma_profile=self.plasma_profile,
             current_drive=self.current_drive,
             plasma_beta=self.plasma_beta,
             plasma_inductance=self.plasma_inductance,
+            plasma_density_limit=self.plasma_density_limit,
         )
         self.physics_detailed = DetailedPhysics(
             plasma_profile=self.plasma_profile,