5 * Revision 1.1.1.1 1995/10/24 10:21:55 cernlib
9 #include "geant321/pilot.h"
10 *CMZ : 3.21/04 23/02/95 14.46.01 by S.Giani
12 SUBROUTINE CEVAP(E,Q,ATAR,CB,EX)
13 C THIS ROUTINE SAMPLES AN EXIT ENERGY FROM AN
14 C EVAPORATION SPECTRUM
15 #include "geant321/minput.inc"
17 C CONVERT THE COULOMB BARRIER (CB) TO UNITS OF EV
19 C CALCULATE THE MAXIMUM ENERGY AVAILABLE
22 IF(EMAX.GT.0.0)GO TO 10
25 IF(EMAX.GT.0.0)GO TO 10
28 IF(EMAX.GT.0.0)GO TO 10
29 WRITE(IOUT,10000)E,EMAX,Q,CBI
30 10000 FORMAT(' MICAP: NEGATIVE MAXIMUM ENERGY CALCULATED IN ROUTINE ',
31 1'EVAP --- INDICATING PROBABLE CROSS SECTION ERROR ALLOWING ',
32 2'THE REACTION TO OCCUR',/,10X,'E,EMAX,Q,CB=',4E13.5)
33 WRITE(6,*) ' CALOR: Fatal ERROR in EVAP ====> STOP '
35 C CALCULATE THE NUCLEAR TEMPERATURE (THETA)
36 10 THETA=4.0161E+03*(SQRT(E+Q-CB)/(ATAR**0.8333333))
37 C SELECT THE EXIT ENERGY FROM AN EVAPORATION SPECTRUM
42 IF(EX.LT.0.0) EX = 0.0
44 C RESAMPLE 75% OF THE TIME IF EX IS GREATER THAN EMAX
47 #if defined(CERNLIB_MDEBUG)
48 IF(EX.GT.20.*EMAX) WRITE(IOUT,10100)EX,EMAX
49 10100 FORMAT(' MICAP: WARNING EX > EMAX=',2E13.5,' IN ROUTINE EVAP')