add optional time series outputted by StagYY

Author: amorison

src/stagpy/dimensions.py

@@ -71,6 +71,11 @@ def density(self) -> float:
         """Density in kg/m3."""
         return self.par.nml["refstate"]["dens_dimensional"]
 
+    @property
+    def mass(self) -> float:
+        """Mass in kg."""
+        return self.density * self.length**3
+
     @property
     def th_cond(self) -> float:
         """Thermal conductivity in W/(m.K)."""

src/stagpy/phyvars.py

@@ -345,6 +345,61 @@
         "H_AH": Vart("Adiabatic heating", "Power density", "W/m3"),
         "H_cool": Vart("Cooling", "Heat flux", "W/m2"),
         "H_melt": Vart("Melting heat flux", "Heat flux", "W/m2"),
+        "erupta_total": Vart("Total erupta", "Mass", "kg"),
+        "Psurf": Vart("Ground pressure", "Pressure", "Pa"),
+        "weathering": Vart("Weathering", "Mass flux", "tons/yr"),
+        "mH2O_total": Vart("Total H2O mass", "Mass", "kg"),
+        "mH2O_mantle": Vart("H2O in mantle", "Mass", "kg"),
+        "mH2O_crust": Vart("H2O in crust", "Mass", "kg"),
+        "dmH2O": Vart("Absorbed H2O", "Mass", "kg"),
+        "cH2O_mean": Vart("H2O concentration", "Concentration", "1"),
+        "F_top_H2O_diff": Vart("H2O diffusion flux at top", "Mass flux", "kg/s"),
+        "F_top_H2O_erupted": Vart("H2O eruption flux at top", "Mass flux", "kg/s"),
+        "outgassed_H2O": Vart("H2O outgassed", "Fraction", "1"),
+        "cCO2_mean": Vart("CO2 concentration", "Concentration", "1"),
+        "outgassed_CO2": Vart("CO2 outgassed", "Fraction", "1"),
+        "cN2_mean": Vart("N2 concentration", "Concentration", "1"),
+        "outgassed_N2": Vart("N2 outgassed", "Fraction", "1"),
+        "c40Ar_mean": Vart("40Ar concentration", "Concentration", "1"),
+        "c36Ar_mean": Vart("36Ar concentration", "Concentration", "1"),
+        "outgassed_40Ar": Vart("40Ar outgassed", "Fraction", "1"),
+        "outgassed_36Ar": Vart("36Ar outgassed", "Fraction", "1"),
+        "c4He_mean": Vart("4He concentration", "Concentration", "1"),
+        "c3He_mean": Vart("3He concentration", "Concentration", "1"),
+        "outgassed_4He": Vart("4He outgassed", "Fraction", "1"),
+        "outgassed_3He": Vart("3He outgassed", "Fraction", "1"),
+        "IngassedH2O": Vart("Ingassed H2O", "Fraction", "1"),
+        "OutputtedNotEruptedH2O": Vart("Outputted H2O", "Fraction", "1"),
+        "SaturationOutgasH2O": Vart("Saturation-outgas H2O", "Fraction", "1"),
+        "EruptedH2O": Vart("Erupted H2O", "Fraction", "1"),
+        "SurfOceanMass3D": Vart("Scaled H2O mass", "Mass", "kg"),
+        "TotH2OMassOceanPlusMantle": Vart("Total H2O mass", "Mass", "kg"),
+        "GS_min": Vart("Min grain size", "Size", "1"),
+        "GS_mean": Vart("Grain size", "Size", "1"),
+        "GS_max": Vart("Max grain size", "Size", "1"),
+        "intruda": Vart("Intruda", "Mantle mass", "1"),
+        "erupta_TTG": Vart("TTG erupta", "Mantle mass", "1"),
+        "intruda_TTG": Vart("TTG intruda", "Mantle mass", "1"),
+        "TTG1": Vart("TTG1", "Mantle mass", "1"),
+        "TTG2": Vart("TTG2", "Mantle mass", "1"),
+        "TTG3": Vart("TTG3", "Mantle mass", "1"),
+        "s_core": Vart("Core size", "Length", "m"),
+        "tc_core": Vart("Temperature at center", "Temperature", "K"),
+        "ts_core": Vart("Temperature at core surface", "Temperature", "K"),
+        "m_Fe_S2TMO": Vart("FeO mass flux solid to TMO", "Mass", "kg"),
+        "m_Fe_TMO2S": Vart("FeO mass flux TMO to solid", "Mass", "kg"),
+        "m_Fe_S2BMO": Vart("FeO mass flux solid to BMO", "Mass", "kg"),
+        "m_Fe_BMO2S": Vart("FeO mass flux BMO to solid", "Mass", "kg"),
+        "c_Fe_TMO": Vart("FeO content TMO", "Concentration", "1"),
+        "c_Fe_BMO": Vart("FeO content BMO", "Concentration", "1"),
+        "c_Fe_sol": Vart("FeO content solid", "Concentration", "1"),
+        "denstramin": Vart("Min tracers density", "Density", "kg/m3"),
+        "denstramean": Vart("Tracers density", "Density", "kg/m3"),
+        "denstramax": Vart("Max tracers density", "Density", "kg/m3"),
+        "H_impacts": Vart("Impact heating", "Power", "W"),
+        "SurfNetRotX": Vart("Surface net rotation around X", "Rotation", "1/s"),
+        "SurfNetRotY": Vart("Surface net rotation around Y", "Rotation", "1/s"),
+        "SurfNetRotZ": Vart("Surface net rotation around Z", "Rotation", "1/s"),
     }
 )
 """Time series output by StagYY."""
@@ -405,18 +460,22 @@
     {
         "m": attrgetter("length"),
         "kg/m3": attrgetter("density"),
+        "kg": attrgetter("mass"),
+        "kg/s": lambda scl: scl.mass / scl.time,
         "K": attrgetter("temperature"),
         "W/m2": attrgetter("heat_flux"),
         "Pa": attrgetter("stress"),
         "Pa.s": attrgetter("dyn_visc"),
         "s": attrgetter("time"),
+        "W": attrgetter("power"),
         "W/(m.K)": attrgetter("th_cond"),
         "m2/s": attrgetter("th_diff"),
         "W/m3": attrgetter("heat_production"),
         "1/s": lambda scl: 1 / scl.time,
         "K/s": lambda scl: scl.temperature / scl.time,
         "m/s": attrgetter("velocity"),
         "m/s2": attrgetter("acceleration"),
+        "tons/yr": lambda scl: scl.mass / scl.time * 3.15e4,
     }
 )
 """Scales to make values dimensional, see [`Scales`][stagpy.dimensions.Scales]."""