# #!/usr/bin/env python
# -*-coding:Utf-8 -*-

#########  Effets temperature et temps d'incubation sur la vitesse enzyme (E. Jaspard - 2018) ############

import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np

R = 8.314472   			# Cte gaz parfaits
k = 1.3806503E-23 			# Cte Boltzmann
h = 6.62606876E-34 			# Cte Planck
Cat = 50E3   				# delta Gcat kJ.mol-1
Inact = 96E3   				# delta Ginact kJ.mol-1
E0 = 1e-5					# concentration molaire enzyme

def VitesseFonctionTemperature(temps,Temperature):
    return (((k*Temperature)/h)*np.exp(-(Cat/(R*Temperature))))*(E0)*np.exp(-(((k*Temperature)/h)*np.exp(-(Inact/(R*Temperature)))*temps))

def plot_VitesseFonctionTemperature():
    temps = np.arange(0, 100, 0.5) 			# Plage de temps etudiee (en secondes) / Syntaxe : arange(3, 15, 2) => array([ 3,  5,  7,  9, 11, 13])
    Temperature = np.arange(280, 350, 0.5)   # Plage de temperature etudiee (en Kelvin)
    temps,Temperature = np.meshgrid(temps,Temperature)
    Vitesse = VitesseFonctionTemperature(temps,Temperature)

    fig = plt.figure()
    ax = plt.axes(projection='3d')
    #ax.contour3D(temps,Temperature,Vitesse, 1000, cmap='binary')     # Representation beaucoup plus belle mais pas assez constrastee
    ax.plot_surface(temps,Temperature,Vitesse, 1500, cmap='binary')
 
    ax.set_xlabel('temps (s)')
    ax.set_ylabel('temperature (K)')
    ax.set_zlabel('vitesse (microM.s-1)')
    ax.view_init(10, 40)
    plt.show()

plot_VitesseFonctionTemperature()