save_profile_x_fe Subroutine

  1. mod_inout_files.f90
  2. inout_files
  3. save_profile_x_fe

public subroutine save_profile_x_fe(fe_f, file_name, lx, zy)

Arguments

Type IntentOptional AttributesName
type(FE_FILM), intent(in) :: fe_f

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_fe~~CallsGraph proc~save_profile_x_fe save_profile_x_fe proc~get_unit get_unit proc~save_profile_x_fe->proc~get_unit

Contents

Source Code

save_profile_x_fe

Source Code

   subroutine save_profile_x_fe(fe_f, file_name, lx, zy)
   implicit none
   type(FE_FILM),    intent(in) :: fe_f         !! [[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, k

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

   return
   endsubroutine save_profile_x_fe