This is the code described in Barceló et al. Projecting the time scale of initial increase in fishery yield after implementation of marine protected areas, submitted to ICES Journal of Marine Science (2020). These scripts build off of R code presented in Kaplan et al. (2019): Setting expectations for the timeline of fished population recovery for the adaptive management of a marine protected area network. Ecological Applications. 

Specific to Barceló et al. ICES JMS (in press 2020): 
Main analysis:

(Code is written in Matlab) 
figure_2.m
figure_3.m
figure_4.m

(Code is written in R)
figure_5_6.R

Specific to Barceló et al. ICES JMS (in press 2020): 
Supplementary materials:

(Code is written in Matlab)
MPA_sims.m
iterate_model.m
get_Leslie.m
define_Params.m
get_FLEP.m
get_F.m
Fig_SI_sim_eq2.m
porous_model.m