We present a two-particle cluster variation theory to treat the inhomogeneous nematic system confined by two parallel planes which align the surface molecules perpendicular to them. Numerical calculations are made for a film composed of 20 molecular layers. The nematic-isotropic transition temperatures for typical values of the strength of the surface field are given and the layer dependences of the order parameter are shown. The Monte Carlo simulation data are needed to demonstrate our results.