A spectral quadrilateral multidomain penalty method model for high Reynolds number incompressible stratified flows. (13th March 2014)