function prezone_Cs(Par_fix,Par,forcing,Temp,WE,landscape_per)

    include("forest_Cs")
    include("grassland_Cs")
    include("wetland_Cs")

    #include("CompleteMixing")
     Ks=Par[4]
     Ss_p=Par_fix
     prec=forcing[:,1]
     Ep=forcing[:,2]
      # O18_pre=forcing[:,3]
      H3_pre=forcing[:,3]#check
     #O18_river=forcing[:,5]
      tmax=length(prec)
     WEForest=WE[:,1]
     WEGrass=WE[:,2]
     WEWetland=WE[:,3]

     Ss=zeros(tmax,1)
     Sf=zeros(tmax,1)
     Qf=zeros(tmax,1)
     Qs=zeros(tmax,1)
     Qtot=zeros(tmax,1)
     SsTrVol=zeros(tmax,1)
     #  PE_IN=zeros(tmax,1)
       Melt=zeros(tmax,1)
       Su=zeros(tmax,1)
     # PE=zeros(tmax,1)
      E=zeros(tmax,1)
       # Ea=zeros(tmax,1)
     # H2_SsTrVol=zeros(tmax,1)
     Con_Qf=zeros(tmax,1)
     Con_Qs=zeros(tmax,1)
     Con_Qtot=zeros(tmax,1)
     # H2Con_Qf=zeros(tmax,1)
     # H2Con_Qs=zeros(tmax,1)
     # H2Con_Qtot=zeros(tmax,1)
    S_forest=zeros(1,2*Emax+4)
    S_grass=zeros(1,2*Emax+4)
    S_wet=zeros(1,Emax+3)
    STrVol_forest=zeros(2*Emax+4)
    STrVol_grass=zeros(2*Emax+4)
    STrVol_wet=zeros(Emax+3)

    #intinal value of the isotope concentration in Ss-p is the average value of 5%low discharge
    #intial O18=-7.76375
    #intial H2=-56.275
     # O18Tracer_intial=(Ss_p+Ss[1])*(-7.76375)
    # H2Tracer_intial=(Ss_p+Ss[1])*(-56.275)
     H3Tracer_intial=(Ss_p+Ss[1])*50.0
     #48.16#from1990
     #57.89# from 1988
     #121.11# from 1980
     #141.51#from1978
     #257.79#from1970
     #50.0#from1961



     # STrVol_wet[1,3]=O18Tracer_intial
    STrVol_wet[3]=H3Tracer_intial
    #SsTrVol_in[1]=(Ss_p+Ss[1])*-7.76375
    for i in 1:tmax
        ForestResult=[]
        GrassResult=[]
        WetlandResult=[]
        ForestResult= forestM_Cs(Par_fix,Par,S_forest,STrVol_forest,forcing[i,:],Temp[i,:],WEForest)
        GrassResult= grasslandM_Cs(Par_fix,Par,S_grass,STrVol_grass,forcing[i,:],Temp[i,:],WEGrass)
        WetlandResult= wetlandM_Cs(Par_fix,Par,S_wet,STrVol_wet,forcing[i,:],Temp[i,:],WEWetland,landscape_per[3])
        S_forest=ForestResult[1]
        S_grass=GrassResult[1]
        S_wet=WetlandResult[1]

        # Pre_forest=ForestResult[3][5]
        # Pre_grass=GrassResult[3][5]
        # Pre_wet=WetlandResult[3][3]
        #
        melt_forest=ForestResult[3][5]
        melt_grass=GrassResult[3][5]
        melt_wet=WetlandResult[3][3]
        #
        Ei_forest=ForestResult[3][1]
        Ei_grass=GrassResult[3][1]

        Ea_forest=ForestResult[3][2]
        Ea_grass=GrassResult[3][2]
        Ea_wet=WetlandResult[3][1]

       Su[i]=S_forest[2]*landscape_per[1]+S_grass[2]*landscape_per[2]+S_wet[1]*landscape_per[3]
       Sf[i]=S_forest[3]*landscape_per[1]+S_grass[3]*landscape_per[2]
       #  PE_IN[i]=Pre_forest*landscape_per[1]+Pre_grass*landscape_per[2]+Pre_wet*landscape_per[3]
      Melt[i]=melt_forest*landscape_per[1]+melt_grass*landscape_per[2]+melt_wet*landscape_per[3]
       # PE[i]=PE_IN[i]+Melt[i]
       Ei=Ei_forest*landscape_per[1]+Ei_grass*landscape_per[2]
       Ea=Ea_forest*landscape_per[1]+Ea_grass*landscape_per[2]+Ea_wet*landscape_per[3]
       E[i]=Ei+Ea

        #Tracer in storage
        STrVol_forest=ForestResult[2].*delay
        STrVol_grass=GrassResult[2].*delay
        STrVol_wet=WetlandResult[2]
    #Qf and Tracer in Qf
        Qf[i]=ForestResult[3][4]*landscape_per[1]+GrassResult[3][4]*landscape_per[2]+WetlandResult[3][2]*landscape_per[3]
        if Qf[i]>0.0
            Con_Qf[i]=(ForestResult[3][4]*ForestResult[4][4]*landscape_per[1]+GrassResult[3][4]*GrassResult[4][4]*landscape_per[2]+
            WetlandResult[3][2]*WetlandResult[4][2]*landscape_per[3])/Qf[i]
        else
            Con_Qf[i]=0.0
        end
        #Ss and tracer in Ss
        Ss_mix=WetlandResult[1][3]+ForestResult[3][3]*landscape_per[1]+GrassResult[3][3]*landscape_per[2]
        SsTrVol[i]=WetlandResult[2][3]+ForestResult[3][3]*ForestResult[4][3]*landscape_per[1]+
        GrassResult[3][3]*GrassResult[4][3]*landscape_per[2]
        # H2_SsTrVol[i]=WetlandResult[3][3]+ForestResult[4][3]*ForestResult[6][3]*landscape_per[1]+
        # GrassResult[4][3]*GrassResult[6][3]*landscape_per[2]
        Qs[i]= Ks*Ss_mix
        Qs[i]=min(Qs[i],Ss_mix)
        Ss[i]=Ss_mix-Qs[i];
        # Tracer slow reservior
        Con_Qs[i]=SsTrVol[i]/(Ss_mix+Ss_p)
        SsTrVol[i]=SsTrVol[i]-Con_Qs[i]*Qs[i]
        S_wet[3]=Ss[i];
        STrVol_wet[3]=SsTrVol[i]
         STrVol_wet=STrVol_wet.*delay
        # H2Con_Qs[i]=H2_SsTrVol[i]/(Ss_mix+Ss_p)
        # H2_SsTrVol[i]=H2_SsTrVol[i]-H2Con_Qs[i]*Qs[i]
        # STrVol_wet[2,3]=H2_SsTrVol[i]
        #Total outflow water
        Qtot[i]=Qs[i]+Qf[i]
        # Total Tracer Concentration
        if Qtot[i]>0.0
                Con_Qtot[i]=(Qs[i]*Con_Qs[i]+Qf[i]*Con_Qf[i])/Qtot[i]
            else
                Con_Qtot[i]=0.0
            end

    end

 return (Su,Sf,Qf,Qs)

 # ,WB,TracerB
end
