kinematicSingleLayer::CourantNumber() decides on local number of facess
CourantNumber() calls surfaceSum (which requires parallel comms) under condition of local number of internal faces.
CourantNumber() calls surfaceSum (which requires parallel comms) under condition of local number of internal faces.