! read a formatted file, see "example.txt"
      program TestReadFile
        implicit none

        character (len=80) :: filename
        integer, parameter :: fileunit=15
        logical :: bExist, bDM1, bDM2
        integer :: iOK, pos1, pos2
        character (len=80) :: line
        real :: var1, var2, var3

        ! name of input file
        filename="example.txt"

        ! check if file exist
        inquire (file=filename, exist=bExist)
        if (bExist) then
        else
          print *, filename, " not found!"
          stop
        endif

        ! open file
        open (unit=fileunit, file=filename, status='old', iostat=iOK)
        if (iOK /= 0) then
          print *, "can't open ", filename
          stop
        endif

        ! read file line by line
        do while (iOK == 0)
          read (fileunit, "(A)", iostat=iOK) line
         
          ! if line starts with "DM1"
          pos1=index(line, "DM1")
          if (pos1 == 1) then
            bDM1=.true.
            continue
          endif

          ! if line starts with "DM2"
          pos2=index(line, "DM2")
          if (pos2 == 1) then
            bDM1=.false.
            bDM2=.true.
            continue
          endif

          ! if in the 'DM1' block, read variables
          if ((bDM1) .and. (pos1 /= 1)) then
            read (line, *)var1, var2, var3
            print *, var1, var2, var3
          endif
         
        enddo

        ! close file
        close (fileunit)

      end program TestReadFile