これはいけました．ドライベットです． 流速0かつ凸または凹の処理は，中央差分の限界かもしれないですね． 何とか一工夫して使えるようにしたいです． なお，gifはここにあげるために強烈に画質を落としてます．

集中格子の場合 H(i-1)=H(i+1)=0.0,H(i)=1.0の場合，dH/dx=0 これは、各水面勾配の重みが0.５となり、打ち消しあって0になるようですね。 物理現象からすると、iの地点においてi+1とi-1の地点にそれぞれ流れが発生するようになるはずですが、 集中格子が原因…