Two array-out-of-bound errors are fixed. They get the code quietly stuck on our (UF) HPC in 2D parallel simulations.
I'll run the test suite tonight and if it passes I will merge!
Just a quick followup: passes the test suite, i'm going to approve and merge. Nice catch on this one, @shuokong !