Discussion :

[Rémirem]

Bonjour,

j'essaie de rentrer un chemin conditionnel selon des valeurs trouvées dans un fichier texte mais cela ne marche pas...

J'obtiens au final la reponse "chemin"...
Voici mon code:

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
 DO 440 I=1,IG 
         IF(PR(I).EQ.'IPE')THEN 
           FICH=CHEMIN//'ipe.txt'
           JP=18  
           TYP=1
         ENDIF
         IF(PR(I).EQ.'HEA')THEN
           FICH=CHEMIN//'hea.txt'
           JP=24  
           TYP=1
         ENDIF
         IF(PR(I).EQ.'HEB')THEN
           FICH=CHEMIN//'heb.txt'
           JP=29  
           TYP=1
         ENDIF
         
WRITE(*,*)FICH
         OPEN(UNIT=40,
     $     FILE=FICH,
     $     STATUS='OLD',
     $     ACCESS='SEQUENTIAL',
     $     FORM='FORMATTED',
     $     IOSTAT=IOS)

Le write en rouge me renvoie la variable CHEMIN.

Si vous pouvez m'aider, merci d'avance !

[DaZumba]

Si FICH est initialise a CHEMIN quelque-part avant les if, c'est que tes conditions ne sont jamais vraies. Sinon, c'est que FICH est un tableau trop petit pour contenir la chaine que tu ajoutes.

[Rémirem]

Merci pour votre réponse!

Le problème est la, FICH n'est pas initialisé à CHEMIN et est défini comme suit:
CHARACTER*52 FICH, CHEMIN

Je traine la dessus depuis ce matin

Ben le temps que j'y suis, je viens de taper un bête:
write(*,*)coucou
et j'ai comme réponse sur ma console:
-1.#QNAN

Je sais que le programme marché il y a plusieurs années sous Linux mais que depuis le passage à WIndows, une routine possédant une lecture direct d'un fichier texte ne marche plus. Je pensais que le problème venais de là, donc j'ai essayé de passer en sequentiel mais ca n'a rien changé...
Si vous avez une idée...?
Moi j'suis sec

voila la routine:

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
 
      SUBROUTINE NOMPR(CHEMIN,NG,IG,PR,JP,CAR,IT,NON)
      CHARACTER*4 PR(NG) 
      CHARACTER*52 FICH,CHEMIN      
      INTEGER TYP
      DOUBLE PRECISION IT(NG),A,C1,C2,C3,NON(NG,3)
      DOUBLE PRECISION CAR(NG,13),LARGEUR,EPSEM       
      DOUBLE PRECISION A1,B1,D1,R1,ITOR1
 
      DO 440 I=1,IG 
         IF(PR(I).EQ.'IPE')THEN 
           FICH=CHEMIN//'ipe.txt'
           JP=18  
           TYP=1
         ENDIF
         IF(PR(I).EQ.'HEA')THEN
           FICH=CHEMIN//'hea.txt'
           JP=24  
           TYP=1
         ENDIF
         IF(PR(I).EQ.'HEB')THEN
           FICH=CHEMIN//'heb.txt'
           JP=29  
           TYP=1
         ENDIF
         IF(PR(I).EQ.'HEM')THEN
           FICH=CHEMIN//'hem.txt'
           JP=24
           TYP=1  
         ENDIF
         IF(PR(I).EQ.'UPN')THEN
           FICH=CHEMIN//'upn.txt'
           JP=18     
           TYP=2
         ENDIF
         IF(PR(I).EQ.'IPEA')THEN
           FICH=CHEMIN//'ipea.txt'
           JP=15 
           TYP=1
         ENDIF
         IF(PR(I).EQ.'IPER')THEN
           FICH=CHEMIN//'iper.txt'
           JP=15 
           TYP=1
         ENDIF
         IF(PR(I).EQ.'HEAA')THEN 
           FICH=CHEMIN//'heaa.txt'
           JP=24    
           TYP=1
         ENDIF
         IF(PR(I).EQ.'UAP')THEN
           FICH=CHEMIN//'uap.txt'
           JP=10
           TYP=2 
         ENDIF
         IF(PR(I).EQ.'IPN')THEN
           FICH=CHEMIN//'ipn.txt'
           JP=21 
           TYP=1 
         ENDIF                 
C ...... CORNIERES A AILES EGALES
         IF(PR(I).EQ.'COR')THEN
           FICH=CHEMIN//'cor.txt'
           JP=16  
           TYP=3 
         ENDIF 
C ...... CORNIERES A AILES EGALES MONTEE EN AILE DE MOULIN
         IF(PR(I).EQ.'CORA')THEN
           FICH=CHEMIN//'cor.txt'
           JP=16     
           TYP=4 
         ENDIF 
C ...... CORNIERES A AILES EGALES MONTEE DOUBLE
         IF(PR(I).EQ.'CORB')THEN
           FICH=CHEMIN//'cor.txt'
           JP=16     
           TYP=5 
         ENDIF 
C ...... CORNIERES A AILES EGALES MONTEE QUADRUPLE
         IF(PR(I).EQ.'CORC')THEN
           FICH=CHEMIN//'cor.txt'
           JP=16     
           TYP=6 
         ENDIF            
 
         IF(PR(I).EQ.'TUBC')TYP=7
         IF(PR(I).EQ.'TUBR')TYP=8
 
C ...... CAS TUBE CYLINDRIQUE OU CYLINDRE
 
         IF(TYP.EQ.7)THEN
           CAR(I,9)=30
           CALL ICYL(NON(I,1),NON(I,2),CAR(I,2),CAR(I,4),CAR(I,6),
     $                                 CAR(I,1),CAR(I,7))
           CAR(I,5)=CAR(I,4)              
           CAR(I,3)=CAR(I,2) 
           GOTO 440                                         
         ENDIF                               
C
C ...... CAS RECTANGLE CREUX OU PLEIN
C
         IF(TYP.EQ.8)THEN
           IF(NON(I,3).EQ.0)CAR(I,9)=40
           IF(NON(I,3).NE.0)CAR(I,9)=50
           CALL IRECT(NON(I,1),NON(I,2),NON(I,3),CAR(I,3),CAR(I,5),
     $                CAR(I,2),CAR(I,4),CAR(I,1),CAR(I,7),CAR(I,6))
           GOTO 440                                         
         ENDIF
C
C     OUVERTURE DU FICHIER ***.TXT DANS OTUA.DIR
C     ET LECTURE DES CARACTERISTIQUES PROFILS
C           
         OPEN(UNIT=40,
     $     FILE='otua\\cor.txt',
     $     STATUS='OLD',
     $     ACCESS='DIRECT',
     $     FORM='FORMATTED',
     $     RECL=12)
 
      WRITE(*,*)COUCOU
         IREC=0   
         DO 200 J = 1,JP
            READ(40,450,REC=1+IREC)A
            IF(A.EQ.IT(I))THEN
              JJJ=14+IREC  
              READ(40,450,REC=JJJ)CAR(I,1)    
              JJJ=18+IREC  
              READ(40,450,REC=JJJ)CAR(I,2)
              JJJ=21+IREC  
              READ(40,450,REC=JJJ)CAR(I,3)  
              JJJ=19+IREC  
              READ(40,450,REC=JJJ)CAR(I,4)
              JJJ=22+IREC  
              READ(40,450,REC=JJJ)CAR(I,5)
              JJJ=26+IREC  
              READ(40,450,REC=JJJ)CAR(I,6) 
C
C... PROFIL SYMETRIQUE OU NON (10: OUI, 20: NON)
C                            
              IF((TYP.EQ.3).OR.(TYP.EQ.4).OR.(TYP.EQ.5).OR.
     $           (TYP.EQ.6))THEN
                CAR(I,9)=20 
                GOTO 1 
              ENDIF
              IF((TYP.EQ.1).OR.(TYP.EQ.2))THEN 
                CAR(I,9)=10     
                JJJ=3+IREC
                READ(40,450,REC=JJJ)LARGEUR 
                JJJ=5+IREC
                READ(40,450,REC=JJJ)EPSEM
                CAR(I,8)=LARGEUR*EPSEM  
                GOTO 1  
              ENDIF    
C                                  
C ... SECTION DES SEMELLES
C
1             IF((TYP.EQ.1).OR.(TYP.EQ.2))THEN
                READ(40,450,REC=5+IREC)CAR(I,10)  
                READ(40,450,REC=4+IREC)CAR(I,11)
                READ(40,450,REC=2+IREC)CAR(I,12)
                READ(40,450,REC=3+IREC)CAR(I,13) 
                CAR(I,12)=CAR(I,12)-2.D00*CAR(I,10)  
              ENDIF  
              IF(TYP.EQ.3)THEN             
                READ(40,450,REC=1+IREC)C1    
                READ(40,450,REC=17+IREC)C2
                C1=(C1/10.D00)-C2            
                CAR(I,7)=DSQRT(C1*C1+C2*C2)  
              ENDIF  
              IF(TYP.EQ.4)THEN 
                C1=CAR(I,1)                
                READ(40,450,REC=17+IREC)C2
                CAR(I,1)=CAR(I,1)*2.D00  
                CAR(I,2)=2.D00*CAR(I,2)+2.D00*C1*C2*C2 
                CAR(I,3)=CAR(I,2)
                READ(40,450,REC=1+IREC)C3
                C3=C3/10.D00  
                CAR(I,4)=CAR(I,2)/C3
                CAR(I,5)=CAR(I,4)
                READ(40,450,REC=2+IREC)B1
                B1=B1/10.D00   
                R1=C3-B1
                CALL TORX(B1,R1,ITOR1)
                CAR(I,6)=ITOR1
                CAR(I,7)=C3    
              ENDIF 
              IF(TYP.EQ.5)THEN 
                C1=CAR(I,1)                         
                CAR(I,1)=CAR(I,1)*2.D00  
                READ(40,450,REC=17+IREC)C2            
                CAR(I,2)=2.D00*CAR(I,2)+2.D00*C1*C2*C2   
                CAR(I,3)=CAR(I,3)*2.D00   
                READ(40,450,REC=1+IREC)C3  
                C3=C3/10.D00  
                CAR(I,4)=CAR(I,2)/C3
                CAR(I,5)=CAR(I,5)*2.D00   
                A1=2.D00*C3
                READ(40,450,REC=2+IREC)B1
                B1=B1/10.D00
                D1=2.D00*B1
                C1=C3-B1
                READ(40,450,REC=3+IREC)R1
                R1=R1/10.D00
                CALL TORT(A1,B1,C1,D1,R1,ITOR1)   
                CAR(I,6)=ITOR1
                CAR(I,7)=DSQRT(C3*C3+C2*C2)  
              ENDIF  
              IF(TYP.EQ.6)THEN             
                C1=CAR(I,1)                
                CAR(I,1)=CAR(I,1)*4.D00 
                READ(40,450,REC=17+IREC)C2  
                CAR(I,2)=4.D00*CAR(I,2)+4.D00*C1*C2*C2                
                CAR(I,3)=CAR(I,2) 
                READ(40,450,REC=1+IREC)C3
                C3=C3/10.D00        
                CAR(I,4)=CAR(I,2)/C3
                CAR(I,5)=CAR(I,4)    
                READ(40,450,REC=2+IREC)B1
                B1=B1/10.D00   
                R1=C3-B1
                B1=B1*2.D00
                CALL TORX(B1,R1,ITOR1)
                CAR(I,6)=ITOR1
                CAR(I,7)=C3      
              ENDIF  
            ENDIF 
            IREC=IREC+30      
200      CONTINUE
         CLOSE(40)       
440   CONTINUE   
450   FORMAT(F11.0)
      RETURN
      END

En tout cas, merci à tous ceux qui jeterront un oeil, même de loin!

[Sylvain Bergeron]

Simple problème de compréhension du comportement des chaînes en Fortran.

Si Fich et Chemin sont déclarées ainsi :

Code :

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

CHARACTER*10 FICH, CHEMIN

alors

Code :

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

1
2
Chemin = 'Hello'
Fich = Chemin // ' World'

retournera toujours le contenu de Chemin seulement.

Pourquoi ?

Parce qu'une variable caractère Fortran a une longueur égale à sa longueur déclarée et non égale à la longueur de ce qui lui est assigné. Lors de l'assignation, la chaîne est allongée par des blancs ou tronquée pour avoir exactement la longueur déclarée de la chaîne.

Dans l'exemple,

Code :

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

1
2
Chemin = 'Hello' ! Chemin contient 'Hello     '
Fich = Chemin // ' World' ! assigne 'Hello     ', soit 'Hello      World' tronqué à 10 carcatères

Alors comment faire ?

En F90 : Fich = trim(Chemin) // ' World'
En F77 :
k = len_trim(Chemin)
Fich = Chemin(1:k) // ' World'

si len_trim() n'existe pas, elle est facile à créer...

[Rémirem]

Merci pour cette réponse très clair!!
Les choses parraissent simple quand elles sont aussi bien expliqué!