Petit probleme fortran95

[invite395ab1e7]

Bonsoir
Ca fais de jours que je me creuse la tête sur ce petit problème mais je vois pas ce qui cloche je vous met mon prog en dessous et l'erreur qui me renvoie en esperant que vous allez pouvoir m'aider

program tp3
implicit none
!Déclaration des variables
real::qs,A0,dz,z,P
integer::i,T1,Cond0,L,h,j,k
!Déclaration des tableaux
integer,parameter::n=1001
real,dimension (n):: T
real,dimension (n):: Cond

!Présentation du programe
write(*,*)'===================================='
write(*,*)'| Programme Equation de la chaleur |'
write(*,*)'|          Version 1.0             |'
write(*,*)'|       Created by K.Freville      |'
write(*,*)'|            11/01/2012            |'
Write(*,*)'===================================='
! Lecture des entrées/sorties
open(10,file="tp3.in")
open(11,file="tp3.out")
!T0 T°initial, qs flux de surface,K0 conductivité initial,A0 production radiogénique,h profondeur caractériqtique, L profondeur
read(10,*)T1
read(10,*)qs
read(10,*)Cond0
read(10,*)A0
read(10,*)h
read(10,*)L
!*********************************************************************************************************************************************************
!*********************************************************************************************************************************************************
!*********************************************************************************************************************************************************
!Coeur du programe
!Décalration de dz
dz=L/(n-1)
		!Calcule T1 et T2 et boucle d'écriture
		Do k=1,n
			z=dz*(k-1)
			T(1)=T1
			T(2)=((qs*dz)/Cond0)+T(1)
			
		!Boucle des T°
			Do i=3,n
				
				!Déclaration des conductivitées initiales
					Cond(i)=Cond0
					Cond(i-1)=Cond0
				
			
					!Boucle des conductivitées
					Do j=1,5
						!Calcule de P
						!Pour le manteau lithosphèrique
						if (z>30000) then
							P=3300*9.81*z
						!Pour la croute
						elseif (z<=30000) then
							P=2900*9.81*z
						endif
							!solution par différence fini
							T(i)=T(i-1)+((-A0*dz*dz)/Cond(i))*exp((-dz*(i-1))/h)+(Cond(i-1)/Cond(i))*(T(i-1)-T(i-2))
							!Calcul des conductivite
						if (z>30000) then
							!Déclaration des conductivitées initiales
								
							Cond(i)=((1/(O.156+5.45E-4*T(i)))+(0.763E-9*T(i)**3))*(1+0.1*P(i))
								
							Cond(i-1)=((1/(O.156+5.45E-4*T(i-1)))+(0.763E-9*T(i-1)**3))*(1+0.1*P(i-1))
						elseif (z>=30000) then
							Cond(i)=Cond0
							Cond(i-1)=Cond0
						endif
					Enddo
			!
			
			Enddo
			write(11,*)-z,T(k),Cond(k)
		Enddo

Close(10)
Close(11)
!*********************************************************************************************************************************************************
!*********************************************************************************************************************************************************
!*********************************************************************************************************************************************************

stop 
end

Et il me renvoie
[kevin@localhost tp3]$ gfortran tp3.f95 -o tp3
tp3.f95:64.21:

Cond(i)=((1/(O.156+5.45E-4*T(i)))+(0.763E-9*T(i)**3))*(1+0.1*P(i))
1
Error: Expected a right parenthesis in expression at (1)
tp3.f95:66.23:

Cond(i-1)=((1/(O.156+5.45E-4*T(i-1)))+(0.763E-9*T(i-1)**3))*(1+0.1*P(i-1
1
Error: Expected a right parenthesis in expression at (1)




Merci d'avance pour votre aide
Kévin
-----

[invite4492c379]

Hello,

heu ... si tu remplace le O par un 0 ? ça fonctionne mieux ?????

[kevin@localhost tp3]$ gfortran tp3.f95 -o tp3
tp3.f95:64.21:

Cond(i)=((1/(O.156+5.45E-4*T(i)))+(0.763E-9*T(i)**3))*(1+0.1*P(i))
1
Error: Expected a right parenthesis in expression at (1)
tp3.f95:66.23:

Cond(i-1)=((1/(O.156+5.45E-4*T(i-1)))+(0.763E-9*T(i-1)**3))*(1+0.1*P(i-1
1
Error: Expected a right parenthesis in expression at (1)

Peut-être aussi déclarer P comme un tableau, je suppose que ça peut aider si tu utilises P comme un tableau.

[Calvert]

Ce que je vois dans un premier temps, c'est que les "0" de 0.156 ne sont pas des "zéros" mais de "O" majuscules.

Edit: grilled !

[invite395ab1e7]

Merci il y avais bien un petit probleme de O a la place du 0 ceci corriger le probleme persiste et pour ce uqi est du P il ne me semble pas neccessaire de l'inclure d'en un tableau vu que je ne m'en sert q'une fois dans la boucle et que je ne réutilise pas la valeur précédente?


Edit: c bon j'ai régler le pb j j'utiliser P comme un tableau alors que non lol merci a vous

[A voir en vidéo sur Futura]

[invite4492c379]

Merci il y avais bien un petit probleme de O a la place du 0 ceci corriger le probleme persiste et pour ce uqi est du P il ne me semble pas neccessaire de l'inclure d'en un tableau vu que je ne m'en sert q'une fois dans la boucle et que je ne réutilise pas la valeur précédente?

Le problème persiste ? Tu peux préciser ? Je ne connais pas fortran alors il faut que tu sois tres précis.

[Calvert]

Si tu as encore des problèmes, reposte ton nouveau code, ainsi que les messages du compilateur.

[invite395ab1e7]

Merci a tous tou fonctionne comme sur des roulettes je passe au programme suivant maintenant et oui les joies de surbooking scolaire