Fix one issue where bc was being used instead of h_bc.
Fixed another (more important) issue where the slope limiting was partially disabled at physical boundares. It turns out that this approach causes under/overshoot when a sharp gradient of something approaches a physical boundary.
Fix one issue where bc was being used instead of h_bc.
Fixed another (more important) issue where the slope limiting was partially disabled at physical boundares. It turns out that this approach causes under/overshoot when a sharp gradient of something approaches a physical boundary.