[en] The authors have developed the GNASH code to include photonuclear reactions for incident energies up to 140 MeV. Photoabsorption is modeled through the giant resonance at the lower energies, and the quasideuteron mechanism at the higher energies, and the angular momentum coupling of the incident photon to the target is properly accounted for. After the initial interaction, primary and multiple preequilibrium emission of fast particles can occur before compound nucleus decay from the equilibrated compound nucleus. The angular distributions from compound nucleus decay are taken as isotropic, and those from preequilibrium emission (which they obtain from a phase-space model which conserves momentum) are forward-peaked. To test the new modeling they apply the code to calculate photonuclear reactions on ²⁰⁸Pb for incident energies up to 140 MeV