sp_ab_gr.RdAutomatically assigns one color per species.
sp_ab_gr(simul, n_patches, sp, time)results array of a given simulation (a multigen() output)
an integer indicating the patch number
integer for the number of species to consider
integer determining the time step of interest