Programme Fortran ne fonctionnant pas

**gibobo** · 24/03/2007, 18h19

ça fait plusieurs jours que je n'arrive pas a faire fonctionner ce programme quelqu'un peut m'aider à trouver les erreurs:

fichier d'entrée:

19-03-2007
m kg/m3
Hspirale P(1) T(1) Gamma rgaz Lignes
0.11648 625000 291.5 1.1311 100 25

UV UC RC UN RN UK RK SURF
vilbrequin u r u r u r 1 poche
rad rad m rad m rad m m2
U C C N N K K S
0.000 1.484 0.335 -1.658 0.310 -4.800 0.280 0.055
-1.047 0.436 0.325 -2.705 0.298 -5.847 0.270 0.053
-2.094 -0.611 0.323 -3.752 0.295 -6.894 0.268 0.051
-3.142 -1.658 0.310 -4.800 0.280 -7.941 0.255 0.050
-4.189 -2.705 0.298 -5.847 0.270 -8.988 0.245 0.047
-5.236 -3.752 0.295 -6.894 0.268 -10.036 0.235 0.046
-6.283 -4.800 0.280 -7.941 0.255 -11.083 0.220 0.045
-7.330 -5.847 0.270 -8.988 0.245 -12.130 0.210 0.043
-8.378 -6.894 0.268 -10.036 0.235 -13.177 0.205 0.041
-9.425 -7.941 0.255 -11.083 0.220 -14.224 0.185 0.039
-10.472 -8.988 0.245 -12.130 0.210 -15.272 0.180 0.037
-11.519 -10.036 0.235 -13.177 0.205 -16.319 0.165 0.035
-12.566 -11.083 0.220 -14.224 0.185 -17.366 0.155 0.032
-13.614 -12.130 0.210 -15.272 0.180 -18.413 0.150 0.031
-14.661 -13.177 0.205 -16.319 0.165 -19.460 0.140 0.030
-15.708 -14.224 0.185 -17.366 0.155 -20.508 0.130 0.028
-16.755 -15.272 0.180 -18.413 0.150 -21.555 0.125 0.026
-17.802 -16.319 0.165 -19.460 0.140 -22.602 0.110 0.024
-18.850 -17.366 0.155 -20.508 0.130 -23.649 0.100 0.022
-19.897 -18.413 0.150 -21.555 0.125 -24.696 0.000 0.001
-20.944 -19.460 0.140 -22.602 0.110 -25.744 0.000 0.001
-21.991 -20.508 0.130 -23.649 0.100 -26.791 0.000 0.001
-23.038 -21.555 0.125 -24.696 0.000 -27.838 0.000 0.001
-24.086 -22.602 0.110 -25.744 0.000 -28.885 0.000 0.001
-25.133 -23.649 0.100 -26.791 0.000 -29.932 0.000 0.001

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
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
 
       PROGRAM   D
       INTEGER NT,N,IL
       REAL Roe,PINT,PEXT,DP,C,D,HT,k,rgaz,a,b,Fx,Fy,F,Ga,m
 
       DIMENSION UV1(3600),UC1(3600),RC1(3600),UN1(3600),RN1(3600),
     &  UK1(3600),RK1(3600),SURF1(3600)
 
 
       DIMENSION UV(3600),UC(3600),RC(3600),UN(3600),RN(3600),
     &  UK(3600),RK(3600),P(3600),SURF(3600),PX(3600),m(3600),T(3600),
     &  Fx1(3600),Fy1(3600),Fx2(3600),Fy2(3600),Fx(3600),Fy(3600),
     &  F1(3600),F2(3600),F(3600),MMT1(3600),MMT2(3600),MMT(3600)
 
       double precision pi
 
       OPEN(UNIT=5,status='OLD',form='formatted',FILE='tab4.txt') 
       Read (5,*) 
       Read (5,*)
       Read (5,*)
       Read (5,*) 
       read (5,*) HT,P(1),T(1),Ga,rgaz,IL
 
       write (*,*)'    HT = ',HT
       write (*,*)'    P(1)=',P(1)
       write (*,*)'    T(1)=',T(1)
       write (*,*)'    Ga = ',Ga
       write (*,*)'    IL = ',IL
 
       Read (5,*) 
       Read (5,*)
       Read (5,*)
       Read (5,*) 
       Read (5,*)
 
 
       DO N=1,IL   
       READ(5,*) UV1(N),UC1(N),RC1(N),UN1(N),RN1(N),UK1(N),RK1(N),
     &           SURF1(N)  
 
 
         UV(N)=UV1(N)
         UC(N)=UC1(N)     
         RC(N)=RC1(N)
         UN(N)=UN1(N) 
         RN(N)=RN1(N)        
         UK(N)=UK1(N)         
         RK(N)=RK1(N)      
         SURF(N)=SURF1(N)
         write (*,*)'UV(N),UC(N),RC(N),UN(N),RN(N),UK(N),RK(N),
     &              SURF(N)'
 
        write (*,*)UV(N),UC(N),RC(N),UN(N),RN(N),UK(N),RK(N),
     &              SURF(N)
 
       END DO
 
 
       write(*,*) 'interpolation'
      !    go to 101 
           call interpol(IL,360,UV1,UC1,UV,UC)
101   continue
 
 
       read(*,*)     
       write(*,*) 'début du traitement'
       read(*,*)
 
C     'Calculer la masse,pression,PX et la température'
       N = 1
       m(N) = P(1)*(HT*SURF(1))/(rgaz*T(1))
       write(*,*) 'premier m(1)=',m(N)
       read(*,*)
       DO N=2,IL
         m(N)=m(N-1)
         P(N)=P(N-1)*(SURF(N-1)/SURF(N))**(1/Ga)
         write(*,*) 'm(N),P(N)',m(N),P(N)
         read(*,*)   
         IF (N.LE.6) THEN
              PX(N)=P(1)
              ELSE
              PX(N) = P(N-6)
         END IF
 
         T(N)=T(N-1)*(SURF(N-1)/SURF(N))**(Ga-1)
         write(*,*) 'N = ',N,' m = ',m(N),' P = ',P(N),' T = ',T(N)
       END DO
 
       read(*,*)
c      Calcul de la force      
 
       DO N=1,IL
 
        CALL FORCE(HT,P(N),PX(N),UC(N),RC(N),UN(N),RN(N),Fx1(N),
     &       Fy1(N),F1(N))
 
        CALL FORCE(HT,P(N),PX(N),UN(N),RN(N),UK(N),RK(N),Fx2(N),
     &       Fy2(N),F2(N))
        Fx(N) = Fx1(N) + Fx2(N)
        Fy(N) = Fy1(N) + Fy2(N)
        F(N) = sqrt(Fx(N)**2 + Fy(N)**2)
      !write(*,*) F(N)
       END DO
C      Calcul du moment      
 
       DO N=1,IL
 
        CALL MOMENT(UC(N),RC(N),UN(N),RN(N),Fx1(N),Fy1(N),MMT1(N))
        CALL MOMENT(UN(N),RN(N),UK(N),RK(N),Fx2(N),Fy2(N),MMT2(N))
        MMT(N) = MMT1(N)+MMT2(N)
      !write(*,*) MMT(N)
       END DO 
 
      end
C     SUBROUTINES
 
      SUBROUTINE FORCE(HT,PINT,PEXT,CC,rr,DD,kk,FFx,FFy,FF)
 
      !          FORCE(HT,P(N),PX(N),UC(N),RC(N),UN(N),RN(N),Fx1(N),
      !  &       Fy1(N),F1(N))
 
 
 
      INTEGER NX
 
      REAL HT,PINT,PEXT,CC,rr,DD,kk,FFx,FFy,FF,DPX,aa,bb
      double precision pi
      pi=3.14159
 
      DPX = PEXT-PINT
 
      aa=-rr*cos(CC)+kk*cos(DD)
      bb=-rr*sin(CC)+kk*sin(DD)
      FFx=aa*H*DPX
      FFy=bb*H*DPX
      FF=sqrt(FFx*FFx+FFy*FFy)
 
      WRITE (*,*) '    DPX = ',DPX,'    r = ',rr,'    k = ',kk,
     f'    a = ',aa,'    b = ',bb,'    Fx = ',FFx,'    Fy = ',FFy,
     &'    F = ',FF
 
      READ (*,*)    
      RETURN    
      END
 
 
      SUBROUTINE MOMENT(CC,rr,DD,kk,Fx,Fy,MMT)
 
 
      INTEGER NX
 
      REAL CC,rr,DD,kk,Fx,Fy,MMT,aa,bb
      double precision pi
      pi=3.14159
 
 
      aa=(rr*cos(CC)+kk*cos(DD))/2
      bb=(rr*sin(CC)+kk*sin(DD))/2
 
      CALL PRODVECT(aa,bb,0,Fx,Fy,0,MX,MY,MZ)
      MMT = MZ
 
      WRITE (*,*) '    DPX = ',DPX,'    r = ',rr,'    k = ',kk,
     f'    a = ',aa,'    b = ',bb,'    Fx = ',Fx,'    Fy = ',Fy,
     g '    MMT = ',MMT
 
      READ (*,*)    
      RETURN    
      END
 
 
      SUBROUTINE PRODVECT(a,b,c,Fx,Fy,Fz,MX,MY,MZ)  
 
      INTEGER NX
 
      REAL a,b,c,Fx,Fy,Fz,MX,MY,MZ
      double precision pi
      pi=3.14159
 
      MX = b*Fz-c*Fz
      MY = -a*Fz-c*Fx
      MZ =a*Fy-b*Fx
 
 
      WRITE (*,*) 'a,b,c,Fx,Fy,Fz,MX,MY,MZ',
     &             a,b,c,Fx,Fy,Fz,MX,MY,MZ
 
      READ (*,*)    
      RETURN    
      END
 
 
      subroutine INTERPOL(IMAX,KMAX,UV1,UC1,UV,UC)
 
C***********************************************************************
C OBJET   : Interpolations des données contenues dans le fichier       *
C	       d'entrées                                               *
C                                    	                               *
C ENTREES : LIGNE   Nombre de lignes composant le fichier d'entrées    * 
C                    des surfaces                                      *                                   
C	    ITMAX   Itération maximale ou encre angle maximal          *
C	    SURF    Tableau de "LIGNE" lignes contenant les differentes*			      		        
C		    surfaces de fuites				       *
C ES      : SURF   Tableau de "ITMAX" lignes  contenant les differentes*
C		    surfaces de fuites                                 * 
C            	                   				       	    		   *                          
C COMMENT :                                                            *
C MODIF   : JJ/MM/AAAA AUT : Creation                                  *
C***********************************************************************	
 
        dimension  UV1(3600),UC1(3600),UV(3600),UC(3600)
    	real  UVINF,UVSUP,UCINF,UCSUP,pi
    	integer IMAX,KMAX,N,IL
 
        pi=3.14159
 
       UV(1) = UV1(1)
	   UV(KMAX) = UV1(IMAX)
       UC(1) = UC1(1)
	   UC(KMAX) = UC1(IMAX)
 
	  WRITE (*,*) 'INTERPOL xxxxxxxxxxxxxxxxx'
	  READ (*,*)	
C 	INTERPOLATION DES DONNEES CONTENUES DANS  TABINT1 VERS TABINT
       I = 1
      UVINF = UV1(I)
	  UVSUP = UV1(I+1)
	  UCINF = UC1(I)
	  UCSUP = UC1(I+1)
 
	  DO 100 K=2,KMAX-1
	    WRITE(*,*) 'K = ',K
 
        UV(K) = -(K*2*pi)/KMAX
		IF( UV(K).GT.UV1(I+1) )  THEN	      
	      I = I + 1
		  UVINF = UV1(I)
	      UVSUP = UV1(I+1)
	      UCINF = UC1(I)
	      UCSUP = UC1(I+1)
	      WRITE(*,*)'UV(K) = ',UV(K)
		  WRITE (*,*) UVINF,' ',UVSUP,'   ',UCINF,' ',UCSUP
           WRITE(*,*) ' I =',I,' UV1 =',UV1(I),' UC1 =', UC1(I)
		  READ (*,*)                       
								  ELSE
 
          ENDIF
 
          UC(K)=UCINF+(UV(K)-UVINF)*(UCSUP-UCINF)/(UVSUP-UVINF) 
		WRITE (*,*) 'UC(K) = ',UC(K)		
100	   CONTINUE
	  RETURN	
	  END

**Sylvain Bergeron** · 26/03/2007, 13h33

Le source tel que présenté ni respecte ni le format fixe, ni le format libre. En reformatant le tout vers le format fixe, j'obtiens les messages d'erreurs suivants :

Compiling file EDIT1.F.
Compiling program unit D at line 1:
1710-S: "EDIT1.F", line 3: D already used as a program name.
Compiling program unit FORCE at line 114:
2005-W: "EDIT1.F", line 133: H is used but never set.
Compiling program unit MOMENT at line 144:
2005-W: "EDIT1.F", line 162: DPX is used but never set.
Compiling program unit PRODVECT at line 169:
Compiling program unit INTERPOL at line 190:
Encountered 1 error, 2 warnings in file EDIT1.F.

En changeant le nom du programme pour DProg, j'obtiens :

...
K = 23
UV(K) = -0.401425391
-23.0380001 -24.0860004 -21.5550003 -22.6019993
I = 23 UV1 = -23.0380001 UC1 = -21.5550003
UC(K) = 1.05994534
K = 24
UV(K) = -0.418878675
-24.0860004 -25.1329994 -22.6019993 -23.6490002
I = 24 UV1 = -24.0860004 UC1 = -22.6019993
UC(K) = 1.06516564
K = 25
The variable (UV1(26)) has an undefined value.
Error occurs at or near line 237 of _interpol_
Called from or near line 60 of _MAIN__

Commentaires :

Tu devrais toujours utiliser la commande implicite none
dimension est un archaisme, surtout qu'elle est incompatible avec implicit none
Tu ne devrais pas utiliser le unit 5 pour lire un fichier. Prends des units > 10. Pour mon compilateur, read(5 et read(*, c'est la même chose...

**genteur slayer** · 26/03/2007, 16h18

alors, avant toute chose, utilise le [CODE] pour tes source, c plus facil à lire après.... de plus, je suis d'accords que ton prog ne respect ni le format libre ni le fixe

moi je prèfère le format libre, mais EXPLICIT !!! c'est çà dire en

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

implicit none

de plus,

dimension est un archaisme, surtout qu'elle est incompatible avec implicit none

je suis archi pas d'accord!!! comment tu défini des tableaux sinon!!!!

une syntaxe correct pour déclarer un tableau (en norme 95):

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

Real(kind=kind(1.d0)), dimension(<taille>) :: <nom>

utilise efectivement des unit entre 10 et 100 (il ne sont pas réservé par le compilo)

un autre truc aussi, utilise les commande -Wall et -Wextra du compilo (gcc,g77,g95) "Warning all" et "Warning Extra" ainsi que l'option "-fbound-check" lors du link de g95 qui te raportera des erreurs dès la compilation... un example de Makefile:

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
FC=i686-pc-cygwin-g95
Cflag=-info -Wall -Wextra -L./ -I./ -fshort-circuit -fzero -O4 -fno-underscoring
Lflag=-info -Wall -Wextra -L./ -I./ -fbound-check -O4
 
all: chutebille2.exe
 
chutebille2.exe: main.o
	$(FC) -o chutebille2.exe main.o $(Lflag)
 
main.o: main.f95
	$(FC) -c main.f95 $(Cflag)
 
clean:
	rm -rf *.o *.mod *~
 
mrproper: clean
	rm -rf chutebille.exe

Sinon, réagit en programmeur: 3lignes de commentaire par ligne de code: ecris d'abord en français ce que DOIT faire la fonction avant de la programmer, ensuite tu compare le français au fortran, ce doit être la Même chose!!!

**Sylvain Bergeron** · 26/03/2007, 17h04

Envoyé par genteur slayer

de plus, je suis archi pas d'accord!!! comment tu défini des tableaux sinon!!!!

une syntaxe correct pour déclarer un tableau (en norme 95):

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

Real(kind=kind(1.d0)), dimension(<taille>) :: <nom>

Je persiste, la commande «dimension» est un archaisme. Je parle de la commande, pas de l'attribut. Gibobo utilise la commande dans son code. Pour utiliser la commande dimension avec un implicit none, il faut faire :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
«type» Variable
dimension Variable(...)

Il est beaucoup plus simple de dimensionner le tableau en même temps que la déclaration du type. En F95, tu as les choix suivants :

«type», dimension(...) :: Tableau, ou
«type» :: Tableau(...), ou
«type» Tableau(...)

En F77, tu dois utiliser la 3ième forme.

En pratique, quand tu insères des implicit none, tu remplace simplement les commandes dimension par le type pertinent. Pour Gibobo, ça veut dire qu'il doit remplacer :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
DIMENSION UV1(3600),UC1(3600),RC1(3600),UN1(3600),RN1(3600),
& UK1(3600),RK1(3600),SURF1(3600)

par :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
real UV1(3600),UC1(3600),RC1(3600),UN1(3600),RN1(3600),
& UK1(3600),RK1(3600),SURF1(3600)

**gibobo** · 27/03/2007, 00h50

Je tient tout d'abord a vous remercier vous m'avez appris pas mal de petit truc que je ne connaissez pas .Malheureusement a la 25 boucle le programme ne fonctionne pas.Je suis bloqué aider avec l'erreur:

Citation:

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
...
K = 23
UV(K) = -0.401425391
-23.0380001 -24.0860004 -21.5550003 -22.6019993
I = 23 UV1 = -23.0380001 UC1 = -21.5550003
UC(K) = 1.05994534
K = 24
UV(K) = -0.418878675
-24.0860004 -25.1329994 -22.6019993 -23.6490002
I = 24 UV1 = -24.0860004 UC1 = -22.6019993
UC(K) = 1.06516564
K = 25
The variable (UV1(26)) has an undefined value.
Error occurs at or near line 237 of _interpol_
Called from or near line 60 of _MAIN__

désolé si j'ecris pas avec le bon format mais j'utilise le bon vieux copier coller
Merci de m'aider à corriger cette erreur

Programme Fortran ne fonctionnant pas

Fortran

Vue hybride

Discussions similaires

Partager

Partager