Abstract: Stability of multiple equilibria in the thermohaline double-diffusive system is studied here with numerical method in a body-fitted curvilinear coordinate system by introducing nonlinear equation of state, temperature-flux boundary condition and variable limit of forced boundary condition into this system. With heat-to-salt buoyancy ratio A larger than, equal to and smaller than 1, the flow is divided into temperature-dominated regime, temperature-salinity equilibrant regime and salinity-dominated regime. F...