#remove history
rm(list = ls(all = TRUE))
#working directory
setwd("C:/Users/Marie_PC/Documents/Projets/EUROCHAMP/eurochamp2020/Parois_WP2/handbook_EUROCHAMP/plots/")  

################################################################################
#Température (K)
Temp<-300
#Densité de l'air (molécules.cm-3)
airm<-2.50E+19
#------------------------
#Reaction2
#O3+NO =>  NO2 	
A2<-1.80E-12
#Ea/R
E2<-1370
#k
k2<-A2*exp(-E2/Temp)
#------------------------
#Reaction3
#NO2+HV => O3+NO   	
A3<-7.70E-03
E3<-0
k3<-A3*exp(-E3/Temp)
#------------------------
#Reaction6
#OH+NO2 => HNO3	
A6<-1.30E-11
E6<-0	
k6<-A6*exp(-E6/Temp)
#------------------------
#Reaction4
#NO+HO2 =>  NO2+OH 	
A4<-3.70E-12	
E4<-240	
k4<-A4*exp(-E4/Temp)
################################################################################
#concentrations intiales (molec.cm-3)
NO_0<-6.00E+11
NO2_0<-2.50E+10
O3_0<-1.00E+12
OH_0<-1.00E+05
HO2_0<-1.00E+05
HNO3_0<-0
################################################################################
#declaration des vecteurs
ti<-0
tf<-6000
dt<-1
Time<-seq(ti,tf,dt)
ntime<-length(Time)
NO<-numeric(ntime)
sources_NO<-numeric(ntime)
puits_NO<-numeric(ntime)
NO2<-numeric(ntime)
sources_NO2<-numeric(ntime)
puits_NO2<-numeric(ntime)
O3<-numeric(ntime)
sources_O3<-numeric(ntime)
puits_O3<-numeric(ntime)
HNO3<-numeric(ntime)
sources_HNO3<-numeric(ntime)
puits_HNO3<-numeric(ntime)
#initialisation
NO[1]<-NO_0
NO2[1]<-NO2_0
O3[1]<-O3_0
HNO3[1]<-HNO3_0
################################################################################
#resolution
for(i in seq(2,ntime)){

}

################################################################################
#plot