Bonjour. Je débute en fortran et j'ai quelques (voire beaucoup de) soucis. L'exercice que je suis en train d'effectuer doit permettre de représenter un faisceau d'isothermes en partant de l'équation de Van Der Waals (en P(V)).
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
PROGRAM vanderwaals
IMPLICIT NONE
 
REAL               :: P, V, T, x, Y, Z
REAL, PARAMETER    :: a=1, b=1, R=0.08206
INTEGER, DIMENSION (0:100, 100)::bob
REAL, EXTERNAL     :: Van (Y, Z)
bob(:)=0
 
Y=V
Z=T
 
OPEN(10, file="resultat.dat")
 
DO T=0.0, 100.0, 10.0
  DO V=0.0, 1.0, 0.1
 
  x=Van(T, V)
  write(10,*)i, bob(V,T)
 
  ENDDO
ENDDO
close(10)
END PROGRAM
 
 
FUNCTION Van (T, V, a, b)
REAL, PARAMETER ::a=1, b=1
REAL            ::V, T
 
 
Van= R*T/(V-b)-a/(V**2)
END FUNCTION
Voilà mon programme. Je ne suis pas du tout sûr de la déclaration de ma fonction et de mon tableau à 2 dimensions (je pense que ce doit être un des problèmes majeurs).

Merci d'avance pour votre aide.