anchor.diff:          if( rres > fstrSOLID%step_ctrl(cstep)%maxres    ) fstrSOLID%NRstat_i(knstDRESN) = 2
           return
         end if
  -
       enddo
       ! ----- end of inner loop
   
  
positive.code:     if( fstrSOLID%step_ctrl(cstep)%solution == stepStatic ) tincr = 0.d0
  
      P = 0.0d0
      rres = 0.0d0
      stepcnt = 0
      fstrSOLID%dunode(:) = 0.0d0
      fstrSOLID%NRstat_i(:) = 0 ! logging newton iteration(init)
      hecMAT%X = 0.0d0
  
      call fstr_ass_load(cstep, hecMESH, hecMAT, fstrSOLID, fstrPARAM)
  
  
negative.code:     integer(kind=kint) :: i, iter
      integer(kind=kint) :: stepcnt
      integer(kind=kint) :: restrt_step_num
      real(kind=kreal)   :: tt0, tt, res, qnrm, rres, tincr, xnrm, dunrm, rxnrm
      real(kind=kreal), pointer :: coord(:), P(:)
      logical :: isLinear = .false.
  
  