profil_theo_trie_1D Subroutine

subroutine profil_theo_trie_1D(tab, lg, x, mx)

Arguments

Type IntentOptional Attributes Name
real(kind=R8), intent(out), dimension(1:lg) :: tab
integer(kind=I4), intent(in) :: lg
real(kind=R8), intent(in), dimension(1:2) :: x
type(moment_stat), intent(out) :: mx

Calls

proc~~profil_theo_trie_1d~~CallsGraph proc~profil_theo_trie_1d profil_theo_trie_1D proc~calc_moments_1d calc_moments_1D proc~profil_theo_trie_1d->proc~calc_moments_1d

Called by

proc~~profil_theo_trie_1d~~CalledByGraph proc~profil_theo_trie_1d profil_theo_trie_1D program~test_algen test_algen program~test_algen->proc~profil_theo_trie_1d

Source Code

   subroutine profil_theo_trie_1D(tab, lg, x, mx)
   implicit none
   integer(kind=I4) , intent(in )                  :: lg
   real(kind=R8)    , intent(out), dimension(1:lg) :: tab
   real(kind=R8)    , intent(in ), dimension(1:2)  :: x
   type(moment_stat), intent(out)                  :: mx

      real(kind=R8)    :: b1, b2, alp, bet
      integer(kind=I4) :: i

      b1 = -PI_R8/2 *(1. -x(1))
      b2 = +PI_R8/2 *(1. -x(2))
      alp = -(b2- lg*b1)/(b2 -b1)
      bet =      (lg- 1)/(b2 -b1)
      do i = 1, lg
         tab(i) = tan( (i +alp)/bet )
      enddo
      call calc_moments_1D(tab, mx, nb_mom=4, lg=lg)
      tab(1:lg) = (tab(1:lg) -mx%mu)/mx%si

      mx%mu = 0._R8
      mx%si = 1._R8

   return
   endsubroutine profil_theo_trie_1D