This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / GEANT321 / block / fdnopt.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1995/10/24 10:19:35  cernlib
6 * Geant
7 *
8 *
9 #include "geant321/pilot.h"
10 *CMZ :  3.21/02 29/03/94  15.41.45  by  S.Giani
11 *-- Author :
12 *$ CREATE BDNOPT.FOR
13 *COPY BDNOPT
14 *
15 *=== bdnopt ===========================================================*
16 *==                                                                    *
17       SUBROUTINE FDNOPT
18  
19 #include "geant321/dblprc.inc"
20 #include "geant321/dimpar.inc"
21 #include "geant321/iounit.inc"
22 *
23 *----------------------------------------------------------------------*
24 *                                                                      *
25 *   Created on  20 september 1989    by  Alfredo Ferrari - Infn Milan  *
26 *                                                                      *
27 *         Last change on 09-mar-93   by  Alfredo Ferrari               *
28 *                                                                      *
29 *----------------------------------------------------------------------*
30 *
31 #include "geant321/fheavy.inc"
32 #include "geant321/finuc.inc"
33 #include "geant321/parevt.inc"
34 #include "geant321/resnuc.inc"
35 *
36       CHARACTER*8 ANHEA1(12)
37       DIMENSION AMHEA1(12),ICHEA1(12),IBHEA1(12)
38       LOGICAL LDIFF1,LPOWE1,LINCT1,LEVPR1,LHEAV1,LDEEX1,LGDHP1
39      &       ,LPREE1,LHLFI1,LPRFI1,LPARW1
40       DIMENSION LDIFF1(39)
41       LOGICAL LRNFS1
42  
43 * /fheavy/
44       DATA AMHEA1 / 12 * 0.D+00 /
45       DATA ANHEA1 / 'NEUTRON ', 'PROTON  ', 'DEUTERON', '3-H     ',
46      &              '3-He    ', '4-He    ', 'H-FRAG-1', 'H-FRAG-2',
47      &              'H-FRAG-3', 'H-FRAG-4', 'H-FRAG-5', 'H-FRAG-6'/
48       DATA ICHEA1 / 0, 1, 1, 1, 2, 2, 6*0 /,
49      &     IBHEA1 / 1, 1, 2, 3, 3, 4, 6*0 /
50       DATA NPHEA1 / 0 /
51  
52 * /finuc/
53       DATA NP1/ 0 /, TV1/ 0.D+00 /, TVCMS1/ 0.D+00 /, TVREC1 / 0.D+00/,
54      &     TVHEA1 / 0.D+00 /, TVBIN1 / 0.D+00 /
55  
56 * /parevt/
57       DATA DPOWE1 /-13.D+00 /, FSPRD1 / 0.6D+00 /, FSHPF1 / 0.0D+00 /,
58      &     RN1GS1 /-1.0D+00 /, RN2GS1 /-1.0D+00 /
59       DATA LDIFF1 /  .TRUE., .FALSE., 6 * .TRUE., .FALSE., 8 * .TRUE.,
60      &              .FALSE., 4 * .TRUE., .TRUE., 2 * .TRUE.,
61      &              8 * .FALSE., .TRUE., .FALSE., .TRUE., .FALSE.,
62      &              .TRUE., .FALSE. /
63       DATA LPOWE1 / .TRUE.  /, LINCT1 / .TRUE.  /, LEVPR1 / .TRUE.  /,
64      &     LHEAV1 / .TRUE.  /, LDEEX1 / .TRUE.  /, LGDHP1 / .TRUE.  /,
65      &     LPREE1 / .TRUE.  /, LHLFI1 / .FALSE. /, LPRFI1 / .FALSE. /,
66      &     LPARW1 / .TRUE. /
67       DATA ILVMO1 / 0 /
68  
69 * /resnuc/
70       DATA IEVPL1 / 0 /, IEVPH1 / 0 /, IEVNE1 / 0 /, IEVPR1 / 0 /,
71      &     IEVTR1 / 0 /, IEVDE1 / 0 /, IEV3H1 / 0 /, IEV4H1 / 0 /,
72      &     IDEEX1 / 0 /
73       DATA LRNFS1 / .FALSE. /
74 *
75       DO 10 J=1,12
76          AMHEAV(J) = AMHEA1(J)
77          ANHEAV(J) = ANHEA1(J)
78          ICHEAV(J) = ICHEA1(J)
79          IBHEAV(J) = IBHEA1(J)
80   10  CONTINUE
81       NPHEAV = NPHEA1
82 *
83       NP=NP1
84       TV=TV1
85       TVCMS=TVCMS1
86       TVRECL=TVREC1
87       TVHEAV=TVHEA1
88       TVBIND=TVBIN1
89 *
90       DPOWER = DPOWE1
91       FSPRD0 = FSPRD1
92       FSHPFN = FSHPF1
93       RN1GSC = RN1GS1
94       RN2GSC = RN2GS1
95       DO 20 J=1,39
96         LDIFFR(J) = LDIFF1(J)
97   20  CONTINUE
98       LPOWER = LPOWE1
99       LINCTV = LINCT1
100       LEVPRT = LEVPR1
101       LHEAVY = LHEAV1
102       LDEEXG = LDEEX1
103       LGDHPR = LGDHP1
104       LPREEX = LPREE1
105       LHLFIX = LHLFI1
106       LPRFIX = LPRFI1
107       LPARWV = LPARW1
108       ILVMOD = ILVMO1
109 *
110       IEVAPL = IEVPL1
111       IEVAPH = IEVPH1
112       IEVNEU = IEVNE1
113       IEVPRO = IEVPR1
114       IEVTRI = IEVTR1
115       IEVDEU = IEVDE1
116       IEV3HE = IEV3H1
117       IEV4HE = IEV4H1
118       IDEEXG = IDEEX1
119       LRNFSS = LRNFS1
120 *=                                               end*block.bdnopt      *
121       END