Nearly got all models (full B, and full U B) working. Just have an MPI bug of some sort. The single core version matches with matlab full quasi-linear DNS code