Modeling free-electron laser (FEL) oscillators requires calculation of both the light-beam interaction within the undulator and the propagation of the light outside the undulator. We present a paraxial Optical Propagation Code (OPC) based on the Spectral Method and Fresnel Diffraction Integral, which in combination with Genesis 1.3 can be used to perform either steady-state or time-dependent FEL oscillator simulations. A flexible scripting interface is used both to describe the optical resonator and to control the codes for propagation and amplification. OPC enables modeling of complex resonator designs that may include hard-edge elements (apertures) or hole-coupled mirrors with arbitrary shapes. Some capabilities of OPC are illustrated using the FELIX system as an example