Q-switched lasers are often used as a non-contact ultrasound source in non-destructive evaluation (NDE) of materials [1]. Q-switched lasers typically have ns pulse durations and generate broadband ultrasound waves, though longer laser pulses, of 100 microseconds or greater, have also been used [2] for NDE. These longer pulses tend to produce somewhat lower center frequencies than do Q-switched pulses, though they are still a broadband source. But it would be desirable in some NDE applications to narrow the signal bandwidth to improve the signal to noise ration (SNR), and also to have direct control over the center frequency of the generated ultrasound. In principle, this may be achieved by temporal [3,4] or spatial modulation [5,6] of the laser pulse, or both [7]. The purpose of this work was to develop a numerical model of a single, temporally modulated laser source of ultrasound in the thermoelastic regime, for isotropic metals.