Let $p\colon \mathbb{R}^2 \to \mathbb{R}$ be a polynomial with a non-vanishing gradient at $p^{-1}(0)$. Then, the implicit function theorem says that $S = \{(x,y) \in \mathbb{R}^2 \mid p(x,y) = 0\}$ is an embedded smooth curve. In more algebraic terms, its coordinate ring $\mathbb{R}[x,y]/(p)$ has Krull dimension $1$ and its localization at every maximal ideal is regular.
Consider now a linear differential operator $T$ in two variables (ie some finite expression involving linear combinations of $\frac{\partial^{i+j}}{\partial x^i\partial y^j}$ with coefficients in $\mathbb{R}[x,y]$), and lets say it is acting formally on the space of all power series $\mathbb{R}[\![x,y]\!]$. I am not a priori concerned with convergence, but if needed we may restrict to all analytic functions defined on some open ball centered at the origin. Is there any meaningful way of saying that the set $$S = \{f \in \mathbb{R}[\![x,y]\!] \mid Tf = 0\}$$ looks locally like a one-variable power series ring (under suitable regularity conditions on $T$)? I'm particularly interested in a more formal, algebraic approach, that is not very dependend on the base field $\mathbb{R}$.