#remove history
rm(list = ls(all = TRUE))
#working directory
setwd("C:/Users/Marie_PC/Documents/Enseignements/M1_Model_Sys_Env/R")  

################################################################################
#Constante de Stefan-Bolzmann (w.m-2.K-4) 
sigma<-5.67E-08

#Constante solaire de Floréal #actuelle (W.m-2)
F0<-1370

#Albedo
as<-0.50 #sol nu
ab<-0.65 #paquerettes blanches
an<-0.35 #paquerettes noires
ag<-0.55 #paquerettes grises

#Croissance des paquerettes en fonction de la temperature : k(T)=1-b(To-T)^2					
To<-295.5
b<-3.27E-03
#Taux de mortalite des paquerettes
kp<-3.00E-01
#valeurs intiales
xn0<-0.001
xb0<-0.001
xg0<-0.0

#phi
phimin<-0.2
phimax<-2.5
dphi<-0.01
phi<-seq(from=phimin,to=phimax,by=dphi)
nphi<-length(phi)

#Temps
tmin<-0
tmax<-500
dt<-0.01
time<-seq(tmin,tmax,dt)
ntime<-length(time)
################################################################################