diff --git a/exploration/GR3D_Rdescription/GR3Dfunction.R b/exploration/GR3D_Rdescription/GR3Dfunction.R index ce75d0625815810b711afeaf14583018823f2851..d39cbe230975e728babd3a875bfedbd61cf7148c 100644 --- a/exploration/GR3D_Rdescription/GR3Dfunction.R +++ b/exploration/GR3D_Rdescription/GR3Dfunction.R @@ -172,8 +172,17 @@ stockRecruitementRelationship <- function(Triver,Tmin, Topt, Tmax, survivalStock # ---------------------------------------------- # Spawner survival after reproduction #---------------------------------------------- -#Logit with temperature effect -spawnerSurvivalPostReproduction <- function(Triver,Tref, coeffa, coeffb){ +#Logit with temperature effect depending on a Triver +spawnerSurvivalPostReproductionTempRiver <- function(Triver, coeffa, coeffb){ + + spRiverPostSpawn = (coeffb/ (1 + exp(coeffa*(Triver)))) + + return(spRiverPostSpawn) + +} + +#Logit with temperature effect depending on the difference between the river temperature and Tref +spawnerSurvivalPostReproductionTempRef <- function(Triver,Tref, coeffa, coeffb){ spRiverPostSpawn = (coeffb/ (1 + exp(coeffa*(Triver-Tref))))