Changeset 687 for trunk/MagicSoft/Simulation/Corsika/Mmcs/main.f

Timestamp:

03/05/01 16:20:45 (24 years ago)

Author:

harald

Message:

Changes of Ciro and Dennis to write all events of one run in only
one file. Be careful, from now on you need also a new reflector.

File:

• trunk/MagicSoft/Simulation/Corsika/Mmcs/main.f (modified) (16 diffs)

trunk/MagicSoft/Simulation/Corsika/Mmcs/main.f

@@ -425 +425 @@
       ILEFTA = 0
 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
- 
-c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+C
+CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+C
+C     Modified by C. Bigongiari 2001 Jan 16 
+C 
+C
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
       print *,'JCIO::========================================'
       print *,'JCIO:: Initializing JCIO system for advanced'
       print *,'JCIO:: saving of data.'
       print *,'JCIO::========================================'
-c- initialize jcio system
+C
+Cc- initialize jcio system
+C      
       call jcinitio(dsn,nrrun)
-c- create file run######
-      call jcstartrun(runh)
-c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-
+Cc- create file run######
+C      call jcstartrun(runh)
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+C
+C- Modified JCSTARTRUN creates cer###### and dat###### files !
+C
+C     ###### is the RUN number !
+C
+
+      call jcstartrun(RUNH)
+
+C
+C- write Run Header on cer and dat files  
+C 
+      CALL TOBUF(RUNH,0)
+      IF ( LCERFI ) CALL TOBUFC(RUNH,0)
+
+C
+CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+C
 C-----------------------------------------------------------------------
 C  LOOP OVER SHOWERS
+
       DO 2  ISHW = 1,NSHOW
+
  
 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
@@ -482 +507 @@
           FPRINT = .FALSE.
         ENDIF
- 
-c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-c   Create cer######,dat######,sta###### files
-c------------------------------------------------------------
-        call jcnewshower
-c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-
+CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+Cc   Create cer######,dat######,sta###### files
+Cc------------------------------------------------------------
+C        call jcnewshower
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+ 
 C  RESET COUNTERS
         DO 447  K = 1,25
@@ -823 +849 @@
         ENDIF
         EVTH(147) = 0.
- 
+
         IF ( FEGS ) THEN
 C  PARAMETER FOR ELECTRON AND PHOTON REJECT (CONVERT ENERGY TO MEV)
@@ -838 +864 @@
         EVTH( 9) = PTOT0 * ST * SIN(PHI)
         EVTH(10) = PTOT0 * COSTHE
+     
+
 c
 c [*] one block from here sent above
 c
+
         EVTH(85) = CERSIZ
  
@@ -856 +885 @@
         IF ( EVTH(3) .GT. 3.0  .OR.  .NOT. FEGS ) THEN
           CALL TOBUF ( EVTH,0 )
-          IF ( LCERFI ) CALL TOBUFC( EVTH,0 )
+          IF ( LCERFI ) CALL TOBUFC( EVTH,0 ) 
         ENDIF
  
@@ -865 +894 @@
      *          '    GEN      LEVEL  E ON STACK'/)
         NOPART = 0
- 
- 
+
         IF ( CURPAR(1) .LE. 3.D0  .OR.
      *      (CURPAR(1) .GE. 5.D0  .AND.  CURPAR(1) .LE. 7.D0) ) THEN
@@ -889 +917 @@
 C  HADRONIC PARTICLES
           FNPRIM = .TRUE.
- 
+
 C  FILL LONGITUDINAL DISTRIBUTION FOR THE PRIMARY PARTICLE
 C  THE PARTICLE IS TRACKED FROM THICK0 DOWN TO THICK0+CHI*PRMPAR(3)
@@ -917 +945 @@
             ENDIF
           ENDIF
- 
+
 C  CHECK OBSERVATION LEVEL PASSAGE AND UPDATE PARTICLE COORDINATES
           HNEW = H
@@ -944 +972 @@
         ENDIF
  
+
 C-----------------------------------------------------------------------
 C  NORMAL CYCLE
@@ -1027 +1056 @@
 C  FINISH SHOWER AND PRINT INFORMATION
         CALL OUTEND
- 
  
 *       IF ( DEBUG ) WRITE(MDEBUG,442) NPARTO
@@ -1221 +1249 @@
         ENDIF
  
-c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+C
+CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+C
+C     Modified by C. Bigongiari 2001 Jan 16 
+C 
+C
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
 c   Saves statistics to sta###### file
-        call jcstadata(EVTH,EVTE,
-     +      NPROTO,NPROTB,NNEUTR,NNEUTB,NPHOTO,NELECT,NPOSIT,
-     +      NNU   ,NMUM  ,NMUP  ,NPI0  ,NPIM  ,NPIP  ,NK0L  ,
-     +      NK0S  ,NKMI  ,NKPL  ,NHYP  ,NDEUT ,NTRIT ,NALPHA,
-     +      NOTHER,IFINNU,IFINPI,IFINET,IFINKA,IFINHY,
-     +      CERELE,CERHAD,PLONG,LPCT1,NSTEP,THSTEP)
-c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-
-c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+C        call jcstadata(EVTH,EVTE,
+C     +      NPROTO,NPROTB,NNEUTR,NNEUTB,NPHOTO,NELECT,NPOSIT,
+C     +      NNU   ,NMUM  ,NMUP  ,NPI0  ,NPIM  ,NPIP  ,NK0L  ,
+C     +      NK0S  ,NKMI  ,NKPL  ,NHYP  ,NDEUT ,NTRIT ,NALPHA,
+C     +      NOTHER,IFINNU,IFINPI,IFINET,IFINKA,IFINHY,
+C     +      CERELE,CERHAD,PLONG,LPCT1,NSTEP,THSTEP)
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
 C  WRITE SHOWER END TO OUTPUT BUFFER
-c        CALL TOBUF( EVTE,0 )
-        CALL TOBUF( EVTE,1 )
-c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+Cc        CALL TOBUF( EVTE,0 )
+C        CALL TOBUF( EVTE,1 )
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+C        IF ( LCERFI ) THEN
+C          CALL OUTND2
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+Cc          CALL TOBUFC( EVTE,0 )
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+C        ENDIF
+
+C
+C  WRITE SHOWER END TO OUTPUT BUFFER
+        CALL TOBUF( EVTE,0 )
         IF ( LCERFI ) THEN
           CALL OUTND2
-c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-c          CALL TOBUFC( EVTE,0 )
-c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-        ENDIF
+          CALL TOBUFC( EVTE,0 )
+        ENDIF
+
+C
+CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+C
  
         IF ( FPRINT .OR. DEBUG ) WRITE(MONIOU,*)
@@ -1259 +1304 @@
           JHBIN(J) = JHBIN(J) + IHBIN(J)
   19    CONTINUE
- 
+
+
   2   CONTINUE
+
 C  END OF SHOWER LOOP
  
@@ -1271 +1318 @@
       RUNE(3) = REAL(ISHW)
 C  WRITE RUN END TO OUTPUT BUFFER AND FINISH OUTPUT
-c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-c      CALL TOBUF ( RUNE,1 )
-      call jcendrun(rune)
-c      IF ( LCERFI ) CALL TOBUFC( RUNE,1 )
+C
+CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+C
+C     Modified by C. Bigongiari 2001 Jan 16 
+C
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+Cc      CALL TOBUF ( RUNE,1 )
+C      call jcendrun(rune)
+Cc      IF ( LCERFI ) CALL TOBUFC( RUNE,1 )
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+C
+C- write Run End      
+C
+
+      CALL TOBUF ( RUNE,1 )
+      IF ( LCERFI ) CALL TOBUFC( RUNE,1)
+
+C
+CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+C   
 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
 C  TIME SINCE BEGINNING
@@ -1308 +1371 @@
  209  FORMAT(//' ENERGY - MULTIPLICITY MATRIX FOR ALL SHOWERS'/
      *       ' ENERGY RUNS VERTICALLY, MULTIPLICITY HORIZONTALLY'//,
-     *       ' ',6X,5I9,3I8,5I7,'   ENERGY RANGE (GEV)'/
+     *       ' ',6X,5I9,3I8,5I7'   ENERGY RANGE (GEV)'/
      *       37(/' ',I4,1X,I10,4I9,3I8,5I7,1X,1P,2E10.1,0P)//
      *       ' MULT. ',5I9,3I8,5I7,4X,'LOWER BIN LIMIT'/
@@ -1580 +1643 @@
       ENDIF
  
-c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-      call jcenddata(runh,rune)
-c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-
+C
+CBC+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+C
+C     Modified by C. Bigongiari 2001 Jan 16 
+C
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+C      call jcenddata(runh,rune)
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+C
+C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+ 
       WRITE(MONIOU,*)' '
       CALL PRTIME(TTIME)
@@ -1593 +1663 @@
       IF ( MDEBUG .NE. 6 ) CLOSE( MDEBUG )
       CLOSE( EXST )
-c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
-c      CLOSE( PATAPE )
-c      IF ( LCERFI ) CLOSE( CETAPE )
-c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+
+C
+CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+C
+C     Modified by C. Bigongiari 2001 Jan 16 
+C
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+Cc      CLOSE( PATAPE )
+Cc      IF ( LCERFI ) CLOSE( CETAPE )
+Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
+C
+C- close files       
+C
+      
+      CLOSE( PATAPE )
+      IF ( LCERFI ) CLOSE( CETAPE )
+C      CALL JCENDRUN( RUNE ) 
+
+C
+CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+C    
  
       STOP