anchor.diff:      hecMAT%NDOF = 1
       hecMAT%Iarray(98) = 1 !Assmebly complete
       hecMAT%X = 0.0d0
  -    current_time = start_time
  +    current_time = 0.0d0
  +    next_time = 0.0d0
   
       if(fstrHEAT%is_steady == 1)then
         fstrHEAT%beta = 1.0d0
  
positive.code:     if(ISTEP == 1)then
        start_time = 0.0d0
      else
        start_time = fstrHEAT%STEP_EETIME(ISTEP - 1)
      endif
  
      delta_time_base = fstrHEAT%STEP_DLTIME(ISTEP)
  
negative.code:       if( fstrHEAT%is_steady == 1 ) then
          is_end = .true.
        else
          if( (end_time - next_time) / end_time < 1.d-12 ) is_end = .true.
        endif
        if( 0.0d0 < DELMIN .and. fstrHEAT%timepoint_id > 0 ) then
          outflag = is_end .or. is_at_timepoints(next_time, 0.0d0, fstrPARAM%timepoints(fstrHEAT%timepoint_id))
        else
          outflag = is_end
        endif
        total_time = start_time + next_time
  
        if( hecMESH%my_rank.eq.0 ) then
          write(IMSG,"(a,i8,a,1pe12.5,a,1pe12.5)") " ** Increment No. :", total_step, ", current time: ", &
          & current_time, ", delta t: ", delta_time
          write(*,   "(a,i8,a,1pe12.5,a,1pe12.5)") " ** Increment No. :", total_step, ", current time: ", &
          & current_time, ", delta t: ", delta_time
        endif
  
        if(delta_time_base < DELMIN .and. (.not. outflag))then
  