anchor.diff:  C
    3000 CONTINUE
   C
  +      WRITE(9,*) 'About to call PRINT_THERM'
  +      CLOSE(9)
  +      OPEN (9, FORM='FORMATTED', FILE='debug.log', POSITION='APPEND')
  +      
         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(9,*) 'Finished everything'
  +      CLOSE(9)
         WRITE(LOUT,*) 'GATPFIT_HAS_FINISHED_ONE_INPUT'
         FLUSH(LOUT)
         GOTO 34501
   
   34511 CONTINUE
  +      CLOSE(9)
   C
         END
   C***********************************************************************
  
positive.code:          CH(I) = CH(I)*1.D3
        ENDDO
  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
        IF (DATATYPE .EQ. 'WILH') GOTO 3000
  C
  
negative.code:  104  FORMAT(F8.1)
  C      CLOSE(LIN)
  C     PAUSE
        WRITE(9,*) 'Read up to and including the No. of Rotors: ', ROTORS
  C
        CLOSE(9)
        OPEN (9, FORM='FORMATTED', FILE='debug.log', POSITION='APPEND')
        WRITE (9,*) 'I am able to close and re-open this log file'
        FLUSH(9)
  C: Keywords input over
  C
  C: INPUT the Cp, H, and S sample dataset
  