Formatted Input/Output
List-directed IO is convenient.
But that frequently results in sprawling output that is difficult to read. Formatted output is frequently required for legibility.
Formatted output in Fortran is similar to other languages (the general layout descends from Fortran, the oldest higher-level programming language).
Edit Descriptors
The edit descriptor modifies how to output the variables. They are combined into forms like
RdF.w
where R is a repeat count, d is the descriptor, F is the total field width including space for +-, and if requested +-e and exponent, and w is the number of digits to the right of the decimal point. If you are willing to let the compiler calculate the number of characters to use, use Rd0.w. You may also omit the number of decimal places with Rd0 and the compiler will use its default for the type.
For floating-point numbers, Rd.0 will print the integer part.
Strings take only RaF and do not usually require the F since the length will be known to the compiler.
Integers can be written as iF and any of the F spaces not needed will be blank filled, with the digits right justified. When written as iF.m they will be printed with at least m digits and the rest of the field zero-filled on the left if all of F is not needed.
If the field width is specified and the requested literal does not fit, the compiler will output a string of asterisks, e.g. ********.
If you allow the compiler to compute the total field width, note that it will not include spaces before or after the item.
Common Edit Descriptors
As usual, they are not case sensitive.
I !integer
F !real (decimal output)
E !real (exponential output)
ES !like E but use scientific notation.
G !general
D !double precision (prints D rather than E for exponent)
A !character (does not require a field width in most cases)
X !space
/ !write an EOL and go to the next line (record) within the format
: !terminate the output if there are no more variables to write
The real descriptors F, E, G, and D all work for both single and double precision. G allows the compiler to choose whether to use decimal or exponential format.
The default exponential format writes in machine normalization, with the leading digit between 0 and 1. ES causes it to write with the leading digit between 1 and 9, which is what most humans can read most easily. ES ignores p on output.
Modifiers
Some modifiers can change the appearance of the output.
p !multiply by 10
kp !multiply by 10k
The p descriptor applies till the next scale factor is encountered.
Fill an integer field with zeros
I4.4
Format Strings
The format string is constructed as a list of how to output the variables. Unlike some other languages, literal strings are never included, but must have their own edit descriptors.
The format string can be placed directly into the write statement or it can be in a separate format statement. In the write it is enclosed in parentheses and quotes.
For most purposes it is best to put the format string into the write statement. The format statement is older and will be in old code, but it is usually harder to see what is happening. It is useful for particularly long strings, however.
Examples
write(*,'(i5,2x,i6)') i1,i2
write(*,'(i5,a,i6)') i1," ",i2
write(*,'(a,f0.6)') "The result is ",res
write(*,'(a,i4,es15.7)') "The answer",i1,dpi
write(*,'(2p,f8.2,0p,f8.2)') rpi, dpi
write(*,'(a,f8.2,/,a,i6)') mess1,res,mess2,i1
write(*,'(a)') ' '
write(*,'("first value ",f8.2,", second value ",i6)') res
write(*,'(a)') ' '
write(*,'("first value ",f8.2,:," second value ",i6)') res
A format string may be a variable.
character(len=32) :: formatstr
code
formt='(f8.3,es15.7)'
write(*,formatstr) A, B
Repetition
Format strings can be repeated for multiple variables. If more than one descriptor is present, the format to be repeated should be enclosed in parentheses.
write(*,'(2L)')is_zero,is_finite
write(*,'(2f8.2)') z !complex
write(*,'4(f0.6)') a(1,:)
write(*,'(4(i2,3x,f8.3))') (j,b(j),j=1,4)
Especially when an array is allocatable, it may be awkward to specify the repeat count if it is unknown at compile time. List-directed I/O allows the compiler to choose to add end-of-line markers to line up columns, so the output can differ between different compilers and may not be what is desired for later processing. In principle a variable format string can be constructed with internal writes, but this can be complicated. The Fortran 2008 standard introduced the * repetition count. The compiler will repeat until it runs out of items.
do i=1,size(arr,1)
write(*,'(*(f12.4))') arr(i,:)
enddo
Format Statements
Format statements are abundant in older code, before the strings could be inserted into writes. FORMAT is non-executable but can appear anywhere in the source. It is the only non-executable statement that can do so.
It can still be useful for a particularly complex format (to keep the write statement short and readable) or for formats that are repeated in many write statements. The second parameter to the write is then an integer statement label. The label marks the format statement.
Example
write(*,100)x,y,z
100 format(3e15.8)
Traditionally the format is placed immediately below the line which refers to it, or else all format statements are grouped together just before the end statement of their program unit.
Formatted Input
In Fortran it is best to avoid formatted input as much as possible, as it can lead to errors. For historical reasons, if a format is specified to be real (floating point) but no decimal point is included in the data, the compiler inserts it based on the format. For example, suppose we had data in the form
112 9876 12
with a format string of
read(infile,'(f8.4,f6.2,i4)'
We intended to read three values, two reals and an integer, but this format results in input values of
112.987602 0.119999997 0
The errors in the real values are the consequence of converting from decimal to binary and back.
Note that the spaces between the values are ignored.
We get the expected result with read(infile,*):
112.000000 9876.00000 12
The compiler has more freedom so the conversion is also more accurate.
Unformatted input is permitted when casting from a character to a real, so there is no need for formatted input at all.
character(len=12) :: quantity_char
real :: quantity
quantity_char="100"
read(quantity_char,*) quantity
Fortran Non-Advancing IO
Unlike most languages, Fortran print and write by default add an end-of-line marker at the end of the output. If we’d like to suppress this so that we can write multiple groups of output on the same line, or we would like to write some output to the console and read something from the console, we can use non-advancing IO.
write(*,'(a)',advance='no') "Enter input value:"
read(*,*) value
Non-advancing IO must be formatted
- ‘yes’ for advance is valid also but is the default.
- The argument to
advancecan be a character variable so you can decide based on conditionals to advance or not. - If you do not want to advance, use
advance='no'
Exercises
- Examine this example code:
A variety of formatting examples
program printit
implicit none
integer, parameter :: rk=kind(1.0)
integer, parameter :: dk=kind(1.d0)
complex :: z
real(rk) :: rpi=4.0*atan(1.0)
real(rk) :: dpi=4.d0*atan(1.d0)
real(rk) :: res
integer :: i,j
integer :: i1=42, i2=99
logical :: is_zero, is_finite
character(len=12) :: mess1, mess2
real(rk), dimension(4) :: b
real(rk), dimension(4,4) :: arr
res=cos(rpi)
do i=1,size(arr,1)
b(i)=real(i)
do j=1,size(arr,2)
arr(i,j)=real(i+j)
enddo
enddo
z=cmplx(res,rpi)
mess1="First result"
mess2="An integer"
is_zero=.false.
is_finite=.true.
write(*,'(i5,2x,i6)') i1,i2
write(*,'(i5,a,i6)') i1," ",i2
write(*,'(a,f0.6)') "The result is ",res
write(*,'(a,i4,es15.7)') "The answer",i1,dpi
write(*,'(4(i2,3x,f8.3))') (j,b(j),j=1,4)
write(*,'(2f8.2)') z !complex
write(*,'(2L)')is_zero,is_finite
write(*,'(2p,f8.2,0p,f8.2)') rpi, dpi
write(*,'(a,f8.2,/,a,i6)') mess1,res,mess2,i1
write(*,'(a)') ' '
write(*,'("first value ",f8.2,", second value ",i6)') res
write(*,'(a)') ' '
write(*,'("first value ",f8.2,:," second value ",i6)') res
write(*,*) '----------------------'
do i=1,size(arr,1)
write(*,'(*(g0))') arr(i,:)
enddo
write(*,*) '----------------------'
do i=1,size(arr,1)
write(*,'(*(f12.4))') arr(i,:)
enddo
! Old style
WRITE ( *, 2 ) i1
WRITE ( *, 2 ) i1,i2
2 FORMAT ( 1X 'RESULT = ', I2, :, 3X, 'AT INDEX = ', I2 )
write(*,*)
write(*,'(f8.2)',advance='no') res
write(*,'(i6)',advance='yes') i2
write(*,*) 'All done'
end program
Make sure you understand all the formats. Correct unattractive outputs. Experiment with changing the formats.
- Write a program that computes pi using a trig identity such as
pi=4*atan(1).- Use kind to switch between real and double precision
- integer, parameter ::rk=kind(1.0) or (1.0d0)
- Using single precision, print pi in
- E format
- Scientific notation
- Scientific notation with 8 decimal places Repeat for double precision.
- Use kind to switch between real and double precision
Solution with variable strings.
program print_pi
implicit none
integer, parameter :: rk=kind(1.0)
!integer, parameter :: rk=kind(1.0d0)
real(rk) :: pi=4.0_rk*atan(1.0_rk)
write(*,'(g12.6)') pi
write(*,'(e12.3)') pi
write(*,'(es12.3)') pi
write(*,'(es15.8)') pi
end program
- In an “infinite” while loop: Request an integer from the user with non-advancing input/output, e.g.
"Please enter an integer:" <then read integer>
If the integer is 1, print “zebra”. If it is 2, print “kangaroo”. If it is anything else other than zero, print “not found”. If it is 0, exit the loop. This requires only a simple change to a previous program.
Solution with variable strings.
program get_animal
implicit none
integer :: ans
character(len=8) :: text
do
write(*,'(a)',advance='no') "Please enter a digit 1 or 2, or 0 to exit:"
read(*,*) ans
if (ans==0) then
exit
else if (ans==1) then
print *, "zebra"
else if (ans==2) then
print *, "kangaroo"
else
print *, "not found"
endif
enddo
end program