source: JEM-EUSO/esaf_lal/tags/v1_r0/esaf/packages/simulation/externals/slast/src/init/depth_init.F @ 117

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

ESAF version compilable on mac OS
File size: 1.7 KB
Line 
1      SUBROUTINE DEPTH_INIT
2#include "detector.inc"
3#include "book_histos.inc"
4      CHARACTER*12 CHTAG1(2)/'ISOTHERMIC',
5     +     'US STANDARD'/
6      CHARACTER*20 CHTAG2(2)/'SPHERICAL ATMOSPHERE','PLANAR ATMOSPHERE'/
7      IF(DEPT.EQ.1) THEN
8         CALL BOOK_HISTOS('DUMP')
9*     ... CHECK ATMOSHERE FUNCTION ...
10         IATMO = ATMO
11         DO L = 1, 2
12            ID = 100*L
13            ATMO = L
14            DO h = 0., 100., 0.1
15               CALL HFILL(ID,h,1.e-4*atm(h),1.) ! in g/cm^3
16            ENDDO
17         ENDDO
18         ATMO = IATMO
19
20*     ... COMPUTE DEPTH'S INTEGRALS AND SAVE THEM IN AUX.HBK
21
22         WRITE(6,*) ' '
23         WRITE(6,98)
24         step_h    = Hmax/Nh       ! km
25         Nint      = Nh*Nth
26         i_portion = Nint/10
27         CURV0     = CURV
28         DO CURV = 1, 2
29            DO L = 1, 2
30               k    = 0
31               ID   = 100*L
32               ATMO = L
33               WRITE(6,99) CHTAG1(L),CHTAG2(CURV)
34               DO i = 0, Nh-1
35                  h = i*step_h + step_h/2
36                  DO j = 0, Nth -1
37                     th = j*thmax/Nth + thmax/Nth/2
38                     k = k + 1
39                     depth = DEPTH_INTEGRAL(th,h)
40                     CALL HFILL(ID+int(CURV),th,h,depth)
41                     ntmp = k/i_portion
42                     ntmp = k - ntmp*i_portion
43                     IF(ntmp.EQ.0) WRITE(6,100)
44                  ENDDO
45               ENDDO
46               WRITE(6,101)
47               WRITE(6,*) ' '
48            ENDDO
49         ENDDO
50         WRITE(6,*) ' '
51         CURV = CURV0
52         ATMO = IATMO
53      ENDIF
54 98   FORMAT(1x,'COMPUTING THE DEPTHS INTEGRALS: ')
55 99   FORMAT($,1x,'FOR ',2A,1x,'(breath deeply)')
56 100  FORMAT($,1x,'.')
57 101  FORMAT($,1x,'[done]')
58      END
Note: See TracBrowser for help on using the repository browser.