source: JEM-EUSO/esaf_lal/tags/v1_r0/esaf/packages/simulation/externals/slast/src/init/fluor_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.5 KB
Line 
1      SUBROUTINE FLUOR_INIT
2#include "event.inc"
3#include "shower.inc"
4#include "detector.inc"
5#include "paw.inc"
6      LOGICAL DEBUG/.FALSE./
7      CHARACTER*22 CHTAG(1)/'FLUORESCENT YIELDS    '/
8      LOGICAL HEXIST
9      IF(IFLU.NE.1) RETURN
10      Ns   = 200
11      Smin = 1.e-5
12      Smax = 2.0
13      step_s = (Smax-Smin)/(Ns-1)
14      Nint = Ns
15      i_portion = Nint/10
16*
17      CALL BOOK_HISTOS('DUMP')
18*
19      WRITE(6,*) ' '
20      WRITE(6,98)
21      WRITE(6,99)
22*
23      l = 0
24      ID = 100000
25      DO j = 1, Ns
26        s = Smin + (j-1)*step_s
27        l = l + 1
28        CALL HFILL(ID,s,FLUOR_ENERGY(s,1),1.)
29        CALL HFILL(ID+1,s,FLUOR_ENERGY(s,2),1.)
30        ntmp = l/i_portion
31        ntmp = l - ntmp*i_portion
32        IF(ntmp.EQ.0) WRITE(6,100)
33      ENDDO
34      WRITE(6,101)
35      WRITE(6,*) ' '
36*
37      WRITE(6,99)
38      Nh   = 500
39      Hmax = 1.e2
40      step_h = Hmax/Nh
41      Ns = 20
42      Smin = 0.1
43      Smax = 2.0
44      step_s = (Smax-Smin)/(Ns-1)
45      i_portion = Nh*(Ns-1)/10
46      l = 0
47*
48      DO j = 1, Ns-1
49        s = Smin + (j-1)*step_s + step_s/2
50        DO i = 0, Nh-1
51          h = i*step_h + step_h/2
52          l = l + 1
53          CALL HFILL(ID+2,s,h,FLUOR_YIELD(s,h))
54          ntmp = l/i_portion
55          ntmp = l - ntmp*i_portion
56          IF(ntmp.EQ.0) WRITE(6,100)
57        ENDDO
58      ENDDO
59      WRITE(6,101)
60      WRITE(6,*) ' '
61*
62      IF(DEBUG) CALL RANDOM_OUTPUT
63 98   FORMAT(1x,'COMPUTING THE FLUORESCENT YIELDS INTEGRALS: ')
64 99   FORMAT($,32x,' ')
65 100  FORMAT($,1x,'.')
66 101  FORMAT($,1x,'[done]')
67      END
68     
Note: See TracBrowser for help on using the repository browser.