The Reynolds stress model proposed in the 1st Report is applied to the development of the boundary layers in pressure gradients. First, predictions are made for two sink flow boundary layers : a laminarescent and a laminarized flow. The numerical solutions are extensively compared with experimental data including the distributions of the streamwise turbulence intensity. It is shown that the predictions are generally in good agreement with the data, reproducing a distinctive feature between the two flows. Next, two adverse gradient cases are tested. For the case of a non-equilibrium boundary layer in a strong pressure gradient, extensive comparisons of the numerical solutions with the data are made including the distribution of all the Reynolds stress components. It is shown that the overall agreement is satisfactory. However, a sound conclusion is not drawn concerning the difference in the logarithmic law intercept, which suggests a need for further consideration.