Programming

s = " 1.0007126 6.8643880e+05 3.2924832e-09 9.4557639e+02  1.6407053 4.3481956e+03"

open(unit=unit_bgnd, file=fin, form='formatted', status='old')
! Compute the number of lines in the file
nfl = 0
nhdr = 0
do
  read(unit_bgnd, '(A)', iostat=ios) s
  if (ios > 0) then
    write(*,*) ' Error reading file: ', fin
    exit
  else if (ios < 0) then
    write(*,'(A,I4)') '  No. of lines read from file: ', nfl
    exit
  else
    if (s(:1) == '#') nhdr = nhdr + 1
    nfl = nfl + 1
  end if
end do
! Move read position to beginning of file
rewind(unit_bgnd, iostat = ios)
nfv = nfl - nhdr
print *, 'nhdr, nfl, nv = ', nhdr, nfl, nfv
allocate(zfv(nfv))
allocate(cfv(nfv))
! Read values from the file containing Model S
ifv = 0
allocate(val(6))
do ifl = 1, nfl
  read(unit_bgnd, '(A)', iostat=ios) s
  print *, 's = ', s, ', s(1:1) = ', s(1:1)
  if (s(1:1) /= '#') then
      ifv = ifv + 1
      read(s, *) (val(i), i = 1, 6)
      print *, val(1), val(2)
      zfv(ifv) = val(1)  ! Gets the ratio radius
      cfv(ifv) = val(2)  ! Gets the sound speed in cm/sec
  endif
end do
deallocate(val)

float a, b;
if (sscanf(s, " %f %f", &a, &b) != 2)
    printf("Error\n");
else
    printf("Values are %f and %f\n", a, b);