Discussion :

[mari2]

bonjour, j'ai un fichier que je doit le lire par fortran , mais j'arrive pas
quelqu'un peut m'aider s'il vous plaît!
je serai ravi
le fichier est de forme :

$MeshFormat
2.2 0 8
$EndMeshFormat
$Nodes
13
1 0 0 0
2 0 1 0
3 1 1 0
4 1 0 0
5 0 0.5 0
6 0.5 0.5 0
7 0.5 0.5 0
8 0.5 0.5 0
9 0.75 0.5 0
10 0.75 0.75 0
11 0.5 0.75 0
12 0.5 0.5 0
13 0.625 0.625 0
$EndNodes
$Elements
27
1 15 2 0 1 1
2 15 2 0 2 2
3 15 2 0 3 3
4 15 2 0 4 4
5 15 2 0 5 5
6 15 2 0 6 6
7 15 2 0 7 7
8 15 2 0 8 8
9 15 2 0 9 9
10 15 2 0 10 10
11 15 2 0 11 11
12 1 2 0 1 1 4
13 1 2 0 2 3 4
14 1 2 0 3 3 2
15 1 2 0 4 1 2
16 1 2 0 8 11 10
17 1 2 0 9 10 9
18 1 2 0 10 9 6
19 1 2 0 11 6 11
20 2 2 0 12 1 12 2
21 2 2 0 12 1 4 12
22 2 2 0 12 2 12 3
23 2 2 0 12 3 12 4
24 2 2 0 14 6 13 9
25 2 2 0 14 6 11 13
26 2 2 0 14 9 13 10
27 2 2 0 14 10 13 11
$EndElements

[Sylvain Bergeron]

Tu donnes très peu d'information...

Voici la stratégie que j'utiliserais :

• Je ferais 3 subroutines, pour lire chacune des sections
• Dans chaque subroutine, je testerais que les mots clés de début et de fin sont bien rencontrés, j'allouerais les tableaux s'il y a lieu, et je lirais les données.

J'aurais donc quelque chose comme :

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
 
...
open(11,...)
call LireMeshFormat(11,...)
call LireNodes(11,...)
call LireElements(11,...)
...
 
subroutine LireMeshFormat(NumFic,...)
   ...
   character (20) :: MotCle
   read(NumFic,'(a)') MotCle
   if (MotCle /= "$MeshFormat") --> erreur !
   read(NumFic, *) données de la ligne
   read(NumFic,'(a)') MotCle
   if (MotCle /= "$EndMeshFormat") --> erreur !
   end
 
subroutine LireNodes(NumFic,...)
   ...
   character (20) :: MotCle
   integer           :: NbrNodes
   integer           :: k
 
   read(NumFic,'(a)') MotCle
   if (MotCle /= "$Nodes") --> erreur !
   read(NumFic, *) NbrNodes
   allocate(...(NbrNodes))
   do k = 1, NbrNodes
       read(NumFic, *) données (k) de la ligne
   enddo
   read(NumFic,'(a)') MotCle
   if (MotCle /= "$EndNodes") --> erreur !
   end
 
subroutine LireElements(NumFic,...)
...

[wafiflicha]

Bonjour,
je pourrais répondre à votre question par l'exemple décrit en dessous, et à vous de généraliser pour le reste de votre 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
25
26
27
28
29
30
$Nodes
13
1 0 0 0
2 0 1 0
3 1 1 0
4 1 0 0
5 0 0.5 0
6 0.5 0.5 0
7 0.5 0.5 0
8 0.5 0.5 0
9 0.75 0.5 0
10 0.75 0.75 0
11 0.5 0.75 0
12 0.5 0.5 0
13 0.625 0.625 0
 
 
character (20) :: node
integer :: nbrnode, i
real, allocate(:,:), table :: n
 
 
 
! ouvrir votre fichier
open (1,'nom_de_fichier' )
read(1,'(a)') node
read(1,*) nbrnode
allocate (n(nbrnode,4))
read(1,*) (n(i,1), n(i,2), n(i,3), n(i,4), i=1,nbrnode)
close(1)

et voilà, tu vas enregistrer vos variables dans un tableau n.