Convective edge condition
The code runs as is, but there are strange results. Do not trust the values yet
Improved numerics, still getting strange results with realistic parameter values