anchor.diff:           ENDIF
         ENDDO
   C
  -
  +      OPEN (9, FORM='FORMATTED', FILE='debug.log', POSITION='APPEND')
  +      WRITE(9,*) 'Read elements and H298'
   C   We have already read H_298 line, so continue with S_298
   C     READ (LIN, *) MARK, H_298
   C     WRITE(*,*) MARK, H_298
  
positive.code: C
   3000 CONTINUE
  C
        CALL PRINT_THERM(T_int,T_int1,SNAM,ENAM,ELNO,MW,H_298,DLTH,
       &                   DATATYPE,LOUT,THERM1,THERM2,THERM3,J_elem,
       &                 T_min,T_max)
  C
  C     CLOSE(LOUT)
        WRITE(LOUT,*) 'GATPFIT_HAS_FINISHED_ONE_INPUT'
        FLUSH(LOUT)
        GOTO 34501
  
  34511 CONTINUE
  C
        END
  C***********************************************************************
  
negative.code:          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
  