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: 
          !res = dsqrt(bsize)/n_node_global
          res = dsqrt(bsize/resb)
          if( hecMESH%my_rank==0 ) then
            write(*,'(a,i5,a,1pe12.4)')"iter: ",iter,", res: ",res
            write(ISTA,'(''iter='',I5,''- Residual'',E15.7)')iter,res
          endif
          if( res<fstrSOLID%step_ctrl(cstep)%converg ) exit
  !C
  !C-- linear solver [A]{X} = {B}
  !C
  
negative.code:       call dynamic_output_monit(hecMESH, fstrPARAM, fstrDYNAMIC, myEIG, fstrSOLID)
      end if
  
      fstrDYNAMIC%VEC3(:) =0.d0
      hecMAT%X(:) =0.d0
  
  !!
  !!    step = 1,2,....,fstrDYNAMIC%n_step
  !!
      do i= restrt_step_num, fstrDYNAMIC%n_step
          if(ndof == 4 .and. hecMESH%my_rank==0) write(*,'(a,i5)')"iter: ",i
  
         fstrDYNAMIC%i_step = i
         fstrDYNAMIC%t_curr = fstrDYNAMIC%t_delta * i
  