anchor.diff:           CH(I) = CH(I)*1.D3
         ENDDO
   C
  -      WRITE (9,*) 'About to call SWILT', CH(M_in)
  -      CLOSE(9)
  -      OPEN (9, FORM='FORMATTED', FILE='debug.log', POSITION='APPEND')
  -C
   C First fit SWILT because it takes into account limits at 0 and infinite T
            CALL SWILT(SNAM,ENAM,ELNO,STRUC_MOL,ATOMS,ROTORS,THERM1,
        &                TEMP,CPT,CH,CS,M_in)
  -C
  -      WRITE (9,*) 'Completed the Wilhoit fitting in SWILT'
  -      CLOSE(9)
  -      OPEN (9, FORM='FORMATTED', FILE='debug.log', POSITION='APPEND')
   C
         IF (DATATYPE .EQ. 'WILH') GOTO 3000
   C
  
positive.code:         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
  
negative.code:          ENDIF
        ENDDO
  C
  
  C   We have already read H_298 line, so continue with S_298
  C     READ (LIN, *) MARK, H_298
  C     WRITE(*,*) MARK, H_298
  