# -*- coding: utf-8 -*-

import numpy as np


#nombre de points
N=26
#champs de température au départ
T=np.linspace(80.0,80.0,N)+40*np.random.rand(N)
T1=np.linspace(80.0,80.0,N)
#déclaration des constantes avec theta air=0°C
e=0.02
cp=390
lamda=390       #on ne paut pas utiliser lambda
mu=8900
h=10
a=e/(N-1)
D=lamda/mu/cp

#calcul des constantes de tamps
tauplaque=e*e/D
tauair=mu*cp*e/2/h
taulim=a*a/2/D
pente=lamda/(h*a+lamda)


#affichage des constantes de temps
print('tauplaque=',tauplaque,'  s')
print('taulim=',taulim*1000,'  ms')
print('tauair= ',tauair/60,'  min')

temps=[]
Temperature=[]
t=0.0
#construction des 600 valeurs
for n in range(1,601):
    
    for k in range (0,5000):
        for p in range(N):
            T1[p]=T[p]
        for p in range(1,N-1):
            T[p]=???
        T[0]=???
        T[N-1]=???
        t=t+taulim
    Temperature.append(sum(T)/N)
    temps.append(t/60.0)
    #print('point n° ',n,'fait')

print('tfinal= ',t/60,'  min')
  
f=open("datater.txt",'w')
for i in range(0,599):
    f.write(str(temps[i])),f.write("\t"), f.write(str(Temperature[i])),f.write("\n")
f.close()

print ("C'est fini !")