Can’t get similar result with sigma coordinate and rectangular coordinate systems in the computation of density stratified flow.  Dose anybody know why?