Error when reading in float in Fortran

Question

This should be quite simple, but I can\'t manage to read in a floating point number in Fortran. My program test.f looks like this:

  PROGRAM TEST
  ope         


        

Answer 1

In your read statement

read(1,'(f3.0)')line

the f3.0 tells tour program to read 3 digits with 0 digits after the decimal (this is what the n.m syntax means). So I presume that the program is just reading 1 from the file (not 1.2), which is an integer. Try replacing that line with something like

read(1,'(f3.1)')line

although, if the number in your file is likely to change and be larger than 9.9 or have more than one decimal place you should increase the field width to something larger than 3.

See the documentation of the read intrinsic and for data edit descriptors for more information on reading and writing in Fortran.

Edit: the format specifier, the second argument in quotes in your read statment, has the form fw.d, where f indicates that the data to read is a floating point number, w is the width of the field including all blanks and decimal points and d specifies the number of digits to the right of the decimal point.