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: C     WRITE(*,*) STRUC_MOL
        READ (LIN, *) ATOMS
  C     WRITE(*,*) 'NO. OF ATOMS', ATOMS
        READ (LIN, *) ROTORS
  C     WRITE(*,*) 'NO. OF ROTORS', ROTORS
   104  FORMAT(F8.1)
  C
  C: Keywords input over
  C
  C: INPUT the Cp, H, and S sample dataset
  