def METHOD_NAME(self):
    """Calculate soil production rate."""
    core = self._grid.core_nodes
    self.grid.at_node["soil_production__rate"][core] = self._w0[core] * np.exp(
        -self.grid.at_node["soil__depth"][core] / self._wstar[core]
    )