save_profile_x_ms Subroutine

  1. mod_inout_files.f90
  2. inout_files
  3. save_profile_x_ms

public subroutine save_profile_x_ms(ms_fe_f, file_name, lx, zy)

Arguments

Type IntentOptional AttributesName
type(MS_FE_FILM), intent(in) :: ms_fe_f

MS_FE_FILM element
character(len=*), intent(in) :: file_name

filename
real(kind=R8), intent(in) :: lx

surface length
real(kind=R8), intent(in) :: zy

distance from a border

Calls

proc~~save_profile_x_ms~~CallsGraph proc~save_profile_x_ms save_profile_x_ms proc~get_unit get_unit proc~save_profile_x_ms->proc~get_unit

Contents

Source Code

save_profile_x_ms

Source Code

   subroutine save_profile_x_ms(ms_fe_f, file_name, lx, zy)
   implicit none
   type(MS_FE_FILM), intent(in) :: ms_fe_f      !! [[MS_FE_FILM]] *element*
   character(len=*), intent(in) :: file_name    !! *filename*
   real(kind=R8),    intent(in) :: lx           !! *surface length*
   real(kind=R8),    intent(in) :: zy           !! *distance from a border*

      real(kind=R8)    :: dr
      integer(kind=I4) :: i, e, k

      dr = lx / ms_fe_f%bs_fe_f(1)%m%n
      dr = dr / 10
      call get_unit(k)
      open(k, file = file_name, status = 'unknown')
         do e = 1, ms_fe_f%ts_fe_f%m%ne
         do i = 1, ms_fe_f%bs_fe_f(e)%m%n
            if (((ms_fe_f%bs_fe_f(e)%m%y(i) - zy)**2) < (dr ** 2)) then
               write (k, *) ms_fe_f%bs_fe_f(e)%m%x(i), ms_fe_f%bs_fe_f(e)%vn(i, H_N), &
                                                       ms_fe_f%bs_fe_f(e)%vn(i, P_N)
            endif
         enddo
         enddo
      close(k)

   return
   endsubroutine save_profile_x_ms