How to handle Fortran global allocatable variables in a module across subroutines

Question

I have the following module with an allocatable variable which is defined in the module, allocated in a subroutine, and then also used in a second subroutine called by the first subroutine. In this situation do I have to pass the variable to the second subroutine and declare INTENT(inout)? Or since it's a global variable it doesn't need to be passed as an argument?

MODULE test

  IMPLICIT NONE
  SAVE

  REAL,ALLOCATABLE,DIMENSION(:,:,:) :: total

CONTAINS

  !--- 1st subroutine
  SUBROUTINE my_subr1(n,m,z)
    IMPLICIT NONE
    INTEGER,INTENT(in) :: n,m,z
    ALLOCATE(total (n,m,z))
    total=.9
    CALL my_subr2(n)

  END SUBROUTINE my_subr1

  !-- 2nd subroutine
  SUBROUTINE my_subr2(n)
    IMPLICIT NONE
    INTEGER,INTENT(in) :: n

    total(n,:,:)=total(n-1,:,:)
  END SUBROUTINE my_subr2
END MODULE test

All subroutine within the module have access to total. Module variables are "save" by default, also. Your code hass issues though (no declaration of dummy arguments in my_subr1 and no end module .... — Pierre de Buyl
– Pierre de Buyl, Commented Feb 25, 2019 at 9:36

Rodrigo Rodrigues · Accepted Answer · 2019-02-25 09:59:17Z

5

do I have to pass the variable to the second subroutine and declare INTENT(inout)?

No, you don't. Any variable decalred in the body of the module has the save attribute by default. You must, though, ensure that the second subroutine is only called after the first was executed, or else the program will fail because total would not be initialized yet.

All functions and subroutines declared in the module will have access to total by host association.

By the way, there are some issues you should address in your code, as mentioned by @PierredeBuyl in the comments:

Variables declared in the module body are saved by default; you should remove the SAVE statement.
Procedures declared in a module inherit the IMPLICIT directive from the module scope, there is no need to redeclare it in the subroutine if you won't change it.
You are missing the declaration of the arguments in my_subr1.

edited Feb 25, 2019 at 9:59

answered Feb 25, 2019 at 9:51

Rodrigo Rodrigues

8,6321 gold badge32 silver badges43 bronze badges

Sign up to request clarification or add additional context in comments.

2 Comments

francescalus Over a year ago

The save attribute isn't terribly relevant here. Even under Fortran standards where there isn't that attribute by default, the call from my_subr1 to my_subr2 means the variable doesn't go out of scope between the two.

Herman Toothrot Over a year ago

thanks for the explanation, I have read about SAVE and IMPLICIT

Collectives™ on Stack Overflow

How to handle Fortran global allocatable variables in a module across subroutines

1 Answer 1

2 Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

2 Comments

Your Answer

Sign up or log in

Post as a guest

Related