Many people perceive Fortran as an archaic and "dead" programming language. However, most scientific and engineering code is written in Fortran. As such, programming in F77 and F90 remains a necessary skill for most technical programmers. Moreover, the latest Fortran standards (2003, 2008, 2015) allow the programmer to write highly efficient code with minimum effort, while utilizing all of the modern language features, such as OOP (object-oriented programming). FORTRAN is an acronym for "FORmula TRANslation", and is best suited for mathematical and numerical applications rather than graphics or database applications. Most fortran codes take text input from a file or command-line rather than from a menu or GUI interface.
Writing and Compiling a Simple Program
Write a “Hello World” program.
This is usually the first program to write in any language, and it just prints “Hello world” to the screen. Write the following code in any text editor and save it as helloworld.f. Pay attention that there must be exactly 6 spaces in front of every line. program helloworld implicit none character*13 hello_string hello_string = “Hello, world!” write (*,*) hello_string end program helloworld
Compile the program.
To do this, type f77 helloworld.f into the command line. If this gives an error, you probably haven’t installed a Fortran compiler like for example gfortran yet.
Run your program.
The compiler has produced a file called a.out. Run this file by typing ./a.out.
Understand what you just wrote.
- program helloworld indicates the start of the program “helloworld”. Similarly, end program helloworld indicates its end.
- By default, if you don’t declare a variable type, Fortran treats a variable with a name that begins with a letter from i to n as integer, and all others as a real number. It is recommended to use implicit none if you don’t need that behaviour.
- character*13 hello_string declares an array of characters which is called hello_string.
- hello_string = “Hello, world!” assigns the value “Hello, world!” to the declared array. Unlike in other languages like C, this can’t be done on the same line as declaring the array.
- write (*,*) hello_string prints the value of hello_string to the standard output. The first * means to write to standard output, as opposed to some file. The second * means not to use any special formatting.
Add a comment.
This isn’t necessary in such a simple program, but it will be useful when you write something more complex, so you should know how to add them. There are two ways to add a comment.
- To add a comment that has an entire line on its own, write a c directly into a new line, without the 6 spaces. After that, write your comment. You should leave a space between the c and your comment for better readability, but this isn’t required. Note that you have to use a ! instead of a c in Fortran 95 and newer.
- To add a comment in the same line as code, add a ! where you want your comment to begin. Again, a space isn’t required, but improves readability.
Using Input and If-Constructions
Understand different data types.
- INTEGER is used for whole numbers, like 1, 3, or -3.
- REAL can also contain a number that isn’t whole, like 2.5.
- COMPLEX is used to store complex numbers. The first number is the real and the second the imaginary part.
- CHARACTER is used for characters, like letters or punctuation.
- LOGICAL can be either .true. or .false.. This is like the boolean type in other programming languages.
Get the user’s input.
In the “Hello world” program that you wrote before, getting user input would be useless. So open a new file and name it compnum.f. When you’ve finished it, it will tell the user whether the number they entered is positive, negative or equal to zero. program compnum real r write (*,*) “Enter a real number:” read (*,*) r end program
- Enter the lines program compnum and end program compnum.
- Then, declare a variable of the type REAL. Make sure that your declaration is between the beginning and the end of the program.
- Explain the user what they’re supposed to do. Write some text with the write function.
- Read the user’s input into the variable you declared with the read function.
Process the user’s input with an if-construction.
Put it between the read (*,*) r and the end program. if (r .gt. 0) then write (*,*) “That number is positive.” else if (r .lt. 0) then write (*,*) “That number is negative.” else write (*,*) “That number is 0.” end if
- Comparison is done with .gt. (greater than), .lt. (less than) and .eq. (equals) in Fortran.
- Fortran supports if, else if, and else
- A Fortran if-construction always ends with end if.
Compile and run your program.
Input some numbers to test it. If you enter a letter, it will raise an error, but that’s okay because the program doesn’t check whether the input is a letter, a number, or something else.
Using Loops and Arrays
Open a new file.
Since this concept is different, you’ll have to write a new program again. Name the file addmany.f. Insert the corresponding program and end program statements, as well as an implicit none. When you’re finished, this program will read 10 numbers and print their sum.
Declare an array of length 10.
This is where you will store the numbers. Since you probably want a sum of real numbers, you should declare the array as real. You declare such an array with real numbers(50) (numbers is the name of the array, not an expression).
Declare some variables.
Declare numSum as a real number. You will use it to store the sum later, but since sum is already taken by a Fortran expression, you have to use a name like numSum. Set it to 0. Declare i as an integer and don’t assign it any value yet. That will be done in the do-loop.
Create a do-loop.
The equivalent of that in other programming languages would be a for-loop.Your do loop should now look like this: do 1 i = 1, 10 numSum = numSum + numbers(i) 1 continue
- A do-loop always starts with do.
- On the same line as the do, separated from it by a space, is the label to which the program will go when it’s finished. For now, just write a 1, you’ll set the label later.
- After that, again only separated by a space, type i = 1,10. This will make the variable i, which you had declared before the loop, go from 1 to 10 in steps of 1. The steps aren’t mentioned in this expression, so Fortran uses the default value of 1. You could also have written i = 1,10,1
- Put some code inside the loop (indent with spaces for better readability). For this program, you should increase the variable numSum with the i-th element of the array numbers. This is done with the expression numSum = numSum + number(i)
- End the loop with a continue statement that has a label. Type only 4 spaces. After that, type a 1. That’s the label which you told the do-loop to go to after it finishes. Then, type a space and continue. The continue expression does nothing, but it gives a good spot to place a label, as well as showing that the do-loop ended.
Print numSum.
Also, it would make sense to give some context, for example “The sum of your numbers is:”. Use the write function for both. Your entire code should now look as follows: program addmany implicit none real numbers(10) real numSum integer i numSum = 0 write (*,*) “Enter 10 numbers:” read (*,*) numbers do 1 i = 1, 10 numSum = numSum + numbers(i) 1 continue write (*,*) “Their sum is:” write (*,*) numSum end program addmany
Compile and run your code.
Don’t forget to test it. You can either press ↵ Enter after each number you enter or enter many numbers on the same line and separate them with a space.
Understanding Advanced Concepts
Have a good idea of what your program will do.
Think about what sort of data is needed as input, how to structure the output, and include some intermediate output so you can monitor the progress of your calculation. This will be very useful if you know your calculation will run for a long time or involves multiple complicated steps.
Find a good Fortran reference.
Fortran has many more functions than explained in this article, and they might be useful for the program you want to write. A reference lists all functions a programming language has. This is one for Fortran 77 and this is one for Fortran 90/95.
Learn about subroutines and functions.
Learn how to read and write from/to files.
Also learn how to format your input/output.
Learn about the new features of Fortran 90/95 and newer.
Skip this step if you know that you’ll only be writing/maintaining Fortran 77 code.
- Remember that Fortran 90 introduced the “Free Form” source code, allowing code to be written without the spaces and without the 72 character limit.
Read or look up some books on Scientific Programming.
For example, the book “Numerical Recipes in Fortran” is both a good text on scientific programming algorithms and a good introduction to how to put together codes. More recent editions include chapters on how to program in a mixed-language environment and parallel programming. Another example is “Modern Fortran in Practice” written by Arjen Markus. The book gives an insight into how to write Fortran programs in twenty-first-century style in accordance with the latest Fortran standards.
Learn how to compile a program spread across multiple files.
Let’s assume that your Fortran program is spread across the files main.f and morestuff.f, and that you want the resulting binary to be named allstuff. Then you’ll have to write following commands into the command line:f77 -c morestuff.ff77 -c main.ff77 -c morestuff.ff77 -o allstuff main.o morestuff.fThen run the file by typing ./allstuff.
Use the optimization your compiler provides.
Most compilers include optimization algorithms that improve the efficiency of your code. These are typically turned on by including a -O , -O2, or -O3 flag when compiling (again depending upon your version of fortran).
- Generally, the lowest level -O or -O2 level is best. Be aware that using the more aggressive optimization option can introduce errors in complex codes and may even slow things down! Test your code.
Tips
- You might find it easier to use an online IDE (integrated development environment) at first. A good option is Coding Ground. You will find a multitude of programming languages there, including Fortran-95. Another option is Ideone.
- Start with small programs. When you are making your own code, try to identify the most essential part of the problem – is it the data input or the calling of the functions, the structure of the loop (these are some very elementary examples) and start from there. Then build upon that in small increments.
- Fortran is not case-sensitive. You could, for example, declare a variable “real Num” and write “num = 1” in the next line to assign a value to it. But that’s a bad style, so avoid it. More importantly, Fortran doesn’t care about the case of functions and statements either. It’s quite common to write functions and statements in UPPERCASE and variables in lowercase.