anchor.diff:           CS(I)  =0.D0
         ENDDO
   C
  -!      CALL DATAINPUT (M_in,TEMP,CPT,CH,CS)
  +C      CALL DATAINPUT (M_in,TEMP,CPT,CH,CS)
         CALL DATAGROUP (M_in,TEMP,CPT,CH,CS,LIN)
   C
         TINT  = 0.D0
         TINT1 = 0.D0
   C
  -! SWILT uses the enthalpy in: cal/mol
  +C SWILT uses the enthalpy in: cal/mol
         DO I=1,M_in
            CH(I) = CH(I)*1.D3
         ENDDO
   C
  -! First fit SWILT because it takes into account limits at 0 and infinite T
  +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
  
positive.code:       ENDDO
  C
  
  !   We have already read H_298 line, so continue with S_298
  !     READ (LIN, *) MARK, H_298
  !     WRITE(*,*) MARK, H_298
        READ (LIN, *) MARK, S_298
  !     WRITE(*,*) MARK, S_298
        READ (LIN, *) MARK, DLTH
  !     WRITE(*,*) 'HEAT OF FORMATION at 298K', DLTH
        READ (LIN, *) MARK, MW
  !     WRITE(*,*) 'MOLAR WEIGHT', MW
        READ (LIN, *) MARK, T_int
  !     WRITE(*,*) MARK, T_int
        READ (LIN, *) MARK, T_min
  !     WRITE(*,*) MARK, T_min
        READ (LIN, *) MARK, T_max
  !     WRITE(*,*) MARK, T_max
        if (T_max .gt. 6000.0) then
        write(*,*) 'Warning!!!'
        write(*,*) 'The maximum T should be no greater than 6000K!'
  
negative.code:  201  FORMAT(1X,E15.6,2X,E15.6,2X,E15.6,2X,E15.6,2X,E15.6,2X,E15.6)
  C
        ELSEIF (DATATYPE .EQ. 'CHEM') THEN
  C      write(*,*) SNAM
  C
         WRITE(LOUT,*) 'The Chemkin polynomial coefficients calculated:'
        if (J_elem .eq. 2) then
  