Triflow is not able to handle 2D simution by itself. We have to help and provide our own routine. It will be slower than Triflow compiler, but enough for our computation. We have to provide:

The derivative F

A template for the fields (which describe their dimension and coordinates)

The Jacobian for implicits schemes. Lazy me, we will drop the jacobian for now.