COSMOS v7.655  COSMOSv7655
(AirShowerMC)
cgetMoliereU.f
Go to the documentation of this file.
1 ! implicit none
2 ! real*8 dep, cosz, rmu
3 ! dep = 9000.
4 ! cosz = 0.5
5 ! call cgetMoliereU(dep, cosz, rmu)
6 ! write(*, *) ' rmu=',rmu
7 ! end
8 ! get moliere unit (in m) at 2 r.l above dep along cosz
9 ! direction.
10  subroutine cgetmoliereu(dep, cosz, rmu)
11  implicit none
12 #include "Zelemagp.h"
13  real*8 dep ! input. depth in kg/m2
14  real*8 cosz ! input. zenith angle
15  real*8 rmu ! output. Moliere length in m
16 !
17  real*8 cthick2den, tmp
18 
19  tmp=max( dep - 2*x0*cosz, 100.d0) ! 2 r%l above
20  rmu = es/ecrit * x0 / cthick2den(tmp)
21  end
22 
subroutine cgetmoliereu(dep, cosz, rmu)
Definition: cgetMoliereU.f:11
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
Definition: cblkEvhnp.h:5