A <- c(3000, 4000, 5000, 2000)
A0 <- 70000
bounds <- c(100, 90, 70, 80)
m <- c(100, 90, 70, 50)
M <- bounds
ucosts <- c(0.7, 15, 2, 3) # unit costs

N_pop507 <- pop507s_ucost[, "N"]
A_pop507 <- N_pop507 * pop507s_ucost[, "S"]
ucosts_pop507 <- pop507s_ucost[, "unit_cost"]