table_dyn.f90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
! Copyright (c) 2019 FrontISTR Commons
! This software is released under the MIT License, see LICENSE.txt
!-------------------------------------------------------------------------------
module m_table_dyn
 
  use m_fstr
  use hecmw
 
  implicit none
 
contains
 
  !C================================================================C
  !C-- subroutine fstr_get_amplitude_dyn
  !C================================================================C
  subroutine fstr_get_amplitude_dyn(hecMESH, fstrSOLID, fstrDYNAMIC, ig0, t_curr, value, flag_u)
    type(hecmwst_local_mesh) :: hecMESH
    type(fstr_solid)         :: fstrSOLID
    type(fstr_dynamic)       :: fstrDYNAMIC
    integer(kind=kint)       :: ig0
    real(kind=kreal)         :: t_curr
    real(kind=kreal)         :: value    
    integer(kind=kint)       :: flag_u
 
    integer(kind=kint) :: jj_n_amp
    real(kind=kreal)   :: t_eval
 
    jj_n_amp = 0
    if( flag_u == 1 ) then
      jj_n_amp = fstrsolid%BOUNDARY_ngrp_amp(ig0)
    else if( flag_u == 2 ) then
      jj_n_amp = fstrsolid%VELOCITY_ngrp_amp(ig0)
    else if( flag_u == 3 ) then
      jj_n_amp = fstrsolid%ACCELERATION_ngrp_amp(ig0)
    else if( flag_u == 0 ) then
      jj_n_amp = fstrsolid%CLOAD_ngrp_amp(ig0)
    else if( flag_u == 10 ) then
      jj_n_amp = fstrsolid%DLOAD_ngrp_amp(ig0)
    end if
 
    if( jj_n_amp <= 0 ) return  ! no amplitude: value unchanged (factor 1.0)
 
    t_eval = t_curr
    if( fstrdynamic%idx_eqa == 11 ) t_eval = t_curr - fstrdynamic%t_delta
 
    call hecmw_get_amplitude_value(hecmesh%amp, jj_n_amp, t_eval, value)
 
  end subroutine fstr_get_amplitude_dyn
 
end module m_table_dyn