anchor.diff:  ! ----- update the strain, stress, and internal force
           call fstr_UpdateNewton( hecMESH, hecMAT, fstrSOLID,fstrDYNAMIC%t_delta,iter,fstrDYNAMIC%strainEnergy )
   
  +        if(ndof == 4) exit
         enddo
   
   ! -----  not convergence
  
positive.code: !C-- RHS LOAD VECTOR CHECK
  !C
          numnp=hecMAT%NP
          bsize=0.0
          do iiii5 = 1,numnp*ndof
            bsize=bsize+hecMAT%B(iiii5)**2
          enddo
  !C-- Gather RHS vector
          call hecmw_allREDUCE_R1( hecMESH,bsize,hecmw_sum )
  
          if(iter == 1)then
            resb = bsize
  
negative.code: 
          !res = dsqrt(bsize)/n_node_global
          res = dsqrt(bsize/resb)
          if( ndof /= 4 ) then
            if(hecMESH%my_rank==0) write(*,'(a,i5,a,1pe12.4)')"iter: ",iter,", res: ",res
            if(hecMESH%my_rank==0) write(ISTA,'(''iter='',I5,''- Residual'',E15.7)')iter,res
            if( res<fstrSOLID%step_ctrl(cstep)%converg ) exit
          endif
  !C
  !C-- linear solver [A]{X} = {B}
  !C
  