20 type(hecmwst_local_mesh) :: hecmesh
21 type(hecmwst_matrix) :: hecMAT
27 integer(kind=kint) :: j
31 if( fstrsolid%VELOCITY_type ==
kbcinitial )
then
32 do j = 1, hecmesh%n_node*hecmesh%n_dof
33 fstrdynamic%ACC(j,1)=(hecmat%B(j)-fstrdynamic%ray_m*fstreig%mass(j)*fstrdynamic%VEL(j,1))/&
36 elseif( fstrsolid%ACCELERATION_type ==
kbcinitial )
then
37 do j = 1, hecmesh%n_node*hecmesh%n_dof
38 fstrdynamic%VEL(j,1)=(hecmat%B(j)-fstreig%mass(j)*fstrdynamic%ACC(j,1))/&
39 (fstrdynamic%ray_m*fstreig%mass(j))