source: JEM-EUSO/esaf_cc_at_lal/packages/simulation/externals/slast/src/init/mie_init.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: 1.1 KB
Line 
1      SUBROUTINE MIE_INIT
2      IMPLICIT NONE
3#include "detector.inc"
4#include "book_histos.inc"
5      REAL MIE_INTEGRAL
6      REAL step_h,h,th
7      INTEGER i_portion,i,j,k,L,ntmp,Nint,ID
8*
9      IF(IMIE.EQ.1) THEN
10         CALL BOOK_HISTOS('DUMP')
11
12*     ... COMPUTE MIE'S INTEGRALS AND SAVE THEM IN AUX.HBK
13
14         WRITE(6,*) ' '
15         WRITE(6,98)
16         WRITE(6,99)
17         step_h = Hmax/Nh       ! km
18         Nint = Nh*Nth
19         i_portion = Nint/10
20         DO L = 1, 1
21            k = 0
22            ID = 10000*L
23            DO i = 0, Nh-1
24               h = i*step_h + step_h/2
25               DO j = 0, Nth -1
26                  th = j*thmax/Nth + thmax/Nth/2
27                  k = k + 1
28                  CALL HFILL(ID,th,h,MIE_INTEGRAL(th,h))
29                  ntmp = k/i_portion
30                  ntmp = k - ntmp*i_portion
31                  IF(ntmp.EQ.0) WRITE(6,100)
32               ENDDO
33            ENDDO
34            WRITE(6,101)
35            WRITE(6,*) ' '
36         ENDDO
37         WRITE(6,*) ' '
38      ENDIF
39 98   FORMAT(1x,'COMPUTING THE MIE INTEGRALS: ')
40 99   FORMAT($,52x,' ')
41 100  FORMAT($,1x,'.')
42 101  FORMAT($,1x,'[done]')
43      END
Note: See TracBrowser for help on using the repository browser.