anchor.diff:           call dynamic_mat_ass_load (hecMESH, hecMAT, fstrSOLID, fstrDYNAMIC)
            do j=1, hecMESH%n_node*  hecMESH%n_dof
              hecMAT%B(j)=hecMAT%B(j)- fstrSOLID%QFORCE(j) + myEIG%mass(j)*( fstrDYNAMIC%VEC1(j)-a3*fstrSOLID%dunode(j)   &
  -		     + fstrDYNAMIC%ray_m* hecMAT%X(j) ) + fstrDYNAMIC%ray_k*fstrDYNAMIC%VEC3(j)
  +         + fstrDYNAMIC%ray_m* hecMAT%X(j) ) + fstrDYNAMIC%ray_k*fstrDYNAMIC%VEC3(j)
            end do
   
   !C ********************************************************************************
  
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:       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
  