source: JEM-EUSO/esaf_cc_at_lal/packages/simulation/externals/slast/src/fluor/estoping_power.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: 746 bytes
Line 
1      REAL FUNCTION ESTOPING_POWER(E)
2      IMPLICIT NONE
3      REAL e,e_data(81),s_data(81),HX
4      INTEGER iopen,i
5      SAVE iopen,e_data,s_data
6      DATA iopen/0/
7*
8 3    CONTINUE
9      IF(iopen.EQ.0) THEN
10        iopen = 1
11        i = 0
12        OPEN(66,file='config/slast/stoping_power.data',STATUS='OLD')
13 1      i = i + 1
14        IF (i.GT.81) GOTO 2
15        READ(66,*,end=2) e_data(i),s_data(i)
16        GOTO 1
17 2      CLOSE(66)
18        GOTO 3
19      ELSE
20        IF(e.LT.e_data(1))  ESTOPING_POWER = s_data(1)
21        IF(e.GT.e_data(81)) ESTOPING_POWER = s_data(81)
22        DO i = 1, 80
23          IF(e_data(i).LE.E.AND.E.LT.e_data(i+1)) THEN
24            ESTOPING_POWER = s_data(i)
25          ENDIF
26        ENDDO
27      ENDIF   
28*
29      END
Note: See TracBrowser for help on using the repository browser.