The compression molding is widely used in the automotive industry to produce products that are large, thin, lightweight and stiff. The molded product is formed by squeezing a fiber-reinforced plastic composites. During a molding process of fiber reinforced thermoplastic composites, control of filling patterns in mold, orientation and distribution of fibers are needed to predict the effects of molding parameters on the flow characteristics. It is the objective of this paper to develop an isothermal compression molding simulation that can handle both thin and thick charges and motion of the flow front, and can predict pressure distributions and accurate velocity gradients. The composites are treated as an incompressible Newtonian fluid. The effects of slip parameter $alpha$ and extensional/shear viscosity ratio $zeta$ on the mold filling parameters are also discussed.