| 1 | DOUBLE PRECISION FUNCTION CGHSIG( PPART,EKIN,LPART )
|
|---|
| 2 |
|
|---|
| 3 | C-----------------------------------------------------------------------
|
|---|
| 4 | C C(ORSIKA) GH(EISHA) SIG(MA)
|
|---|
| 5 | C
|
|---|
| 6 | C CALCULATION OF THE PROBABILITIES FOR (IN)ELASTIC INTERACTIONS ***
|
|---|
| 7 | C THIS FUNCTION IS CALLED FROM BOX2
|
|---|
| 8 | C
|
|---|
| 9 | C ORIGIN : F.CARMINATI, H.FESEFELDT (ROUTINE GHESIG)
|
|---|
| 10 | C REDESIGN: P. GABRIEL IK1 FZK KARLSRUHE
|
|---|
| 11 | C-----------------------------------------------------------------------
|
|---|
| 12 | C *** IPART DENOTES THE GHEISHA PARTICLE INDEX ***
|
|---|
| 13 | C
|
|---|
| 14 | C CONVENTION :
|
|---|
| 15 | C
|
|---|
| 16 | C PARTICLE IPART
|
|---|
| 17 | C ------------------------------
|
|---|
| 18 | C GAMMA 1
|
|---|
| 19 | C NEUTRINO 2
|
|---|
| 20 | C POSITRON 3
|
|---|
| 21 | C ELECTRON 4
|
|---|
| 22 | C MUON + 5
|
|---|
| 23 | C MUON - 6
|
|---|
| 24 | C PION + 7
|
|---|
| 25 | C PION 0 8
|
|---|
| 26 | C PION - 9
|
|---|
| 27 | C KAON + 10
|
|---|
| 28 | C KAON 0 S 11
|
|---|
| 29 | C KAON 0 L 12
|
|---|
| 30 | C KAON - 13
|
|---|
| 31 | C PROTON 14
|
|---|
| 32 | C PROTON BAR 15
|
|---|
| 33 | C NEUTRON 16
|
|---|
| 34 | C NEUTRON BAR 17
|
|---|
| 35 | C LAMBDA 18
|
|---|
| 36 | C LAMBDA BAR 19
|
|---|
| 37 | C SIGMA + 20
|
|---|
| 38 | C SIGMA 0 21
|
|---|
| 39 | C SIGMA - 22
|
|---|
| 40 | C SIGMA + BAR 23
|
|---|
| 41 | C SIGMA 0 BAR 24
|
|---|
| 42 | C SIGMA - BAR 25
|
|---|
| 43 | C XSI 0 26
|
|---|
| 44 | C XSI - 27
|
|---|
| 45 | C XSI 0 BAR 28
|
|---|
| 46 | C XSI - BAR 29
|
|---|
| 47 | C DEUTERON 30
|
|---|
| 48 | C TRITON 31
|
|---|
| 49 | C ALPHA 32
|
|---|
| 50 | C OMEGA - 33
|
|---|
| 51 | C OMEGA - BAR 34
|
|---|
| 52 | C NEW PARTICLES 35
|
|---|
| 53 | C
|
|---|
| 54 | C-----------------------------------------------------------------------
|
|---|
| 55 |
|
|---|
| 56 | *KEEP,CGCOMP.
|
|---|
| 57 | PARAMETER (KK=3)
|
|---|
| 58 | COMMON/CGCOMP/ ACOMP,ZCOMP,WCOMP
|
|---|
| 59 | REAL ACOMP(KK),ZCOMP(KK),WCOMP(KK)
|
|---|
| 60 | *KEEP,RUNPAR.
|
|---|
| 61 | COMMON /RUNPAR/ FIXHEI,THICK0,HILOECM,HILOELB,
|
|---|
| 62 | * STEPFC,NRRUN,NSHOW,PATAPE,MONIIN,
|
|---|
| 63 | * MONIOU,MDEBUG,NUCNUC,
|
|---|
| 64 | * CETAPE,
|
|---|
| 65 | * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
|
|---|
| 66 | * N1STTR,MDBASE,
|
|---|
| 67 | * DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
|
|---|
| 68 | * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
|
|---|
| 69 | * ,GHEISH,GHESIG
|
|---|
| 70 | COMMON /RUNPAC/ DSN,HOST,USER
|
|---|
| 71 | DOUBLE PRECISION FIXHEI,THICK0,HILOECM,HILOELB
|
|---|
| 72 | REAL STEPFC
|
|---|
| 73 | INTEGER NRRUN,NSHOW,PATAPE,MONIIN,MONIOU,MDEBUG,NUCNUC,
|
|---|
| 74 | * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
|
|---|
| 75 | * N1STTR,MDBASE
|
|---|
| 76 | INTEGER CETAPE
|
|---|
| 77 | CHARACTER*79 DSN
|
|---|
| 78 | CHARACTER*20 HOST,USER
|
|---|
| 79 |
|
|---|
| 80 | LOGICAL DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
|
|---|
| 81 | * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
|
|---|
| 82 | * ,GHEISH,GHESIG
|
|---|
| 83 | *KEND.
|
|---|
| 84 |
|
|---|
| 85 | COMMON/GSECTI/ AIEL(20),AIIN(20),AIFI(20),AICA(20),ALAM,K0FLAG
|
|---|
| 86 | INTEGER K0FLAG
|
|---|
| 87 | REAL AIEL,AIIN,AIFI,AICA,ALAM
|
|---|
| 88 |
|
|---|
| 89 | C --- GHEISHA COMMONS ---
|
|---|
| 90 | COMMON/RESULT/XEND,YEND,ZEND,RCA,RCE,AMAS,NCH,TOF,PX,PY,PZ,
|
|---|
| 91 | * USERW,INTCT,P,EN,EK,AMASQ,DELTN,ITK,NTK,IPART,IND,
|
|---|
| 92 | * LCALO,ICEL,SINL,COSL,SINP,COSP,
|
|---|
| 93 | * XOLD,YOLD,ZOLD,POLD,PXOLD,PYOLD,PZOLD,
|
|---|
| 94 | * XSCAT,YSCAT,ZSCAT,PSCAT,PXSCAT,PYSCAT,PZSCAT
|
|---|
| 95 | REAL NCH,INTCT
|
|---|
| 96 |
|
|---|
| 97 | COMMON/PRNTFL/INBCD,NEWBCD,INBIN,NEWBIN,NPEVT,NEVTP,LPRT,NPRT(10)
|
|---|
| 98 | LOGICAL LPRT,NPRT
|
|---|
| 99 |
|
|---|
| 100 |
|
|---|
| 101 | DIMENSION ALPHA(35),ALPHAC(41),IPART2(7),CSA(4)
|
|---|
| 102 | DIMENSION PARTEL(35),PARTIN(35),ICORR(35),INTRC(35)
|
|---|
| 103 |
|
|---|
| 104 | C --- DIMENSION STATEMENTS FOR CROSS SECTION DATA ---
|
|---|
| 105 | DIMENSION PLAB(41),CSEL(35,41),CSIN(35,41),CSPIEL(3,41),
|
|---|
| 106 | $ CSPIIN(3,41),CSPNEL(3,41),CSPNIN(3,41),
|
|---|
| 107 | $ ELAB(17),CNLWAT(15),CNLWEL(15,17),CNLWIN(15,17),
|
|---|
| 108 | $ CSCAP(100)
|
|---|
| 109 |
|
|---|
| 110 | C --- DIMENSION STMTS. FOR GEANT/GHEISHA PARTICLE CODE CONVERSIONS ---
|
|---|
| 111 | C --- KIPART(I)=GHEISHA CODE CORRESPONDING TO GEANT CODE I ---
|
|---|
| 112 | C --- IKPART(I)=GEANT CODE CORRESPONDING TO GHEISHA CODE I ---
|
|---|
| 113 |
|
|---|
| 114 | DIMENSION KIPART(48),IKPART(35)
|
|---|
| 115 |
|
|---|
| 116 | C --- CROSS-SECTION DATA BY "PCSDAT" 01-FEB-1989 ---
|
|---|
| 117 | DATA PLAB /
|
|---|
| 118 | $ 0.00000E+00, 0.10000 , 0.15000 , 0.20000 , 0.25000 ,
|
|---|
| 119 | $ 0.30000 , 0.35000 , 0.40000 , 0.45000 , 0.50000 ,
|
|---|
| 120 | $ 0.55000 , 0.60000 , 0.65000 , 0.70000 , 0.75000 ,
|
|---|
| 121 | $ 0.80000 , 0.85000 , 0.90000 , 0.95000 , 1.0000 ,
|
|---|
| 122 | $ 1.1000 , 1.2000 , 1.3000 , 1.4000 , 1.5000 ,
|
|---|
| 123 | $ 1.6000 , 1.8000 , 2.0000 , 2.2000 , 2.4000 ,
|
|---|
| 124 | $ 2.6000 , 2.8000 , 3.0000 , 4.0000 , 5.0000 ,
|
|---|
| 125 | $ 6.0000 , 8.0000 , 10.000 , 20.000 , 100.00 ,
|
|---|
| 126 | $ 1000.0 /
|
|---|
| 127 |
|
|---|
| 128 | DATA ((CSEL(I,J),I=1,6),J=1,41) / 246 * 0. /
|
|---|
| 129 | DATA (CSEL( 8,J),J=1,41) / 41 * 0. /
|
|---|
| 130 | DATA (CSEL(21,J),J=1,41) / 41 * 0. /
|
|---|
| 131 | DATA (CSEL(24,J),J=1,41) / 41 * 0. /
|
|---|
| 132 | DATA ((CSEL(I,J),I=30,32),J=1,41) / 123 * 0. /
|
|---|
| 133 | DATA (CSEL(35,J),J=1,41) / 41 * 0. /
|
|---|
| 134 |
|
|---|
| 135 | DATA (CSEL( 7,J),J=1,41) /
|
|---|
| 136 | $ 0.00000E+00, 6.0000 , 20.000 , 71.000 , 155.00 ,
|
|---|
| 137 | $ 195.00 , 130.00 , 78.000 , 60.000 , 32.000 ,
|
|---|
| 138 | $ 23.500 , 18.500 , 15.000 , 12.500 , 10.000 ,
|
|---|
| 139 | $ 9.1000 , 8.6000 , 8.8000 , 9.5000 , 10.600 ,
|
|---|
| 140 | $ 13.000 , 15.500 , 17.100 , 17.200 , 16.200 ,
|
|---|
| 141 | $ 15.000 , 12.300 , 10.200 , 9.0000 , 8.0000 ,
|
|---|
| 142 | $ 7.3000 , 6.8000 , 6.5000 , 5.8000 , 5.4000 ,
|
|---|
| 143 | $ 5.2000 , 5.0000 , 4.9000 , 3.8000 , 3.2000 ,
|
|---|
| 144 | $ 3.5000 /
|
|---|
| 145 | DATA (CSEL( 9,J),J=1,41) /
|
|---|
| 146 | $ 0.00000E+00, 1.0000 , 3.0000 , 8.0000 , 18.000 ,
|
|---|
| 147 | $ 25.000 , 27.500 , 12.300 , 10.600 , 11.000 ,
|
|---|
| 148 | $ 12.500 , 14.500 , 17.000 , 19.400 , 19.800 ,
|
|---|
| 149 | $ 16.800 , 14.000 , 14.800 , 20.000 , 26.100 ,
|
|---|
| 150 | $ 19.500 , 15.000 , 12.800 , 11.500 , 10.500 ,
|
|---|
| 151 | $ 9.8000 , 8.8000 , 8.2000 , 7.8000 , 7.5000 ,
|
|---|
| 152 | $ 7.2000 , 7.0000 , 6.8000 , 6.1000 , 5.7000 ,
|
|---|
| 153 | $ 5.4000 , 4.9000 , 4.6000 , 4.0000 , 3.3000 ,
|
|---|
| 154 | $ 3.5000 /
|
|---|
| 155 | DATA (CSEL(10,J),J=1,41) /
|
|---|
| 156 | $ 10.000 , 11.200 , 11.300 , 11.400 , 11.500 ,
|
|---|
| 157 | $ 11.600 , 11.800 , 12.000 , 12.100 , 12.200 ,
|
|---|
| 158 | $ 12.300 , 12.400 , 12.500 , 12.500 , 12.500 ,
|
|---|
| 159 | $ 12.400 , 12.300 , 12.200 , 12.000 , 11.800 ,
|
|---|
| 160 | $ 11.200 , 11.500 , 9.9000 , 9.4000 , 8.8000 ,
|
|---|
| 161 | $ 8.4000 , 7.5000 , 6.9000 , 6.3000 , 5.9000 ,
|
|---|
| 162 | $ 5.5000 , 5.2000 , 5.0000 , 4.0000 , 3.5000 ,
|
|---|
| 163 | $ 3.3000 , 3.1000 , 3.1000 , 3.0000 , 2.5000 ,
|
|---|
| 164 | $ 3.0000 /
|
|---|
| 165 | DATA (CSEL(11,J),J=1,41) /
|
|---|
| 166 | $ 10.000 , 11.200 , 11.300 , 11.400 , 11.500 ,
|
|---|
| 167 | $ 11.600 , 11.800 , 12.000 , 12.100 , 12.200 ,
|
|---|
| 168 | $ 12.300 , 12.400 , 12.500 , 12.500 , 12.500 ,
|
|---|
| 169 | $ 12.400 , 12.300 , 12.200 , 12.000 , 11.800 ,
|
|---|
| 170 | $ 11.200 , 11.500 , 9.9000 , 9.4000 , 8.8000 ,
|
|---|
| 171 | $ 8.4000 , 7.5000 , 6.9000 , 6.3000 , 5.9000 ,
|
|---|
| 172 | $ 5.5000 , 5.2000 , 5.0000 , 4.0000 , 3.5000 ,
|
|---|
| 173 | $ 3.3000 , 3.1000 , 3.1000 , 3.0000 , 2.5000 ,
|
|---|
| 174 | $ 3.0000 /
|
|---|
| 175 | DATA (CSEL(12,J),J=1,41) /
|
|---|
| 176 | $ 160.83 , 82.800 , 58.575 , 43.683 , 34.792 ,
|
|---|
| 177 | $ 28.650 , 24.367 , 20.917 , 18.192 , 16.300 ,
|
|---|
| 178 | $ 14.608 , 13.017 , 12.250 , 11.700 , 12.017 ,
|
|---|
| 179 | $ 14.075 , 15.842 , 16.433 , 16.042 , 15.008 ,
|
|---|
| 180 | $ 12.575 , 10.708 , 9.2000 , 8.0167 , 7.2833 ,
|
|---|
| 181 | $ 7.0750 , 6.6333 , 6.1250 , 5.6583 , 5.2750 ,
|
|---|
| 182 | $ 4.9333 , 4.6250 , 4.4583 , 3.7333 , 3.3833 ,
|
|---|
| 183 | $ 3.1833 , 2.9833 , 2.7500 , 2.3667 , 2.2000 ,
|
|---|
| 184 | $ 2.6000 /
|
|---|
| 185 | DATA (CSEL(13,J),J=1,41) /
|
|---|
| 186 | $ 300.00 , 140.00 , 97.000 , 70.000 , 55.000 ,
|
|---|
| 187 | $ 45.000 , 37.000 , 31.000 , 26.000 , 23.000 ,
|
|---|
| 188 | $ 20.000 , 17.000 , 15.500 , 14.500 , 14.700 ,
|
|---|
| 189 | $ 18.500 , 22.000 , 23.000 , 22.500 , 20.700 ,
|
|---|
| 190 | $ 16.500 , 14.000 , 11.500 , 9.6000 , 8.6000 ,
|
|---|
| 191 | $ 8.5000 , 8.3000 , 7.6000 , 7.0000 , 6.4000 ,
|
|---|
| 192 | $ 5.9000 , 5.5000 , 5.3000 , 4.4000 , 4.1000 ,
|
|---|
| 193 | $ 3.9000 , 3.7000 , 3.3000 , 2.6000 , 2.5000 ,
|
|---|
| 194 | $ 3.0000 /
|
|---|
| 195 | DATA (CSEL(14,J),J=1,41) /
|
|---|
| 196 | $ 1100.0 , 115.00 , 105.00 , 100.00 , 56.000 ,
|
|---|
| 197 | $ 40.000 , 27.000 , 22.000 , 21.000 , 20.000 ,
|
|---|
| 198 | $ 20.000 , 20.000 , 20.500 , 21.000 , 22.000 ,
|
|---|
| 199 | $ 23.000 , 24.000 , 24.000 , 24.400 , 24.500 ,
|
|---|
| 200 | $ 25.000 , 25.500 , 26.000 , 26.500 , 27.000 ,
|
|---|
| 201 | $ 27.000 , 26.000 , 23.000 , 21.500 , 20.000 ,
|
|---|
| 202 | $ 19.000 , 18.000 , 17.000 , 13.000 , 11.500 ,
|
|---|
| 203 | $ 10.300 , 9.4000 , 9.0000 , 8.8000 , 7.0000 ,
|
|---|
| 204 | $ 7.5000 /
|
|---|
| 205 | DATA (CSEL(15,J),J=1,41) /
|
|---|
| 206 | $ 200.00 , 163.00 , 141.00 , 120.00 , 111.00 ,
|
|---|
| 207 | $ 99.500 , 92.500 , 86.500 , 82.000 , 78.000 ,
|
|---|
| 208 | $ 74.000 , 71.000 , 67.500 , 65.000 , 62.500 ,
|
|---|
| 209 | $ 59.700 , 58.100 , 56.300 , 54.700 , 52.700 ,
|
|---|
| 210 | $ 50.000 , 48.400 , 47.000 , 46.000 , 45.200 ,
|
|---|
| 211 | $ 42.800 , 39.200 , 36.300 , 32.800 , 30.400 ,
|
|---|
| 212 | $ 28.100 , 26.300 , 24.500 , 19.250 , 16.840 ,
|
|---|
| 213 | $ 14.600 , 12.340 , 11.210 , 8.8500 , 7.5000 ,
|
|---|
| 214 | $ 7.5000 /
|
|---|
| 215 | DATA (CSEL(16,J),J=1,41) /
|
|---|
| 216 | $ 4200.0 , 440.00 , 420.00 , 400.00 , 230.00 ,
|
|---|
| 217 | $ 160.00 , 105.00 , 80.000 , 62.000 , 50.000 ,
|
|---|
| 218 | $ 45.000 , 41.000 , 38.000 , 36.000 , 35.000 ,
|
|---|
| 219 | $ 34.000 , 33.000 , 32.000 , 31.500 , 31.000 ,
|
|---|
| 220 | $ 30.500 , 30.000 , 29.500 , 29.000 , 28.500 ,
|
|---|
| 221 | $ 28.000 , 26.000 , 23.000 , 21.500 , 20.000 ,
|
|---|
| 222 | $ 19.000 , 18.000 , 17.000 , 13.000 , 11.500 ,
|
|---|
| 223 | $ 10.300 , 9.4000 , 9.0000 , 8.8000 , 7.0000 ,
|
|---|
| 224 | $ 7.5000 /
|
|---|
| 225 | DATA (CSEL(17,J),J=1,41) /
|
|---|
| 226 | $ 185.88 , 133.23 , 119.37 , 102.86 , 93.102 ,
|
|---|
| 227 | $ 82.752 , 76.205 , 71.008 , 67.366 , 64.096 ,
|
|---|
| 228 | $ 60.891 , 58.501 , 55.735 , 53.773 , 51.839 ,
|
|---|
| 229 | $ 49.671 , 48.485 , 47.045 , 45.803 , 44.306 ,
|
|---|
| 230 | $ 42.623 , 41.786 , 41.115 , 40.630 , 40.129 ,
|
|---|
| 231 | $ 38.242 , 35.233 , 32.662 , 29.639 , 27.573 ,
|
|---|
| 232 | $ 25.536 , 23.948 , 22.356 , 17.723 , 15.614 ,
|
|---|
| 233 | $ 13.653 , 11.675 , 10.653 , 8.6198 , 7.4464 ,
|
|---|
| 234 | $ 7.4821 /
|
|---|
| 235 | DATA (CSEL(18,J),J=1,41) /
|
|---|
| 236 | $ 1100.0 , 115.00 , 105.00 , 100.00 , 56.000 ,
|
|---|
| 237 | $ 40.000 , 27.000 , 22.000 , 21.000 , 20.000 ,
|
|---|
| 238 | $ 20.000 , 19.067 , 19.333 , 19.500 , 19.833 ,
|
|---|
| 239 | $ 20.567 , 21.800 , 22.900 , 23.869 , 23.809 ,
|
|---|
| 240 | $ 22.161 , 21.488 , 19.732 , 19.433 , 19.345 ,
|
|---|
| 241 | $ 19.029 , 18.121 , 16.280 , 15.258 , 14.280 ,
|
|---|
| 242 | $ 13.644 , 12.963 , 12.316 , 9.5333 , 8.4333 ,
|
|---|
| 243 | $ 7.5728 , 6.9696 , 6.7518 , 6.6175 , 5.6000 ,
|
|---|
| 244 | $ 6.1145 /
|
|---|
| 245 | DATA (CSEL(19,J),J=1,41) /
|
|---|
| 246 | $ 157.65 , 73.701 , 76.096 , 68.571 , 57.305 ,
|
|---|
| 247 | $ 49.257 , 43.616 , 40.024 , 38.098 , 36.287 ,
|
|---|
| 248 | $ 34.674 , 33.105 , 31.712 , 30.685 , 29.613 ,
|
|---|
| 249 | $ 28.602 , 28.336 , 28.075 , 27.786 , 27.215 ,
|
|---|
| 250 | $ 26.380 , 26.146 , 25.108 , 24.783 , 24.360 ,
|
|---|
| 251 | $ 23.219 , 21.431 , 20.095 , 18.382 , 17.267 ,
|
|---|
| 252 | $ 16.100 , 15.175 , 14.271 , 11.573 , 10.305 ,
|
|---|
| 253 | $ 9.1471 , 8.0149 , 7.4349 , 6.2499 , 5.8928 ,
|
|---|
| 254 | $ 6.0774 /
|
|---|
| 255 | DATA (CSEL(20,J),J=1,41) /
|
|---|
| 256 | $ 1100.0 , 115.00 , 105.00 , 100.00 , 56.000 ,
|
|---|
| 257 | $ 40.000 , 27.000 , 22.000 , 21.000 , 20.000 ,
|
|---|
| 258 | $ 20.000 , 19.067 , 19.333 , 19.500 , 19.833 ,
|
|---|
| 259 | $ 20.567 , 21.800 , 22.900 , 23.869 , 23.809 ,
|
|---|
| 260 | $ 22.161 , 21.488 , 19.732 , 19.433 , 19.345 ,
|
|---|
| 261 | $ 19.029 , 18.121 , 16.280 , 15.258 , 14.280 ,
|
|---|
| 262 | $ 13.644 , 12.963 , 12.316 , 9.5333 , 8.4333 ,
|
|---|
| 263 | $ 7.5728 , 6.9696 , 6.7518 , 6.6175 , 5.6000 ,
|
|---|
| 264 | $ 6.1145 /
|
|---|
| 265 | DATA (CSEL(22,J),J=1,41) /
|
|---|
| 266 | $ 1100.0 , 115.00 , 105.00 , 100.00 , 56.000 ,
|
|---|
| 267 | $ 40.000 , 27.000 , 22.000 , 21.000 , 20.000 ,
|
|---|
| 268 | $ 20.000 , 19.067 , 19.333 , 19.500 , 19.833 ,
|
|---|
| 269 | $ 20.567 , 21.800 , 22.900 , 23.869 , 23.809 ,
|
|---|
| 270 | $ 22.161 , 21.488 , 19.732 , 19.433 , 19.345 ,
|
|---|
| 271 | $ 19.029 , 18.121 , 16.280 , 15.258 , 14.280 ,
|
|---|
| 272 | $ 13.644 , 12.963 , 12.316 , 9.5333 , 8.4333 ,
|
|---|
| 273 | $ 7.5728 , 6.9696 , 6.7518 , 6.6175 , 5.6000 ,
|
|---|
| 274 | $ 6.1145 /
|
|---|
| 275 | DATA (CSEL(23,J),J=1,41) /
|
|---|
| 276 | $ 185.88 , 133.23 , 119.37 , 102.86 , 93.102 ,
|
|---|
| 277 | $ 82.752 , 76.205 , 71.008 , 67.366 , 64.096 ,
|
|---|
| 278 | $ 60.891 , 58.104 , 55.241 , 53.140 , 50.934 ,
|
|---|
| 279 | $ 48.660 , 47.566 , 46.585 , 45.581 , 44.003 ,
|
|---|
| 280 | $ 41.134 , 39.374 , 36.878 , 35.523 , 34.503 ,
|
|---|
| 281 | $ 32.334 , 29.365 , 27.370 , 24.705 , 22.921 ,
|
|---|
| 282 | $ 21.229 , 19.879 , 18.559 , 14.625 , 12.758 ,
|
|---|
| 283 | $ 11.041 , 9.3440 , 8.5484 , 6.7104 , 6.0000 ,
|
|---|
| 284 | $ 6.1131 /
|
|---|
| 285 | DATA (CSEL(25,J),J=1,41) /
|
|---|
| 286 | $ 157.65 , 73.701 , 76.096 , 68.571 , 57.305 ,
|
|---|
| 287 | $ 49.257 , 43.616 , 40.024 , 38.098 , 36.287 ,
|
|---|
| 288 | $ 34.674 , 33.105 , 31.712 , 30.685 , 29.613 ,
|
|---|
| 289 | $ 28.602 , 28.336 , 28.075 , 27.786 , 27.215 ,
|
|---|
| 290 | $ 26.380 , 26.146 , 25.108 , 24.783 , 24.360 ,
|
|---|
| 291 | $ 23.219 , 21.431 , 20.095 , 18.382 , 17.267 ,
|
|---|
| 292 | $ 16.100 , 15.175 , 14.271 , 11.573 , 10.305 ,
|
|---|
| 293 | $ 9.1471 , 8.0149 , 7.4349 , 6.2499 , 5.8928 ,
|
|---|
| 294 | $ 6.0774 /
|
|---|
| 295 | DATA (CSEL(26,J),J=1,41) /
|
|---|
| 296 | $ 1100.0 , 115.00 , 105.00 , 100.00 , 56.000 ,
|
|---|
| 297 | $ 40.000 , 27.000 , 22.000 , 21.000 , 20.000 ,
|
|---|
| 298 | $ 20.000 , 18.133 , 18.167 , 18.000 , 17.667 ,
|
|---|
| 299 | $ 18.133 , 19.600 , 21.800 , 23.338 , 23.118 ,
|
|---|
| 300 | $ 19.323 , 17.476 , 13.464 , 12.367 , 11.691 ,
|
|---|
| 301 | $ 11.057 , 10.242 , 9.5593 , 9.0151 , 8.5591 ,
|
|---|
| 302 | $ 8.2884 , 7.9253 , 7.6311 , 6.0667 , 5.3667 ,
|
|---|
| 303 | $ 4.8456 , 4.5392 , 4.5036 , 4.4351 , 4.2000 ,
|
|---|
| 304 | $ 4.7289 /
|
|---|
| 305 | DATA (CSEL(27,J),J=1,41) /
|
|---|
| 306 | $ 1100.0 , 115.00 , 105.00 , 100.00 , 56.000 ,
|
|---|
| 307 | $ 40.000 , 27.000 , 22.000 , 21.000 , 20.000 ,
|
|---|
| 308 | $ 20.000 , 18.133 , 18.167 , 18.000 , 17.667 ,
|
|---|
| 309 | $ 18.133 , 19.600 , 21.800 , 23.338 , 23.118 ,
|
|---|
| 310 | $ 19.323 , 17.476 , 13.464 , 12.367 , 11.691 ,
|
|---|
| 311 | $ 11.057 , 10.242 , 9.5593 , 9.0151 , 8.5591 ,
|
|---|
| 312 | $ 8.2884 , 7.9253 , 7.6311 , 6.0667 , 5.3667 ,
|
|---|
| 313 | $ 4.8456 , 4.5392 , 4.5036 , 4.4351 , 4.2000 ,
|
|---|
| 314 | $ 4.7289 /
|
|---|
| 315 | DATA (CSEL(28,J),J=1,41) /
|
|---|
| 316 | $ 157.65 , 73.701 , 76.096 , 68.571 , 57.305 ,
|
|---|
| 317 | $ 49.257 , 43.616 , 40.024 , 38.098 , 36.287 ,
|
|---|
| 318 | $ 34.674 , 32.708 , 31.218 , 30.052 , 28.707 ,
|
|---|
| 319 | $ 27.591 , 27.417 , 27.615 , 27.564 , 26.913 ,
|
|---|
| 320 | $ 24.891 , 23.734 , 20.871 , 19.677 , 18.734 ,
|
|---|
| 321 | $ 17.311 , 15.563 , 14.803 , 13.448 , 12.615 ,
|
|---|
| 322 | $ 11.794 , 11.106 , 10.474 , 8.4745 , 7.4498 ,
|
|---|
| 323 | $ 6.5350 , 5.6835 , 5.3300 , 4.3406 , 4.4464 ,
|
|---|
| 324 | $ 4.7083 /
|
|---|
| 325 | DATA (CSEL(29,J),J=1,41) /
|
|---|
| 326 | $ 143.53 , 43.935 , 54.462 , 51.429 , 39.407 ,
|
|---|
| 327 | $ 32.510 , 27.321 , 24.532 , 23.465 , 22.383 ,
|
|---|
| 328 | $ 21.566 , 20.209 , 19.453 , 18.825 , 18.046 ,
|
|---|
| 329 | $ 17.562 , 17.802 , 18.360 , 18.667 , 18.519 ,
|
|---|
| 330 | $ 17.514 , 17.120 , 14.985 , 14.306 , 13.663 ,
|
|---|
| 331 | $ 12.753 , 11.596 , 11.165 , 10.287 , 9.7882 ,
|
|---|
| 332 | $ 9.2294 , 8.7539 , 8.3300 , 6.9480 , 6.2234 ,
|
|---|
| 333 | $ 5.5881 , 5.0189 , 4.7733 , 4.1104 , 4.3929 ,
|
|---|
| 334 | $ 4.6905 /
|
|---|
| 335 | DATA (CSEL(33,J),J=1,41) /
|
|---|
| 336 | $ 1100.0 , 115.00 , 105.00 , 100.00 , 56.000 ,
|
|---|
| 337 | $ 40.000 , 27.000 , 22.000 , 21.000 , 20.000 ,
|
|---|
| 338 | $ 20.000 , 18.133 , 18.167 , 18.000 , 17.667 ,
|
|---|
| 339 | $ 18.133 , 19.600 , 21.800 , 23.338 , 23.118 ,
|
|---|
| 340 | $ 19.323 , 17.476 , 13.464 , 12.367 , 11.691 ,
|
|---|
| 341 | $ 11.057 , 10.242 , 9.5593 , 9.0151 , 8.5591 ,
|
|---|
| 342 | $ 8.2884 , 7.9253 , 7.6311 , 6.0667 , 5.3667 ,
|
|---|
| 343 | $ 4.8456 , 4.5392 , 4.5036 , 4.4351 , 4.2000 ,
|
|---|
| 344 | $ 4.7289 /
|
|---|
| 345 | DATA (CSEL(34,J),J=1,41) /
|
|---|
| 346 | $ 143.53 , 43.935 , 54.462 , 51.429 , 39.407 ,
|
|---|
| 347 | $ 32.510 , 27.321 , 24.532 , 23.465 , 22.383 ,
|
|---|
| 348 | $ 21.566 , 20.209 , 19.453 , 18.825 , 18.046 ,
|
|---|
| 349 | $ 17.562 , 17.802 , 18.360 , 18.667 , 18.519 ,
|
|---|
| 350 | $ 17.514 , 17.120 , 14.985 , 14.306 , 13.663 ,
|
|---|
| 351 | $ 12.753 , 11.596 , 11.165 , 10.287 , 9.7882 ,
|
|---|
| 352 | $ 9.2294 , 8.7539 , 8.3300 , 6.9480 , 6.2234 ,
|
|---|
| 353 | $ 5.5881 , 5.0189 , 4.7733 , 4.1104 , 4.3929 ,
|
|---|
| 354 | $ 4.6905 /
|
|---|
| 355 |
|
|---|
| 356 | DATA ((CSIN(I,J),I=1,6),J=1,41) / 246 * 0. /
|
|---|
| 357 | DATA (CSIN( 8,J),J=1,41) / 41 * 0. /
|
|---|
| 358 | DATA (CSIN(21,J),J=1,41) / 41 * 0. /
|
|---|
| 359 | DATA (CSIN(24,J),J=1,41) / 41 * 0. /
|
|---|
| 360 | DATA ((CSIN(I,J),I=30,32),J=1,41) / 123 * 0. /
|
|---|
| 361 | DATA (CSIN(35,J),J=1,41) / 41 * 0. /
|
|---|
| 362 |
|
|---|
| 363 | DATA (CSIN( 7,J),J=1,41) /
|
|---|
| 364 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 365 | $ 0.00000E+00, 0.00000E+00, 0.50000 , 1.2000 , 1.7000 ,
|
|---|
| 366 | $ 2.2500 , 3.0000 , 3.6000 , 4.5000 , 5.4000 ,
|
|---|
| 367 | $ 6.3000 , 8.6000 , 9.0000 , 10.000 , 11.500 ,
|
|---|
| 368 | $ 14.000 , 17.000 , 19.500 , 22.000 , 24.000 ,
|
|---|
| 369 | $ 21.500 , 18.500 , 19.000 , 20.500 , 22.200 ,
|
|---|
| 370 | $ 23.000 , 23.300 , 23.000 , 21.000 , 20.500 ,
|
|---|
| 371 | $ 20.200 , 20.100 , 20.000 , 20.000 , 20.000 ,
|
|---|
| 372 | $ 21.000 /
|
|---|
| 373 | DATA (CSIN( 9,J),J=1,41) /
|
|---|
| 374 | $ 0.00000E+00, 3.0000 , 9.2000 , 20.500 , 36.500 ,
|
|---|
| 375 | $ 45.000 , 28.000 , 19.500 , 15.500 , 14.200 ,
|
|---|
| 376 | $ 15.500 , 17.500 , 20.000 , 23.000 , 26.000 ,
|
|---|
| 377 | $ 20.000 , 23.000 , 26.500 , 32.000 , 35.000 ,
|
|---|
| 378 | $ 28.500 , 22.000 , 22.500 , 23.500 , 24.000 ,
|
|---|
| 379 | $ 24.500 , 26.000 , 27.500 , 27.500 , 27.000 ,
|
|---|
| 380 | $ 26.500 , 25.500 , 25.000 , 23.000 , 22.500 ,
|
|---|
| 381 | $ 22.200 , 22.000 , 22.000 , 21.200 , 20.700 ,
|
|---|
| 382 | $ 21.000 /
|
|---|
| 383 | DATA (CSIN(10,J),J=1,41) /
|
|---|
| 384 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 385 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 386 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 387 | $ 0.50000 , 1.5000 , 2.7000 , 3.8000 , 4.8000 ,
|
|---|
| 388 | $ 6.5000 , 7.6000 , 8.4000 , 9.0000 , 9.4000 ,
|
|---|
| 389 | $ 9.8000 , 10.500 , 11.000 , 11.500 , 11.800 ,
|
|---|
| 390 | $ 12.200 , 12.400 , 12.600 , 13.200 , 13.500 ,
|
|---|
| 391 | $ 13.700 , 14.000 , 14.200 , 14.500 , 16.400 ,
|
|---|
| 392 | $ 17.000 /
|
|---|
| 393 | DATA (CSIN(11,J),J=1,41) /
|
|---|
| 394 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 395 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 396 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 397 | $ 0.50000 , 1.5000 , 2.7000 , 3.8000 , 4.8000 ,
|
|---|
| 398 | $ 6.5000 , 7.6000 , 8.4000 , 9.0000 , 9.4000 ,
|
|---|
| 399 | $ 9.8000 , 10.500 , 11.000 , 11.500 , 11.800 ,
|
|---|
| 400 | $ 12.200 , 12.400 , 12.600 , 13.200 , 13.500 ,
|
|---|
| 401 | $ 13.700 , 14.000 , 14.200 , 14.500 , 16.400 ,
|
|---|
| 402 | $ 17.000 /
|
|---|
| 403 | DATA (CSIN(12,J),J=1,41) /
|
|---|
| 404 | $ 266.67 , 133.33 , 83.333 , 57.083 , 44.500 ,
|
|---|
| 405 | $ 33.250 , 24.583 , 20.833 , 18.333 , 16.083 ,
|
|---|
| 406 | $ 15.625 , 15.083 , 14.833 , 15.083 , 15.833 ,
|
|---|
| 407 | $ 17.042 , 18.958 , 20.758 , 22.533 , 22.825 ,
|
|---|
| 408 | $ 21.250 , 18.567 , 17.767 , 18.100 , 19.933 ,
|
|---|
| 409 | $ 20.783 , 21.225 , 21.000 , 20.558 , 20.258 ,
|
|---|
| 410 | $ 20.017 , 19.767 , 19.600 , 19.183 , 18.850 ,
|
|---|
| 411 | $ 18.575 , 18.350 , 18.175 , 17.808 , 17.558 ,
|
|---|
| 412 | $ 19.250 /
|
|---|
| 413 | DATA (CSIN(13,J),J=1,41) /
|
|---|
| 414 | $ 400.00 , 200.00 , 120.00 , 81.000 , 62.000 ,
|
|---|
| 415 | $ 47.000 , 35.000 , 28.000 , 24.000 , 21.000 ,
|
|---|
| 416 | $ 19.500 , 19.000 , 18.800 , 19.000 , 20.000 ,
|
|---|
| 417 | $ 21.000 , 23.000 , 25.000 , 27.000 , 27.500 ,
|
|---|
| 418 | $ 25.500 , 22.000 , 20.800 , 21.000 , 23.000 ,
|
|---|
| 419 | $ 24.000 , 24.000 , 23.800 , 23.000 , 22.500 ,
|
|---|
| 420 | $ 22.000 , 21.600 , 21.400 , 21.000 , 20.500 ,
|
|---|
| 421 | $ 20.200 , 19.800 , 19.500 , 18.600 , 17.500 ,
|
|---|
| 422 | $ 20.000 /
|
|---|
| 423 | DATA (CSIN(14,J),J=1,41) /
|
|---|
| 424 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 425 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 426 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 427 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.10000 , 1.5000 ,
|
|---|
| 428 | $ 7.0000 , 12.000 , 17.000 , 19.500 , 20.500 ,
|
|---|
| 429 | $ 22.000 , 23.500 , 24.800 , 25.800 , 26.500 ,
|
|---|
| 430 | $ 27.000 , 27.500 , 28.000 , 30.000 , 31.000 ,
|
|---|
| 431 | $ 32.000 , 32.500 , 32.500 , 33.000 , 33.500 ,
|
|---|
| 432 | $ 34.000 /
|
|---|
| 433 | DATA (CSIN(15,J),J=1,41) /
|
|---|
| 434 | $ 1500.0 , 1160.0 , 310.00 , 230.00 , 178.00 ,
|
|---|
| 435 | $ 153.00 , 134.00 , 124.00 , 113.00 , 106.00 ,
|
|---|
| 436 | $ 101.00 , 96.000 , 92.000 , 89.000 , 87.000 ,
|
|---|
| 437 | $ 84.000 , 81.000 , 78.500 , 76.500 , 75.000 ,
|
|---|
| 438 | $ 72.000 , 70.000 , 68.000 , 64.500 , 63.000 ,
|
|---|
| 439 | $ 62.000 , 61.000 , 59.500 , 58.500 , 56.500 ,
|
|---|
| 440 | $ 56.500 , 56.000 , 55.500 , 52.000 , 50.000 ,
|
|---|
| 441 | $ 48.000 , 45.000 , 44.000 , 39.200 , 34.500 ,
|
|---|
| 442 | $ 34.500 /
|
|---|
| 443 | DATA (CSIN(16,J),J=1,41) /
|
|---|
| 444 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 445 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 446 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 447 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.10000 , 1.5000 ,
|
|---|
| 448 | $ 7.0000 , 12.000 , 17.000 , 19.500 , 20.500 ,
|
|---|
| 449 | $ 22.000 , 23.500 , 24.800 , 25.800 , 26.500 ,
|
|---|
| 450 | $ 27.000 , 27.500 , 28.000 , 30.000 , 31.000 ,
|
|---|
| 451 | $ 32.000 , 32.500 , 32.500 , 33.000 , 33.500 ,
|
|---|
| 452 | $ 34.000 /
|
|---|
| 453 | DATA (CSIN(17,J),J=1,41) /
|
|---|
| 454 | $ 1394.1 , 948.17 , 262.43 , 197.14 , 149.30 ,
|
|---|
| 455 | $ 127.25 , 110.39 , 101.79 , 92.834 , 87.104 ,
|
|---|
| 456 | $ 83.109 , 79.099 , 75.965 , 73.627 , 72.161 ,
|
|---|
| 457 | $ 69.889 , 67.595 , 65.595 , 64.057 , 63.054 ,
|
|---|
| 458 | $ 61.377 , 60.434 , 59.485 , 56.970 , 55.931 ,
|
|---|
| 459 | $ 55.398 , 54.827 , 53.538 , 52.861 , 51.247 ,
|
|---|
| 460 | $ 51.344 , 50.992 , 50.644 , 47.876 , 46.358 ,
|
|---|
| 461 | $ 44.887 , 42.577 , 41.815 , 38.180 , 34.254 ,
|
|---|
| 462 | $ 34.418 /
|
|---|
| 463 | DATA (CSIN(18,J),J=1,41) /
|
|---|
| 464 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 465 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 466 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 467 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.97815E-01, 1.4577 ,
|
|---|
| 468 | $ 6.2052 , 10.112 , 12.902 , 14.300 , 14.688 ,
|
|---|
| 469 | $ 15.505 , 16.379 , 17.554 , 18.309 , 18.920 ,
|
|---|
| 470 | $ 19.389 , 19.804 , 20.284 , 22.000 , 22.733 ,
|
|---|
| 471 | $ 23.527 , 24.097 , 24.382 , 24.816 , 26.800 ,
|
|---|
| 472 | $ 27.719 /
|
|---|
| 473 | DATA (CSIN(19,J),J=1,41) /
|
|---|
| 474 | $ 1182.4 , 524.50 , 167.30 , 131.43 , 91.895 ,
|
|---|
| 475 | $ 75.743 , 63.184 , 57.376 , 52.502 , 49.313 ,
|
|---|
| 476 | $ 47.326 , 44.762 , 43.222 , 42.015 , 41.221 ,
|
|---|
| 477 | $ 40.244 , 39.504 , 39.145 , 38.860 , 38.731 ,
|
|---|
| 478 | $ 37.987 , 37.814 , 36.326 , 34.750 , 33.953 ,
|
|---|
| 479 | $ 33.635 , 33.349 , 32.938 , 32.785 , 32.092 ,
|
|---|
| 480 | $ 32.373 , 32.312 , 32.329 , 31.261 , 30.597 ,
|
|---|
| 481 | $ 30.073 , 29.228 , 29.182 , 27.683 , 27.107 ,
|
|---|
| 482 | $ 27.956 /
|
|---|
| 483 | DATA (CSIN(20,J),J=1,41) /
|
|---|
| 484 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 485 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 486 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 487 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.97815E-01, 1.4577 ,
|
|---|
| 488 | $ 6.2052 , 10.112 , 12.902 , 14.300 , 14.688 ,
|
|---|
| 489 | $ 15.505 , 16.379 , 17.554 , 18.309 , 18.920 ,
|
|---|
| 490 | $ 19.389 , 19.804 , 20.284 , 22.000 , 22.733 ,
|
|---|
| 491 | $ 23.527 , 24.097 , 24.382 , 24.816 , 26.800 ,
|
|---|
| 492 | $ 27.719 /
|
|---|
| 493 | DATA (CSIN(22,J),J=1,41) /
|
|---|
| 494 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 495 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 496 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 497 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.97815E-01, 1.4577 ,
|
|---|
| 498 | $ 6.2052 , 10.112 , 12.902 , 14.300 , 14.688 ,
|
|---|
| 499 | $ 15.505 , 16.379 , 17.554 , 18.309 , 18.920 ,
|
|---|
| 500 | $ 19.389 , 19.804 , 20.284 , 22.000 , 22.733 ,
|
|---|
| 501 | $ 23.527 , 24.097 , 24.382 , 24.816 , 26.800 ,
|
|---|
| 502 | $ 27.719 /
|
|---|
| 503 | DATA (CSIN(23,J),J=1,41) /
|
|---|
| 504 | $ 1394.1 , 948.17 , 262.43 , 197.14 , 149.30 ,
|
|---|
| 505 | $ 127.25 , 110.39 , 101.79 , 92.834 , 87.104 ,
|
|---|
| 506 | $ 83.109 , 78.563 , 75.292 , 72.760 , 70.900 ,
|
|---|
| 507 | $ 68.467 , 66.314 , 64.955 , 63.746 , 62.623 ,
|
|---|
| 508 | $ 59.233 , 56.946 , 53.355 , 49.810 , 48.090 ,
|
|---|
| 509 | $ 46.839 , 45.695 , 44.863 , 44.062 , 42.599 ,
|
|---|
| 510 | $ 42.684 , 42.328 , 42.041 , 39.508 , 37.880 ,
|
|---|
| 511 | $ 36.299 , 34.075 , 33.553 , 29.723 , 27.600 ,
|
|---|
| 512 | $ 28.120 /
|
|---|
| 513 | DATA (CSIN(25,J),J=1,41) /
|
|---|
| 514 | $ 1182.4 , 524.50 , 167.30 , 131.43 , 91.895 ,
|
|---|
| 515 | $ 75.743 , 63.184 , 57.376 , 52.502 , 49.313 ,
|
|---|
| 516 | $ 47.326 , 44.762 , 43.222 , 42.015 , 41.221 ,
|
|---|
| 517 | $ 40.244 , 39.504 , 39.145 , 38.860 , 38.731 ,
|
|---|
| 518 | $ 37.987 , 37.814 , 36.326 , 34.750 , 33.953 ,
|
|---|
| 519 | $ 33.635 , 33.349 , 32.938 , 32.785 , 32.092 ,
|
|---|
| 520 | $ 32.373 , 32.312 , 32.329 , 31.261 , 30.597 ,
|
|---|
| 521 | $ 30.073 , 29.228 , 29.182 , 27.683 , 27.107 ,
|
|---|
| 522 | $ 27.956 /
|
|---|
| 523 | DATA (CSIN(26,J),J=1,41) /
|
|---|
| 524 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 525 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 526 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 527 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.95639E-01, 1.4154 ,
|
|---|
| 528 | $ 5.4104 , 8.2240 , 8.8031 , 9.1000 , 8.8761 ,
|
|---|
| 529 | $ 9.0095 , 9.2576 , 10.307 , 10.818 , 11.341 ,
|
|---|
| 530 | $ 11.778 , 12.108 , 12.569 , 14.000 , 14.467 ,
|
|---|
| 531 | $ 15.054 , 15.694 , 16.263 , 16.632 , 20.100 ,
|
|---|
| 532 | $ 21.438 /
|
|---|
| 533 | DATA (CSIN(27,J),J=1,41) /
|
|---|
| 534 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 535 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 536 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 537 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.95639E-01, 1.4154 ,
|
|---|
| 538 | $ 5.4104 , 8.2240 , 8.8031 , 9.1000 , 8.8761 ,
|
|---|
| 539 | $ 9.0095 , 9.2576 , 10.307 , 10.818 , 11.341 ,
|
|---|
| 540 | $ 11.778 , 12.108 , 12.569 , 14.000 , 14.467 ,
|
|---|
| 541 | $ 15.054 , 15.694 , 16.263 , 16.632 , 20.100 ,
|
|---|
| 542 | $ 21.438 /
|
|---|
| 543 | DATA (CSIN(28,J),J=1,41) /
|
|---|
| 544 | $ 1182.4 , 524.50 , 167.30 , 131.43 , 91.895 ,
|
|---|
| 545 | $ 75.743 , 63.184 , 57.376 , 52.502 , 49.313 ,
|
|---|
| 546 | $ 47.326 , 44.225 , 42.549 , 41.148 , 39.960 ,
|
|---|
| 547 | $ 38.822 , 38.223 , 38.505 , 38.549 , 38.301 ,
|
|---|
| 548 | $ 35.843 , 34.326 , 30.196 , 27.590 , 26.112 ,
|
|---|
| 549 | $ 25.076 , 24.217 , 24.264 , 23.985 , 23.445 ,
|
|---|
| 550 | $ 23.713 , 23.647 , 23.726 , 22.892 , 22.119 ,
|
|---|
| 551 | $ 21.485 , 20.726 , 20.921 , 19.226 , 20.454 ,
|
|---|
| 552 | $ 21.658 /
|
|---|
| 553 | DATA (CSIN(29,J),J=1,41) /
|
|---|
| 554 | $ 1076.5 , 312.66 , 119.74 , 98.571 , 63.193 ,
|
|---|
| 555 | $ 49.990 , 39.579 , 35.168 , 32.335 , 30.417 ,
|
|---|
| 556 | $ 29.434 , 27.325 , 26.514 , 25.775 , 25.120 ,
|
|---|
| 557 | $ 24.711 , 24.818 , 25.600 , 26.106 , 26.355 ,
|
|---|
| 558 | $ 25.220 , 24.760 , 21.681 , 20.060 , 19.044 ,
|
|---|
| 559 | $ 18.474 , 18.044 , 18.301 , 18.347 , 18.192 ,
|
|---|
| 560 | $ 18.557 , 18.639 , 18.870 , 18.769 , 18.478 ,
|
|---|
| 561 | $ 18.372 , 18.302 , 18.735 , 18.206 , 20.207 ,
|
|---|
| 562 | $ 21.576 /
|
|---|
| 563 | DATA (CSIN(33,J),J=1,41) /
|
|---|
| 564 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 565 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 566 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
|
|---|
| 567 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.95639E-01, 1.4154 ,
|
|---|
| 568 | $ 5.4104 , 8.2240 , 8.8031 , 9.1000 , 8.8761 ,
|
|---|
| 569 | $ 9.0095 , 9.2576 , 10.307 , 10.818 , 11.341 ,
|
|---|
| 570 | $ 11.778 , 12.108 , 12.569 , 14.000 , 14.467 ,
|
|---|
| 571 | $ 15.054 , 15.694 , 16.263 , 16.632 , 20.100 ,
|
|---|
| 572 | $ 21.438 /
|
|---|
| 573 | DATA (CSIN(34,J),J=1,41) /
|
|---|
| 574 | $ 1076.5 , 312.66 , 119.74 , 98.571 , 63.193 ,
|
|---|
| 575 | $ 49.990 , 39.579 , 35.168 , 32.335 , 30.417 ,
|
|---|
| 576 | $ 29.434 , 27.325 , 26.514 , 25.775 , 25.120 ,
|
|---|
| 577 | $ 24.711 , 24.818 , 25.600 , 26.106 , 26.355 ,
|
|---|
| 578 | $ 25.220 , 24.760 , 21.681 , 20.060 , 19.044 ,
|
|---|
| 579 | $ 18.474 , 18.044 , 18.301 , 18.347 , 18.192 ,
|
|---|
| 580 | $ 18.557 , 18.639 , 18.870 , 18.769 , 18.478 ,
|
|---|
| 581 | $ 18.372 , 18.302 , 18.735 , 18.206 , 20.207 ,
|
|---|
| 582 | $ 21.576 /
|
|---|
| 583 |
|
|---|
| 584 | DATA (CSPIEL( 1,J),J=1,41) /
|
|---|
| 585 | $ 0.00000E+00, 350.00 , 580.00 , 600.00 , 550.00 ,
|
|---|
| 586 | $ 450.00 , 410.00 , 370.00 , 340.00 , 230.00 ,
|
|---|
| 587 | $ 220.00 , 205.00 , 180.00 , 155.00 , 145.00 ,
|
|---|
| 588 | $ 140.00 , 160.00 , 195.00 , 235.00 , 250.00 ,
|
|---|
| 589 | $ 270.00 , 280.00 , 300.00 , 300.00 , 290.00 ,
|
|---|
| 590 | $ 285.00 , 265.00 , 240.00 , 230.00 , 222.00 ,
|
|---|
| 591 | $ 204.00 , 196.00 , 190.00 , 170.00 , 170.00 ,
|
|---|
| 592 | $ 160.00 , 150.00 , 140.00 , 120.00 , 80.000 ,
|
|---|
| 593 | $ 80.000 /
|
|---|
| 594 | DATA (CSPIEL( 2,J),J=1,41) /
|
|---|
| 595 | $ 0.00000E+00, 700.00 , 1000.0 , 1200.0 , 1300.0 ,
|
|---|
| 596 | $ 1300.0 , 1250.0 , 1250.0 , 1100.0 , 1000.0 ,
|
|---|
| 597 | $ 940.00 , 740.00 , 700.00 , 670.00 , 660.00 ,
|
|---|
| 598 | $ 670.00 , 680.00 , 700.00 , 735.00 , 800.00 ,
|
|---|
| 599 | $ 810.00 , 820.00 , 820.00 , 810.00 , 800.00 ,
|
|---|
| 600 | $ 800.00 , 700.00 , 600.00 , 500.00 , 470.00 ,
|
|---|
| 601 | $ 440.00 , 410.00 , 380.00 , 330.00 , 330.00 ,
|
|---|
| 602 | $ 330.00 , 330.00 , 330.00 , 285.00 , 240.00 ,
|
|---|
| 603 | $ 240.00 /
|
|---|
| 604 | DATA (CSPIEL( 3,J),J=1,41) /
|
|---|
| 605 | $ 0.00000E+00, 1700.0 , 2200.0 , 2200.0 , 1800.0 ,
|
|---|
| 606 | $ 1300.0 , 1200.0 , 900.00 , 900.00 , 1000.0 ,
|
|---|
| 607 | $ 1100.0 , 1300.0 , 1400.0 , 1420.0 , 1490.0 ,
|
|---|
| 608 | $ 1560.0 , 1580.0 , 1690.0 , 1795.0 , 2000.0 ,
|
|---|
| 609 | $ 2070.0 , 2140.0 , 2050.0 , 2010.0 , 1970.0 ,
|
|---|
| 610 | $ 1880.0 , 1690.0 , 1500.0 , 1420.0 , 1390.0 ,
|
|---|
| 611 | $ 1350.0 , 1360.0 , 1370.0 , 1280.0 , 1290.0 ,
|
|---|
| 612 | $ 1295.0 , 1250.0 , 1200.0 , 1050.0 , 900.00 ,
|
|---|
| 613 | $ 900.00 /
|
|---|
| 614 | DATA (CSPIIN( 1,J),J=1,41) /
|
|---|
| 615 | $ 0.00000E+00, 200.00 , 320.00 , 500.00 , 600.00 ,
|
|---|
| 616 | $ 600.00 , 590.00 , 530.00 , 510.00 , 470.00 ,
|
|---|
| 617 | $ 430.00 , 425.00 , 420.00 , 425.00 , 425.00 ,
|
|---|
| 618 | $ 430.00 , 430.00 , 435.00 , 435.00 , 440.00 ,
|
|---|
| 619 | $ 430.00 , 430.00 , 420.00 , 420.00 , 420.00 ,
|
|---|
| 620 | $ 415.00 , 415.00 , 410.00 , 410.00 , 408.00 ,
|
|---|
| 621 | $ 406.00 , 404.00 , 400.00 , 380.00 , 340.00 ,
|
|---|
| 622 | $ 340.00 , 340.00 , 340.00 , 340.00 , 340.00 ,
|
|---|
| 623 | $ 340.00 /
|
|---|
| 624 | DATA (CSPIIN( 2,J),J=1,41) /
|
|---|
| 625 | $ 0.00000E+00, 400.00 , 800.00 , 1000.0 , 1100.0 ,
|
|---|
| 626 | $ 1200.0 , 1150.0 , 1050.0 , 1000.0 , 900.00 ,
|
|---|
| 627 | $ 860.00 , 860.00 , 850.00 , 850.00 , 840.00 ,
|
|---|
| 628 | $ 830.00 , 820.00 , 810.00 , 805.00 , 800.00 ,
|
|---|
| 629 | $ 800.00 , 800.00 , 800.00 , 800.00 , 800.00 ,
|
|---|
| 630 | $ 800.00 , 800.00 , 800.00 , 800.00 , 780.00 ,
|
|---|
| 631 | $ 760.00 , 740.00 , 720.00 , 720.00 , 700.00 ,
|
|---|
| 632 | $ 690.00 , 680.00 , 670.00 , 665.00 , 660.00 ,
|
|---|
| 633 | $ 660.00 /
|
|---|
| 634 | DATA (CSPIIN( 3,J),J=1,41) /
|
|---|
| 635 | $ 0.00000E+00, 1000.0 , 1900.0 , 2600.0 , 2900.0 ,
|
|---|
| 636 | $ 3000.0 , 2800.0 , 2600.0 , 2500.0 , 2300.0 ,
|
|---|
| 637 | $ 2200.0 , 2000.0 , 1900.0 , 1880.0 , 1860.0 ,
|
|---|
| 638 | $ 1840.0 , 1820.0 , 1810.0 , 1805.0 , 1800.0 ,
|
|---|
| 639 | $ 1780.0 , 1760.0 , 1750.0 , 1740.0 , 1730.0 ,
|
|---|
| 640 | $ 1720.0 , 1710.0 , 1700.0 , 1680.0 , 1660.0 ,
|
|---|
| 641 | $ 1650.0 , 1640.0 , 1630.0 , 1620.0 , 1610.0 ,
|
|---|
| 642 | $ 1605.0 , 1600.0 , 1600.0 , 1550.0 , 1500.0 ,
|
|---|
| 643 | $ 1500.0 /
|
|---|
| 644 | DATA (CSPNEL( 1,J),J=1,41) /
|
|---|
| 645 | $ 2100.0 , 1800.0 , 1500.0 , 1050.0 , 900.00 ,
|
|---|
| 646 | $ 950.00 , 800.00 , 650.00 , 570.00 , 390.00 ,
|
|---|
| 647 | $ 300.00 , 240.00 , 230.00 , 230.00 , 220.00 ,
|
|---|
| 648 | $ 220.00 , 225.00 , 225.00 , 240.00 , 240.00 ,
|
|---|
| 649 | $ 290.00 , 330.00 , 335.00 , 350.00 , 355.00 ,
|
|---|
| 650 | $ 370.00 , 350.00 , 330.00 , 310.00 , 290.00 ,
|
|---|
| 651 | $ 270.00 , 265.00 , 260.00 , 230.00 , 210.00 ,
|
|---|
| 652 | $ 210.00 , 200.00 , 200.00 , 190.00 , 180.00 ,
|
|---|
| 653 | $ 180.00 /
|
|---|
| 654 | DATA (CSPNEL( 2,J),J=1,41) /
|
|---|
| 655 | $ 3800.0 , 2900.0 , 1850.0 , 1550.0 , 1450.0 ,
|
|---|
| 656 | $ 1520.0 , 1460.0 , 1300.0 , 1140.0 , 880.00 ,
|
|---|
| 657 | $ 700.00 , 620.00 , 540.00 , 560.00 , 460.00 ,
|
|---|
| 658 | $ 460.00 , 470.00 , 470.00 , 480.00 , 480.00 ,
|
|---|
| 659 | $ 580.00 , 600.00 , 610.00 , 620.00 , 620.00 ,
|
|---|
| 660 | $ 620.00 , 590.00 , 580.00 , 460.00 , 440.00 ,
|
|---|
| 661 | $ 420.00 , 400.00 , 480.00 , 430.00 , 380.00 ,
|
|---|
| 662 | $ 380.00 , 380.00 , 380.00 , 380.00 , 380.00 ,
|
|---|
| 663 | $ 380.00 /
|
|---|
| 664 | DATA (CSPNEL( 3,J),J=1,41) /
|
|---|
| 665 | $ 7000.0 , 6000.0 , 4500.0 , 3350.0 , 2700.0 ,
|
|---|
| 666 | $ 3000.0 , 3550.0 , 3970.0 , 3280.0 , 2490.0 ,
|
|---|
| 667 | $ 2100.0 , 1510.0 , 1440.0 , 1370.0 , 1370.0 ,
|
|---|
| 668 | $ 1370.0 , 1400.0 , 1400.0 , 1420.0 , 1420.0 ,
|
|---|
| 669 | $ 1440.0 , 1460.0 , 1460.0 , 1450.0 , 1450.0 ,
|
|---|
| 670 | $ 1470.0 , 1400.0 , 1400.0 , 1380.0 , 1370.0 ,
|
|---|
| 671 | $ 1360.0 , 1350.0 , 1340.0 , 1330.0 , 1320.0 ,
|
|---|
| 672 | $ 1310.0 , 1305.0 , 1300.0 , 1300.0 , 1300.0 ,
|
|---|
| 673 | $ 1300.0 /
|
|---|
| 674 | DATA (CSPNIN( 1,J),J=1,41) /
|
|---|
| 675 | $ 0.00000E+00, 200.00 , 400.00 , 800.00 , 800.00 ,
|
|---|
| 676 | $ 550.00 , 500.00 , 450.00 , 430.00 , 410.00 ,
|
|---|
| 677 | $ 400.00 , 390.00 , 380.00 , 370.00 , 370.00 ,
|
|---|
| 678 | $ 370.00 , 365.00 , 365.00 , 360.00 , 360.00 ,
|
|---|
| 679 | $ 360.00 , 360.00 , 365.00 , 370.00 , 375.00 ,
|
|---|
| 680 | $ 380.00 , 400.00 , 410.00 , 420.00 , 430.00 ,
|
|---|
| 681 | $ 440.00 , 440.00 , 440.00 , 440.00 , 440.00 ,
|
|---|
| 682 | $ 440.00 , 440.00 , 440.00 , 440.00 , 440.00 ,
|
|---|
| 683 | $ 440.00 /
|
|---|
| 684 | DATA (CSPNIN( 2,J),J=1,41) /
|
|---|
| 685 | $ 0.00000E+00, 400.00 , 950.00 , 1050.0 , 1050.0 ,
|
|---|
| 686 | $ 980.00 , 940.00 , 900.00 , 860.00 , 820.00 ,
|
|---|
| 687 | $ 800.00 , 780.00 , 760.00 , 740.00 , 740.00 ,
|
|---|
| 688 | $ 740.00 , 730.00 , 730.00 , 720.00 , 720.00 ,
|
|---|
| 689 | $ 720.00 , 720.00 , 730.00 , 740.00 , 750.00 ,
|
|---|
| 690 | $ 760.00 , 800.00 , 820.00 , 820.00 , 820.00 ,
|
|---|
| 691 | $ 820.00 , 820.00 , 820.00 , 820.00 , 820.00 ,
|
|---|
| 692 | $ 820.00 , 820.00 , 820.00 , 820.00 , 820.00 ,
|
|---|
| 693 | $ 820.00 /
|
|---|
| 694 | DATA (CSPNIN( 3,J),J=1,41) /
|
|---|
| 695 | $ 0.00000E+00, 0.00000E+00, 500.00 , 1450.0 , 1700.0 ,
|
|---|
| 696 | $ 1800.0 , 1750.0 , 1730.0 , 1720.0 , 1710.0 ,
|
|---|
| 697 | $ 1700.0 , 1690.0 , 1660.0 , 1630.0 , 1630.0 ,
|
|---|
| 698 | $ 1630.0 , 1600.0 , 1600.0 , 1580.0 , 1580.0 ,
|
|---|
| 699 | $ 1580.0 , 1580.0 , 1600.0 , 1630.0 , 1650.0 ,
|
|---|
| 700 | $ 1670.0 , 1760.0 , 1800.0 , 1800.0 , 1800.0 ,
|
|---|
| 701 | $ 1800.0 , 1800.0 , 1800.0 , 1800.0 , 1800.0 ,
|
|---|
| 702 | $ 1800.0 , 1800.0 , 1800.0 , 1800.0 , 1800.0 ,
|
|---|
| 703 | $ 1800.0 /
|
|---|
| 704 | DATA ELAB /
|
|---|
| 705 | $ 0.10000E-03, 0.20000E-03, 0.30000E-03, 0.40000E-03, 0.50000E-03,
|
|---|
| 706 | $ 0.70000E-03, 0.10000E-02, 0.20000E-02, 0.30000E-02, 0.40000E-02,
|
|---|
| 707 | $ 0.50000E-02, 0.70000E-02, 0.10000E-01, 0.15000E-01, 0.20000E-01,
|
|---|
| 708 | $ 0.25000E-01, 0.32700E-01/
|
|---|
| 709 | DATA CNLWAT /
|
|---|
| 710 | $ 1.0000 , 16.000 , 27.000 , 56.000 , 59.000 ,
|
|---|
| 711 | $ 64.000 , 91.000 , 112.00 , 119.00 , 127.00 ,
|
|---|
| 712 | $ 137.00 , 181.00 , 207.00 , 209.00 , 238.00 /
|
|---|
| 713 | DATA (CNLWEL( 1,J),J=1,17) /
|
|---|
| 714 | $ 6000.0 , 5500.0 , 5200.0 , 4900.0 , 4800.0 ,
|
|---|
| 715 | $ 4400.0 , 4000.0 , 2900.0 , 2200.0 , 1800.0 ,
|
|---|
| 716 | $ 1400.0 , 1100.0 , 900.00 , 700.00 , 600.00 ,
|
|---|
| 717 | $ 560.00 , 520.00 /
|
|---|
| 718 | DATA (CNLWEL( 2,J),J=1,17) /
|
|---|
| 719 | $ 5400.0 , 5050.0 , 4800.0 , 4600.0 , 4399.0 ,
|
|---|
| 720 | $ 4090.0 , 3700.0 , 2600.0 , 1950.0 , 1600.0 ,
|
|---|
| 721 | $ 1300.0 , 900.00 , 700.00 , 800.00 , 1050.0 ,
|
|---|
| 722 | $ 1250.0 , 1320.0 /
|
|---|
| 723 | DATA (CNLWEL( 3,J),J=1,17) /
|
|---|
| 724 | $ 5500.0 , 5150.0 , 4900.0 , 4699.0 , 4490.0 ,
|
|---|
| 725 | $ 4150.0 , 3750.0 , 2790.0 , 2100.0 , 1650.0 ,
|
|---|
| 726 | $ 1300.0 , 950.00 , 800.00 , 860.00 , 1000.0 ,
|
|---|
| 727 | $ 1090.0 , 1080.0 /
|
|---|
| 728 | DATA (CNLWEL( 4,J),J=1,17) /
|
|---|
| 729 | $ 5499.0 , 4970.0 , 4450.0 , 4080.0 , 3750.0 ,
|
|---|
| 730 | $ 3380.0 , 2900.0 , 2400.0 , 2380.0 , 2350.0 ,
|
|---|
| 731 | $ 2300.0 , 2100.0 , 1720.0 , 1370.0 , 1200.0 ,
|
|---|
| 732 | $ 1060.0 , 870.00 /
|
|---|
| 733 | DATA (CNLWEL( 5,J),J=1,17) /
|
|---|
| 734 | $ 5399.0 , 4710.0 , 4180.0 , 3760.0 , 3460.0 ,
|
|---|
| 735 | $ 3150.0 , 2730.0 , 2270.0 , 1850.0 , 1850.0 ,
|
|---|
| 736 | $ 2130.0 , 2330.0 , 2120.0 , 1640.0 , 1310.0 ,
|
|---|
| 737 | $ 1100.0 , 1050.0 /
|
|---|
| 738 | DATA (CNLWEL( 6,J),J=1,17) /
|
|---|
| 739 | $ 5099.0 , 4405.0 , 3825.0 , 3455.0 , 3125.0 ,
|
|---|
| 740 | $ 2695.0 , 2350.0 , 1850.0 , 1580.0 , 1820.0 ,
|
|---|
| 741 | $ 2050.0 , 2210.0 , 2000.0 , 1590.0 , 1310.0 ,
|
|---|
| 742 | $ 1120.0 , 1040.0 /
|
|---|
| 743 | DATA (CNLWEL( 7,J),J=1,17) /
|
|---|
| 744 | $ 6290.0 , 5960.0 , 5640.0 , 5370.0 , 5150.0 ,
|
|---|
| 745 | $ 4800.0 , 4250.0 , 3150.0 , 2470.0 , 2100.0 ,
|
|---|
| 746 | $ 2230.0 , 2420.0 , 2450.0 , 2050.0 , 1760.0 ,
|
|---|
| 747 | $ 1550.0 , 1330.0 /
|
|---|
| 748 | DATA (CNLWEL( 8,J),J=1,17) /
|
|---|
| 749 | $ 6885.0 , 6650.0 , 6350.0 , 6150.0 , 6000.0 ,
|
|---|
| 750 | $ 5700.0 , 5360.0 , 4250.0 , 2800.0 , 1870.0 ,
|
|---|
| 751 | $ 1810.0 , 1820.0 , 2170.0 , 2450.0 , 2150.0 ,
|
|---|
| 752 | $ 1700.0 , 1390.0 /
|
|---|
| 753 | DATA (CNLWEL( 9,J),J=1,17) /
|
|---|
| 754 | $ 6600.0 , 6500.0 , 6400.0 , 6249.0 , 6190.0 ,
|
|---|
| 755 | $ 5950.0 , 5520.0 , 4250.0 , 2750.0 , 1900.0 ,
|
|---|
| 756 | $ 1850.0 , 1950.0 , 2340.0 , 2800.0 , 2540.0 ,
|
|---|
| 757 | $ 2100.0 , 1760.0 /
|
|---|
| 758 | DATA (CNLWEL(10,J),J=1,17) /
|
|---|
| 759 | $ 7400.0 , 7200.0 , 6999.0 , 6840.0 , 6655.0 ,
|
|---|
| 760 | $ 6320.0 , 5820.0 , 4400.0 , 2850.0 , 2000.0 ,
|
|---|
| 761 | $ 1800.0 , 1800.0 , 2150.0 , 2600.0 , 2350.0 ,
|
|---|
| 762 | $ 1950.0 , 2100.0 /
|
|---|
| 763 | DATA (CNLWEL(11,J),J=1,17) /
|
|---|
| 764 | $ 7900.0 , 7700.0 , 7499.0 , 7390.0 , 7202.0 ,
|
|---|
| 765 | $ 6810.0 , 6360.0 , 4920.0 , 3450.0 , 2600.0 ,
|
|---|
| 766 | $ 2200.0 , 1950.0 , 2300.0 , 2800.0 , 2650.0 ,
|
|---|
| 767 | $ 2250.0 , 2050.0 /
|
|---|
| 768 | DATA (CNLWEL(12,J),J=1,17) /
|
|---|
| 769 | $ 7900.0 , 7750.0 , 7699.0 , 7590.0 , 7450.0 ,
|
|---|
| 770 | $ 7200.0 , 6850.0 , 5650.0 , 4400.0 , 3700.0 ,
|
|---|
| 771 | $ 3400.0 , 2800.0 , 2700.0 , 3100.0 , 3250.0 ,
|
|---|
| 772 | $ 3100.0 , 2750.0 /
|
|---|
| 773 | DATA (CNLWEL(13,J),J=1,17) /
|
|---|
| 774 | $ 6100.0 , 5950.0 , 5750.0 , 5599.0 , 5440.0 ,
|
|---|
| 775 | $ 5200.0 , 4800.0 , 4300.0 , 5800.0 , 5750.0 ,
|
|---|
| 776 | $ 4800.0 , 3420.0 , 2650.0 , 3200.0 , 3650.0 ,
|
|---|
| 777 | $ 3500.0 , 2980.0 /
|
|---|
| 778 | DATA (CNLWEL(14,J),J=1,17) /
|
|---|
| 779 | $ 6100.0 , 5950.0 , 5750.0 , 5599.0 , 5440.0 ,
|
|---|
| 780 | $ 5200.0 , 4800.0 , 4300.0 , 5800.0 , 5750.0 ,
|
|---|
| 781 | $ 4800.0 , 3420.0 , 2650.0 , 3200.0 , 3650.0 ,
|
|---|
| 782 | $ 3500.0 , 2980.0 /
|
|---|
| 783 | DATA (CNLWEL(15,J),J=1,17) /
|
|---|
| 784 | $ 6600.0 , 6350.0 , 6100.0 , 5899.0 , 5690.0 ,
|
|---|
| 785 | $ 5300.0 , 4850.0 , 4450.0 , 5650.0 , 5700.0 ,
|
|---|
| 786 | $ 4950.0 , 3850.0 , 3050.0 , 3050.0 , 3460.0 ,
|
|---|
| 787 | $ 3650.0 , 3340.0 /
|
|---|
| 788 | DATA (CNLWIN( 1,J),J=1,17) / 17*0.0E+00 /
|
|---|
| 789 | DATA (CNLWIN( 2,J),J=1,17) /
|
|---|
| 790 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 ,
|
|---|
| 791 | $ 10.000 , 50.000 , 100.00 , 200.00 , 300.00 ,
|
|---|
| 792 | $ 400.00 , 600.00 , 700.00 , 750.00 , 700.00 ,
|
|---|
| 793 | $ 700.00 , 680.00 /
|
|---|
| 794 | DATA (CNLWIN( 3,J),J=1,17) /
|
|---|
| 795 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 ,
|
|---|
| 796 | $ 50.000 , 100.00 , 260.00 , 450.00 , 600.00 ,
|
|---|
| 797 | $ 700.00 , 800.00 , 900.00 , 940.00 , 900.00 ,
|
|---|
| 798 | $ 860.00 , 820.00 /
|
|---|
| 799 | DATA (CNLWIN( 4,J),J=1,17) /
|
|---|
| 800 | $ 1.0000 , 80.000 , 200.00 , 320.00 , 400.00 ,
|
|---|
| 801 | $ 520.00 , 700.00 , 1000.0 , 1120.0 , 1200.0 ,
|
|---|
| 802 | $ 1200.0 , 1200.0 , 1180.0 , 1130.0 , 1100.0 ,
|
|---|
| 803 | $ 1090.0 , 1080.0 /
|
|---|
| 804 | DATA (CNLWIN( 5,J),J=1,17) /
|
|---|
| 805 | $ 1.0000 , 90.000 , 220.00 , 340.00 , 420.00 ,
|
|---|
| 806 | $ 550.00 , 720.00 , 1080.0 , 1300.0 , 1400.0 ,
|
|---|
| 807 | $ 1420.0 , 1420.0 , 1380.0 , 1260.0 , 1190.0 ,
|
|---|
| 808 | $ 1150.0 , 1100.0 /
|
|---|
| 809 | DATA (CNLWIN( 6,J),J=1,17) /
|
|---|
| 810 | $ 1.0000 , 95.000 , 225.00 , 345.00 , 425.00 ,
|
|---|
| 811 | $ 555.00 , 750.00 , 1150.0 , 1500.0 , 1680.0 ,
|
|---|
| 812 | $ 1700.0 , 1690.0 , 1550.0 , 1360.0 , 1240.0 ,
|
|---|
| 813 | $ 1180.0 , 1120.0 /
|
|---|
| 814 | DATA (CNLWIN( 7,J),J=1,17) /
|
|---|
| 815 | $ 10.000 , 140.00 , 260.00 , 380.00 , 450.00 ,
|
|---|
| 816 | $ 600.00 , 750.00 , 1200.0 , 1580.0 , 1800.0 ,
|
|---|
| 817 | $ 1820.0 , 1830.0 , 1800.0 , 1750.0 , 1690.0 ,
|
|---|
| 818 | $ 1650.0 , 1620.0 /
|
|---|
| 819 | DATA (CNLWIN( 8,J),J=1,17) /
|
|---|
| 820 | $ 15.000 , 150.00 , 300.00 , 400.00 , 500.00 ,
|
|---|
| 821 | $ 650.00 , 840.00 , 1500.0 , 2100.0 , 2130.0 ,
|
|---|
| 822 | $ 2140.0 , 2130.0 , 2080.0 , 2000.0 , 1950.0 ,
|
|---|
| 823 | $ 1900.0 , 1860.0 /
|
|---|
| 824 | DATA (CNLWIN( 9,J),J=1,17) /
|
|---|
| 825 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 ,
|
|---|
| 826 | $ 150.00 , 380.00 , 1000.0 , 1650.0 , 2100.0 ,
|
|---|
| 827 | $ 2100.0 , 2100.0 , 2060.0 , 1950.0 , 1860.0 ,
|
|---|
| 828 | $ 1800.0 , 1740.0 /
|
|---|
| 829 | DATA (CNLWIN(10,J),J=1,17) /
|
|---|
| 830 | $ 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 , 45.000 ,
|
|---|
| 831 | $ 180.00 , 380.00 , 1050.0 , 1900.0 , 2300.0 ,
|
|---|
| 832 | $ 2300.0 , 2200.0 , 2150.0 , 2000.0 , 1900.0 ,
|
|---|
| 833 | $ 1800.0 , 1750.0 /
|
|---|
| 834 | DATA (CNLWIN(11,J),J=1,17) /
|
|---|
| 835 | $ 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 , 48.000 ,
|
|---|
| 836 | $ 190.00 , 390.00 , 1080.0 , 2000.0 , 2400.0 ,
|
|---|
| 837 | $ 2400.0 , 2300.0 , 2200.0 , 2100.0 , 1950.0 ,
|
|---|
| 838 | $ 1850.0 , 1800.0 /
|
|---|
| 839 | DATA (CNLWIN(12,J),J=1,17) /
|
|---|
| 840 | $ 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 , 50.000 ,
|
|---|
| 841 | $ 200.00 , 400.00 , 1100.0 , 2100.0 , 2500.0 ,
|
|---|
| 842 | $ 2500.0 , 2450.0 , 2300.0 , 2100.0 , 2000.0 ,
|
|---|
| 843 | $ 1900.0 , 1850.0 /
|
|---|
| 844 | DATA (CNLWIN(13,J),J=1,17) /
|
|---|
| 845 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 ,
|
|---|
| 846 | $ 100.00 , 350.00 , 900.00 , 1400.0 , 2000.0 ,
|
|---|
| 847 | $ 2300.0 , 2380.0 , 2400.0 , 2300.0 , 2250.0 ,
|
|---|
| 848 | $ 2200.0 , 2120.0 /
|
|---|
| 849 | DATA (CNLWIN(14,J),J=1,17) /
|
|---|
| 850 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 ,
|
|---|
| 851 | $ 100.00 , 350.00 , 900.00 , 1400.0 , 2000.0 ,
|
|---|
| 852 | $ 2300.0 , 2380.0 , 2400.0 , 2300.0 , 2250.0 ,
|
|---|
| 853 | $ 2200.0 , 2120.0 /
|
|---|
| 854 | DATA (CNLWIN(15,J),J=1,17) /
|
|---|
| 855 | $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 ,
|
|---|
| 856 | $ 100.00 , 400.00 , 950.00 , 1600.0 , 2200.0 ,
|
|---|
| 857 | $ 2550.0 , 2750.0 , 2700.0 , 2600.0 , 2540.0 ,
|
|---|
| 858 | $ 2450.0 , 2360.0 /
|
|---|
| 859 | DATA (CSCAP(J),J=1,50) /
|
|---|
| 860 | $ 6.0000 , 5.7000 , 5.5000 , 5.3000 , 5.2000 ,
|
|---|
| 861 | $ 5.1000 , 5.0000 , 4.9000 , 4.8000 , 4.8000 ,
|
|---|
| 862 | $ 4.8000 , 4.8000 , 4.8000 , 4.8000 , 4.8000 ,
|
|---|
| 863 | $ 4.8000 , 4.9000 , 5.0000 , 5.2000 , 5.5000 ,
|
|---|
| 864 | $ 6.0000 , 6.7000 , 7.5000 , 8.5000 , 10.000 ,
|
|---|
| 865 | $ 12.000 , 14.500 , 19.000 , 26.500 , 40.000 ,
|
|---|
| 866 | $ 75.000 , 120.00 , 180.00 , 260.00 , 360.00 ,
|
|---|
| 867 | $ 330.00 , 60.000 , 7.0000 , 9.5000 , 20.000 ,
|
|---|
| 868 | $ 75.000 , 140.00 , 250.00 , 360.00 , 480.00 ,
|
|---|
| 869 | $ 580.00 , 590.00 , 500.00 , 300.00 , 100.00 /
|
|---|
| 870 | DATA (CSCAP(J),J=51,100) /
|
|---|
| 871 | $ 200.00 , 300.00 , 400.00 , 470.00 , 500.00 ,
|
|---|
| 872 | $ 430.00 , 100.00 , 20.000 , 22.000 , 40.000 ,
|
|---|
| 873 | $ 560.00 , 950.00 , 1000.0 , 1000.0 , 1000.0 ,
|
|---|
| 874 | $ 990.00 , 920.00 , 860.00 , 790.00 , 740.00 ,
|
|---|
| 875 | $ 650.00 , 600.00 , 540.00 , 470.00 , 440.00 ,
|
|---|
| 876 | $ 390.00 , 360.00 , 340.00 , 320.00 , 310.00 ,
|
|---|
| 877 | $ 280.00 , 2.0000 , 2.5000 , 6.0000 , 13.000 ,
|
|---|
| 878 | $ 38.000 , 65.000 , 140.00 , 280.00 , 300.00 ,
|
|---|
| 879 | $ 430.00 , 580.00 , 650.00 , 800.00 , 920.00 ,
|
|---|
| 880 | $ 1100.0 , 1250.0 , 1400.0 , 1550.0 , 1700.0 /
|
|---|
| 881 | C --- END OF CROSS SECTION DATA STATEMENTS ---
|
|---|
| 882 |
|
|---|
| 883 | C --- DATA STMTS. FOR GEANT/GHEISHA PARTICLE CODE CONVERSIONS ---
|
|---|
| 884 | C --- KIPART(I)=GHEISHA CODE CORRESPONDING TO GEANT CODE I ---
|
|---|
| 885 | C --- IKPART(I)=GEANT CODE CORRESPONDING TO GHEISHA CODE I ---
|
|---|
| 886 |
|
|---|
| 887 | DATA KIPART/
|
|---|
| 888 | $ 1, 3, 4, 2, 5, 6, 8, 7,
|
|---|
| 889 | $ 9, 12, 10, 13, 16, 14, 15, 11,
|
|---|
| 890 | $ 35, 18, 20, 21, 22, 26, 27, 33,
|
|---|
| 891 | $ 17, 19, 23, 24, 25, 28, 29, 34,
|
|---|
| 892 | $ 35, 35, 35, 35, 35, 35, 35, 35,
|
|---|
| 893 | $ 35, 35, 35, 35, 30, 31, 32, 35/
|
|---|
| 894 |
|
|---|
| 895 | DATA IKPART/
|
|---|
| 896 | $ 1, 4, 2, 3, 5, 6, 8, 7,
|
|---|
| 897 | $ 9, 11, 16, 10, 12, 14, 15, 13,
|
|---|
| 898 | $ 25, 18, 26, 19, 20, 21, 27, 28,
|
|---|
| 899 | $ 29, 22, 23, 30, 31, 45, 46, 47,
|
|---|
| 900 | $ 24, 32, 48/
|
|---|
| 901 |
|
|---|
| 902 |
|
|---|
| 903 | C PARAMETER (ONETHR=1./3.)
|
|---|
| 904 | DATA ONETHR / .33333333/
|
|---|
| 905 | DATA ALPHA / 6*0.7,
|
|---|
| 906 | + 0.75 ,0.75 ,0.75 ,
|
|---|
| 907 | + 0.76,0.76 ,0.76 ,0.76 ,
|
|---|
| 908 | + 0.685,0.63 ,0.685,0.63,0.685,0.63,
|
|---|
| 909 | + 3*0.685,3*0.63,2*0.685,2*0.63,
|
|---|
| 910 | + 3*0.7,0.685,0.63,0.7/
|
|---|
| 911 | DATA ALPHAC /1.2,1.2,1.2,1.15,0.90,0.91,0.98,1.06,1.10,1.11,
|
|---|
| 912 | + 1.10,1.08,1.05,1.01,0.985,0.962,0.945,0.932,
|
|---|
| 913 | + 0.925,0.920,0.920,0.921,0.922,0.923,0.928,0.931,
|
|---|
| 914 | + 0.940,0.945,0.950,0.955,0.958,0.962,0.965,0.976,
|
|---|
| 915 | + 0.982,0.988,0.992,1.010,1.020,1.030,1.040/
|
|---|
| 916 | DATA PARTEL/6*0.,29*1./
|
|---|
| 917 | DATA PARTIN/6*0.,1.00,0.00,1.05,1.20,1.35,1.30,1.20,1.00,1.30,
|
|---|
| 918 | + 1.00,1.30,1.00,1.30,1.00,1.00,1.00,1.30,1.30,1.30,
|
|---|
| 919 | + 1.00,1.00,1.30,1.30,1.00,1.,1.,1.,1.3,1./
|
|---|
| 920 | DATA ICORR /14*1, 0, 1, 0, 1, 0, 3*1, 3*0, 2*1, 2*0, 4*1, 2*0/
|
|---|
| 921 | C-- SET INTRC TO 0 FOR IPART = 26-29, 33, 34 ( XI'S AND OMEGA'S )
|
|---|
| 922 | C-DH- DATA INTRC /6*0, 1, 0, 12*1, 0, 2*2, 0, 1, 4*0, 3*1, 3*0 /
|
|---|
| 923 | C-- RESET INTRC FOR IPART = 26-29, 33, 34 ( XI'S AND OMEGA'S )
|
|---|
| 924 | DATA INTRC /6*0, 1, 0, 12*1, 0, 2*2, 0, 10*1, 0/
|
|---|
| 925 |
|
|---|
| 926 | C CROSS SECTIONS ON NUCLEUS ARE KNOWN ONLY FOR PIONS AND PROTONS.
|
|---|
| 927 | C THE GENERAL LAW SIGMA(A)=1.25*SIGMA(TOT,PROTON)*A**ALPHA IS VALID
|
|---|
| 928 | C ONLY FOR MOMENTA > 2 GEV.THE PARAMETRIZATION DONE HERE GIVES ONLY
|
|---|
| 929 | C A BEHAVIOUR AVERAGED OVER MOMENTA AND PARTICLE TYPES.
|
|---|
| 930 | C FOR A DETECTOR WITH ONLY A FEW MATERIALS IT'S OF COURSE MUCHBETTER
|
|---|
| 931 | C TO USE TABLES OF THE MEASURED CROSS SECTIONS .
|
|---|
| 932 | C FOR ELEMENTS WITH THE FOLLOWING ATOMIC NUMBERS MEASURED CROSS
|
|---|
| 933 | C SECTIONS ARE AVAILABLE (SEE "PCSDATA").
|
|---|
| 934 |
|
|---|
| 935 | C H AL CU PB
|
|---|
| 936 | DATA CSA /1. ,27.00 ,63.54 ,207.19 /
|
|---|
| 937 | DATA IPART2/9,8,7,11,10,13,12/
|
|---|
| 938 | SAVE ALPHA,ALPHAC,PARTEL,PARTIN,CSA,IPART2,ICORR,INTRC
|
|---|
| 939 | C-----------------------------------------------------------------------
|
|---|
| 940 |
|
|---|
| 941 | IF ( DEBUG ) THEN
|
|---|
| 942 | WRITE(MDEBUG,*) 'CGHSIG:'
|
|---|
| 943 | NPRT(4)=.TRUE.
|
|---|
| 944 | NPRT(9)=.TRUE.
|
|---|
| 945 | ENDIF
|
|---|
| 946 |
|
|---|
| 947 | C --- INITIALISE CGHSIG AND SWITCH TO GHEISHA PARTICLE CODE ---
|
|---|
| 948 | CGHSIG=0.0
|
|---|
| 949 | IF ( LPART .LE. 48 ) THEN
|
|---|
| 950 | IPART = KIPART(LPART)
|
|---|
| 951 | ELSEIF ( LPART .EQ. 201 ) THEN
|
|---|
| 952 | IPART = 30
|
|---|
| 953 | ELSEIF ( LPART .EQ. 301 ) THEN
|
|---|
| 954 | IPART = 31
|
|---|
| 955 | ELSEIF ( LPART .EQ. 402 ) THEN
|
|---|
| 956 | IPART = 32
|
|---|
| 957 | ELSE
|
|---|
| 958 | GOTO 160
|
|---|
| 959 | ENDIF
|
|---|
| 960 |
|
|---|
| 961 | C --- NO INTERACTION FOR GAMMAS, NEUTRINOS, ELECTRONS, POSITRONS, MUONS,
|
|---|
| 962 | C --- NEUTRAL PIONS, NEUTRAL SIGMAS AND ANTISIGMAS AND NEW PARTICLES.
|
|---|
| 963 | IF(INTRC(IPART).EQ. 0) GO TO 160
|
|---|
| 964 | P=PPART
|
|---|
| 965 | EK=EKIN
|
|---|
| 966 |
|
|---|
| 967 | C --- INITIALISE THE CROSS-SECTIONS WITH 0.0 ---
|
|---|
| 968 | DO 10 K=1,KK
|
|---|
| 969 | AIEL(K)=0.0
|
|---|
| 970 | AIIN(K)=0.0
|
|---|
| 971 | AICA(K)=0.0
|
|---|
| 972 | 10 CONTINUE
|
|---|
| 973 | C
|
|---|
| 974 | IF ((IPART .GE. 30) .AND. (IPART .LE. 32)) THEN
|
|---|
| 975 |
|
|---|
| 976 | C --- TAKE GEOMETRICAL CROSS SECTIONS FOR INELASTIC SCATTERING ---
|
|---|
| 977 | C --- OF DEUTERONS, TRITONS AND ALPHAS ---
|
|---|
| 978 | IF ( IPART .EQ. 30 ) THEN
|
|---|
| 979 | APART=2.0**ONETHR
|
|---|
| 980 | ELSEIF ( IPART .EQ. 31 ) THEN
|
|---|
| 981 | APART=3.0**ONETHR
|
|---|
| 982 | ELSEIF ( IPART .EQ. 32 ) THEN
|
|---|
| 983 | APART=4.0**ONETHR
|
|---|
| 984 | ENDIF
|
|---|
| 985 | DO 20 K=1,KK
|
|---|
| 986 | AIIN(K)=49.0*(APART+ACOMP(K)**ONETHR)**2
|
|---|
| 987 | 20 CONTINUE
|
|---|
| 988 | IF (NPRT(9)) WRITE(MDEBUG,10000)
|
|---|
| 989 |
|
|---|
| 990 | ELSEIF ((IPART .EQ. 16) .AND. (EK .LE. 0.0327)) THEN
|
|---|
| 991 |
|
|---|
| 992 | C --- USE TABLES FOR LOW ENERGY NEUTRONS ---
|
|---|
| 993 | C --- GET ENERGY BIN ---
|
|---|
| 994 | JE2=17
|
|---|
| 995 | DO 30 J=2,17
|
|---|
| 996 | IF (EK .LT. ELAB(J)) THEN
|
|---|
| 997 | JE2=J
|
|---|
| 998 | GO TO 40
|
|---|
| 999 | ENDIF
|
|---|
| 1000 | 30 CONTINUE
|
|---|
| 1001 |
|
|---|
| 1002 | 40 JE1=JE2-1
|
|---|
| 1003 | EKX=MAX(EK,1.0E-9)
|
|---|
| 1004 | DELAB=ELAB(JE2)-ELAB(JE1)
|
|---|
| 1005 | DO 70 K=1,KK
|
|---|
| 1006 |
|
|---|
| 1007 | C --- GET A BIN ---
|
|---|
| 1008 | JA2=15
|
|---|
| 1009 | DO 50 J=2,15
|
|---|
| 1010 | IF (ACOMP(K) .LT. CNLWAT(J)) THEN
|
|---|
| 1011 | JA2=J
|
|---|
| 1012 | GO TO 60
|
|---|
| 1013 | ENDIF
|
|---|
| 1014 | 50 CONTINUE
|
|---|
| 1015 |
|
|---|
| 1016 | 60 JA1=JA2-1
|
|---|
| 1017 | DNLWAT=CNLWAT(JA2)-CNLWAT(JA1)
|
|---|
| 1018 |
|
|---|
| 1019 | C --- USE LINEAR INTERPOLATION OR EXTRAPOLATION BY Y=RCE*X+RCA*X+B ---
|
|---|
| 1020 |
|
|---|
| 1021 | C --- ELASTIC CROSS SECTION ---
|
|---|
| 1022 | C --- E INTERPOLATION OR EXTRAPOLATION AT JA1 ---
|
|---|
| 1023 | DY=CNLWEL(JA1,JE2)-CNLWEL(JA1,JE1)
|
|---|
| 1024 | RCE=DY/DELAB
|
|---|
| 1025 | C --- A INTERPOLATION OR EXTRAPOLATION AT JE1 ---
|
|---|
| 1026 | DY=CNLWEL(JA2,JE1)-CNLWEL(JA1,JE1)
|
|---|
| 1027 | RCA=DY/DNLWAT
|
|---|
| 1028 | B=CNLWEL(JA1,JE1)-RCE*ELAB(JE1)-RCA*CNLWAT(JA1)
|
|---|
| 1029 | AIEL(K)=RCE*EK+RCA*ACOMP(K)+B
|
|---|
| 1030 |
|
|---|
| 1031 | C --- INELASTIC CROSS SECTION ---
|
|---|
| 1032 | C --- E INTERPOLATION OR EXTRAPOLATION AT JA1 ---
|
|---|
| 1033 | DY=CNLWIN(JA1,JE2)-CNLWIN(JA1,JE1)
|
|---|
| 1034 | RCE=DY/DELAB
|
|---|
| 1035 | C --- A INTERPOLATION OR EXTRAPOLATION AT JE1 ---
|
|---|
| 1036 | DY=CNLWIN(JA2,JE1)-CNLWIN(JA1,JE1)
|
|---|
| 1037 | RCA=DY/DNLWAT
|
|---|
| 1038 | B=CNLWIN(JA1,JE1)-RCE*ELAB(JE1)-RCA*CNLWAT(JA1)
|
|---|
| 1039 | AIIN(K)=RCE*EK+RCA*ACOMP(K)+B
|
|---|
| 1040 |
|
|---|
| 1041 | IZNO=ZCOMP(K)+0.01
|
|---|
| 1042 | AICA(K)=11.12*CSCAP(IZNO)/(EKX*1.0E6)**0.577
|
|---|
| 1043 | 70 CONTINUE
|
|---|
| 1044 | IF (NPRT(9)) WRITE(MDEBUG,10100)
|
|---|
| 1045 | ELSE
|
|---|
| 1046 |
|
|---|
| 1047 | C --- USE PARAMETRIZATION OF CROSS SECTION DATA FOR ALL OTHER CASES ---
|
|---|
| 1048 |
|
|---|
| 1049 | IF (NPRT(9)) WRITE(MDEBUG,10200)
|
|---|
| 1050 |
|
|---|
| 1051 | C --- GET MOMENTUM BIN ---
|
|---|
| 1052 | J=40
|
|---|
| 1053 | DO 80 I=2,41
|
|---|
| 1054 | IF (P .LT. PLAB(I)) THEN
|
|---|
| 1055 | J=I-1
|
|---|
| 1056 | GO TO 90
|
|---|
| 1057 | ENDIF
|
|---|
| 1058 | 80 CONTINUE
|
|---|
| 1059 |
|
|---|
| 1060 | C --- START WITH CROSS SECTIONS FOR SCATTERING ON FREE PROTONS ---
|
|---|
| 1061 | C --- USE LINEAR INTERPOLATION OR EXTRAPOLATION BY Y=RC*X+B ---
|
|---|
| 1062 | 90 DX=PLAB(J+1)-PLAB(J)
|
|---|
| 1063 | C --- ELASTIC CROSS SECTION ---
|
|---|
| 1064 | DY=CSEL(IPART,J+1)-CSEL(IPART,J)
|
|---|
| 1065 | RC=DY/DX
|
|---|
| 1066 | B=CSEL(IPART,J)-RC*PLAB(J)
|
|---|
| 1067 | AIELIN=RC*P+B
|
|---|
| 1068 | C --- INELASTIC CROSS SECTION ---
|
|---|
| 1069 | DY=CSIN(IPART,J+1)-CSIN(IPART,J)
|
|---|
| 1070 | RC=DY/DX
|
|---|
| 1071 | B=CSIN(IPART,J)-RC*PLAB(J)
|
|---|
| 1072 | AIININ=RC*P+B
|
|---|
| 1073 | ALPH=ALPHA(IPART)
|
|---|
| 1074 | IF(IPART.LT.14) THEN
|
|---|
| 1075 | DY=ALPHAC(J+1)-ALPHAC(J)
|
|---|
| 1076 | RC=DY/DX
|
|---|
| 1077 | B=ALPHAC(J)-RC*PLAB(J)
|
|---|
| 1078 | CORFAC=RC*P+B
|
|---|
| 1079 | ALPH=ALPH*CORFAC
|
|---|
| 1080 |
|
|---|
| 1081 | IPART3=IPART2(IPART-6)
|
|---|
| 1082 |
|
|---|
| 1083 | C --- ELASTIC CROSS SECTION ---
|
|---|
| 1084 | DY=CSEL(IPART3,J+1)-CSEL(IPART3,J)
|
|---|
| 1085 | RC=DY/DX
|
|---|
| 1086 | B=CSEL(IPART3,J)-RC*PLAB(J)
|
|---|
| 1087 | XSECEL=RC*P+B
|
|---|
| 1088 | C --- INELASTIC CROSS SECTION ---
|
|---|
| 1089 | DY=CSIN(IPART3,J+1)-CSIN(IPART3,J)
|
|---|
| 1090 | RC=DY/DX
|
|---|
| 1091 | B=CSIN(IPART3,J)-RC*PLAB(J)
|
|---|
| 1092 | XSECIN=RC*P+B
|
|---|
| 1093 |
|
|---|
| 1094 | ENDIF
|
|---|
| 1095 |
|
|---|
| 1096 | DO 100 K=1,KK
|
|---|
| 1097 | AIEL(K)=AIELIN
|
|---|
| 1098 | AIIN(K)=AIININ
|
|---|
| 1099 |
|
|---|
| 1100 | IF (ACOMP(K) .GE. 1.5) THEN
|
|---|
| 1101 |
|
|---|
| 1102 | C --- A-DEPENDENCE FROM PARAMETRIZATION ---
|
|---|
| 1103 | CREL=1.0
|
|---|
| 1104 | CRIN=1.0
|
|---|
| 1105 | C --- GET MEDIUM BIN 1=HYDR. 2=AL 3=CU 4=PB ---
|
|---|
| 1106 | I=3
|
|---|
| 1107 | IF (ACOMP(K) .LT. 50.0) I=2
|
|---|
| 1108 | IF (ACOMP(K) .GT. 100.0) I=4
|
|---|
| 1109 | IF ((IPART .EQ. 14) .OR. (IPART .EQ. 16)) THEN
|
|---|
| 1110 |
|
|---|
| 1111 | C --- PROTONS AND NEUTRONS ---
|
|---|
| 1112 |
|
|---|
| 1113 | C --- ELASTIC CROSS SECTION ---
|
|---|
| 1114 | DY=CSPNEL(I-1,J+1)-CSPNEL(I-1,J)
|
|---|
| 1115 | RC=DY/DX
|
|---|
| 1116 | B=CSPNEL(I-1,J)-RC*PLAB(J)
|
|---|
| 1117 | XSECEL=RC*P+B
|
|---|
| 1118 | C --- INELASTIC CROSS SECTION ---
|
|---|
| 1119 | DY=CSPNIN(I-1,J+1)-CSPNIN(I-1,J)
|
|---|
| 1120 | RC=DY/DX
|
|---|
| 1121 | B=CSPNIN(I-1,J)-RC*PLAB(J)
|
|---|
| 1122 | XSECIN=RC*P+B
|
|---|
| 1123 | IF (AIEL(K) .GE. 0.001) CREL=XSECEL/(0.36*AIEL(K)*
|
|---|
| 1124 | + CSA(I)**1.17)
|
|---|
| 1125 | AITOT=AIEL(K)+AIIN(K)
|
|---|
| 1126 | IF (AITOT .GE. 0.001) CRIN=XSECIN/(AITOT*CSA(I)**
|
|---|
| 1127 | + ALPH)
|
|---|
| 1128 |
|
|---|
| 1129 | ELSEIF (IPART .LT. 15) THEN
|
|---|
| 1130 |
|
|---|
| 1131 | C --- CALCULATE CORRECTION FACTORS FROM VALUES ON AL,CU,PB FOR ALL ---
|
|---|
| 1132 | C --- MESONS USE LINEAR INTERPOLATION OR EXTRAPOLATION BY Y=RC*X+B ---
|
|---|
| 1133 | C --- NOTE THAT DATA IS ONLY AVAILABLE FOR PIONS AND PROTONS
|
|---|
| 1134 | WGCH=0.5
|
|---|
| 1135 | IF (ACOMP(K) .LT. 20.0) WGCH=0.5+0.5*EXP(-(ACOMP(K)-1.0))
|
|---|
| 1136 | AIEL(K)=WGCH*AIEL(K)+(1.0-WGCH)*XSECEL
|
|---|
| 1137 | AIIN(K)=WGCH*AIIN(K)+(1.0-WGCH)*XSECIN
|
|---|
| 1138 |
|
|---|
| 1139 | C --- THIS SECTION NOT FOR KAONS ---
|
|---|
| 1140 | IF (IPART .LT. 10) THEN
|
|---|
| 1141 |
|
|---|
| 1142 | C --- ELASTIC CROSS SECTION ---
|
|---|
| 1143 | DY=CSPIEL(I-1,J+1)-CSPIEL(I-1,J)
|
|---|
| 1144 | RC=DY/DX
|
|---|
| 1145 | B=CSPIEL(I-1,J)-RC*PLAB(J)
|
|---|
| 1146 | XSPIEL=RC*P+B
|
|---|
| 1147 | C --- INELASTIC CROSS SECTION ---
|
|---|
| 1148 | DY=CSPIIN(I-1,J+1)-CSPIIN(I-1,J)
|
|---|
| 1149 | RC=DY/DX
|
|---|
| 1150 | B=CSPIIN(I-1,J)-RC*PLAB(J)
|
|---|
| 1151 | XSPIIN=RC*P+B
|
|---|
| 1152 |
|
|---|
| 1153 | IF (AIEL(K) .GE. 0.001) CREL=XSPIEL/(0.36* AIEL(K)
|
|---|
| 1154 | + *CSA(I)**1.17)
|
|---|
| 1155 | AITOT=AIEL(K)+AIIN(K)
|
|---|
| 1156 | IF (AITOT .GE. 0.001) CRIN=XSPIIN/(AITOT*CSA(I)
|
|---|
| 1157 | + **ALPH)
|
|---|
| 1158 | ENDIF
|
|---|
| 1159 | ENDIF
|
|---|
| 1160 | AIIN(K)=CRIN*(AIIN(K)+AIEL(K))*ACOMP(K)**ALPH
|
|---|
| 1161 | AIEL(K)=CREL*0.36*AIEL(K)*ACOMP(K)**1.17
|
|---|
| 1162 | AIEL(K)=AIEL(K)*PARTEL(IPART)
|
|---|
| 1163 | AIIN(K)=AIIN(K)*PARTIN(IPART)
|
|---|
| 1164 | ENDIF
|
|---|
| 1165 | 100 CONTINUE
|
|---|
| 1166 |
|
|---|
| 1167 | ENDIF
|
|---|
| 1168 |
|
|---|
| 1169 | C --- CALCULATE INTERACTION PROBABILITY ---
|
|---|
| 1170 |
|
|---|
| 1171 | ALAM=0.0
|
|---|
| 1172 | DO 150 K=1,KK
|
|---|
| 1173 | AIEL(K) = AIEL(K)*WCOMP(K)
|
|---|
| 1174 | AIIN(K) = AIIN(K)*WCOMP(K)
|
|---|
| 1175 | AICA(K) = AICA(K)*WCOMP(K)
|
|---|
| 1176 | ALAM = ALAM + AIEL(K) + AIIN(K) + AICA(K)
|
|---|
| 1177 | 150 CONTINUE
|
|---|
| 1178 |
|
|---|
| 1179 | C --- PASS THE CROSS SECTION (MBARN) TO CORSIKA ---
|
|---|
| 1180 | CGHSIG=ALAM
|
|---|
| 1181 |
|
|---|
| 1182 | GO TO 999
|
|---|
| 1183 |
|
|---|
| 1184 | C --- PRINTOUT OF SKIPPED PARTICLES IN CASE OF INTERFACE DEBUG ---
|
|---|
| 1185 | 160 IF (NPRT(9)) WRITE(MDEBUG,10300) IPART
|
|---|
| 1186 | 10000 FORMAT(' *CGHSIG* GEOM X-SECT. FOR INEL. SCAT. OF D,T AND ALPHA')
|
|---|
| 1187 | 10100 FORMAT(' *CGHSIG* X-SECT. FROM LOW ENERGY NEUTRON TABLES')
|
|---|
| 1188 | 10200 FORMAT(' *CGHSIG* X-SECT. FROM PARAMETRIZATION OF DATA')
|
|---|
| 1189 | 10300 FORMAT(' *CGHSIG* GHEISHA PARTICLE ',I3,' SKIPPED')
|
|---|
| 1190 | 999 RETURN
|
|---|
| 1191 | END
|
|---|
