source: JEM-EUSO/esaf_cc_at_lal/packages/simulation/externals/slast/src/atm/depth.F @ 114

Last change on this file since 114 was 114, checked in by moretto, 11 years ago

actual version of ESAF at CCin2p3
File size: 881 bytes
Line 
1*
2*     ... THIS FUNCTION RETURNS THE DEPTH IN THE ATMOSPHERE ...
3*     ... BETWEEN THE POINT WITH RADIUS_VECTOR R(3) AND     ...
4*     ... INFINITY  ALONG THE INCOMING DIRECTION            ...
5*     ... DATE: 09.06.2002
6*     ... AUTHOR: D.V. NAUMOV
7*
8      REAL FUNCTION DEPTH(R)
9      IMPLICIT NONE
10#include "detector.inc"
11#include "event.inc"
12      INTEGER i_aux_hbk, ID
13      REAL R(3), PI, HALF_PI, Hstar, H_INT,THETA_INT,H_R,THETA_R
14      REAL DEPTH_R, DEPTH_INT, DEPTH_INTEGRAL
15*
16      PI      = ACOS(-1.)
17      HALF_PI = PI/2
18*     
19      CALL IMPACT(R,Hstar,H_INT,THETA_INT,H_R,THETA_R)
20      IF(THETA_R.LE.HALF_PI) THEN
21        DEPTH_R = DEPTH_INTEGRAL(THETA_R,H_R)
22      ELSE       
23         THETA_R = PI - THETA_R
24         DEPTH_R = 2*DEPTH_INTEGRAL(HALF_PI,Hstar)-
25     &        DEPTH_INTEGRAL(H_R,THETA_R)
26      ENDIF
27      DEPTH = DEPTH_R
28 1000 CONTINUE
29      END
Note: See TracBrowser for help on using the repository browser.