I looked at it and it looks very messy to include the evaporation term. But maybe I'm missing something. So instead, I worked an analogous problem by filling an inverted right triangle (like the cross-section of a trough made with two pieces of wood) with area in the same way: We can fill the area at a rate of [itex]k[/itex] and let the evaporation rate be proportional to the length of the base (on top of it) with proportionality constant [itex]c[/itex]. So, I would write the change in area as:
[tex]\frac{dA}{dt}=f(A,t)=k-cL[/tex]
but I can get the length of the base as a function of A. For a right triangle, that's [itex]L(A)=\sqrt{4 A}[/itex] so then we have including evaporation:
[tex]\frac{dA}{dt}=k-c\sqrt{4A}[/tex]
Integrating:
[tex]\int_0^A \frac{dy}{k-c\sqrt{4y}}=\int_0^t dt[/tex]
And even that simple case is still very messy and the function A(t) is in terms of a Lambert function.
Now, if you used the same principle, in your case, expressing area of the pond surface as a function of volume beneath it, the DE I think would be:
[tex]\frac{dV}{dt}=\pi-0.01 A(V)[/tex]
but the function A(V) looks way too messy to do it symbolically.
Maybe someone else can help.