A k-ε turbulence model is developed to calculate wall turbulent shear flows under various pressure gradient conditions. In the present model, we set the dissipation rate of turbulent energy at a wall equal to zero, though the wall limiting behavior of velocity fluctuations is reproduced exactly. Thus, the model assures computational expediency and convergence. The proposed model is constructed to properly take into account the effect of pressure gradients on shear layers. It was found by Nagano, Tagawa and Tsuji (Turbulent Shear Flows 8, (1992), 7) that in adverse pressure gradient flows the Van Driest constant decreased with increasing dimensionless pressure gradient parameter P⁺. Therefore, the present model incorporates the modified Van Driest dumping function which is a function of P⁺. The validity of the proposed model was tested by application to a turbulent channel flow and boundary layers with P⁺ < 0, P⁺ = 0 and P⁺ > 0. The model predictions indicate that agreement with the experiment and the direct simulation data is very good over a wide range of pressure variations.