Nays2DHのソースコードの芦田・江頭・劉の式の差分式等をみていると(特にip1,im1,jp1,jm1を用いて条件分岐をしているあたりについて)，計算格子に対して，半mesh分ずれているのではないかと思いました．
この点について詳細に確認したいのですが，もし何か差分式等を確認できる資料があればご提示いただけないでしょうか？