anchor.diff:          CALL DATAFIND(THERM1,TEMF(I),CPF(I),CHF(I),CSF(I))
           CHF(I) = CHF(I)/1.D3
         ENDDO
  +      WRITE (9,*) CHF(10)
   C
         DO I=12,101
   C     90 points between 298 and 6000
  
positive.code:  104  FORMAT(F8.1)
  C      CLOSE(LIN)
  C     PAUSE
  C
  C: Keywords input over
  C
  C: INPUT the Cp, H, and S sample dataset
  
negative.code:          ENDIF
        ENDDO
  C
        OPEN (9, FORM='FORMATTED', FILE='debug.log', POSITION='APPEND')
        WRITE(9,*) 'Read elements and H298'
  C   We have already read H_298 line, so continue with S_298
  C     READ (LIN, *) MARK, H_298
  C     WRITE(*,*) MARK, H_298
  