anchor.diff:      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)   :: start_time, delta_time_base, delta_time, current_time, next_time, total_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
  
positive.code:         call hecmw_abort(hecmw_comm_get_comm())
        endif
  
        call heat_solve_main(hecMESH, hecMAT, hecMATmpc, fstrPARAM, fstrHEAT, ISTEP, iterALL, next_time, delta_time)
  
        if(0.0d0 < DELMIN)then
          tmpmax = 0.0d0
  
negative.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
  