anchor.diff:   102  FORMAT (A5,1X,A2,1X,A4)
    103  FORMAT (A5,1X,A8)
         READ (LIN, *) DATATYPE
  -!     WRITE(*,*) 'DATA FORMAT:  ', DATATYPE
  +C     WRITE(*,*) 'DATA FORMAT:  ', DATATYPE
   C
         IF (DATATYPE .EQ. 'NASA') THEN
            READ (LIN, *) MARK, T_int1
  -!        WRITE(*,*) MARK, T_int1
  +C        WRITE(*,*) MARK, T_int1
         ELSE
            READ (LIN, 101)
            T_int1=0.D0
         ENDIF
   C
         READ (LIN, *) STRUC_MOL
  -!     WRITE(*,*) STRUC_MOL
  +C     WRITE(*,*) STRUC_MOL
         READ (LIN, *) ATOMS
  -!     WRITE(*,*) 'NO. OF ATOMS', ATOMS
  +C     WRITE(*,*) 'NO. OF ATOMS', ATOMS
         READ (LIN, *) ROTORS
  -!     WRITE(*,*) 'NO. OF ROTORS', ROTORS
  +C     WRITE(*,*) 'NO. OF ROTORS', ROTORS
    104  FORMAT(F8.1)
  -!      CLOSE(LIN)
  -!     PAUSE
  +C      CLOSE(LIN)
  +C     PAUSE
   C
   C: Keywords input over
   C
  
positive.code:        WRITE(LOUT,303) THERM2(4),THERM2(5),THERM2(6),THERM2(7)
  C
        ELSE IF (DATATYPE .EQ. 'NASA') THEN
  !      write(*,*) 'THE COEFFICIENTS OF NASA FORM'
  !      write(*,*) THERM1(1:9)
  !      write(*,*) THERM2(1:9)
  !      write(*,*) THERM3(1:9)
  C
         WRITE(LOUT,*) 'The NASA polynomial coefficients calculated:'
         WRITE(LOUT,399) SNAM
  ! HF_298 in 'J/mol' in NASA format: 1 cal = 4.184 J
  !      WRITE(*,*) MW
  !      WRITE(*,*) DLTH*4.184D0
         WRITE(LOUT,400) 3,ENAM(1),ELNO(1),'.00',ENAM(2),ELNO(2),'.00',
       &                   ENAM(3),ELNO(3),'.00',ENAM(4),ELNO(4),'.00',
       &                 ENAM(5),ELNO(5),'.00',MW,DLTH*4.184D0
  
negative.code: C
  C: 'M' is the size of sample dataset
        INTEGER M_in
  C     PARAMETER (M_in=59)
        PARAMETER (M_in=7)
        DOUBLE PRECISION TEMP,CPT,CH,CS,TEMF,CPF,CHF,CSF
        DIMENSION TEMP(M_in),CPT(M_in),CH(M_in),CS(M_in)
  