anchor.diff:  
       fstrSOLID%dunode(:) =0.d0
   
  -    a1 = 1.d0/fstrDYNAMIC%t_delta**2
  -    a2 = 1.d0/(2.d0*fstrDYNAMIC%t_delta)
  +    a1 = 1.d0/fstrDYN%t_delta**2
  +    a2 = 1.d0/(2.d0*fstrDYN%t_delta)
   
       call setMASS(fstrSOLID,hecMESH,hecMAT,fstrEIG)
       call hecmw_mpc_trans_mass(hecMESH, hecMAT, fstrEIG%mass)
  
positive.code:           endif
            if( fstrPARAM%fg_couple_window > 0 ) then
              j = i - restrt_step_num + 1
              kk = fstrDYNAMIC%n_step - restrt_step_num + 1
              bsize = 0.5*(1.0-cos(2.0*PI*dfloat(j)/dfloat(kk)))
              do kkk0 = 1, fstrCPL%coupled_node_n
                kkk1 = 3 * kkk0
  
negative.code:     endif
  
      call hecmw_mpc_mat_finalize(hecMESH, hecMAT, hecMATmpc)
  
    end subroutine fstr_solve_dynamic_nlexplicit
  
  