Simple numerical experiments on two-dimensional coastal upwelling are made with emphasis on the role of non-geostrophic solenoidal field of density in the formation of double-celled circulation and multi-celled density front. Geometry of shelf and slope is not taken into account. Existence of poleward undercurrent presumably caused by the longshore variation of the large scale pressure field is also suppressed for the sake of simplicity.