How to create functions in Fortran?

Question

I'm sure the solution to this is extremely basic, but I'm having a hard time figuring out how to use functions in Fortran. I have the following simple program:

  PROGRAM main
    IMPLICIT NONE
    INTEGER :: a,b
    a = 3
    b = 5
    PRINT *,funct(a,b)
  END PROGRAM

  FUNCTION funct(a,b)
    IMPLICIT NONE
    INTEGER :: funct
    INTEGER :: a,b

    funct = a + b
  END FUNCTION

I've tried several variations of this, including assigning a data type before FUNCTION, assigning the result of funct to another variable in the main program and printing that variable, and moving the FUNCTION block above the PROGRAM block. None of these worked. With the current program I get an error on line 6 (the line with the PRINT statement):

Error: Return type mismatch of function 'funct' (UNKNOWN/INTEGER(4))
Error: Function 'funct' has no IMPLICIT type

From all of the guides I've tried, I seem to be doing it right; at least one of the variations, or a combination of some of them, should have worked. How do I need to change this code to use the function?

the problem is that you used implicit none while you did not give funct a declaration. you must add integer ,external :: funct. — M. Chinoune
– M. Chinoune, Commented Aug 29, 2017 at 8:20

Pierre de Buyl · Accepted Answer · 2017-08-29 07:03:15Z

5

Simply putting the function in the file will not make it accessible to the main program.

Traditionally, you could simply declare a function as external and the compiler would simply expect to find a suitable declaration at compile-time.

Modern Fortran organizes code and data in "modules". For your purpose, however, it is simpler to "contain" the function within the scope of the main program as follows:

PROGRAM main
  IMPLICIT NONE
  INTEGER :: a,b
  a = 3
  b = 5
  PRINT *,funct(a,b)

CONTAINS

  FUNCTION funct(a,b)
    IMPLICIT NONE
    INTEGER :: funct
    INTEGER :: a,b

    funct = a + b
  END FUNCTION funct
END PROGRAM main

answered Aug 29, 2017 at 7:03

Pierre de Buyl

7,3232 gold badges18 silver badges23 bronze badges

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

Comments

veryreverie · Accepted Answer · 2021-03-25 10:06:41Z

0

A simpler solution can be the following code

 PROGRAM main
    IMPLICIT NONE
    INTEGER :: a,b, funct
    a = 3
    b = 5
    PRINT *,funct(a,b)
  END PROGRAM

  FUNCTION funct(a,b)
    IMPLICIT NONE
    INTEGER :: funct
    INTEGER :: a,b

    funct = a + b
  END FUNCTION

where the only difference is in the third line, where I have declared funct as an integer. It compiles and it prints 8 as result.

edited Mar 25, 2021 at 10:06

veryreverie

3,0012 gold badges16 silver badges26 bronze badges

answered Mar 25, 2021 at 9:32

Jonny_92

1552 silver badges15 bronze badges

2 Comments

francescalus Over a year ago

You've (correctly) declared funct's return type as integer not real. For general circumstances, though, this is quite bad advice.

Jonny_92 Over a year ago

Thanks a lot, I am quite new to Fortran.

Collectives™ on Stack Overflow

How to create functions in Fortran?

2 Answers 2

Comments

2 Comments

Your Answer

Linked

Hot Network Questions

Collectives™ on Stack Overflow

2 Answers 2

Comments

2 Comments

Your Answer

Sign up or log in

Post as a guest

Linked

Related