anchor.diff:        COMMON /INDICTOR/INDC_ELNO
   C
         IF (DATATYPE .EQ. 'WILH') THEN
  -!       write(*,*) 'THE COEFFICIENTS OF WILHOIT FORM'
  -!       write(*,*) THERM1(1:6)
  +C       write(*,*) 'THE COEFFICIENTS OF WILHOIT FORM'
  +C       write(*,*) THERM1(1:6)
           WRITE(LOUT,*) 'The Wilhoit polynomial coefficients calculated:'
           WRITE(LOUT,200)
           WRITE(LOUT,201) THERM1(1:6)
  
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:       DOUBLE PRECISION TEMP,CPT,CH,CS
        DIMENSION TEMP(*),CPT(*),CH(*),CS(*)
        DOUBLE PRECISION TEMP1,C1,CH1,CS1
  C     DIMENSION TEMP1(50),C1(50),CH1(50),CS1(50)
  C
  C     DATA TEMP1/290.262D0, 393.258D0, 496.255D0, 589.888D0, 795.88D0,
  C     &           1001.87D0, 1488.76D0/
  C
  C     DATA  C1/17.7725D0, 22.4645D0, 27.0142D0, 30.7109D0, 36.9668D0,
  C     &            41.8009D0, 49.3365D0/
  C
  C     DATA CH1/0.D0,0.D0,0.D0,0.D0,0.D0,0.D0,0.D0/
  C
  C     DATA CS1/0.D0,0.D0,0.D0,0.D0,0.D0,0.D0,0.D0/
  C
        MARK=''
        TEXT=''
  