Dear Firedrakers,

I'm trying to implement system (11)-(14) in the attached .pdf, but with Hcm being a Dirichlet boundary condition having a constant value 0.07 at y = 0.

The solution rises uniformly, where it should dip near y = 0, and it seems to do so regardless of the BC value I choose. Which leads to me to believe

that I'm not implementing it correctly.

How could I fix this?

Regards,
Will