cgetMoliereU.f

Go to the documentation of this file.

!      implicit none
!      real*8 dep, cosz, rmu
!      dep = 9000.
!      cosz = 0.5
!      call cgetMoliereU(dep, cosz, rmu)
!      write(*, *) ' rmu=',rmu
!      end
!            get moliere unit (in m) at 2 r.l above dep along cosz
!            direction.
       subroutine cgetmoliereu(dep, cosz, rmu)
       implicit none
#include "Zelemagp.h"
       real*8 dep  ! input.  depth in kg/m2
       real*8 cosz ! input.  zenith angle
       real*8 rmu  ! output. Moliere length in m
!
       real*8  cthick2den, tmp

       tmp=max( dep - 2*x0*cosz, 100.d0)    ! 2 r%l above
       rmu = es/ecrit * x0 / cthick2den(tmp)
       end