1- Change ref pop size min = 0 by min = 30
Thus modified
  N_ref[N_ref == 0] <- NaN
   by
  N_ref[N_ref < 30] <- NaN

2- added warning message (if too many iterations with N_ref < 30)

3- updated description