Discussion :

[wasabi1024]

Bonjour,

Malgré des recherches sur plusieurs forums je n'ai pas pu trouver la réponse à mon problème. J'espère la trouver ici

Mon problème est le suivant : j'essaie de passer en argument a une subroutine un vecteur allocatable mais ça ne marche pas. Voici un bout de code ou le problème apparait.

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
 
      MODULE allocvars
      implicit none
      real,dimension(:,:),allocatable :: field
      END MODULE
 
      PROGRAM main
      use allocvars
      implicit none
      allocate(field(5,5))
      field=4.0
      write(*,*) 'field in main : ',field(2,2)
      call sub(field)
      END
 
      SUBROUTINE sub(fieldbyreference)
      use allocvars
      implicit none
      real,dimension(:,:),allocatable :: fieldbyreference
      write(*,*) 'field in sub (from module):',field(2,2)
      write(*,*) 'field in sub (from reference):',fieldbyreference(2,2)
      END

Et voici le résultat :

Code :

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

1
2
3
4
5
6
 
[wasabi1024@computer fortranpb]$ gfortran main.f 
[wasabi1024@computer fortranpb]$ ./a.out 
 field in main :    4.0000000    
 field in sub (from module):   4.0000000    
Segmentation fault

Quelqu'un comprend pourquoi le passage par référence ne fonctionne pas?

[Marlan]

Salut,

Tu n'as pas spéciafié le type entrée/sortie de ton champ fieldbyreference.
Qui plus est, il n'est pas nécessaire de préciser le status ALLOCATABLE sur un tableau en entrée de subroutine.
Le code corrigé serait :

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
      MODULE allocvars
      implicit none
      real,dimension(:,:),allocatable :: field
      END MODULE

      PROGRAM main
      use allocvars
      implicit none
     INTERFACE
               SUBROUTINE sub(fieldbyreference)
               use allocvars
               implicit none
               real,dimension(:,:),INTENT(IN) :: fieldbyreference
               END SUBROUTINE sub
      END INTERFACE
      allocate(field(5,5))
      field=4.0
      write(*,*) 'field in main : ',field(2,2)
      call sub(field)
      END

      SUBROUTINE sub(fieldbyreference)
      use allocvars
      implicit none
      real,dimension(:,:),INTENT(IN) :: fieldbyreference
      write(*,*) 'field in sub (from module):',field(2,2)
      write(*,*) 'field in sub (from reference):',fieldbyreference(2,2)
      END

[Sylvain Bergeron]

La version Marlan ne va pas marcher non plus : il faut une interface explicite pour un « dummy argument » qui est de type « assumed shape », (real,dimension(:,:),INTENT(IN) :: fieldbyreference)

Tu peux faire une recherche de « interface explicite » dans le forum pour plus d'info.

[Marlan]

A, oui, bien vu. Désolé de mon erreur, je corrige le post.