The present numerical work combines gain switching and external optical feedback to generate pulses with a high level of coherence across a large array of semiconductor lasers. The influence of several parameters such as modulation frequency, feedback strength and frequency detuning is analyzed and this dataset gathers the temporal traces generated with our code, that is adapted from the Lang-Kobayashi model to take into account decayed coupling in an array of lasers.