Discussion :

[forcetran90]

Bonjour,

Je viens requérir votre aide car dans le cadre de mon cursus scolaire, on me demande de réaliser un programme en FORTRAN90. Une bonne partie de mon programme est écrit mais j'ai un problème. En fait mon programme est constitué du dit programme et d'un sous programme. Ma question est de savoir comment utiliser un tableau crée dans mon programme principal dans mon sous programme?

Merci d'avance de vos réponses.

[Ladgalen]

Bonjour,

En fortran chaque programme ou sous programme est indépendant, donc tu dois faire communiquer tes sous programmes en passant les variables de l'un à l'autre. Tu as plusieurs possibilités. Dans ton cas (1 seul sous programme et 1 tableau à transférer) l'utilisation d'arguments à ta subroutine semble le plus simple.

à voir dans la FAQ : ici

[forcetran90]

Tout d'abord, merci d'avance de votre réponse.

En fait je n'ais pas vraiment bien décris mon problème. Je vais donc essayer de l'expliciter en vous écrivant une partie de mon code. De plus je précise que lors de la compilation je n'ai aucun problème mais lorsque je lance mon programme après avoir entré mes paramètres j ai une " ERREUR SEGMENTATION " qui s'affiche. Je me suis renseigné et il semblerais que ce soit du au fait que mon tableau accède à une partie de la mémoire qui ne lui ait pas autorisé....

voici une partie de mon code:

PROGRAM PG

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
IMPLICIT NONE
INTEGER(KIND=8)   :: i, n, j
REAL(KIND=8), ALLOCATABLE, DIMENSION( : )   :: x, y, m
REAL(KIND=8) :: tmax, t
READ*, tmax, t
READ*, n
ALLOCATE(x(1:n), y(1:n), m(1:n))
DO i=1, n
READ*, x(i), y(i), m(i)
END DO
DO WHILE(t<=tmax)
DO j=1, n
call TESTE( x(j), y(j), vx(j), vy(j), m(j), N, t, j)
END DO
t=t+h
PRINT*,'&'
END DO
DEALLOCATE(x, y, m)
END PG
 
SUBROUTINE TESTE(x, y, m, n, t, j)
IMPLICIT NONE
INTEGER(KIND=8)   :: i, n, j,
REAL(KIND=8)   :: t
REAL(KIND=8), DIMENSION( : )   :: x, y, m
REAL(KIND=8), ALLOCATABLE, DIMENSION( : )   :: rx, ry
REAL(KIND=8), ALLOCATABLE, DIMENSION( : )   :: r
 
 
ALLOCATE(rx(1:n), ry(1:n))
DO i=1, n
    rx(i)=x(i)-x(j)
    ry(i)=y(i)-y(j)
 
END DO
END SUBROUTINE TESTE

Merci d'avance de vos réponses.

NB: J'ai supprimé une bonne partie de mon programme qui n'est pas en cause dans mon problème.

[Ladgalen]

Bonjour

Ton code semble pas mal (pense à l'indenter c'est plus facile à lire). Par contre il y a une erreur lié au tableau qui doit être l'origine du problème. Tu ne peux pas déclarer des tableaux allocatable dans une subroutine (en fait oui j'expliquerai plus bas). Dans ton programme principal tu alloue tes tableaux et dans les sous programmes il faut qu'il soit statique. Donc tu dois transférer à ta subroutine les dimensions des tableaux et tu déclare alors les tableaux avec la bonne dimension. Par exemple:

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
 
program exemple
  implicit none
 
  integer :: dim, i
  real(kind=8), dimension(:), allocatable :: tab      ! dynamique
 
  write(*,"('dim ?')")
  read(*,*) dim
  write(*,"('dim = ', i12)") dim
 
  allocate( tab(1:dim) )
 
  do i=1,dim
    tab(i) = dble(i**2)
  end do
 
  call print( dim, tab)
 
end program exemple
 
!* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
 
subroutine print( dim, tab)
  implicit none
 
  ! arguments :
  integer, intent(in) :: dim            ! intent(in) est nécessaire pour utiliser dim
                                        ! comme dimension
  real(kind=8), dimension(dim), intent(in) :: tab       ! statique
 
  ! local :
  integer :: i
 
  do i=1,dim
    write(*,*)i,tab(i)
  end do
 
end subroutine print

Il est très utile de précisé dans la déclaration des arguments s'il s'agit d'arguement d'entré intent(in), de sortie intent(ou) ou les deux (inout). Le compilateur peut vérifier ce que tu fais avec tes arguments et ça simplifie la lecture. Ici, il faut dire que dim est un argument d'entré (donc nom modifiable) pour avoir le droit de l'utiliser pour déclarer la dimension.

Pour faire un allocate dans une subroutine il faut utiliser les modules (F90 et suivant uniquement). Si un tableau est déclaré dans un module alors il peut être alloué dans une subroutine qui est contenu dans ce module.

[forcetran90]

Merci beaucoup, j'ai comme vous me l'avez dit donné les dimensions de mon tableau dans ma subroutine, et depuis plus aucun problème.

Encore merci car j'ai passé énormément de temps à rechercher comment résoudre cette erreur!!!