anchor.diff:      Nch(l) = Nch(l) - 1
       Nch(1) = -1
       do ic = 2, Nstop - 1
  -      call sum(ic,Xlnzr,Colno,Zln,Diag,Nch,Parent,Temp,Indx)
  +      call sum(ic,Xlnzr,Colno,Zln,Diag,Nch,Parent,Temp,Indx,ISEm)
       enddo
       !
       ! phase II
  
positive.code:   !          if(iv(22).eq.0)    normal type
    !          if(iv(22).gt.0)    code generation type
    !======================================================================!
    subroutine NUFCT0(Ir)
      use HECMW_UTIL
      implicit none
      !------
      integer, intent(out):: Ir
      !------
      integer:: ndeg2
      integer:: ndegl
  
      if ( STAge/=20 ) then
        print *, '*********Setting Stage 40!*********'
  
negative.code:   !======================================================================!
    !> @brief SUM
    !======================================================================!
    subroutine sum(Ic,Xlnzr,Colno,Zln,Diag,Nch,Par,Temp,Indx,ISEm)
      implicit none
      !------
      integer, intent(in):: Ic
  