Fortran 90 kind parameter

by

The KIND of a variable is an integer label which tells the compiler which of its supported kinds it should use.

Beware that although it is common for the KIND parameter to be the same as the number of bytes stored in a variable of that KIND, it is not required by the Fortran standard.

That is, on a lot of systems, 

REAl(KIND=4) :: xs   ! 4 byte ieee float
REAl(KIND=8) :: xd   ! 8 byte ieee float
REAl(KIND=16) :: xq   ! 16 byte ieee float

but there may be compilers for example with:

REAL(KIND=1) :: XS   ! 4 BYTE FLOAT
REAL(KIND=2) :: XD   ! 8 BYTE FLOAT
REAL(KIND=3) :: XQ   ! 16 BYTE FLOAT

Similarly for integer and logical types.

(If I went digging, I could probably find examples. Search the usenet group comp.lang.fortran for kind to find examples. The most informed discussion of Fortran occurs there, with some highly experienced people contributing.)

So, if you can’t count on a particular kind value giving you the same data representation on different platforms, what do you do? That’s what the intrinsic functions SELECTED_REAL_KIND and SELECTED_INT_KIND are for. Basically, you tell the function what sort of numbers you need to be able to represent, and it will return the kind you need to use.

I usually use these kinds, as they usually give me 4 byte and 8 byte reals:

!--! specific precisions, usually same as real and double precision
integer, parameter :: r6 = selected_real_kind(6) 
integer, parameter :: r15 = selected_real_kind(15)

So I might subsequently declare a variable as:

real(kind=r15) :: xd

Note that this may cause problems where you use mixed language programs, and you need to absolutely specify the number of bytes that variables occupy. If you need to make sure, there are enquiry intrinsics that will tell you about each kind, from which you can deduce the memory footprint of a variable, its precision, exponent range and so on. Or, you can revert to the non-standard but commonplace real*4, real*8 etc declaration style.

When you start with a new compiler, it’s worth looking at the compiler specific kind values so you know what you’re dealing with. Search the net for kindfinder.f90 for a handy program that will tell you about the kinds available for a compiler.

More Related Contents:

  1. Fortran running error
  2. Fortran: integer*4 vs integer(4) vs integer(kind=4)
  3. Why Segmentation fault is happening in this openmp code?
  4. Unexpected “padding” in a Fortran unformatted file
  5. What is an undefined reference/unresolved external symbol error and how do I fix it in Fortran?
  6. Is There a Better Double-Precision Assignment in Fortran 90?
  7. MPI_Recv overwrites parts of memory it should not access
  8. Why should I use interfaces?
  9. Sending 2D arrays in Fortran with MPI_Gather
  10. Where to put `implicit none` in Fortran
  11. Print to standard output from a function defined in an Fortran module
  12. Why is this a function declared inside the module and then used somewhere else in the same module not seen by the linker?
  13. Does Fortran have undefined behavior?
  14. Reading data file in Fortran with known number of lines but unknown number of entries in each line
  15. What are the ways to pass a set of variable values through the subroutine to a function without common block?
  16. What does “real*8” mean?
  17. Line truncated, Syntax error in argument list
  18. Fortran IF statement with numbers/labels rather than another statement
  19. Using a deferred-length character string to read user input
  20. Function in fortran, passing array in, receiving array out
  21. Providing an argument that has not the TARGET attribute to a procedure with a dummy argument that has the TARGET attribute
  22. Fortran OpenMP with subroutines and functions
  23. Difference between intent(out) and intent(inout)
  24. Error: Nonnegative width required in format string at (1)
  25. Fortran intent(inout) versus omitting intent
  26. Fortran array cannot be returned in function: not a DUMMY variable
  27. Is Fortran unable to do the addition between 865398.78 and -865398.78? Why the answer is -0.03?
  28. Using do loop in a Fortran 90 program to read different number of lines for n frames?
  29. How to pass subroutine names as arguments in Fortran?
  30. Confusing double precision real in Fortran

Leave a Comment