rm ciover, cfover, advth extra rprof

Those are only relevant for very specific problems and therefore
should be in user scripts rather than directly exposed in StagPy.

Author: amorison

src/stagpy/phyvars.py

@@ -277,9 +277,6 @@
         "advds": processing.advds_prof,
         "advas": processing.advas_prof,
         "energy": processing.energy_prof,
-        "ciover": processing.init_c_overturn,
-        "cfover": processing.c_overturned,
-        "advth": processing.advth,
     }
 )
 """Radial profiles that StagPy can compute."""

src/stagpy/processing.py

@@ -240,88 +240,6 @@ def energy_prof(step: Step) -> Rprof:
     )
 
 
-def advth(step: Step) -> Rprof:
-    """Compute theoretical advection.
-
-    This compute the theoretical profile of total advection as function of
-    radius.
-
-    Args:
-        step: a `Step` of a `StagyyData` instance.
-
-    Returns:
-        the theoretical advection and radius.
-    """
-    rbot, rtop = step.rprofs.bounds
-    rmean = 0.5 * (rbot + rtop)
-    rad = step.rprofs.centers
-    radio = step.timeinfo["H_int"]
-    if rbot != 0:  # spherical
-        th_adv = -(rtop**3 - rad**3) / rmean**2 / 3
-    else:
-        th_adv = rad - rtop
-    th_adv *= radio
-    th_adv += step.timeinfo["Nutop"]
-    return Rprof(th_adv, rad, Varr("Theoretical advection", "Heat flux", "W/m2"))
-
-
-def init_c_overturn(step: Step) -> Rprof:
-    """Compute concentration before overturn.
-
-    This compute the resulting composition profile if fractional
-    crystallization of a SMO is assumed.
-
-    Args:
-        step: a `Step` of a `StagyyData` instance.
-
-    Returns:
-        the composition and radius.
-    """
-    rbot, rtop = step.rprofs.bounds
-    par = step.sdat.par
-    xieut = par.nml["tracersin"]["fe_eut"]
-    k_fe = par.nml["tracersin"]["k_fe"]
-    xi0l = par.nml["tracersin"]["fe_cont"]
-    xi0s = k_fe * xi0l
-    xired = xi0l / xieut
-    rsup = (rtop**3 - xired ** (1 / (1 - k_fe)) * (rtop**3 - rbot**3)) ** (1 / 3)
-
-    def initprof(rpos: float) -> float:
-        """Theoretical initial profile."""
-        if rpos < rsup:
-            return xi0s * ((rtop**3 - rbot**3) / (rtop**3 - rpos**3)) ** (1 - k_fe)
-        return xieut
-
-    rad = np.linspace(rbot, rtop, 500)
-    initprof_vec = np.vectorize(initprof)
-    return Rprof(
-        initprof_vec(rad),
-        rad,
-        Varr("Concentration before overturn", "Concentration", "1"),
-    )
-
-
-def c_overturned(step: Step) -> Rprof:
-    """Compute theoretical overturned concentration.
-
-    This compute the resulting composition profile if fractional
-    crystallization of a SMO is assumed and then a purely radial
-    overturn happens.
-
-    Args:
-        step: a `Step` of a `StagyyData` instance.
-
-    Returns:
-        the composition and radius.
-    """
-    rbot, rtop = step.rprofs.bounds
-    cinit = init_c_overturn(step)
-    radf = (rtop**3 + rbot**3 - cinit.rad**3) ** (1 / 3)
-    return Rprof(
-        cinit.values, radf, Varr("Overturned concentration", "Concentration", "1")
-    )
-
-
 def stream_function(step: Step) -> Field:
     """Compute the stream function in 2D geometry.