Du DEBUGAGE, YOUPI !

Bonjour à tous,

J'essaye de faire des subroutines afin de généraliser ma procédure pour lire et écrire avec des fichiers et j'ai un ptit bug ...
le programme test se présente comme ça :

Code:

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
PROGRAM TESTLECTURE
 
IMPLICIT NONE
 
TYPE FICHIER
 CHARACTER(LEN=4) 				:: nom
 INTEGER 					:: unite	
END TYPE
INTEGER 					:: n	
DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE	:: t
TYPE (fichier) 					:: sat1,sat2
 
sat1%nom	= "sat1"
sat2%nom	= "sat2"
sat1%unite	= 15
sat2%unite	= 16
 
CALL ouverturecompte(sat1%nom,sat1%unite,n)
ALLOCATE(t(n))
write(*,*) n
CALL lecture(sat1%nom,sat1%unite,n,t)
CALL ecriture(sat2%nom,sat2%unite,n,t)
 
 
CONTAINS 
 
SUBROUTINE OUVERTURECOMPTE(nomfichier,unitee,nlc)
	CHARACTER(LEN=*),INTENT(IN) 			:: nomfichier
	INTEGER 	,INTENT(IN) 			:: unitee
	INTEGER		,INTENT(OUT)			:: nlc
	INTEGER		 				:: ios=0, ios2
 	open(unit=unitee,file = nomfichier,form='formatted',status='old',action='read',iostat=ios)
	if (ios .ne. 0) stop "Pb avec comptage"
 	ios2 = 0
  	nlc = 0
  	do while(ios2 .eq. 0)
      		read(unit=unitee,FMT=*,iostat=ios2)
      		nlc = nlc+1
   	enddo
 	close(unit=unitee)
	nlc=nlc-1
END SUBROUTINE OUVERTURECOMPTE
 
 
SUBROUTINE LECTURE(nomfichier,unitee,nlc,tcheck)
	CHARACTER(LEN=*) ,INTENT(IN)			:: nomfichier
	INTEGER 	 ,INTENT(IN) 			:: unitee
	INTEGER		 ,INTENT(IN) 			:: nlc
	INTEGER						:: k
	DOUBLE PRECISION, DIMENSION(nlc),INTENT(OUT) 	:: tcheck
	INTEGER		 				:: ios=0
	write(*,*)unitee
	open(unit=unitee,file = nomfichier,form='formatted',status='old',action='read',iostat=ios)
	if (ios .ne. 0) stop "Pb avec lecture!"
 	do k=1,nlc
 		read(unit=unitee,FMT=*) tcheck(k)
  	enddo
	close(unit=unitee)
	RETURN
END SUBROUTINE LECTURE
 
 
SUBROUTINE ECRITURE(nomfichier,unitee,nlc,tcheck)
	CHARACTER(LEN=*) ,INTENT(IN)			:: nomfichier
	INTEGER 	 ,INTENT(IN) 			:: unitee
	INTEGER		 ,INTENT(IN) 			:: nlc
	INTEGER						:: k
	DOUBLE PRECISION, DIMENSION(nlc),INTENT(OUT) 	:: tcheck
	INTEGER		 				:: ios=0
	open(unit=unitee,file=nomfichier,form='formatted',status='replace',action='readwrite' ,iostat=ios)
		if (ios .ne. 0) stop "Pb ecriture"
 	do k=1,nlc
   		write(unit=unitee,FMT='(e18.11)') tcheck(k)
	enddo 	
	close(unit=unitee)
	RETURN
END SUBROUTINE ECRITURE
 
 
END PROGRAM TESTLECTURE

J'ai une erreur :

Code:

1
2
3
 
At line 56 of file testlecture.f90 (unit = 15, file = 'sat1')
Fortran runtime error: End of file

Bon ça doit être mon FMT = * qui pose pb mais je sais pas comment le changer ...
Un ptit coups de pouce svp ?

Merci !

Salut.
Le nombre de lignes du fichier lu par la subroutine OUVERTURECOMPTE n'est pas correct. Tu peux la corriger en imposant un format de lecture (apparemment il y a un nombre réel par ligne, si le gabarit est connu), et en affectant la valeur lue à une variable.
J'ai observé qu'une ligne supplémentaire était comptée dans ce cas, lorsque la dernière ligne de données était terminée par un retour à la ligne (donc une ligne vide en plus).

Hello,
Il y avait bien un retour à la ligne en trop .. grr !
J'aime pas donner des formats de lecture et d'écriture.
Ce que je fais, c'est que je déclare tout en double précision, et j'ajoute des options de compilations pour que le compilateur prenne bien cette double précision en compte (et oui faut lui dire 2x à Gfortran!) ...

Bonjour,

Citation:

Envoyé par DonDouego

Ce que je fais, c'est que je déclare tout en double précision, et j'ajoute des options de compilations pour que le compilateur prenne bien cette double précision en compte (et oui faut lui dire 2x à Gfortran!) ...

Non, il suffit de le préciser dans le code ou à la compilation, mais pas besoin des deux, normalement.
Tu as un exemple minimal sous la main qui illustre ce que tu constates?