Add one more fonction to compute survival after spawning

30f3ba4a · Poulet Camille · a2fd77d0 · 30f3ba4a
Commit 30f3ba4a authored 4 years ago by Poulet Camille
Hide whitespace changes
Inline Side-by-side

Showing

with 11 additions and 2 deletions
+11 -2
--- 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))))