Syntaxe Vieux FORTRAN

Bonjour A Tous

J'ai un probleme avec un "vieux" fichier FORTRAN (dévellopé sur une plateforme VAX avec VAX/FORTRAN).....

A priori, je suis parvenu a le convertir pour le compilateur gfortran...

mais il me reste une ligne de commande qui est totalement incomprehensible pour moi, (qui suis totalement novice en FORTRAN) et pour le compilateur ....
Code:

1 2 3 4 5 6 7 8 9 10 ! WRITE(3,'(A8,I8,8X,6I8,''+'',<int((ndof1*2-6)/8.)> ! (/''+ '',8I8,''+''), ! /''+ '',<ndof1*2-6-int((ndof1*2-6)/8.)*8>I8, ! T73,''+'')')'GENEL ',iel,((dof1(j,k),j=1,2),k=1,ndof1) ! WRITE(3,'(''TEST'')') ! PRINT*,A8,I8,8X,6I8,((int((ndof1*2-6)/8.)) ) ! &(/''+ '',8I8,''+''),& ! &/''+ '',(ndof1*2-6-int((ndof1*2-6)/8.)*8)I8,& ! &T73,''+'')')'GENEL ',iel,((dof1(j,k),j=1,2),k=1,ndof1)
Je pense que le probleme est juste une erreur de syntaxe, mais si quelq'un pouvai me donner un petit coup de main, ....

Merci par avance

Le problème vient d'une extension (toujours présente d'ailleurs dans le compilateur Intel). Les bouts "<expression>" dans les formats sont remplacés par leur valeur numérique à l'exécution. Pour contourner le problème, tu dois construire le format dans une variable caractère, puis utiliser le format construit pour écrire.

Le premier write :
Code:

1 2 3 4 WRITE(3,'(A8,I8,8X,6I8,''+'',<int((ndof1*2-6)/8.)> (/''+ '',8I8,''+''), /''+ '',<ndof1*2-6-int((ndof1*2-6)/8.)*8>I8, T73,''+'')')'GENEL ',iel,((dof1(j,k),j=1,2),k=1,ndof1)
devient donc :
Code:

1 2 3 4 5 6 7 8 9 10 11 12 character (256) :: MonFmt ... ! Construction du format par un write interne write(MonFmt,'(a,i0,a,i0,a)') '(A8,I8,8X,6I8,''+''', & int((ndof1*2-6)/8.), & ' (/''+ '',8I8,''+''),/''+ '',', & ndof1*2-6-int((ndof1*2-6)/8.)*8, & 'I8,T73,''+'')' ! Utilisation du format WRITE(3,MonFmt) 'GENEL ',iel,((dof1(j,k),j=1,2),k=1,ndof1)

Merci beaucoup pour la réponse,

a priori, le compilateur a bien "digérer" la solution, par contre j'ai essaye de reproduire la solution sur la deuxième (et dernière) ligne de code qui me résiste, ... mais sans succès,

voici la ligne de code initiale
Code:

1 2 3 4 5 6 ! WRITE(3,'(''+* K '',3E16.6,''+'',/, & ! &<int((nodf1**2/2.+ndof1/2.-3)/4.)>(''+* '',4E16.6,''+'',/), & ! &''+* '',(<int(ndof1**2/2.+ndof1/2.-3-int((ndof1**2/2.+ & ! &ndof1/2.-3)/4.)*4)>E16.6)') & ! &(mat1(k),k=1,int(ndof1**2/2.+ndof1/2.))
et ma tentative de traduction ....
Code:

1 2 3 4 5 6 7 8 9 10 11 character (3000) :: MonFmt2 ... ! Mon deuxieme formatage write (MonFmt2,'a,i0,a,a,i0') '(''+* K '',3E16.6,''+'',/,& &',int((nodf1**2/2.+ndof1/2.-3)/4.),'& &(''+* '',4E16.6,''+'',/), & &''+* '',(',int(ndof1**2/2.+ndof1/2.-3-int((ndof1**2/2.+ & &ndof1/2.-3)/4.)*4),'E16.6)' ! a utiliser avec : WRITE(3,MonFmt2) (mat1(k),k=1,int(ndof1**2/2.+ndof1/2.))
Mais à priori ... cela ne plais pas a mon deuxième compilateur préferé:cry:

A mon avis l'erreur est dans (MonFmt2,'a,i0,a,a,i0') je pense qu'il faut mettre (MonFmt2,'a,i0,a,i0,a')

Je n'ai pas vérifié et je n'ai pas le temps, mais à coup sûr, le format doit être entre parenthèses : (MonFmt2,'(a,i0,a,i0,a)')

La même règle s'applique pour le format construit : il doit contenir des parenthèses au début et à la fin.

Je progresse ....
c'étais bien un probleme de parenthese, le compilateur accepte bien la modification.
Pour memoire le code source
Code:

1 2 3 4 5 6 7 8 9 10 ! Mon deuxieme formatage write (MonFmt2,'(a,i0,a,i0,a)') '(''+* K '',3E16.6,''+'',/,& &',int((nodf1**2/2.+ndof1/2.-3)/4.),'& &(''+* '',4E16.6,''+'',/), & &''+* '',(',int(ndof1**2/2.+ndof1/2.-3-int((ndof1**2/2.+ & &ndof1/2.-3)/4.)*4),'E16.6)' ! a utiliser avec : WRITE(3,MonFmt2) (mat1(k),k=1,int(ndof1**2/2.+ndof1/2.)
par contre je n'ai pas du bien traduire le code, lors de l'execution de mon ".exe", j'ai le probleme suivant :
Code:

1 2 3 4 5 6 7 At line 164 of file genelm.f90 (Unit 3 "Ferrure_A-8_2D_Geo.bdf.GEN") Traceback: not available, compile with -ftrace=frame or -ftrace=full Fortran runtime error: Expected P edit descriptor in format ('+* K ',3E16.6,'+',/,-244668376('+* ',4E16.6,'+',/), '+* ',(2E16.6) ^
j'ai essayer de compiler avec l'option specifié, mais pas de solution ....
pour info, la ligne 164 correspond a la commande write....
Quelq"un a une idee ??? et encore merci beaucoup pour toutes vos reponses

Le problème est le facteur de répétition négatif dans le format (le -244668376). Ce nombre est le résultat de l'expression int((nodf1**2/2.+ndof1/2.-3)/4.).

Désolé, mais je perd pieds ....

je ne comprend pas la dernière réponse , (je suis une truffe en Fortran, je le rappel :? )et cela même si j'ai récupéré et parcouru les docs VAX/FORTRAN a propos du sujet qui me préoccupe " Variable Format Expression" :aie: cela ne m'a pas aidé....

A priori, j'ai tendance à pensée que le problème viens de mon deuxieme formatage temporaire,en effet, le code que j'ai recopié (avec précision...) à été publier il y a plus de 15 ans, mais pour ce problème, je suis sec....

Pour info, la sortie de cette commande Write doit ressemblé à quelque chose comme ça
Code:

1 2 3 4 5 6 7 8 9 10 11 12 * K 3.249404E-05 2.742304E-07 1.260558E-05 * -2.535259E-08 7.395782E-07 3.715230E-09 9.747304E-05 * 2.298511E-06 -1.882033E-06 -2.254303E-08 7.378169E-07 * 1.121394E-09 5.761329E-05 1.070783E-06 2.997293E-05 * -2.436260E-08 7.386966E-07 2.811055E-09 1.327319E-05 * 2.142209E-07 -2.820588E-07 1.044184E-08 -7.886342E-08 * 4.123580E-06 3.650244E-05 -1.044423E-05 -2.805073E-07 * 9.661857E-09 -7.999187E-08 1.765229E-06 2.473880E-05 * -3.597681E-06 -2.814634E-07 1.007685E-08 -7.925026E-08 * 4.237562E-05 -5.551566E-08 5.366323E-07 5.386769E-09 ....
Si quelqu'un a une idée ....

De toute évidence, je dois mieux expliquer le problème...

Ton programme arrête en erreur sur :
Code:

1 2 Fortran runtime error: Expected P edit descriptor in format ('+* K ',3E16.6,'+',/,-244668376('+* ',4E16.6,'+',/), '+* ',(2E16.6)
Le message est bizarre a priori, mais un bout de format est affiché. Si tu observe attentivement ce format, tu remarque que c'est le résultat de :
Code:

1 2 3 4 5 write (MonFmt2,'(a,i0,a,i0,a)') '(''+* K '',3E16.6,''+'',/,& &',int((nodf1**2/2.+ndof1/2.-3)/4.),'& &(''+* '',4E16.6,''+'',/), & &''+* '',(',int(ndof1**2/2.+ndof1/2.-3-int((ndof1**2/2.+ & &ndof1/2.-3)/4.)*4),'E16.6)'
dans lequel l'expression int((nodf1**2/2.+ndof1/2.-3)/4.) a été remplacée par -244668376.

Le problème, c'est que le format devrait contenir un nombre positif indiquant le nombre de fois à répéter la section "('+* ',4E16.6,'+',/)", pas un nombre négatif (le message relatif au P edit descriptor vient du fait que le masque P est le seul à supporter un nombre négatif). On peut donc conclure que l'expression int((nodf1**2/2.+ndof1/2.-3)/4.) est mal interprétée.

Je suggère donc d'ajouter :

Code:

write(*,*) 'ndof1=',ndof1,'; int((nodf1**2/2.+ndof1/2.-3)/4.)',int((nodf1**2/2.+ndof1/2.-3)/4.)

avant le write(MonFmt,... pour voir ce qui se passe.

Bonjour,
j'ai ajouté le code specifié, voici le resultat....

ndof1= 54 ; int((nodf1**2/2.+ndof1/2.-3)/4.) 134304605

Nota : pour info, ndof1 =54 corespond bien à la taille de la matrice que j'essaye d'ecrire (diagonale inferieur), par contre j'ai demandé la visualisation de nodf1**2, qui sauf erreur :mrgreen: de ma part correspond a 54² ... et j'ai le resultat suivant :ndof1**2 = 1074436793 ????:calim2:

A priori, l'erreur venait d'ne erreur de Syntaxe dans mes variable ndof1 ecrit a la place de ndof1..:oops:, une fois corrrigé, j'ai le message suivant .....

ndof1= 54 ; int((nodf1**2/2.+ndof1/2.-3)/4.) 370

At line 169 of file genelm.f90 (Unit 3 "Ferrure_A-8_2D_Geo.bdf.GEN")
Traceback: not available, compile with -ftrace=frame or -ftrace=full
Fortran runtime error: Unexpected end of format string
^

je rappel ci dessous le code,
la ligne 169 correspont à la ligne de code :

WRITE(3,MonFmt2) (mat1(k),k=1,int(ndof1**2/2.+ndof1/2.))
Code:

1 2 3 4 5 6 7 8 9 10 11 12 write (MonFmt2,'(a,i0,a,i0,a)') '(''+* K '',3E16.6,''+'',/,& &',int((ndof1**2/2.+ndof1/2.-3)/4.),'& &(''+* '',4E16.6,''+'',/), & &''+* '',(',int(ndof1**2/2.+ndof1/2.-3-int((ndof1**2/2.+ & &ndof1/2.-3)/4.)*4),'E16.6)' write(*,*) 'ndof1=',ndof1,'; int((ndof1**2/2.+ndof1/2.-3)/4.)',int((ndof1**2/2.+ndof1/2.-3)/4.) write(*,*) ndof1**2 ! a utiliser avec : WRITE(3,MonFmt2) (mat1(k),k=1,int(ndof1**2/2.+ndof1/2.))
et encore merci pour votre aide ...

Ajoute une ligne :

write (*,'(a)') trim(MonFmt2)

juste avant la ligne fatale. Tu pourra ainsi voir en quoi le format est incomplet...

YEEESSSSS !!!!

Merci encore beaucoup pour votre aide !!!!!

j'avais en effet une parenthèse en trop dans mon code...
:ccool: