use formulae rho fun for initial air density in Pyrcel example. closes …

…#1108 (#1110)
open-atmos · Aug 1, 2023 · dc8a402 · dc8a402
1 parent 8add5ec
commit dc8a402
Showing 1 changed file with 4 additions and 5 deletions.
diff --git a/examples/PySDM_examples/Pyrcel/settings.py b/examples/PySDM_examples/Pyrcel/settings.py
@@ -42,15 +42,14 @@ def __init__(
 
     @property
     def initial_air_density(self):
-        const = self.formulae.constants
-        dry_air_density = (
-            self.formulae.trivia.p_d(
+        return self.formulae.state_variable_triplet.rho_of_rhod_qv(
+            rhod=self.formulae.trivia.p_d(
                 self.initial_pressure, self.initial_vapour_mixing_ratio
             )
             / self.initial_temperature
-            / const.Rd
+            / self.formulae.constants.Rd,
+            qv=self.initial_vapour_mixing_ratio,
         )
-        return dry_air_density * (1 + self.initial_vapour_mixing_ratio)
 
     @property
     def nt(self) -> int: