anchor.diff:              write(*,*) ' *** EXCEEDED TOLERANCE OF VARIATION IN TEMPERATUTE.'
               write(*,*) ' : NODE NUMBER  = ', mnod, ' : DELTA TEMP   = ', tmpmax
             endif
  -          delta_time = 0.5d0*delta_time
  +          delta_time_base = 0.5d0*delta_time_base
             cycle tr_loop
           endif
   
  -        if(iterALL <= miniter) delta_time = delta_time*1.5d0
  +        if(iterALL <= miniter) delta_time_base = delta_time_base*1.5d0
         else
           if(fstrHEAT%is_iter_max_limit) call hecmw_abort( hecmw_comm_get_comm() )
         endif
  
positive.code:       enddo
  
        call heat_output_log(hecMESH, fstrPARAM, fstrHEAT, total_step, next_time)
        call heat_output_result(hecMESH, fstrHEAT, total_step, next_time, is_end)
        call heat_output_visual(hecMESH, fstrRESULT, fstrHEAT, total_step, next_time, is_end)
        call heat_output_restart(hecMESH, fstrHEAT, total_step, is_end, next_time)
  
        total_step = total_step + 1
        current_time = next_time
  
negative.code:     implicit none
      integer(kind=kint) :: ISTEP, ITM, incr, iterALL, i, inod, mnod, max_step, tstep, interval, bup_n_dof
      integer(kind=kint) :: total_step
      real(kind=kreal)   :: start_time, delta_time_base, delta_time, current_time, next_time, end_time, DELMAX, DELMIN
      real(kind=kreal)   :: val, CHK, tmpmax, dltmp, tmpmax_myrank, remain_time
      type(hecmwST_local_mesh)  :: hecMESH
      type(hecmwST_matrix)      :: hecMAT
      type(hecmwST_result_data) :: fstrRESULT
  