modify_h_pocket Subroutine

  1. mod_test_musst.f90
  2. test_musst
  3. modify_h_pocket

private subroutine modify_h_pocket(fe_f, lx)

Arguments

Type IntentOptional AttributesName
type(FE_FILM), intent(inout) :: fe_f
real(kind=R8), intent(in) :: lx

Called by

proc~~modify_h_pocket~~CalledByGraph proc~modify_h_pocket modify_h_pocket proc~test_pocket_fe test_pocket_fe proc~test_pocket_fe->proc~modify_h_pocket proc~run_test run_test proc~run_test->proc~test_pocket_fe program~main main program~main->proc~run_test

Contents

Source Code

modify_h_pocket

Source Code

   subroutine modify_h_pocket(fe_f, lx)
   implicit none
   type(FE_FILM), intent(inout) :: fe_f
   real(kind=R8), intent(in)    :: lx

      integer(kind=I4) :: e, i
      logical(kind=I4) :: groove

      fe_f%vn(:,H2_N) = fe_f%data_f%h_0
      fe_f%vn(:,H1_N) = 0._R8
      do e = 1, fe_f%m%ne
         groove = .true.
         do i = 1, 4
            if (fe_f%m%x(fe_f%m%con(e,i)) > (lx / 2._R8)) groove = .false.
         enddo
         if (groove) then
            fe_f%vc(e, HG_C) = fe_f%data_f%h_g
         endif
      enddo
      fe_f%vn(:, H_N) = fe_f%vn(:, H2_N) - fe_f%vn(:, H1_N)

   return
   endsubroutine modify_h_pocket