DOUBLE PRECISION FUNCTION CGHSIG( PPART,EKIN,LPART ) C----------------------------------------------------------------------- C C(ORSIKA) GH(EISHA) SIG(MA) C C CALCULATION OF THE PROBABILITIES FOR (IN)ELASTIC INTERACTIONS *** C THIS FUNCTION IS CALLED FROM BOX2 C C ORIGIN : F.CARMINATI, H.FESEFELDT (ROUTINE GHESIG) C REDESIGN: P. GABRIEL IK1 FZK KARLSRUHE C----------------------------------------------------------------------- C *** IPART DENOTES THE GHEISHA PARTICLE INDEX *** C C CONVENTION : C C PARTICLE IPART C ------------------------------ C GAMMA 1 C NEUTRINO 2 C POSITRON 3 C ELECTRON 4 C MUON + 5 C MUON - 6 C PION + 7 C PION 0 8 C PION - 9 C KAON + 10 C KAON 0 S 11 C KAON 0 L 12 C KAON - 13 C PROTON 14 C PROTON BAR 15 C NEUTRON 16 C NEUTRON BAR 17 C LAMBDA 18 C LAMBDA BAR 19 C SIGMA + 20 C SIGMA 0 21 C SIGMA - 22 C SIGMA + BAR 23 C SIGMA 0 BAR 24 C SIGMA - BAR 25 C XSI 0 26 C XSI - 27 C XSI 0 BAR 28 C XSI - BAR 29 C DEUTERON 30 C TRITON 31 C ALPHA 32 C OMEGA - 33 C OMEGA - BAR 34 C NEW PARTICLES 35 C C----------------------------------------------------------------------- *KEEP,CGCOMP. PARAMETER (KK=3) COMMON/CGCOMP/ ACOMP,ZCOMP,WCOMP REAL ACOMP(KK),ZCOMP(KK),WCOMP(KK) *KEEP,RUNPAR. COMMON /RUNPAR/ FIXHEI,THICK0,HILOECM,HILOELB, * STEPFC,NRRUN,NSHOW,PATAPE,MONIIN, * MONIOU,MDEBUG,NUCNUC, * CETAPE, * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL, * N1STTR,MDBASE, * DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR, * FIX1I,FMUADD,FNKG,FPRINT,FDBASE * ,GHEISH,GHESIG COMMON /RUNPAC/ DSN,HOST,USER DOUBLE PRECISION FIXHEI,THICK0,HILOECM,HILOELB REAL STEPFC INTEGER NRRUN,NSHOW,PATAPE,MONIIN,MONIOU,MDEBUG,NUCNUC, * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL, * N1STTR,MDBASE INTEGER CETAPE CHARACTER*79 DSN CHARACTER*20 HOST,USER LOGICAL DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR, * FIX1I,FMUADD,FNKG,FPRINT,FDBASE * ,GHEISH,GHESIG *KEND. COMMON/GSECTI/ AIEL(20),AIIN(20),AIFI(20),AICA(20),ALAM,K0FLAG INTEGER K0FLAG REAL AIEL,AIIN,AIFI,AICA,ALAM C --- GHEISHA COMMONS --- COMMON/RESULT/XEND,YEND,ZEND,RCA,RCE,AMAS,NCH,TOF,PX,PY,PZ, * USERW,INTCT,P,EN,EK,AMASQ,DELTN,ITK,NTK,IPART,IND, * LCALO,ICEL,SINL,COSL,SINP,COSP, * XOLD,YOLD,ZOLD,POLD,PXOLD,PYOLD,PZOLD, * XSCAT,YSCAT,ZSCAT,PSCAT,PXSCAT,PYSCAT,PZSCAT REAL NCH,INTCT COMMON/PRNTFL/INBCD,NEWBCD,INBIN,NEWBIN,NPEVT,NEVTP,LPRT,NPRT(10) LOGICAL LPRT,NPRT DIMENSION ALPHA(35),ALPHAC(41),IPART2(7),CSA(4) DIMENSION PARTEL(35),PARTIN(35),ICORR(35),INTRC(35) C --- DIMENSION STATEMENTS FOR CROSS SECTION DATA --- DIMENSION PLAB(41),CSEL(35,41),CSIN(35,41),CSPIEL(3,41), $ CSPIIN(3,41),CSPNEL(3,41),CSPNIN(3,41), $ ELAB(17),CNLWAT(15),CNLWEL(15,17),CNLWIN(15,17), $ CSCAP(100) C --- DIMENSION STMTS. FOR GEANT/GHEISHA PARTICLE CODE CONVERSIONS --- C --- KIPART(I)=GHEISHA CODE CORRESPONDING TO GEANT CODE I --- C --- IKPART(I)=GEANT CODE CORRESPONDING TO GHEISHA CODE I --- DIMENSION KIPART(48),IKPART(35) C --- CROSS-SECTION DATA BY "PCSDAT" 01-FEB-1989 --- DATA PLAB / $ 0.00000E+00, 0.10000 , 0.15000 , 0.20000 , 0.25000 , $ 0.30000 , 0.35000 , 0.40000 , 0.45000 , 0.50000 , $ 0.55000 , 0.60000 , 0.65000 , 0.70000 , 0.75000 , $ 0.80000 , 0.85000 , 0.90000 , 0.95000 , 1.0000 , $ 1.1000 , 1.2000 , 1.3000 , 1.4000 , 1.5000 , $ 1.6000 , 1.8000 , 2.0000 , 2.2000 , 2.4000 , $ 2.6000 , 2.8000 , 3.0000 , 4.0000 , 5.0000 , $ 6.0000 , 8.0000 , 10.000 , 20.000 , 100.00 , $ 1000.0 / DATA ((CSEL(I,J),I=1,6),J=1,41) / 246 * 0. / DATA (CSEL( 8,J),J=1,41) / 41 * 0. / DATA (CSEL(21,J),J=1,41) / 41 * 0. / DATA (CSEL(24,J),J=1,41) / 41 * 0. / DATA ((CSEL(I,J),I=30,32),J=1,41) / 123 * 0. / DATA (CSEL(35,J),J=1,41) / 41 * 0. / DATA (CSEL( 7,J),J=1,41) / $ 0.00000E+00, 6.0000 , 20.000 , 71.000 , 155.00 , $ 195.00 , 130.00 , 78.000 , 60.000 , 32.000 , $ 23.500 , 18.500 , 15.000 , 12.500 , 10.000 , $ 9.1000 , 8.6000 , 8.8000 , 9.5000 , 10.600 , $ 13.000 , 15.500 , 17.100 , 17.200 , 16.200 , $ 15.000 , 12.300 , 10.200 , 9.0000 , 8.0000 , $ 7.3000 , 6.8000 , 6.5000 , 5.8000 , 5.4000 , $ 5.2000 , 5.0000 , 4.9000 , 3.8000 , 3.2000 , $ 3.5000 / DATA (CSEL( 9,J),J=1,41) / $ 0.00000E+00, 1.0000 , 3.0000 , 8.0000 , 18.000 , $ 25.000 , 27.500 , 12.300 , 10.600 , 11.000 , $ 12.500 , 14.500 , 17.000 , 19.400 , 19.800 , $ 16.800 , 14.000 , 14.800 , 20.000 , 26.100 , $ 19.500 , 15.000 , 12.800 , 11.500 , 10.500 , $ 9.8000 , 8.8000 , 8.2000 , 7.8000 , 7.5000 , $ 7.2000 , 7.0000 , 6.8000 , 6.1000 , 5.7000 , $ 5.4000 , 4.9000 , 4.6000 , 4.0000 , 3.3000 , $ 3.5000 / DATA (CSEL(10,J),J=1,41) / $ 10.000 , 11.200 , 11.300 , 11.400 , 11.500 , $ 11.600 , 11.800 , 12.000 , 12.100 , 12.200 , $ 12.300 , 12.400 , 12.500 , 12.500 , 12.500 , $ 12.400 , 12.300 , 12.200 , 12.000 , 11.800 , $ 11.200 , 11.500 , 9.9000 , 9.4000 , 8.8000 , $ 8.4000 , 7.5000 , 6.9000 , 6.3000 , 5.9000 , $ 5.5000 , 5.2000 , 5.0000 , 4.0000 , 3.5000 , $ 3.3000 , 3.1000 , 3.1000 , 3.0000 , 2.5000 , $ 3.0000 / DATA (CSEL(11,J),J=1,41) / $ 10.000 , 11.200 , 11.300 , 11.400 , 11.500 , $ 11.600 , 11.800 , 12.000 , 12.100 , 12.200 , $ 12.300 , 12.400 , 12.500 , 12.500 , 12.500 , $ 12.400 , 12.300 , 12.200 , 12.000 , 11.800 , $ 11.200 , 11.500 , 9.9000 , 9.4000 , 8.8000 , $ 8.4000 , 7.5000 , 6.9000 , 6.3000 , 5.9000 , $ 5.5000 , 5.2000 , 5.0000 , 4.0000 , 3.5000 , $ 3.3000 , 3.1000 , 3.1000 , 3.0000 , 2.5000 , $ 3.0000 / DATA (CSEL(12,J),J=1,41) / $ 160.83 , 82.800 , 58.575 , 43.683 , 34.792 , $ 28.650 , 24.367 , 20.917 , 18.192 , 16.300 , $ 14.608 , 13.017 , 12.250 , 11.700 , 12.017 , $ 14.075 , 15.842 , 16.433 , 16.042 , 15.008 , $ 12.575 , 10.708 , 9.2000 , 8.0167 , 7.2833 , $ 7.0750 , 6.6333 , 6.1250 , 5.6583 , 5.2750 , $ 4.9333 , 4.6250 , 4.4583 , 3.7333 , 3.3833 , $ 3.1833 , 2.9833 , 2.7500 , 2.3667 , 2.2000 , $ 2.6000 / DATA (CSEL(13,J),J=1,41) / $ 300.00 , 140.00 , 97.000 , 70.000 , 55.000 , $ 45.000 , 37.000 , 31.000 , 26.000 , 23.000 , $ 20.000 , 17.000 , 15.500 , 14.500 , 14.700 , $ 18.500 , 22.000 , 23.000 , 22.500 , 20.700 , $ 16.500 , 14.000 , 11.500 , 9.6000 , 8.6000 , $ 8.5000 , 8.3000 , 7.6000 , 7.0000 , 6.4000 , $ 5.9000 , 5.5000 , 5.3000 , 4.4000 , 4.1000 , $ 3.9000 , 3.7000 , 3.3000 , 2.6000 , 2.5000 , $ 3.0000 / DATA (CSEL(14,J),J=1,41) / $ 1100.0 , 115.00 , 105.00 , 100.00 , 56.000 , $ 40.000 , 27.000 , 22.000 , 21.000 , 20.000 , $ 20.000 , 20.000 , 20.500 , 21.000 , 22.000 , $ 23.000 , 24.000 , 24.000 , 24.400 , 24.500 , $ 25.000 , 25.500 , 26.000 , 26.500 , 27.000 , $ 27.000 , 26.000 , 23.000 , 21.500 , 20.000 , $ 19.000 , 18.000 , 17.000 , 13.000 , 11.500 , $ 10.300 , 9.4000 , 9.0000 , 8.8000 , 7.0000 , $ 7.5000 / DATA (CSEL(15,J),J=1,41) / $ 200.00 , 163.00 , 141.00 , 120.00 , 111.00 , $ 99.500 , 92.500 , 86.500 , 82.000 , 78.000 , $ 74.000 , 71.000 , 67.500 , 65.000 , 62.500 , $ 59.700 , 58.100 , 56.300 , 54.700 , 52.700 , $ 50.000 , 48.400 , 47.000 , 46.000 , 45.200 , $ 42.800 , 39.200 , 36.300 , 32.800 , 30.400 , $ 28.100 , 26.300 , 24.500 , 19.250 , 16.840 , $ 14.600 , 12.340 , 11.210 , 8.8500 , 7.5000 , $ 7.5000 / DATA (CSEL(16,J),J=1,41) / $ 4200.0 , 440.00 , 420.00 , 400.00 , 230.00 , $ 160.00 , 105.00 , 80.000 , 62.000 , 50.000 , $ 45.000 , 41.000 , 38.000 , 36.000 , 35.000 , $ 34.000 , 33.000 , 32.000 , 31.500 , 31.000 , $ 30.500 , 30.000 , 29.500 , 29.000 , 28.500 , $ 28.000 , 26.000 , 23.000 , 21.500 , 20.000 , $ 19.000 , 18.000 , 17.000 , 13.000 , 11.500 , $ 10.300 , 9.4000 , 9.0000 , 8.8000 , 7.0000 , $ 7.5000 / DATA (CSEL(17,J),J=1,41) / $ 185.88 , 133.23 , 119.37 , 102.86 , 93.102 , $ 82.752 , 76.205 , 71.008 , 67.366 , 64.096 , $ 60.891 , 58.501 , 55.735 , 53.773 , 51.839 , $ 49.671 , 48.485 , 47.045 , 45.803 , 44.306 , $ 42.623 , 41.786 , 41.115 , 40.630 , 40.129 , $ 38.242 , 35.233 , 32.662 , 29.639 , 27.573 , $ 25.536 , 23.948 , 22.356 , 17.723 , 15.614 , $ 13.653 , 11.675 , 10.653 , 8.6198 , 7.4464 , $ 7.4821 / DATA (CSEL(18,J),J=1,41) / $ 1100.0 , 115.00 , 105.00 , 100.00 , 56.000 , $ 40.000 , 27.000 , 22.000 , 21.000 , 20.000 , $ 20.000 , 19.067 , 19.333 , 19.500 , 19.833 , $ 20.567 , 21.800 , 22.900 , 23.869 , 23.809 , $ 22.161 , 21.488 , 19.732 , 19.433 , 19.345 , $ 19.029 , 18.121 , 16.280 , 15.258 , 14.280 , $ 13.644 , 12.963 , 12.316 , 9.5333 , 8.4333 , $ 7.5728 , 6.9696 , 6.7518 , 6.6175 , 5.6000 , $ 6.1145 / DATA (CSEL(19,J),J=1,41) / $ 157.65 , 73.701 , 76.096 , 68.571 , 57.305 , $ 49.257 , 43.616 , 40.024 , 38.098 , 36.287 , $ 34.674 , 33.105 , 31.712 , 30.685 , 29.613 , $ 28.602 , 28.336 , 28.075 , 27.786 , 27.215 , $ 26.380 , 26.146 , 25.108 , 24.783 , 24.360 , $ 23.219 , 21.431 , 20.095 , 18.382 , 17.267 , $ 16.100 , 15.175 , 14.271 , 11.573 , 10.305 , $ 9.1471 , 8.0149 , 7.4349 , 6.2499 , 5.8928 , $ 6.0774 / DATA (CSEL(20,J),J=1,41) / $ 1100.0 , 115.00 , 105.00 , 100.00 , 56.000 , $ 40.000 , 27.000 , 22.000 , 21.000 , 20.000 , $ 20.000 , 19.067 , 19.333 , 19.500 , 19.833 , $ 20.567 , 21.800 , 22.900 , 23.869 , 23.809 , $ 22.161 , 21.488 , 19.732 , 19.433 , 19.345 , $ 19.029 , 18.121 , 16.280 , 15.258 , 14.280 , $ 13.644 , 12.963 , 12.316 , 9.5333 , 8.4333 , $ 7.5728 , 6.9696 , 6.7518 , 6.6175 , 5.6000 , $ 6.1145 / DATA (CSEL(22,J),J=1,41) / $ 1100.0 , 115.00 , 105.00 , 100.00 , 56.000 , $ 40.000 , 27.000 , 22.000 , 21.000 , 20.000 , $ 20.000 , 19.067 , 19.333 , 19.500 , 19.833 , $ 20.567 , 21.800 , 22.900 , 23.869 , 23.809 , $ 22.161 , 21.488 , 19.732 , 19.433 , 19.345 , $ 19.029 , 18.121 , 16.280 , 15.258 , 14.280 , $ 13.644 , 12.963 , 12.316 , 9.5333 , 8.4333 , $ 7.5728 , 6.9696 , 6.7518 , 6.6175 , 5.6000 , $ 6.1145 / DATA (CSEL(23,J),J=1,41) / $ 185.88 , 133.23 , 119.37 , 102.86 , 93.102 , $ 82.752 , 76.205 , 71.008 , 67.366 , 64.096 , $ 60.891 , 58.104 , 55.241 , 53.140 , 50.934 , $ 48.660 , 47.566 , 46.585 , 45.581 , 44.003 , $ 41.134 , 39.374 , 36.878 , 35.523 , 34.503 , $ 32.334 , 29.365 , 27.370 , 24.705 , 22.921 , $ 21.229 , 19.879 , 18.559 , 14.625 , 12.758 , $ 11.041 , 9.3440 , 8.5484 , 6.7104 , 6.0000 , $ 6.1131 / DATA (CSEL(25,J),J=1,41) / $ 157.65 , 73.701 , 76.096 , 68.571 , 57.305 , $ 49.257 , 43.616 , 40.024 , 38.098 , 36.287 , $ 34.674 , 33.105 , 31.712 , 30.685 , 29.613 , $ 28.602 , 28.336 , 28.075 , 27.786 , 27.215 , $ 26.380 , 26.146 , 25.108 , 24.783 , 24.360 , $ 23.219 , 21.431 , 20.095 , 18.382 , 17.267 , $ 16.100 , 15.175 , 14.271 , 11.573 , 10.305 , $ 9.1471 , 8.0149 , 7.4349 , 6.2499 , 5.8928 , $ 6.0774 / DATA (CSEL(26,J),J=1,41) / $ 1100.0 , 115.00 , 105.00 , 100.00 , 56.000 , $ 40.000 , 27.000 , 22.000 , 21.000 , 20.000 , $ 20.000 , 18.133 , 18.167 , 18.000 , 17.667 , $ 18.133 , 19.600 , 21.800 , 23.338 , 23.118 , $ 19.323 , 17.476 , 13.464 , 12.367 , 11.691 , $ 11.057 , 10.242 , 9.5593 , 9.0151 , 8.5591 , $ 8.2884 , 7.9253 , 7.6311 , 6.0667 , 5.3667 , $ 4.8456 , 4.5392 , 4.5036 , 4.4351 , 4.2000 , $ 4.7289 / DATA (CSEL(27,J),J=1,41) / $ 1100.0 , 115.00 , 105.00 , 100.00 , 56.000 , $ 40.000 , 27.000 , 22.000 , 21.000 , 20.000 , $ 20.000 , 18.133 , 18.167 , 18.000 , 17.667 , $ 18.133 , 19.600 , 21.800 , 23.338 , 23.118 , $ 19.323 , 17.476 , 13.464 , 12.367 , 11.691 , $ 11.057 , 10.242 , 9.5593 , 9.0151 , 8.5591 , $ 8.2884 , 7.9253 , 7.6311 , 6.0667 , 5.3667 , $ 4.8456 , 4.5392 , 4.5036 , 4.4351 , 4.2000 , $ 4.7289 / DATA (CSEL(28,J),J=1,41) / $ 157.65 , 73.701 , 76.096 , 68.571 , 57.305 , $ 49.257 , 43.616 , 40.024 , 38.098 , 36.287 , $ 34.674 , 32.708 , 31.218 , 30.052 , 28.707 , $ 27.591 , 27.417 , 27.615 , 27.564 , 26.913 , $ 24.891 , 23.734 , 20.871 , 19.677 , 18.734 , $ 17.311 , 15.563 , 14.803 , 13.448 , 12.615 , $ 11.794 , 11.106 , 10.474 , 8.4745 , 7.4498 , $ 6.5350 , 5.6835 , 5.3300 , 4.3406 , 4.4464 , $ 4.7083 / DATA (CSEL(29,J),J=1,41) / $ 143.53 , 43.935 , 54.462 , 51.429 , 39.407 , $ 32.510 , 27.321 , 24.532 , 23.465 , 22.383 , $ 21.566 , 20.209 , 19.453 , 18.825 , 18.046 , $ 17.562 , 17.802 , 18.360 , 18.667 , 18.519 , $ 17.514 , 17.120 , 14.985 , 14.306 , 13.663 , $ 12.753 , 11.596 , 11.165 , 10.287 , 9.7882 , $ 9.2294 , 8.7539 , 8.3300 , 6.9480 , 6.2234 , $ 5.5881 , 5.0189 , 4.7733 , 4.1104 , 4.3929 , $ 4.6905 / DATA (CSEL(33,J),J=1,41) / $ 1100.0 , 115.00 , 105.00 , 100.00 , 56.000 , $ 40.000 , 27.000 , 22.000 , 21.000 , 20.000 , $ 20.000 , 18.133 , 18.167 , 18.000 , 17.667 , $ 18.133 , 19.600 , 21.800 , 23.338 , 23.118 , $ 19.323 , 17.476 , 13.464 , 12.367 , 11.691 , $ 11.057 , 10.242 , 9.5593 , 9.0151 , 8.5591 , $ 8.2884 , 7.9253 , 7.6311 , 6.0667 , 5.3667 , $ 4.8456 , 4.5392 , 4.5036 , 4.4351 , 4.2000 , $ 4.7289 / DATA (CSEL(34,J),J=1,41) / $ 143.53 , 43.935 , 54.462 , 51.429 , 39.407 , $ 32.510 , 27.321 , 24.532 , 23.465 , 22.383 , $ 21.566 , 20.209 , 19.453 , 18.825 , 18.046 , $ 17.562 , 17.802 , 18.360 , 18.667 , 18.519 , $ 17.514 , 17.120 , 14.985 , 14.306 , 13.663 , $ 12.753 , 11.596 , 11.165 , 10.287 , 9.7882 , $ 9.2294 , 8.7539 , 8.3300 , 6.9480 , 6.2234 , $ 5.5881 , 5.0189 , 4.7733 , 4.1104 , 4.3929 , $ 4.6905 / DATA ((CSIN(I,J),I=1,6),J=1,41) / 246 * 0. / DATA (CSIN( 8,J),J=1,41) / 41 * 0. / DATA (CSIN(21,J),J=1,41) / 41 * 0. / DATA (CSIN(24,J),J=1,41) / 41 * 0. / DATA ((CSIN(I,J),I=30,32),J=1,41) / 123 * 0. / DATA (CSIN(35,J),J=1,41) / 41 * 0. / DATA (CSIN( 7,J),J=1,41) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.50000 , 1.2000 , 1.7000 , $ 2.2500 , 3.0000 , 3.6000 , 4.5000 , 5.4000 , $ 6.3000 , 8.6000 , 9.0000 , 10.000 , 11.500 , $ 14.000 , 17.000 , 19.500 , 22.000 , 24.000 , $ 21.500 , 18.500 , 19.000 , 20.500 , 22.200 , $ 23.000 , 23.300 , 23.000 , 21.000 , 20.500 , $ 20.200 , 20.100 , 20.000 , 20.000 , 20.000 , $ 21.000 / DATA (CSIN( 9,J),J=1,41) / $ 0.00000E+00, 3.0000 , 9.2000 , 20.500 , 36.500 , $ 45.000 , 28.000 , 19.500 , 15.500 , 14.200 , $ 15.500 , 17.500 , 20.000 , 23.000 , 26.000 , $ 20.000 , 23.000 , 26.500 , 32.000 , 35.000 , $ 28.500 , 22.000 , 22.500 , 23.500 , 24.000 , $ 24.500 , 26.000 , 27.500 , 27.500 , 27.000 , $ 26.500 , 25.500 , 25.000 , 23.000 , 22.500 , $ 22.200 , 22.000 , 22.000 , 21.200 , 20.700 , $ 21.000 / DATA (CSIN(10,J),J=1,41) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.50000 , 1.5000 , 2.7000 , 3.8000 , 4.8000 , $ 6.5000 , 7.6000 , 8.4000 , 9.0000 , 9.4000 , $ 9.8000 , 10.500 , 11.000 , 11.500 , 11.800 , $ 12.200 , 12.400 , 12.600 , 13.200 , 13.500 , $ 13.700 , 14.000 , 14.200 , 14.500 , 16.400 , $ 17.000 / DATA (CSIN(11,J),J=1,41) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.50000 , 1.5000 , 2.7000 , 3.8000 , 4.8000 , $ 6.5000 , 7.6000 , 8.4000 , 9.0000 , 9.4000 , $ 9.8000 , 10.500 , 11.000 , 11.500 , 11.800 , $ 12.200 , 12.400 , 12.600 , 13.200 , 13.500 , $ 13.700 , 14.000 , 14.200 , 14.500 , 16.400 , $ 17.000 / DATA (CSIN(12,J),J=1,41) / $ 266.67 , 133.33 , 83.333 , 57.083 , 44.500 , $ 33.250 , 24.583 , 20.833 , 18.333 , 16.083 , $ 15.625 , 15.083 , 14.833 , 15.083 , 15.833 , $ 17.042 , 18.958 , 20.758 , 22.533 , 22.825 , $ 21.250 , 18.567 , 17.767 , 18.100 , 19.933 , $ 20.783 , 21.225 , 21.000 , 20.558 , 20.258 , $ 20.017 , 19.767 , 19.600 , 19.183 , 18.850 , $ 18.575 , 18.350 , 18.175 , 17.808 , 17.558 , $ 19.250 / DATA (CSIN(13,J),J=1,41) / $ 400.00 , 200.00 , 120.00 , 81.000 , 62.000 , $ 47.000 , 35.000 , 28.000 , 24.000 , 21.000 , $ 19.500 , 19.000 , 18.800 , 19.000 , 20.000 , $ 21.000 , 23.000 , 25.000 , 27.000 , 27.500 , $ 25.500 , 22.000 , 20.800 , 21.000 , 23.000 , $ 24.000 , 24.000 , 23.800 , 23.000 , 22.500 , $ 22.000 , 21.600 , 21.400 , 21.000 , 20.500 , $ 20.200 , 19.800 , 19.500 , 18.600 , 17.500 , $ 20.000 / DATA (CSIN(14,J),J=1,41) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.10000 , 1.5000 , $ 7.0000 , 12.000 , 17.000 , 19.500 , 20.500 , $ 22.000 , 23.500 , 24.800 , 25.800 , 26.500 , $ 27.000 , 27.500 , 28.000 , 30.000 , 31.000 , $ 32.000 , 32.500 , 32.500 , 33.000 , 33.500 , $ 34.000 / DATA (CSIN(15,J),J=1,41) / $ 1500.0 , 1160.0 , 310.00 , 230.00 , 178.00 , $ 153.00 , 134.00 , 124.00 , 113.00 , 106.00 , $ 101.00 , 96.000 , 92.000 , 89.000 , 87.000 , $ 84.000 , 81.000 , 78.500 , 76.500 , 75.000 , $ 72.000 , 70.000 , 68.000 , 64.500 , 63.000 , $ 62.000 , 61.000 , 59.500 , 58.500 , 56.500 , $ 56.500 , 56.000 , 55.500 , 52.000 , 50.000 , $ 48.000 , 45.000 , 44.000 , 39.200 , 34.500 , $ 34.500 / DATA (CSIN(16,J),J=1,41) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.10000 , 1.5000 , $ 7.0000 , 12.000 , 17.000 , 19.500 , 20.500 , $ 22.000 , 23.500 , 24.800 , 25.800 , 26.500 , $ 27.000 , 27.500 , 28.000 , 30.000 , 31.000 , $ 32.000 , 32.500 , 32.500 , 33.000 , 33.500 , $ 34.000 / DATA (CSIN(17,J),J=1,41) / $ 1394.1 , 948.17 , 262.43 , 197.14 , 149.30 , $ 127.25 , 110.39 , 101.79 , 92.834 , 87.104 , $ 83.109 , 79.099 , 75.965 , 73.627 , 72.161 , $ 69.889 , 67.595 , 65.595 , 64.057 , 63.054 , $ 61.377 , 60.434 , 59.485 , 56.970 , 55.931 , $ 55.398 , 54.827 , 53.538 , 52.861 , 51.247 , $ 51.344 , 50.992 , 50.644 , 47.876 , 46.358 , $ 44.887 , 42.577 , 41.815 , 38.180 , 34.254 , $ 34.418 / DATA (CSIN(18,J),J=1,41) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.97815E-01, 1.4577 , $ 6.2052 , 10.112 , 12.902 , 14.300 , 14.688 , $ 15.505 , 16.379 , 17.554 , 18.309 , 18.920 , $ 19.389 , 19.804 , 20.284 , 22.000 , 22.733 , $ 23.527 , 24.097 , 24.382 , 24.816 , 26.800 , $ 27.719 / DATA (CSIN(19,J),J=1,41) / $ 1182.4 , 524.50 , 167.30 , 131.43 , 91.895 , $ 75.743 , 63.184 , 57.376 , 52.502 , 49.313 , $ 47.326 , 44.762 , 43.222 , 42.015 , 41.221 , $ 40.244 , 39.504 , 39.145 , 38.860 , 38.731 , $ 37.987 , 37.814 , 36.326 , 34.750 , 33.953 , $ 33.635 , 33.349 , 32.938 , 32.785 , 32.092 , $ 32.373 , 32.312 , 32.329 , 31.261 , 30.597 , $ 30.073 , 29.228 , 29.182 , 27.683 , 27.107 , $ 27.956 / DATA (CSIN(20,J),J=1,41) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.97815E-01, 1.4577 , $ 6.2052 , 10.112 , 12.902 , 14.300 , 14.688 , $ 15.505 , 16.379 , 17.554 , 18.309 , 18.920 , $ 19.389 , 19.804 , 20.284 , 22.000 , 22.733 , $ 23.527 , 24.097 , 24.382 , 24.816 , 26.800 , $ 27.719 / DATA (CSIN(22,J),J=1,41) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.97815E-01, 1.4577 , $ 6.2052 , 10.112 , 12.902 , 14.300 , 14.688 , $ 15.505 , 16.379 , 17.554 , 18.309 , 18.920 , $ 19.389 , 19.804 , 20.284 , 22.000 , 22.733 , $ 23.527 , 24.097 , 24.382 , 24.816 , 26.800 , $ 27.719 / DATA (CSIN(23,J),J=1,41) / $ 1394.1 , 948.17 , 262.43 , 197.14 , 149.30 , $ 127.25 , 110.39 , 101.79 , 92.834 , 87.104 , $ 83.109 , 78.563 , 75.292 , 72.760 , 70.900 , $ 68.467 , 66.314 , 64.955 , 63.746 , 62.623 , $ 59.233 , 56.946 , 53.355 , 49.810 , 48.090 , $ 46.839 , 45.695 , 44.863 , 44.062 , 42.599 , $ 42.684 , 42.328 , 42.041 , 39.508 , 37.880 , $ 36.299 , 34.075 , 33.553 , 29.723 , 27.600 , $ 28.120 / DATA (CSIN(25,J),J=1,41) / $ 1182.4 , 524.50 , 167.30 , 131.43 , 91.895 , $ 75.743 , 63.184 , 57.376 , 52.502 , 49.313 , $ 47.326 , 44.762 , 43.222 , 42.015 , 41.221 , $ 40.244 , 39.504 , 39.145 , 38.860 , 38.731 , $ 37.987 , 37.814 , 36.326 , 34.750 , 33.953 , $ 33.635 , 33.349 , 32.938 , 32.785 , 32.092 , $ 32.373 , 32.312 , 32.329 , 31.261 , 30.597 , $ 30.073 , 29.228 , 29.182 , 27.683 , 27.107 , $ 27.956 / DATA (CSIN(26,J),J=1,41) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.95639E-01, 1.4154 , $ 5.4104 , 8.2240 , 8.8031 , 9.1000 , 8.8761 , $ 9.0095 , 9.2576 , 10.307 , 10.818 , 11.341 , $ 11.778 , 12.108 , 12.569 , 14.000 , 14.467 , $ 15.054 , 15.694 , 16.263 , 16.632 , 20.100 , $ 21.438 / DATA (CSIN(27,J),J=1,41) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.95639E-01, 1.4154 , $ 5.4104 , 8.2240 , 8.8031 , 9.1000 , 8.8761 , $ 9.0095 , 9.2576 , 10.307 , 10.818 , 11.341 , $ 11.778 , 12.108 , 12.569 , 14.000 , 14.467 , $ 15.054 , 15.694 , 16.263 , 16.632 , 20.100 , $ 21.438 / DATA (CSIN(28,J),J=1,41) / $ 1182.4 , 524.50 , 167.30 , 131.43 , 91.895 , $ 75.743 , 63.184 , 57.376 , 52.502 , 49.313 , $ 47.326 , 44.225 , 42.549 , 41.148 , 39.960 , $ 38.822 , 38.223 , 38.505 , 38.549 , 38.301 , $ 35.843 , 34.326 , 30.196 , 27.590 , 26.112 , $ 25.076 , 24.217 , 24.264 , 23.985 , 23.445 , $ 23.713 , 23.647 , 23.726 , 22.892 , 22.119 , $ 21.485 , 20.726 , 20.921 , 19.226 , 20.454 , $ 21.658 / DATA (CSIN(29,J),J=1,41) / $ 1076.5 , 312.66 , 119.74 , 98.571 , 63.193 , $ 49.990 , 39.579 , 35.168 , 32.335 , 30.417 , $ 29.434 , 27.325 , 26.514 , 25.775 , 25.120 , $ 24.711 , 24.818 , 25.600 , 26.106 , 26.355 , $ 25.220 , 24.760 , 21.681 , 20.060 , 19.044 , $ 18.474 , 18.044 , 18.301 , 18.347 , 18.192 , $ 18.557 , 18.639 , 18.870 , 18.769 , 18.478 , $ 18.372 , 18.302 , 18.735 , 18.206 , 20.207 , $ 21.576 / DATA (CSIN(33,J),J=1,41) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.95639E-01, 1.4154 , $ 5.4104 , 8.2240 , 8.8031 , 9.1000 , 8.8761 , $ 9.0095 , 9.2576 , 10.307 , 10.818 , 11.341 , $ 11.778 , 12.108 , 12.569 , 14.000 , 14.467 , $ 15.054 , 15.694 , 16.263 , 16.632 , 20.100 , $ 21.438 / DATA (CSIN(34,J),J=1,41) / $ 1076.5 , 312.66 , 119.74 , 98.571 , 63.193 , $ 49.990 , 39.579 , 35.168 , 32.335 , 30.417 , $ 29.434 , 27.325 , 26.514 , 25.775 , 25.120 , $ 24.711 , 24.818 , 25.600 , 26.106 , 26.355 , $ 25.220 , 24.760 , 21.681 , 20.060 , 19.044 , $ 18.474 , 18.044 , 18.301 , 18.347 , 18.192 , $ 18.557 , 18.639 , 18.870 , 18.769 , 18.478 , $ 18.372 , 18.302 , 18.735 , 18.206 , 20.207 , $ 21.576 / DATA (CSPIEL( 1,J),J=1,41) / $ 0.00000E+00, 350.00 , 580.00 , 600.00 , 550.00 , $ 450.00 , 410.00 , 370.00 , 340.00 , 230.00 , $ 220.00 , 205.00 , 180.00 , 155.00 , 145.00 , $ 140.00 , 160.00 , 195.00 , 235.00 , 250.00 , $ 270.00 , 280.00 , 300.00 , 300.00 , 290.00 , $ 285.00 , 265.00 , 240.00 , 230.00 , 222.00 , $ 204.00 , 196.00 , 190.00 , 170.00 , 170.00 , $ 160.00 , 150.00 , 140.00 , 120.00 , 80.000 , $ 80.000 / DATA (CSPIEL( 2,J),J=1,41) / $ 0.00000E+00, 700.00 , 1000.0 , 1200.0 , 1300.0 , $ 1300.0 , 1250.0 , 1250.0 , 1100.0 , 1000.0 , $ 940.00 , 740.00 , 700.00 , 670.00 , 660.00 , $ 670.00 , 680.00 , 700.00 , 735.00 , 800.00 , $ 810.00 , 820.00 , 820.00 , 810.00 , 800.00 , $ 800.00 , 700.00 , 600.00 , 500.00 , 470.00 , $ 440.00 , 410.00 , 380.00 , 330.00 , 330.00 , $ 330.00 , 330.00 , 330.00 , 285.00 , 240.00 , $ 240.00 / DATA (CSPIEL( 3,J),J=1,41) / $ 0.00000E+00, 1700.0 , 2200.0 , 2200.0 , 1800.0 , $ 1300.0 , 1200.0 , 900.00 , 900.00 , 1000.0 , $ 1100.0 , 1300.0 , 1400.0 , 1420.0 , 1490.0 , $ 1560.0 , 1580.0 , 1690.0 , 1795.0 , 2000.0 , $ 2070.0 , 2140.0 , 2050.0 , 2010.0 , 1970.0 , $ 1880.0 , 1690.0 , 1500.0 , 1420.0 , 1390.0 , $ 1350.0 , 1360.0 , 1370.0 , 1280.0 , 1290.0 , $ 1295.0 , 1250.0 , 1200.0 , 1050.0 , 900.00 , $ 900.00 / DATA (CSPIIN( 1,J),J=1,41) / $ 0.00000E+00, 200.00 , 320.00 , 500.00 , 600.00 , $ 600.00 , 590.00 , 530.00 , 510.00 , 470.00 , $ 430.00 , 425.00 , 420.00 , 425.00 , 425.00 , $ 430.00 , 430.00 , 435.00 , 435.00 , 440.00 , $ 430.00 , 430.00 , 420.00 , 420.00 , 420.00 , $ 415.00 , 415.00 , 410.00 , 410.00 , 408.00 , $ 406.00 , 404.00 , 400.00 , 380.00 , 340.00 , $ 340.00 , 340.00 , 340.00 , 340.00 , 340.00 , $ 340.00 / DATA (CSPIIN( 2,J),J=1,41) / $ 0.00000E+00, 400.00 , 800.00 , 1000.0 , 1100.0 , $ 1200.0 , 1150.0 , 1050.0 , 1000.0 , 900.00 , $ 860.00 , 860.00 , 850.00 , 850.00 , 840.00 , $ 830.00 , 820.00 , 810.00 , 805.00 , 800.00 , $ 800.00 , 800.00 , 800.00 , 800.00 , 800.00 , $ 800.00 , 800.00 , 800.00 , 800.00 , 780.00 , $ 760.00 , 740.00 , 720.00 , 720.00 , 700.00 , $ 690.00 , 680.00 , 670.00 , 665.00 , 660.00 , $ 660.00 / DATA (CSPIIN( 3,J),J=1,41) / $ 0.00000E+00, 1000.0 , 1900.0 , 2600.0 , 2900.0 , $ 3000.0 , 2800.0 , 2600.0 , 2500.0 , 2300.0 , $ 2200.0 , 2000.0 , 1900.0 , 1880.0 , 1860.0 , $ 1840.0 , 1820.0 , 1810.0 , 1805.0 , 1800.0 , $ 1780.0 , 1760.0 , 1750.0 , 1740.0 , 1730.0 , $ 1720.0 , 1710.0 , 1700.0 , 1680.0 , 1660.0 , $ 1650.0 , 1640.0 , 1630.0 , 1620.0 , 1610.0 , $ 1605.0 , 1600.0 , 1600.0 , 1550.0 , 1500.0 , $ 1500.0 / DATA (CSPNEL( 1,J),J=1,41) / $ 2100.0 , 1800.0 , 1500.0 , 1050.0 , 900.00 , $ 950.00 , 800.00 , 650.00 , 570.00 , 390.00 , $ 300.00 , 240.00 , 230.00 , 230.00 , 220.00 , $ 220.00 , 225.00 , 225.00 , 240.00 , 240.00 , $ 290.00 , 330.00 , 335.00 , 350.00 , 355.00 , $ 370.00 , 350.00 , 330.00 , 310.00 , 290.00 , $ 270.00 , 265.00 , 260.00 , 230.00 , 210.00 , $ 210.00 , 200.00 , 200.00 , 190.00 , 180.00 , $ 180.00 / DATA (CSPNEL( 2,J),J=1,41) / $ 3800.0 , 2900.0 , 1850.0 , 1550.0 , 1450.0 , $ 1520.0 , 1460.0 , 1300.0 , 1140.0 , 880.00 , $ 700.00 , 620.00 , 540.00 , 560.00 , 460.00 , $ 460.00 , 470.00 , 470.00 , 480.00 , 480.00 , $ 580.00 , 600.00 , 610.00 , 620.00 , 620.00 , $ 620.00 , 590.00 , 580.00 , 460.00 , 440.00 , $ 420.00 , 400.00 , 480.00 , 430.00 , 380.00 , $ 380.00 , 380.00 , 380.00 , 380.00 , 380.00 , $ 380.00 / DATA (CSPNEL( 3,J),J=1,41) / $ 7000.0 , 6000.0 , 4500.0 , 3350.0 , 2700.0 , $ 3000.0 , 3550.0 , 3970.0 , 3280.0 , 2490.0 , $ 2100.0 , 1510.0 , 1440.0 , 1370.0 , 1370.0 , $ 1370.0 , 1400.0 , 1400.0 , 1420.0 , 1420.0 , $ 1440.0 , 1460.0 , 1460.0 , 1450.0 , 1450.0 , $ 1470.0 , 1400.0 , 1400.0 , 1380.0 , 1370.0 , $ 1360.0 , 1350.0 , 1340.0 , 1330.0 , 1320.0 , $ 1310.0 , 1305.0 , 1300.0 , 1300.0 , 1300.0 , $ 1300.0 / DATA (CSPNIN( 1,J),J=1,41) / $ 0.00000E+00, 200.00 , 400.00 , 800.00 , 800.00 , $ 550.00 , 500.00 , 450.00 , 430.00 , 410.00 , $ 400.00 , 390.00 , 380.00 , 370.00 , 370.00 , $ 370.00 , 365.00 , 365.00 , 360.00 , 360.00 , $ 360.00 , 360.00 , 365.00 , 370.00 , 375.00 , $ 380.00 , 400.00 , 410.00 , 420.00 , 430.00 , $ 440.00 , 440.00 , 440.00 , 440.00 , 440.00 , $ 440.00 , 440.00 , 440.00 , 440.00 , 440.00 , $ 440.00 / DATA (CSPNIN( 2,J),J=1,41) / $ 0.00000E+00, 400.00 , 950.00 , 1050.0 , 1050.0 , $ 980.00 , 940.00 , 900.00 , 860.00 , 820.00 , $ 800.00 , 780.00 , 760.00 , 740.00 , 740.00 , $ 740.00 , 730.00 , 730.00 , 720.00 , 720.00 , $ 720.00 , 720.00 , 730.00 , 740.00 , 750.00 , $ 760.00 , 800.00 , 820.00 , 820.00 , 820.00 , $ 820.00 , 820.00 , 820.00 , 820.00 , 820.00 , $ 820.00 , 820.00 , 820.00 , 820.00 , 820.00 , $ 820.00 / DATA (CSPNIN( 3,J),J=1,41) / $ 0.00000E+00, 0.00000E+00, 500.00 , 1450.0 , 1700.0 , $ 1800.0 , 1750.0 , 1730.0 , 1720.0 , 1710.0 , $ 1700.0 , 1690.0 , 1660.0 , 1630.0 , 1630.0 , $ 1630.0 , 1600.0 , 1600.0 , 1580.0 , 1580.0 , $ 1580.0 , 1580.0 , 1600.0 , 1630.0 , 1650.0 , $ 1670.0 , 1760.0 , 1800.0 , 1800.0 , 1800.0 , $ 1800.0 , 1800.0 , 1800.0 , 1800.0 , 1800.0 , $ 1800.0 , 1800.0 , 1800.0 , 1800.0 , 1800.0 , $ 1800.0 / DATA ELAB / $ 0.10000E-03, 0.20000E-03, 0.30000E-03, 0.40000E-03, 0.50000E-03, $ 0.70000E-03, 0.10000E-02, 0.20000E-02, 0.30000E-02, 0.40000E-02, $ 0.50000E-02, 0.70000E-02, 0.10000E-01, 0.15000E-01, 0.20000E-01, $ 0.25000E-01, 0.32700E-01/ DATA CNLWAT / $ 1.0000 , 16.000 , 27.000 , 56.000 , 59.000 , $ 64.000 , 91.000 , 112.00 , 119.00 , 127.00 , $ 137.00 , 181.00 , 207.00 , 209.00 , 238.00 / DATA (CNLWEL( 1,J),J=1,17) / $ 6000.0 , 5500.0 , 5200.0 , 4900.0 , 4800.0 , $ 4400.0 , 4000.0 , 2900.0 , 2200.0 , 1800.0 , $ 1400.0 , 1100.0 , 900.00 , 700.00 , 600.00 , $ 560.00 , 520.00 / DATA (CNLWEL( 2,J),J=1,17) / $ 5400.0 , 5050.0 , 4800.0 , 4600.0 , 4399.0 , $ 4090.0 , 3700.0 , 2600.0 , 1950.0 , 1600.0 , $ 1300.0 , 900.00 , 700.00 , 800.00 , 1050.0 , $ 1250.0 , 1320.0 / DATA (CNLWEL( 3,J),J=1,17) / $ 5500.0 , 5150.0 , 4900.0 , 4699.0 , 4490.0 , $ 4150.0 , 3750.0 , 2790.0 , 2100.0 , 1650.0 , $ 1300.0 , 950.00 , 800.00 , 860.00 , 1000.0 , $ 1090.0 , 1080.0 / DATA (CNLWEL( 4,J),J=1,17) / $ 5499.0 , 4970.0 , 4450.0 , 4080.0 , 3750.0 , $ 3380.0 , 2900.0 , 2400.0 , 2380.0 , 2350.0 , $ 2300.0 , 2100.0 , 1720.0 , 1370.0 , 1200.0 , $ 1060.0 , 870.00 / DATA (CNLWEL( 5,J),J=1,17) / $ 5399.0 , 4710.0 , 4180.0 , 3760.0 , 3460.0 , $ 3150.0 , 2730.0 , 2270.0 , 1850.0 , 1850.0 , $ 2130.0 , 2330.0 , 2120.0 , 1640.0 , 1310.0 , $ 1100.0 , 1050.0 / DATA (CNLWEL( 6,J),J=1,17) / $ 5099.0 , 4405.0 , 3825.0 , 3455.0 , 3125.0 , $ 2695.0 , 2350.0 , 1850.0 , 1580.0 , 1820.0 , $ 2050.0 , 2210.0 , 2000.0 , 1590.0 , 1310.0 , $ 1120.0 , 1040.0 / DATA (CNLWEL( 7,J),J=1,17) / $ 6290.0 , 5960.0 , 5640.0 , 5370.0 , 5150.0 , $ 4800.0 , 4250.0 , 3150.0 , 2470.0 , 2100.0 , $ 2230.0 , 2420.0 , 2450.0 , 2050.0 , 1760.0 , $ 1550.0 , 1330.0 / DATA (CNLWEL( 8,J),J=1,17) / $ 6885.0 , 6650.0 , 6350.0 , 6150.0 , 6000.0 , $ 5700.0 , 5360.0 , 4250.0 , 2800.0 , 1870.0 , $ 1810.0 , 1820.0 , 2170.0 , 2450.0 , 2150.0 , $ 1700.0 , 1390.0 / DATA (CNLWEL( 9,J),J=1,17) / $ 6600.0 , 6500.0 , 6400.0 , 6249.0 , 6190.0 , $ 5950.0 , 5520.0 , 4250.0 , 2750.0 , 1900.0 , $ 1850.0 , 1950.0 , 2340.0 , 2800.0 , 2540.0 , $ 2100.0 , 1760.0 / DATA (CNLWEL(10,J),J=1,17) / $ 7400.0 , 7200.0 , 6999.0 , 6840.0 , 6655.0 , $ 6320.0 , 5820.0 , 4400.0 , 2850.0 , 2000.0 , $ 1800.0 , 1800.0 , 2150.0 , 2600.0 , 2350.0 , $ 1950.0 , 2100.0 / DATA (CNLWEL(11,J),J=1,17) / $ 7900.0 , 7700.0 , 7499.0 , 7390.0 , 7202.0 , $ 6810.0 , 6360.0 , 4920.0 , 3450.0 , 2600.0 , $ 2200.0 , 1950.0 , 2300.0 , 2800.0 , 2650.0 , $ 2250.0 , 2050.0 / DATA (CNLWEL(12,J),J=1,17) / $ 7900.0 , 7750.0 , 7699.0 , 7590.0 , 7450.0 , $ 7200.0 , 6850.0 , 5650.0 , 4400.0 , 3700.0 , $ 3400.0 , 2800.0 , 2700.0 , 3100.0 , 3250.0 , $ 3100.0 , 2750.0 / DATA (CNLWEL(13,J),J=1,17) / $ 6100.0 , 5950.0 , 5750.0 , 5599.0 , 5440.0 , $ 5200.0 , 4800.0 , 4300.0 , 5800.0 , 5750.0 , $ 4800.0 , 3420.0 , 2650.0 , 3200.0 , 3650.0 , $ 3500.0 , 2980.0 / DATA (CNLWEL(14,J),J=1,17) / $ 6100.0 , 5950.0 , 5750.0 , 5599.0 , 5440.0 , $ 5200.0 , 4800.0 , 4300.0 , 5800.0 , 5750.0 , $ 4800.0 , 3420.0 , 2650.0 , 3200.0 , 3650.0 , $ 3500.0 , 2980.0 / DATA (CNLWEL(15,J),J=1,17) / $ 6600.0 , 6350.0 , 6100.0 , 5899.0 , 5690.0 , $ 5300.0 , 4850.0 , 4450.0 , 5650.0 , 5700.0 , $ 4950.0 , 3850.0 , 3050.0 , 3050.0 , 3460.0 , $ 3650.0 , 3340.0 / DATA (CNLWIN( 1,J),J=1,17) / 17*0.0E+00 / DATA (CNLWIN( 2,J),J=1,17) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , $ 10.000 , 50.000 , 100.00 , 200.00 , 300.00 , $ 400.00 , 600.00 , 700.00 , 750.00 , 700.00 , $ 700.00 , 680.00 / DATA (CNLWIN( 3,J),J=1,17) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 , $ 50.000 , 100.00 , 260.00 , 450.00 , 600.00 , $ 700.00 , 800.00 , 900.00 , 940.00 , 900.00 , $ 860.00 , 820.00 / DATA (CNLWIN( 4,J),J=1,17) / $ 1.0000 , 80.000 , 200.00 , 320.00 , 400.00 , $ 520.00 , 700.00 , 1000.0 , 1120.0 , 1200.0 , $ 1200.0 , 1200.0 , 1180.0 , 1130.0 , 1100.0 , $ 1090.0 , 1080.0 / DATA (CNLWIN( 5,J),J=1,17) / $ 1.0000 , 90.000 , 220.00 , 340.00 , 420.00 , $ 550.00 , 720.00 , 1080.0 , 1300.0 , 1400.0 , $ 1420.0 , 1420.0 , 1380.0 , 1260.0 , 1190.0 , $ 1150.0 , 1100.0 / DATA (CNLWIN( 6,J),J=1,17) / $ 1.0000 , 95.000 , 225.00 , 345.00 , 425.00 , $ 555.00 , 750.00 , 1150.0 , 1500.0 , 1680.0 , $ 1700.0 , 1690.0 , 1550.0 , 1360.0 , 1240.0 , $ 1180.0 , 1120.0 / DATA (CNLWIN( 7,J),J=1,17) / $ 10.000 , 140.00 , 260.00 , 380.00 , 450.00 , $ 600.00 , 750.00 , 1200.0 , 1580.0 , 1800.0 , $ 1820.0 , 1830.0 , 1800.0 , 1750.0 , 1690.0 , $ 1650.0 , 1620.0 / DATA (CNLWIN( 8,J),J=1,17) / $ 15.000 , 150.00 , 300.00 , 400.00 , 500.00 , $ 650.00 , 840.00 , 1500.0 , 2100.0 , 2130.0 , $ 2140.0 , 2130.0 , 2080.0 , 2000.0 , 1950.0 , $ 1900.0 , 1860.0 / DATA (CNLWIN( 9,J),J=1,17) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 , $ 150.00 , 380.00 , 1000.0 , 1650.0 , 2100.0 , $ 2100.0 , 2100.0 , 2060.0 , 1950.0 , 1860.0 , $ 1800.0 , 1740.0 / DATA (CNLWIN(10,J),J=1,17) / $ 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 , 45.000 , $ 180.00 , 380.00 , 1050.0 , 1900.0 , 2300.0 , $ 2300.0 , 2200.0 , 2150.0 , 2000.0 , 1900.0 , $ 1800.0 , 1750.0 / DATA (CNLWIN(11,J),J=1,17) / $ 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 , 48.000 , $ 190.00 , 390.00 , 1080.0 , 2000.0 , 2400.0 , $ 2400.0 , 2300.0 , 2200.0 , 2100.0 , 1950.0 , $ 1850.0 , 1800.0 / DATA (CNLWIN(12,J),J=1,17) / $ 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 , 50.000 , $ 200.00 , 400.00 , 1100.0 , 2100.0 , 2500.0 , $ 2500.0 , 2450.0 , 2300.0 , 2100.0 , 2000.0 , $ 1900.0 , 1850.0 / DATA (CNLWIN(13,J),J=1,17) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 , $ 100.00 , 350.00 , 900.00 , 1400.0 , 2000.0 , $ 2300.0 , 2380.0 , 2400.0 , 2300.0 , 2250.0 , $ 2200.0 , 2120.0 / DATA (CNLWIN(14,J),J=1,17) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 , $ 100.00 , 350.00 , 900.00 , 1400.0 , 2000.0 , $ 2300.0 , 2380.0 , 2400.0 , 2300.0 , 2250.0 , $ 2200.0 , 2120.0 / DATA (CNLWIN(15,J),J=1,17) / $ 0.00000E+00, 0.00000E+00, 0.00000E+00, 1.0000 , 10.000 , $ 100.00 , 400.00 , 950.00 , 1600.0 , 2200.0 , $ 2550.0 , 2750.0 , 2700.0 , 2600.0 , 2540.0 , $ 2450.0 , 2360.0 / DATA (CSCAP(J),J=1,50) / $ 6.0000 , 5.7000 , 5.5000 , 5.3000 , 5.2000 , $ 5.1000 , 5.0000 , 4.9000 , 4.8000 , 4.8000 , $ 4.8000 , 4.8000 , 4.8000 , 4.8000 , 4.8000 , $ 4.8000 , 4.9000 , 5.0000 , 5.2000 , 5.5000 , $ 6.0000 , 6.7000 , 7.5000 , 8.5000 , 10.000 , $ 12.000 , 14.500 , 19.000 , 26.500 , 40.000 , $ 75.000 , 120.00 , 180.00 , 260.00 , 360.00 , $ 330.00 , 60.000 , 7.0000 , 9.5000 , 20.000 , $ 75.000 , 140.00 , 250.00 , 360.00 , 480.00 , $ 580.00 , 590.00 , 500.00 , 300.00 , 100.00 / DATA (CSCAP(J),J=51,100) / $ 200.00 , 300.00 , 400.00 , 470.00 , 500.00 , $ 430.00 , 100.00 , 20.000 , 22.000 , 40.000 , $ 560.00 , 950.00 , 1000.0 , 1000.0 , 1000.0 , $ 990.00 , 920.00 , 860.00 , 790.00 , 740.00 , $ 650.00 , 600.00 , 540.00 , 470.00 , 440.00 , $ 390.00 , 360.00 , 340.00 , 320.00 , 310.00 , $ 280.00 , 2.0000 , 2.5000 , 6.0000 , 13.000 , $ 38.000 , 65.000 , 140.00 , 280.00 , 300.00 , $ 430.00 , 580.00 , 650.00 , 800.00 , 920.00 , $ 1100.0 , 1250.0 , 1400.0 , 1550.0 , 1700.0 / C --- END OF CROSS SECTION DATA STATEMENTS --- C --- DATA STMTS. FOR GEANT/GHEISHA PARTICLE CODE CONVERSIONS --- C --- KIPART(I)=GHEISHA CODE CORRESPONDING TO GEANT CODE I --- C --- IKPART(I)=GEANT CODE CORRESPONDING TO GHEISHA CODE I --- DATA KIPART/ $ 1, 3, 4, 2, 5, 6, 8, 7, $ 9, 12, 10, 13, 16, 14, 15, 11, $ 35, 18, 20, 21, 22, 26, 27, 33, $ 17, 19, 23, 24, 25, 28, 29, 34, $ 35, 35, 35, 35, 35, 35, 35, 35, $ 35, 35, 35, 35, 30, 31, 32, 35/ DATA IKPART/ $ 1, 4, 2, 3, 5, 6, 8, 7, $ 9, 11, 16, 10, 12, 14, 15, 13, $ 25, 18, 26, 19, 20, 21, 27, 28, $ 29, 22, 23, 30, 31, 45, 46, 47, $ 24, 32, 48/ C PARAMETER (ONETHR=1./3.) DATA ONETHR / .33333333/ DATA ALPHA / 6*0.7, + 0.75 ,0.75 ,0.75 , + 0.76,0.76 ,0.76 ,0.76 , + 0.685,0.63 ,0.685,0.63,0.685,0.63, + 3*0.685,3*0.63,2*0.685,2*0.63, + 3*0.7,0.685,0.63,0.7/ DATA ALPHAC /1.2,1.2,1.2,1.15,0.90,0.91,0.98,1.06,1.10,1.11, + 1.10,1.08,1.05,1.01,0.985,0.962,0.945,0.932, + 0.925,0.920,0.920,0.921,0.922,0.923,0.928,0.931, + 0.940,0.945,0.950,0.955,0.958,0.962,0.965,0.976, + 0.982,0.988,0.992,1.010,1.020,1.030,1.040/ DATA PARTEL/6*0.,29*1./ DATA PARTIN/6*0.,1.00,0.00,1.05,1.20,1.35,1.30,1.20,1.00,1.30, + 1.00,1.30,1.00,1.30,1.00,1.00,1.00,1.30,1.30,1.30, + 1.00,1.00,1.30,1.30,1.00,1.,1.,1.,1.3,1./ DATA ICORR /14*1, 0, 1, 0, 1, 0, 3*1, 3*0, 2*1, 2*0, 4*1, 2*0/ C-- SET INTRC TO 0 FOR IPART = 26-29, 33, 34 ( XI'S AND OMEGA'S ) C-DH- DATA INTRC /6*0, 1, 0, 12*1, 0, 2*2, 0, 1, 4*0, 3*1, 3*0 / C-- RESET INTRC FOR IPART = 26-29, 33, 34 ( XI'S AND OMEGA'S ) DATA INTRC /6*0, 1, 0, 12*1, 0, 2*2, 0, 10*1, 0/ C CROSS SECTIONS ON NUCLEUS ARE KNOWN ONLY FOR PIONS AND PROTONS. C THE GENERAL LAW SIGMA(A)=1.25*SIGMA(TOT,PROTON)*A**ALPHA IS VALID C ONLY FOR MOMENTA > 2 GEV.THE PARAMETRIZATION DONE HERE GIVES ONLY C A BEHAVIOUR AVERAGED OVER MOMENTA AND PARTICLE TYPES. C FOR A DETECTOR WITH ONLY A FEW MATERIALS IT'S OF COURSE MUCHBETTER C TO USE TABLES OF THE MEASURED CROSS SECTIONS . C FOR ELEMENTS WITH THE FOLLOWING ATOMIC NUMBERS MEASURED CROSS C SECTIONS ARE AVAILABLE (SEE "PCSDATA"). C H AL CU PB DATA CSA /1. ,27.00 ,63.54 ,207.19 / DATA IPART2/9,8,7,11,10,13,12/ SAVE ALPHA,ALPHAC,PARTEL,PARTIN,CSA,IPART2,ICORR,INTRC C----------------------------------------------------------------------- IF ( DEBUG ) THEN WRITE(MDEBUG,*) 'CGHSIG:' NPRT(4)=.TRUE. NPRT(9)=.TRUE. ENDIF C --- INITIALISE CGHSIG AND SWITCH TO GHEISHA PARTICLE CODE --- CGHSIG=0.0 IF ( LPART .LE. 48 ) THEN IPART = KIPART(LPART) ELSEIF ( LPART .EQ. 201 ) THEN IPART = 30 ELSEIF ( LPART .EQ. 301 ) THEN IPART = 31 ELSEIF ( LPART .EQ. 402 ) THEN IPART = 32 ELSE GOTO 160 ENDIF C --- NO INTERACTION FOR GAMMAS, NEUTRINOS, ELECTRONS, POSITRONS, MUONS, C --- NEUTRAL PIONS, NEUTRAL SIGMAS AND ANTISIGMAS AND NEW PARTICLES. IF(INTRC(IPART).EQ. 0) GO TO 160 P=PPART EK=EKIN C --- INITIALISE THE CROSS-SECTIONS WITH 0.0 --- DO 10 K=1,KK AIEL(K)=0.0 AIIN(K)=0.0 AICA(K)=0.0 10 CONTINUE C IF ((IPART .GE. 30) .AND. (IPART .LE. 32)) THEN C --- TAKE GEOMETRICAL CROSS SECTIONS FOR INELASTIC SCATTERING --- C --- OF DEUTERONS, TRITONS AND ALPHAS --- IF ( IPART .EQ. 30 ) THEN APART=2.0**ONETHR ELSEIF ( IPART .EQ. 31 ) THEN APART=3.0**ONETHR ELSEIF ( IPART .EQ. 32 ) THEN APART=4.0**ONETHR ENDIF DO 20 K=1,KK AIIN(K)=49.0*(APART+ACOMP(K)**ONETHR)**2 20 CONTINUE IF (NPRT(9)) WRITE(MDEBUG,10000) ELSEIF ((IPART .EQ. 16) .AND. (EK .LE. 0.0327)) THEN C --- USE TABLES FOR LOW ENERGY NEUTRONS --- C --- GET ENERGY BIN --- JE2=17 DO 30 J=2,17 IF (EK .LT. ELAB(J)) THEN JE2=J GO TO 40 ENDIF 30 CONTINUE 40 JE1=JE2-1 EKX=MAX(EK,1.0E-9) DELAB=ELAB(JE2)-ELAB(JE1) DO 70 K=1,KK C --- GET A BIN --- JA2=15 DO 50 J=2,15 IF (ACOMP(K) .LT. CNLWAT(J)) THEN JA2=J GO TO 60 ENDIF 50 CONTINUE 60 JA1=JA2-1 DNLWAT=CNLWAT(JA2)-CNLWAT(JA1) C --- USE LINEAR INTERPOLATION OR EXTRAPOLATION BY Y=RCE*X+RCA*X+B --- C --- ELASTIC CROSS SECTION --- C --- E INTERPOLATION OR EXTRAPOLATION AT JA1 --- DY=CNLWEL(JA1,JE2)-CNLWEL(JA1,JE1) RCE=DY/DELAB C --- A INTERPOLATION OR EXTRAPOLATION AT JE1 --- DY=CNLWEL(JA2,JE1)-CNLWEL(JA1,JE1) RCA=DY/DNLWAT B=CNLWEL(JA1,JE1)-RCE*ELAB(JE1)-RCA*CNLWAT(JA1) AIEL(K)=RCE*EK+RCA*ACOMP(K)+B C --- INELASTIC CROSS SECTION --- C --- E INTERPOLATION OR EXTRAPOLATION AT JA1 --- DY=CNLWIN(JA1,JE2)-CNLWIN(JA1,JE1) RCE=DY/DELAB C --- A INTERPOLATION OR EXTRAPOLATION AT JE1 --- DY=CNLWIN(JA2,JE1)-CNLWIN(JA1,JE1) RCA=DY/DNLWAT B=CNLWIN(JA1,JE1)-RCE*ELAB(JE1)-RCA*CNLWAT(JA1) AIIN(K)=RCE*EK+RCA*ACOMP(K)+B IZNO=ZCOMP(K)+0.01 AICA(K)=11.12*CSCAP(IZNO)/(EKX*1.0E6)**0.577 70 CONTINUE IF (NPRT(9)) WRITE(MDEBUG,10100) ELSE C --- USE PARAMETRIZATION OF CROSS SECTION DATA FOR ALL OTHER CASES --- IF (NPRT(9)) WRITE(MDEBUG,10200) C --- GET MOMENTUM BIN --- J=40 DO 80 I=2,41 IF (P .LT. PLAB(I)) THEN J=I-1 GO TO 90 ENDIF 80 CONTINUE C --- START WITH CROSS SECTIONS FOR SCATTERING ON FREE PROTONS --- C --- USE LINEAR INTERPOLATION OR EXTRAPOLATION BY Y=RC*X+B --- 90 DX=PLAB(J+1)-PLAB(J) C --- ELASTIC CROSS SECTION --- DY=CSEL(IPART,J+1)-CSEL(IPART,J) RC=DY/DX B=CSEL(IPART,J)-RC*PLAB(J) AIELIN=RC*P+B C --- INELASTIC CROSS SECTION --- DY=CSIN(IPART,J+1)-CSIN(IPART,J) RC=DY/DX B=CSIN(IPART,J)-RC*PLAB(J) AIININ=RC*P+B ALPH=ALPHA(IPART) IF(IPART.LT.14) THEN DY=ALPHAC(J+1)-ALPHAC(J) RC=DY/DX B=ALPHAC(J)-RC*PLAB(J) CORFAC=RC*P+B ALPH=ALPH*CORFAC IPART3=IPART2(IPART-6) C --- ELASTIC CROSS SECTION --- DY=CSEL(IPART3,J+1)-CSEL(IPART3,J) RC=DY/DX B=CSEL(IPART3,J)-RC*PLAB(J) XSECEL=RC*P+B C --- INELASTIC CROSS SECTION --- DY=CSIN(IPART3,J+1)-CSIN(IPART3,J) RC=DY/DX B=CSIN(IPART3,J)-RC*PLAB(J) XSECIN=RC*P+B ENDIF DO 100 K=1,KK AIEL(K)=AIELIN AIIN(K)=AIININ IF (ACOMP(K) .GE. 1.5) THEN C --- A-DEPENDENCE FROM PARAMETRIZATION --- CREL=1.0 CRIN=1.0 C --- GET MEDIUM BIN 1=HYDR. 2=AL 3=CU 4=PB --- I=3 IF (ACOMP(K) .LT. 50.0) I=2 IF (ACOMP(K) .GT. 100.0) I=4 IF ((IPART .EQ. 14) .OR. (IPART .EQ. 16)) THEN C --- PROTONS AND NEUTRONS --- C --- ELASTIC CROSS SECTION --- DY=CSPNEL(I-1,J+1)-CSPNEL(I-1,J) RC=DY/DX B=CSPNEL(I-1,J)-RC*PLAB(J) XSECEL=RC*P+B C --- INELASTIC CROSS SECTION --- DY=CSPNIN(I-1,J+1)-CSPNIN(I-1,J) RC=DY/DX B=CSPNIN(I-1,J)-RC*PLAB(J) XSECIN=RC*P+B IF (AIEL(K) .GE. 0.001) CREL=XSECEL/(0.36*AIEL(K)* + CSA(I)**1.17) AITOT=AIEL(K)+AIIN(K) IF (AITOT .GE. 0.001) CRIN=XSECIN/(AITOT*CSA(I)** + ALPH) ELSEIF (IPART .LT. 15) THEN C --- CALCULATE CORRECTION FACTORS FROM VALUES ON AL,CU,PB FOR ALL --- C --- MESONS USE LINEAR INTERPOLATION OR EXTRAPOLATION BY Y=RC*X+B --- C --- NOTE THAT DATA IS ONLY AVAILABLE FOR PIONS AND PROTONS WGCH=0.5 IF (ACOMP(K) .LT. 20.0) WGCH=0.5+0.5*EXP(-(ACOMP(K)-1.0)) AIEL(K)=WGCH*AIEL(K)+(1.0-WGCH)*XSECEL AIIN(K)=WGCH*AIIN(K)+(1.0-WGCH)*XSECIN C --- THIS SECTION NOT FOR KAONS --- IF (IPART .LT. 10) THEN C --- ELASTIC CROSS SECTION --- DY=CSPIEL(I-1,J+1)-CSPIEL(I-1,J) RC=DY/DX B=CSPIEL(I-1,J)-RC*PLAB(J) XSPIEL=RC*P+B C --- INELASTIC CROSS SECTION --- DY=CSPIIN(I-1,J+1)-CSPIIN(I-1,J) RC=DY/DX B=CSPIIN(I-1,J)-RC*PLAB(J) XSPIIN=RC*P+B IF (AIEL(K) .GE. 0.001) CREL=XSPIEL/(0.36* AIEL(K) + *CSA(I)**1.17) AITOT=AIEL(K)+AIIN(K) IF (AITOT .GE. 0.001) CRIN=XSPIIN/(AITOT*CSA(I) + **ALPH) ENDIF ENDIF AIIN(K)=CRIN*(AIIN(K)+AIEL(K))*ACOMP(K)**ALPH AIEL(K)=CREL*0.36*AIEL(K)*ACOMP(K)**1.17 AIEL(K)=AIEL(K)*PARTEL(IPART) AIIN(K)=AIIN(K)*PARTIN(IPART) ENDIF 100 CONTINUE ENDIF C --- CALCULATE INTERACTION PROBABILITY --- ALAM=0.0 DO 150 K=1,KK AIEL(K) = AIEL(K)*WCOMP(K) AIIN(K) = AIIN(K)*WCOMP(K) AICA(K) = AICA(K)*WCOMP(K) ALAM = ALAM + AIEL(K) + AIIN(K) + AICA(K) 150 CONTINUE C --- PASS THE CROSS SECTION (MBARN) TO CORSIKA --- CGHSIG=ALAM GO TO 999 C --- PRINTOUT OF SKIPPED PARTICLES IN CASE OF INTERFACE DEBUG --- 160 IF (NPRT(9)) WRITE(MDEBUG,10300) IPART 10000 FORMAT(' *CGHSIG* GEOM X-SECT. FOR INEL. SCAT. OF D,T AND ALPHA') 10100 FORMAT(' *CGHSIG* X-SECT. FROM LOW ENERGY NEUTRON TABLES') 10200 FORMAT(' *CGHSIG* X-SECT. FROM PARAMETRIZATION OF DATA') 10300 FORMAT(' *CGHSIG* GHEISHA PARTICLE ',I3,' SKIPPED') 999 RETURN END