anchor.diff:        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)
  +      call heat_output_result(hecMESH, fstrHEAT, total_step, next_time, outflag)
  +      call heat_output_visual(hecMESH, fstrRESULT, fstrHEAT, total_step, next_time, outflag)
  +      call heat_output_restart(hecMESH, fstrHEAT, total_step, outflag, next_time)
   
         total_step = total_step + 1
         current_time = next_time
  
positive.code:     hecMAT%NDOF = 1
      hecMAT%Iarray(98) = 1 !Assmebly complete
      hecMAT%X = 0.0d0
      current_time = 0.0d0
      next_time = 0.0d0
  
      if(fstrHEAT%is_steady == 1)then
        fstrHEAT%beta = 1.0d0
  
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
  