added model.R which runs one model for all species/forest combinations (no error catching)