Changeset 687 for trunk/MagicSoft


Ignore:
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.
Location:
trunk/MagicSoft/Simulation/Corsika/Mmcs
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Simulation/Corsika/Mmcs/cerenkov.f

    r286 r687  
    344344 102      CARTIM = ((ETADSN*(THCKOB(NOBSLV)-THICK(DBLE(ZEMIS)))
    345345     *         /WEMIS+PATHCR-VCOR/BETAE)/C(25)+TIME(NP))* 1.E9
    346           
     346         
    347347c          CALL OUTPT2(IQ(NP),I)
    348348          CALL OUTPT2(NCT,I)
    349 
    350 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    351           call jctime(cartim)
    352 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     349C
     350CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     351C
     352C     Modified by C. Bigongiari 2001 Jan 16
     353C
     354C
     355Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     356C          call jctime(cartim)
     357Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     358C
     359C
     360CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     361C   
    353362
    354363          GOTO 1000
    … …  
    401410      CHARACTER*4      CRUNH,CRUNE,CEVTH,CEVTE
    402411      EQUIVALENCE      (RUNH(1),CRUNH), (RUNE(1),CRUNE)
    403       EQUIVALENCE      (EVTH(1),CEVTH), (EVTE(1),CEVTE)
     412      EQUIVALENCE      (EVTH(1),CEVTH), (EVTE(1),CEVTE)   
    404413*KEEP,CONST.
    405414      COMMON /CONST/   PI,PI2,OB3,TB3,ENEPER
    … …  
    680689          CALL OUTPT2(NCT,I)
    681690
    682 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    683           call jctime(cartim)
    684 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     691C
     692CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     693C
     694C     Modified by C. Bigongiari 2001 Jan 16
     695C
     696C
     697Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     698C          call jctime(cartim)
     699Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     700C
     701C
     702CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     703C   
    685704
    686705          GOTO 1000

  • trunk/MagicSoft/Simulation/Corsika/Mmcs/jcio.c

    r286 r687  
    11/*----------------------------------------------------------------------
    2   jcio.c :
    3 
    4 Rutinas para la gestion de ficheros en CORSIKA 5.20 Las posibilidades
    5 que contempla son:
     2  jcio.c :   
     3
     4Rutinas para la gestion de ficheros en CORSIKA 5.20
     5Author J.C.Gonzalez
     6
     7Las posibilidades que contempla son:
    68
    79- Creacion de un fichero de inicio de run: run######
    … …  
    1416- La grabacion se realiza en C, por lo que no existen marcas de inicio
    1517  ni final de bloque en los ficheros.
     18
     19Modified by C.Bigongiari 2000 Dec 19.
     20No more RUN###### and END###### files created.
     21Three files, CER######, DAT###### and STA###### files, per run instead
     22of three files per shower.
     23Now ###### is the RUN NUMBER not the shower number !!!
    1624----------------------------------------------------------------------*/
    1725
    … …  
    3038FILE *patape;
    3139FILE *sttape;
     40/* CBC
    3241FILE *runtape;
    3342FILE *endtape;
     43CBC*/
    3444char tpl[100];
    3545char certpl[100];
    3646char dattpl[100];
    37 char statpl[100];
     47/* CBC char statpl[100]; */
    3848char cerfile[100];
    3949char datfile[100];
     50/* CBC
    4051char stafile[100];
    4152char runfile[100];
    4253char endfile[100];
    43 
    44 /*--------------------------------------------------
    45   jcinitio:
    46   graba un bloque al fichero de particulas
     54CBC */
     55
     56/*--------------------------------------------------
     57  jcinitio: 
     58
     59Define file names
    4760--------------------------------------------------*/
    4861#ifdef JC_UNDERSCORES
    … …  
    6578  strcpy (dattpl, tpl);
    6679  strcat (dattpl, "dat%06d");
     80  /* CBC
    6781  strcpy (statpl, tpl);
    6882  strcat (statpl, "sta%06d");
    69 
     83  */
     84  /* CBC
    7085  sprintf (runfile, "%srun%06d", tpl, *runnum);
    7186  sprintf (endfile, "%send%06d", tpl, *runnum);
     87  CBC */
     88
     89  sprintf (cerfile, certpl,  *runnum);         /* Added by CB */
     90  sprintf (datfile, dattpl,  *runnum);         /* Added by CB */
    7291}
    7392
    … …  
    107126  jcstartrun:
    108127  abre el fichero run###### y graba el contenido
    109 --------------------------------------------------*/
    110 #ifdef JC_UNDERSCORES
    111 void
    112 jcstartrun_ (float *runh)
    113 #else /* JC_NO_UNDERSCORES */
    114 void
    115 jcstartrun (float *runh)
    116 #endif /* JC_UNDERSCORES */
    117 {
     128
     129NOW open CER and DAT files 
     130--------------------------------------------------*/
     131#ifdef JC_UNDERSCORES
     132void
     133jcstartrun_ (void)
     134#else /* JC_NO_UNDERSCORES */
     135void
     136jcstartrun (void)
     137#endif /* JC_UNDERSCORES */
     138{
     139  /* CBC
     140
    118141  if ((runtape = fopen (runfile, "wb")) == NULL) {
    119142    printf ("JCIO:: Cannot open RUN file %s. Exiting.\n", runfile);
    120143    exit (1);
    121144  }
     145  CBC */
     146
    122147  /* puts("JCIO:: saving runheader buffer..."); */
     148
     149  /* CBC
    123150  fwrite (runh, sizeof (float) * MAXBUF, 1, runtape);
    124 }
    125 
    126 /*--------------------------------------------------
    127   jcendrun:
    128   abre el fichero run###### y graba el contenido
    129 --------------------------------------------------*/
    130 #ifdef JC_UNDERSCORES
    131 void
    132 jcendrun_ (float *rune)
    133 #else /* JC_NO_UNDERSCORES */
    134 void
    135 jcendrun (float *rune)
    136 #endif /* JC_UNDERSCORES */
    137 {
    138   fwrite (rune, sizeof (float) * MAXBUF, 1, runtape);
    139   fclose (runtape);
    140 }
    141 
    142 /*--------------------------------------------------
    143   jcnewcerfile:
    144   abre un nuevo fichero Cherenkov
    145 --------------------------------------------------*/
    146 #ifdef JC_UNDERSCORES
    147 void
    148 jcnewcerfile_ ()
    149 #else /* JC_NO_UNDERSCORES */
    150 void
    151 jcnewcerfile ()
    152 #endif /* JC_UNDERSCORES */
    153 {
    154   sprintf (cerfile, certpl, nshow);
     151
     152  CBC */
     153
     154    /* Following lines up to CBC added by CB */
     155
     156
    155157  if ((cetape = fopen (cerfile, "wb")) == NULL) {
    156158    printf ("JCIO:: Cannot open CER file %s. Exiting.\n", cerfile);
    157159    exit (1);
    158160  }
    159   timefirst = 9.0e10;
    160   timelast = -9.0e10;
    161 }
    162 
    163 /*--------------------------------------------------
    164   jcnewdatfile:
    165   abre un nuevo fichero de particulas
    166 --------------------------------------------------*/
    167 #ifdef JC_UNDERSCORES
    168 void
    169 jcnewdatfile_ (void)
    170 #else /* JC_NO_UNDERSCORES */
    171 void
    172 jcnewdatfile (void)
    173 #endif /* JC_UNDERSCORES */
    174 {
    175   sprintf (datfile, dattpl, nshow);
    176161  if ((patape = fopen (datfile, "wb")) == NULL) {
    177162    printf ("JCIO:: Cannot open DAT file %s. Exiting.\n", datfile);
    178163    exit (1);
    179164  }
    180 }
    181 
    182 /*--------------------------------------------------
    183   jcnewstafile:
    184   abre un nuevo fichero de estadisticas
    185 --------------------------------------------------*/
    186 #ifdef JC_UNDERSCORES
    187 void
    188 jcnewstafile_ (void)
    189 #else /* JC_NO_UNDERSCORES */
    190 void
    191 jcnewstafile (void)
    192 #endif /* JC_UNDERSCORES */
    193 {
    194   sprintf (stafile, statpl, nshow);
    195   if ((sttape = fopen (stafile, "wb")) == NULL) {
    196     printf ("JCIO:: Cannot open STA file %s. Exiting.\n", stafile);
    197     exit (1);
    198   }
    199 }
    200 
    201 /*--------------------------------------------------
    202   jcnewshower:
    203   abre nuevos ficheros para la nueva cascada
    204 --------------------------------------------------*/
    205 #ifdef JC_UNDERSCORES
    206 void
    207 jcnewshower_ (void)
    208 #else /* JC_NO_UNDERSCORES */
    209 void
    210 jcnewshower (void)
    211 #endif /* JC_UNDERSCORES */
    212 {
    213   if (nshow > 0) {
    214     fclose (patape);
    215     fclose (cetape);
    216     fclose (sttape);
    217   }
    218 
    219   nshow++;
    220 #ifdef JC_UNDERSCORES
    221   jcnewdatfile_ ();
    222   jcnewcerfile_ ();
    223   jcnewstafile_ ();
    224 #else /* JC_NO_UNDERSCORES */
    225   jcnewdatfile ();
    226   jcnewcerfile ();
    227   jcnewstafile ();
    228 #endif /* JC_UNDERSCORES */
    229 }
    230 
    231 /*--------------------------------------------------
    232   jcenddata:
    233   abre el fichero run###### y graba el contenido
    234 --------------------------------------------------*/
    235 #ifdef JC_UNDERSCORES
    236 void
    237 jcenddata_ (float *runh, float *rune)
    238 #else /* JC_NO_UNDERSCORES */
    239 void
    240 jcenddata (float *runh, float *rune)
    241 #endif /* JC_UNDERSCORES */
    242 {
    243   if ((endtape = fopen (endfile, "wb")) == NULL) {
    244     printf ("JCIO:: Cannot open END file %s. Exiting.\n", endfile);
    245     exit (1);
    246   }
    247   fwrite (runh, sizeof (float) * MAXBUF, 1, endtape);
    248   fwrite (rune, sizeof (float) * MAXBUF, 1, endtape);
    249   fclose (endtape);
    250 }
     165
     166  /* CBC */
     167}
     168
     169/* WHAT FOLLOWS IS NO MORE USED !!!!!! CB */
     170
     171/*--------------------------------------------------
     172  jcstartshower:
     173 
     174Function added by CB 
     175Sets timefirst and timelast variables to their initial value at the beginnng
     176of new shower   
     177--------------------------------------------------*/
     178#ifdef JC_UNDERSCORES
     179void
     180jcstartshower_ (float *evth)
     181#else /* JC_NO_UNDERSCORES */
     182void
     183jcstartshower (float *evth)
     184#endif /* JC_UNDERSCORES */
     185{
     186  timefirst = 9.0e10;
     187  timelast = -9.0e10;
     188}
     189
    251190
    252191/*--------------------------------------------------
    … …  
    357296
    358297}
     298
     299/*--------------------------------------------------
     300  jcendrun:
     301 
     302  Closes STA file
     303--------------------------------------------------*/
     304#ifdef JC_UNDERSCORES
     305void
     306jcendrun_ (float *rune)
     307#else /* JC_NO_UNDERSCORES */
     308void
     309jcendrun (float *rune)
     310#endif /* JC_UNDERSCORES */
     311{
     312  /* fwrite (rune, sizeof (float) * MAXBUF, 1, sttape); */
     313  fclose (sttape);
     314}

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

    r287 r687  
    425425      ILEFTA = 0
    426426c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    427  
    428 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     427C
     428CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     429C
     430C     Modified by C. Bigongiari 2001 Jan 16
     431C
     432C
     433Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    429434      print *,'JCIO::========================================'
    430435      print *,'JCIO:: Initializing JCIO system for advanced'
    431436      print *,'JCIO:: saving of data.'
    432437      print *,'JCIO::========================================'
    433 c- initialize jcio system
     438C
     439Cc- initialize jcio system
     440C     
    434441      call jcinitio(dsn,nrrun)
    435 c- create file run######
    436       call jcstartrun(runh)
    437 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    438 
     442Cc- create file run######
     443C      call jcstartrun(runh)
     444Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     445C
     446C- Modified JCSTARTRUN creates cer###### and dat###### files !
     447C
     448C     ###### is the RUN number !
     449C
     450
     451      call jcstartrun(RUNH)
     452
     453C
     454C- write Run Header on cer and dat files 
     455C
     456      CALL TOBUF(RUNH,0)
     457      IF ( LCERFI ) CALL TOBUFC(RUNH,0)
     458
     459C
     460CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     461C
    439462C-----------------------------------------------------------------------
    440463C  LOOP OVER SHOWERS
     464
    441465      DO 2  ISHW = 1,NSHOW
     466
    442467 
    443468c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    … …  
    482507          FPRINT = .FALSE.
    483508        ENDIF
    484  
    485 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    486 c   Create cer######,dat######,sta###### files
    487 c------------------------------------------------------------
    488         call jcnewshower
    489 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    490 
     509CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     510Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     511Cc   Create cer######,dat######,sta###### files
     512Cc------------------------------------------------------------
     513C        call jcnewshower
     514Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     515CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     516 
    491517C  RESET COUNTERS
    492518        DO 447  K = 1,25
    … …  
    823849        ENDIF
    824850        EVTH(147) = 0.
    825  
     851
    826852        IF ( FEGS ) THEN
    827853C  PARAMETER FOR ELECTRON AND PHOTON REJECT (CONVERT ENERGY TO MEV)
    … …  
    838864        EVTH( 9) = PTOT0 * ST * SIN(PHI)
    839865        EVTH(10) = PTOT0 * COSTHE
     866     
     867
    840868c
    841869c [*] one block from here sent above
    842870c
     871
    843872        EVTH(85) = CERSIZ
    844873 
    … …  
    856885        IF ( EVTH(3) .GT. 3.0  .OR.  .NOT. FEGS ) THEN
    857886          CALL TOBUF ( EVTH,0 )
    858           IF ( LCERFI ) CALL TOBUFC( EVTH,0 )
     887          IF ( LCERFI ) CALL TOBUFC( EVTH,0 ) 
    859888        ENDIF
    860889 
    … …  
    865894     *          '    GEN      LEVEL  E ON STACK'/)
    866895        NOPART = 0
    867  
    868  
     896
    869897        IF ( CURPAR(1) .LE. 3.D0  .OR.
    870898     *      (CURPAR(1) .GE. 5.D0  .AND.  CURPAR(1) .LE. 7.D0) ) THEN
    … …  
    889917C  HADRONIC PARTICLES
    890918          FNPRIM = .TRUE.
    891  
     919
    892920C  FILL LONGITUDINAL DISTRIBUTION FOR THE PRIMARY PARTICLE
    893921C  THE PARTICLE IS TRACKED FROM THICK0 DOWN TO THICK0+CHI*PRMPAR(3)
    … …  
    917945            ENDIF
    918946          ENDIF
    919  
     947
    920948C  CHECK OBSERVATION LEVEL PASSAGE AND UPDATE PARTICLE COORDINATES
    921949          HNEW = H
    … …  
    944972        ENDIF
    945973 
     974
    946975C-----------------------------------------------------------------------
    947976C  NORMAL CYCLE
    … …  
    10271056C  FINISH SHOWER AND PRINT INFORMATION
    10281057        CALL OUTEND
    1029  
    10301058 
    10311059*       IF ( DEBUG ) WRITE(MDEBUG,442) NPARTO
    … …  
    12211249        ENDIF
    12221250 
    1223 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     1251C
     1252CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     1253C
     1254C     Modified by C. Bigongiari 2001 Jan 16
     1255C
     1256C
     1257Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    12241258c   Saves statistics to sta###### file
    1225         call jcstadata(EVTH,EVTE,
    1226      +      NPROTO,NPROTB,NNEUTR,NNEUTB,NPHOTO,NELECT,NPOSIT,
    1227      +      NNU   ,NMUM  ,NMUP  ,NPI0  ,NPIM  ,NPIP  ,NK0L  ,
    1228      +      NK0S  ,NKMI  ,NKPL  ,NHYP  ,NDEUT ,NTRIT ,NALPHA,
    1229      +      NOTHER,IFINNU,IFINPI,IFINET,IFINKA,IFINHY,
    1230      +      CERELE,CERHAD,PLONG,LPCT1,NSTEP,THSTEP)
    1231 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    1232 
    1233 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     1259C        call jcstadata(EVTH,EVTE,
     1260C     +      NPROTO,NPROTB,NNEUTR,NNEUTB,NPHOTO,NELECT,NPOSIT,
     1261C     +      NNU   ,NMUM  ,NMUP  ,NPI0  ,NPIM  ,NPIP  ,NK0L  ,
     1262C     +      NK0S  ,NKMI  ,NKPL  ,NHYP  ,NDEUT ,NTRIT ,NALPHA,
     1263C     +      NOTHER,IFINNU,IFINPI,IFINET,IFINKA,IFINHY,
     1264C     +      CERELE,CERHAD,PLONG,LPCT1,NSTEP,THSTEP)
     1265Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     1266Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    12341267C  WRITE SHOWER END TO OUTPUT BUFFER
    1235 c        CALL TOBUF( EVTE,0 )
    1236         CALL TOBUF( EVTE,1 )
    1237 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     1268Cc        CALL TOBUF( EVTE,0 )
     1269C        CALL TOBUF( EVTE,1 )
     1270Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     1271C        IF ( LCERFI ) THEN
     1272C          CALL OUTND2
     1273Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     1274Cc          CALL TOBUFC( EVTE,0 )
     1275Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     1276C        ENDIF
     1277
     1278C
     1279C  WRITE SHOWER END TO OUTPUT BUFFER
     1280        CALL TOBUF( EVTE,0 )
    12381281        IF ( LCERFI ) THEN
    12391282          CALL OUTND2
    1240 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    1241 c          CALL TOBUFC( EVTE,0 )
    1242 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    1243         ENDIF
     1283          CALL TOBUFC( EVTE,0 )
     1284        ENDIF
     1285
     1286C
     1287CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     1288C
    12441289 
    12451290        IF ( FPRINT .OR. DEBUG ) WRITE(MONIOU,*)
    … …  
    12591304          JHBIN(J) = JHBIN(J) + IHBIN(J)
    12601305  19    CONTINUE
    1261  
     1306
     1307
    12621308  2   CONTINUE
     1309
    12631310C  END OF SHOWER LOOP
    12641311 
    … …  
    12711318      RUNE(3) = REAL(ISHW)
    12721319C  WRITE RUN END TO OUTPUT BUFFER AND FINISH OUTPUT
    1273 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    1274 c      CALL TOBUF ( RUNE,1 )
    1275       call jcendrun(rune)
    1276 c      IF ( LCERFI ) CALL TOBUFC( RUNE,1 )
     1320C
     1321CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     1322C
     1323C     Modified by C. Bigongiari 2001 Jan 16
     1324C
     1325Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     1326Cc      CALL TOBUF ( RUNE,1 )
     1327C      call jcendrun(rune)
     1328Cc      IF ( LCERFI ) CALL TOBUFC( RUNE,1 )
     1329Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     1330C
     1331C- write Run End     
     1332C
     1333
     1334      CALL TOBUF ( RUNE,1 )
     1335      IF ( LCERFI ) CALL TOBUFC( RUNE,1)
     1336
     1337C
     1338CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     1339C   
    12771340c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    12781341C  TIME SINCE BEGINNING
    … …  
    13081371 209  FORMAT(//' ENERGY - MULTIPLICITY MATRIX FOR ALL SHOWERS'/
    13091372     *       ' ENERGY RUNS VERTICALLY, MULTIPLICITY HORIZONTALLY'//,
    1310      *       ' ',6X,5I9,3I8,5I7,'   ENERGY RANGE (GEV)'/
     1373     *       ' ',6X,5I9,3I8,5I7'   ENERGY RANGE (GEV)'/
    13111374     *       37(/' ',I4,1X,I10,4I9,3I8,5I7,1X,1P,2E10.1,0P)//
    13121375     *       ' MULT. ',5I9,3I8,5I7,4X,'LOWER BIN LIMIT'/
    … …  
    15801643      ENDIF
    15811644 
    1582 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    1583       call jcenddata(runh,rune)
    1584 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    1585 
     1645C
     1646CBC+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     1647C
     1648C     Modified by C. Bigongiari 2001 Jan 16
     1649C
     1650Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     1651C      call jcenddata(runh,rune)
     1652Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     1653C
     1654C+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     1655 
    15861656      WRITE(MONIOU,*)' '
    15871657      CALL PRTIME(TTIME)
    … …  
    15931663      IF ( MDEBUG .NE. 6 ) CLOSE( MDEBUG )
    15941664      CLOSE( EXST )
    1595 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
    1596 c      CLOSE( PATAPE )
    1597 c      IF ( LCERFI ) CLOSE( CETAPE )
    1598 c>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     1665
     1666C
     1667CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     1668C
     1669C     Modified by C. Bigongiari 2001 Jan 16
     1670C
     1671Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     1672Cc      CLOSE( PATAPE )
     1673Cc      IF ( LCERFI ) CLOSE( CETAPE )
     1674Cc>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
     1675C
     1676C- close files       
     1677C
     1678     
     1679      CLOSE( PATAPE )
     1680      IF ( LCERFI ) CLOSE( CETAPE )
     1681C      CALL JCENDRUN( RUNE )
     1682
     1683C
     1684CBC++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     1685C   
    15991686 
    16001687      STOP
Note: See TracChangeset for help on using the changeset viewer.