Subroutine to define the bearing geometry
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(FE_FILM), | intent(inout) | :: | fe_f | |||
real(kind=R8), | intent(in) | :: | ly |
subroutine modify_h_bearing_y(fe_f, ly)
implicit none
type(FE_FILM), intent(inout) :: fe_f
real(kind=R8), intent(in) :: ly
integer(kind=I4) :: i
do i = 1, fe_f%m%n
fe_f%vn(i, H2_N) = fe_f%data_f%h_0 + 0.5_R8 * fe_f%data_f%h_0 * cos(2 * PI_R8 * fe_f%m%y(i) / ly)
enddo
fe_f%vn(:, H_N) = fe_f%vn(:, H2_N) - fe_f%vn(:, H1_N)
return
endsubroutine modify_h_bearing_y