From 30f3ba4a6d7c8a4031fa5d4049255ae10b64e6bb Mon Sep 17 00:00:00 2001
From: "camille.poulet" <camille.poulet@irstea.fr>
Date: Wed, 14 Apr 2021 02:52:47 +0200
Subject: [PATCH] Add one more fonction to compute survival after spawning

---
 exploration/GR3D_Rdescription/GR3Dfunction.R | 13 +++++++++++--
 1 file changed, 11 insertions(+), 2 deletions(-)

diff --git a/exploration/GR3D_Rdescription/GR3Dfunction.R b/exploration/GR3D_Rdescription/GR3Dfunction.R
index ce75d06..d39cbe2 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))))
   
-- 
GitLab