is the same as figure \[fig:sol\_diag\_lowpre\] but it includes the pressure range of $P_\mathrm{c}>P_\mathrm{max}$. Figure \[fig:sol\_diag\] shows that most of the regions are covered by the C-c solution for $P_\mathrm{c}>P_\mathrm{max}$. The area of C-U solutions spreads to the pressure range of $P_\mathrm{max}<P_\mathrm{c}<2P_\mathrm{max}$. This is because solutions in this pressure range can connect with thermal equilibrium states before reaching at the critical point at a certain range of $j$. Figure \[fig:large\_flux\_hp\] shows trajectories of the solutions for $j/j_\mathrm{crit}=1$, 2, 5, and 20. The CNM pressure is set to $2.65\times10^{4}$ K cm$^{-3}$. One can see that the solution terminates at the critical point whose pressure is equal to $P_\mathrm{c}/2$ and where density increases with the mass flux.

![ The same as figure \[fig:sol\_diag\_lowpre\] but including the pressure range of $P_\mathrm{c}>P_\mathrm{max}$. The dashed line corresponds to $T_\mathrm{crit}=10^4$ K. []{data-label="fig:sol_diag"}](fig7.eps){width="8.0cm"}

![ Trajectories of solutions for $j/j_\mathrm{crit}=1$, 2, 5, and 20. The dashed line corresponds to $P=P_\mathrm{crit}/2$. []{data-label="fig:large_flux_hp"}](fig8.eps){width="8.0cm"}

Discussion {#sec:discuss}
==========

Implications
------------

@KI06 have investigated multi-dimensional nonlinear development of the TI without external forcing. They found that the turbulence is self-sustained by the TI. The velocity dispersion for sufficiently large simulation box is as large as $\sim0.2-0.4$ km s$^{-1}$. They suggested that this velocity 