Create reef zones under different conditions
KdPARs <- seq(0.05, 0.2, by = 0.01) %>% # Choose variable and values to compare map_df(function(x) { # For each value data.frame(Light = seq(0.001, 1, length = 500)) %>% # Make a dataframe of light levels mutate("Bright depths" = depth(Light, KdPAR = x), # Add max depths given KdPAR "Shaded depths" = depth(Light, KdPAR = x, shade = 0.25), # Max depths when shaded "Whole reef" = reef(Light), # Calculate default reef zonation "KdPAR" = x) }) # Include KdPAR value
How about the boundary depths?
KdPAR_boundaries <- seq(0.05, 0.2, by = 0.01) %>% # Choose variable and values to compare map_df(function(x) { # For each value data.frame("Upper_bound" = depth(reef_boundary(), # Max depth for default reef boundary KdPAR = x), # Given a KdPAR value "Lower_bound" = depth(reef_boundary(), shade = 0.25, # Shaded depth for default reef boundary KdPAR = x), # Given a KdPAR value "KdPAR" = x)}) # Include changing variable
