anchor.diff:  c 
         DO I=1,8
           IJ=IJ+1
  -        READ (LIN, 100) MARK, TEXT
  +        READ (LIN, 100, END=34511) MARK, TEXT
           IF (MARK .EQ. 'SPEC') THEN
   !        WRITE(*,*) MARK
   !        WRITE(*,*) TEXT
  
positive.code:  250    CONTINUE
        RETURN
  C
        END
negative.code: C: M_in: the number of input thermo data, use 'DATAGROUP' to read the input data
  C***********************************************************************************
        PROGRAM LSAP
        use, intrinsic :: iso_fortran_env, only : 
       1                  LIN=>input_unit, LOUT=>output_unit
  C
        implicit none
  C
  C     INTEGER LIN,LOUT
        INTEGER I,J,J_elem,IJ
  C max atoms number allowed in a species is 999
  C max number of different atoms is 5 
        integer indc_elno(5,3)
  !      PARAMETER (LIN=5,LOUT=6) ! 5 is STDIN and 6 is STDOUT
  C For input file
        CHARACTER(LEN=4) MARK, DATATYPE
        CHARACTER(LEN=48) TEXT, SNAM ! beware truncation!
  