source: trunk/MagicSoft/Simulation/Corsika/Mmcs/stankg.f

Last change on this file was 286, checked in by harald, 25 years ago
This is the start point for further developments of the Magic Monte Carlo Simulation written by Jose Carlos Gonzales. Now it is under control of one CVS repository for the whole collaboration. Everyone should use this CVS repository for further developments.
File size: 2.0 KB
Line 
1 SUBROUTINE STANKG
2
3C-----------------------------------------------------------------------
4C STA(RT) NKG
5C
6C INITIALIZE ARRAYS FOR SINGLE SHOWERS NKG CALCULATED VARIABLES
7C THIS SUBROUTINE IS CALLED FROM MAIN
8C-----------------------------------------------------------------------
9
10 IMPLICIT NONE
11*KEEP,NKGS.
12 COMMON /NKGS/ CZX,CZY,CZXY,CZYX,SAH,SL,ZNE
13 DOUBLE PRECISION CZX(-10:10,2),CZY(-10:10,2),CZXY(-10:10,2),
14 * CZYX(-10:10,2),SAH(10),SL(10),ZNE(10)
15*KEEP,RUNPAR.
16 COMMON /RUNPAR/ FIXHEI,THICK0,HILOECM,HILOELB,
17 * STEPFC,NRRUN,NSHOW,PATAPE,MONIIN,
18 * MONIOU,MDEBUG,NUCNUC,
19 * CETAPE,
20 * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
21 * N1STTR,MDBASE,
22 * DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
23 * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
24 * ,GHEISH,GHESIG
25 COMMON /RUNPAC/ DSN,HOST,USER
26 DOUBLE PRECISION FIXHEI,THICK0,HILOECM,HILOELB
27 REAL STEPFC
28 INTEGER NRRUN,NSHOW,PATAPE,MONIIN,MONIOU,MDEBUG,NUCNUC,
29 * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
30 * N1STTR,MDBASE
31 INTEGER CETAPE
32 CHARACTER*79 DSN
33 CHARACTER*20 HOST,USER
34
35 LOGICAL DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
36 * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
37 * ,GHEISH,GHESIG
38*KEND.
39
40 INTEGER I,K
41C-----------------------------------------------------------------------
42
43 IF ( DEBUG ) WRITE(MDEBUG,*) 'STANKG:'
44
45C CLEAR ARRAYS FOR AGE PARAMETER CALCULATION FOR EACH SHOWER
46 DO 17 I = 1,10
47 SAH(I) = 0.D0
48 SL (I) = 0.D0
49 ZNE(I) = 0.D0
50 17 CONTINUE
51C CLEAR LATERAL ELECTRON DISTRIBUTION COUNTERS FOR EACH SHOWER
52 DO 45 K = 1,2
53 DO 45 I = -10,10
54 CZX (I,K) = 0.D0
55 CZY (I,K) = 0.D0
56 CZXY(I,K) = 0.D0
57 CZYX(I,K) = 0.D0
58 45 CONTINUE
59
60 RETURN
61 END
Note: See TracBrowser for help on using the repository browser.