Fixed bug calculating flux at south boundary faces