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:  102  FORMAT (A5,1X,A2,1X,A4)
   103  FORMAT (A5,1X,A8)
        READ (LIN, *) DATATYPE
  !     WRITE(*,*) 'DATA FORMAT:  ', DATATYPE
  C
        IF (DATATYPE .EQ. 'NASA') THEN
           READ (LIN, *) MARK, T_int1
  !        WRITE(*,*) MARK, T_int1
        ELSE
           READ (LIN, 101)
           T_int1=0.D0
        ENDIF
  C
        READ (LIN, *) STRUC_MOL
  !     WRITE(*,*) STRUC_MOL
        READ (LIN, *) ATOMS
  !     WRITE(*,*) 'NO. OF ATOMS', ATOMS
        READ (LIN, *) ROTORS
  !     WRITE(*,*) 'NO. OF ROTORS', ROTORS
   104  FORMAT(F8.1)
  !      CLOSE(LIN)
  !     PAUSE
  C
  C: Keywords input over
  C
  
negative.code:  303  FORMAT(ES15.8,ES15.8,ES15.8,ES15.8,15X,4X,'4')
  C
        ELSEIF (DATATYPE .EQ. 'NIST') THEN
  C      write(*,*) 'THE COEFFICIENTS OF NIST FORM'
  C      write(*,*) THERM1(1:7)
  C      write(*,*) THERM2(1:7)
  C
         WRITE(LOUT,*) 'The NIST polynomial coefficients calculated:'
         WRITE(LOUT,*) 'H0-H0_298=A*t+B*t*t/2+C*t*t*t/3+D*t*t*t*t/4'
  