Output formats for DECAPODES

Here is an example of what output from DECAPODES could be for the SW simulations. This is just the advection diffusion equation.

```julia
AdvDiff = quote
C::Form0{X}
Ċ::Form0{X}
V::Form1{X}
ϕ::Form1{X}
ϕ₁::Form1{X}
ϕ₂::Form1{X}

# Fick's first law
ϕ₁ == (k ∘ d₀)(C)
ϕ₂ == ∧₀₁(C,V)
ϕ == ϕ₁ + ϕ₂
# Diffusion equation
∂ₜ(C) == Ċ
Ċ == ∘(⋆₀⁻¹, dual_d₁, ⋆₁)(ϕ)
end

advdiff = parse_decapode(AdvDiff)
advdiffdp = NamedDecapode(advdiff)
gensim(advdiffdp, [:C, :V])
```

Output

```julia
quote
function simulate(mesh)
begin
k = k(mesh)
d₀ = d₀(mesh)
(⋆₀⁻¹) = (⋆₀⁻¹)(mesh)
dual_d₁ = dual_d₁(mesh)
(⋆₁) = (⋆₁)(mesh)
(∧₀₁) = (∧₀₁)(mesh)
end
return begin
f((C, V)) = begin
var"•_1_1" = k(C)
ϕ₁ = d₀(var"•_1_1")
ϕ₂ = C ∧₀₁ V
ϕ = ϕ₁ + ϕ₂
var"•_3_1" = (⋆₀⁻¹)(ϕ)
var"•_3_2" = dual_d₁(var"•_3_1")
Ċ = (⋆₁)(var"•_3_2")
return Ċ
end
end
end
end
```