Changes between Version 25 and Version 26 of InstallingCorsika

Timestamp:

10/20/18 14:33:11 (6 years ago)

Author:

tbretz

Comment:

--

InstallingCorsika

@@ -123 +123 @@
 > **IACTEXT** The interface to theTELOUTfunction is extended by parameters describing theemitting particle.   This extended information is stored as an additional photon bunch(after the normal one) with mass, charge, energy, and emission time replacing thecx,cy,photons, andzemfields, respectively, and are identified by a wavelength of 9999.The compact output format is disabled for making that possible. In addition, all particlesarriving at the CORSIKA observation level are included in the eventio format output file,in a photon-bunch like block identified by array and detector numbers 999.
 
+== The 700nm problem ==
+
+This is were the data is read from the files
+{{{#!
+C=======================================================================
+
+      SUBROUTINE TPDINI( OBLECE )
+
+C-----------------------------------------------------------------------
+[...]
+C  READ THE TABLE OF ATMOSPHERIC EXTINCTION AND STORE IT IN THE
+C  ATMABS MATRIX.
+[...]
+        DO  I = 1, 105
+          READ(MCERABS,21) WLT
+ 21       FORMAT(I4)
+          READ(MCERABS,*,ERR=995,END=995) (ATMABS(I,J), J=0,50)
+        ENDDO
+[...]
+C  CALCULATE THE ATM. EXTINCTION FOR OBS. LEVEL INTERPOLATING THE TABLE
+        DO  IWL = 1, 105
+          FX0 = ATMABS(IWL,X0)
+          FX1 = ATMABS(IWL,X1)
+          ATEOBS(IWL) = LINEAR(X,X0,X1,FX0,FX1)
+        ENDDO
+      ENDIF
+[...]
+C  READ THE QUANTUM EFFECIENCY OF THE PMT
+        OPEN(UNIT=MCERQEF,FILE=DATDIR(1:INDEX(DATDIR,' ')-1)//
+     *          'quanteff.dat',STATUS='OLD',FORM='FORMATTED',ERR=996)
+        READ(MCERQEF,20) TEXTQEF
+        READ(MCERQEF,23,ERR=997,END=997) (QUAEFF(I),I = 1,105)
+ 23     FORMAT(8F6.3)
+        CLOSE( MCERQEF )
+        WRITE(MONIOU,22) TEXTQEF
+      ENDIF
+[...]
+C  READ THE MIRROR REFLECTIVITY
+        OPEN(UNIT=MCERMIR,FILE=DATDIR(1:INDEX(DATDIR,' ')-1)//
+     *          'mirreff.dat',STATUS='OLD',FORM='FORMATTED',ERR=998)
+        READ(MCERMIR,20) TEXTREF
+        READ(MCERMIR,23,END=999,ERR=999) (MIRREF(I), I=1,105)
+        CLOSE( MCERMIR )
+        WRITE(MONIOU,22) TEXTREF
+      ENDIF
+}}}
+
+
+This is where the index or the y-value and the corresponding x-values are calculated for mirror and pde efficiency.
+{{{#!fortran
+C=====================================================================
+
+      SUBROUTINE TELEFF( ABSORB )
+
+C-----------------------------------------------------------------------
+[...]
+C  CALCULATE THE REFERENCE WL AND INDEX OF WL FOR THE INTERPOLATIONS
+      RIWL = 1 + INT( (WL-180.D0)/5.D0 )
+      WLI0 = RIWL*5 + 175
+      WLI1 = RIWL*5 + 175 + 5
+}}}
+
+This is the atmospheric absorption
+{{{#!fortran
+*-- Author : V. de Souza Filho, Uni. Campinas   22/06/1999
+C=======================================================================
+
+      SUBROUTINE ATABSO( ABSORB )
+
+C-----------------------------------------------------------------------
+[...]
+C  CALCULATE THE REFERENCE WL AND INDEX OF WL FOR THE INTERPOLATIONS
+      RIWL = 1 + INT( (WL-180.D0)/5.D0 )
+      WLI0 = RIWL*5 + 175
+      WLI1 = RIWL*5 + 175 + 5
+}}}
+
+
 == Conclusion ==