The description of parameters frequently used:  &Param data
----------------------------------------------------------

<<ASDepthList>>   real*8 
 This is DepthList for AS observation. Used only if Generate contains 
 "as". See DepthList. 

<<ASHeightList>>   real*8 
 This is HeightList for Air Shower observ. Used only if Generate contains 
 "as". See HeightList. 

<<BaseTime>>   real*8 
 Rough cpu time needed for completing one event (say, for protons, or 
 gamma rays) with energy BaseErg. The cpu time estimation is based on 
 A * ( E1ry par nucleon )**BasePower / BaseErg * BaseTime, where A is mass number 
 (for nucleus; otherwise 1). 

<<Cont>>   logical 
 If T, continuation from a previous job is assumed. Contfile content is used. 

<<ContFile>>   character*70 
 Job continuation information file (full or relative path). 
 default is ../Data/.contInfo. This will be created when job 
 is finished normally. 

<<CosZenith>>   complex*16 
 Range of cos(zenith angle). Say, (0.5, 1.0) 

<<CutOffFile>>   character*70 
 Geomagnetic cut-off file 

<<Ddelta>>   real*8 
 SourceDec +- Ddelta is the region for 'aps' (deg). 

<<DeadLine>>   character*8 
 The dead line before which the job should terminate. 
 Should be given like '10.11.15' which means the nearest 10th, 11 O'clock, 
 15 min. Not used if Within has non zero value. 

<<DepthList>>   real*8 
 Depth List of Observation level in kg/m2. If < 0, HeightList has priority. 
 (See HeightList) 

<<DestEventNo>>   integer 
 Final event no. to be generated. 

<<DtGMT>>   real*8 
 Difference of the local time of the observation place from GMT (hour). 

<<Freec>>   logical 
 if F, the first interaction point is made to be the injection point else randomly 
 sampled. 

<<Generate>>   character*16 
 specify the system should generate 1) EM-cascade(em), 2) one dimensional AS(as/qas) 
 and/or AS Lateral distribution(lat). Generate='em,as' or 'em/as' etc (order/case/separator 
 insensitive) is to generate EM-cascade and AS. Generate='as' will generate AS with some 
 adequate EM cascade (i.e.,the EM cascade is not followed down to KEminObs). 
 If 'qas' is given, quick generation of AS for heavy primaries is tried. See chookASbyH.f 

<<HeightList>>   real*8 
 Height of observation levels in m. This is made from DepthList internally. 
 I.e., this one is usually not an input. However, if the DepthList values are 
 negative, this is used as input and corresponding DepthList is computed internally. 

<<HeightOfInj>>   real*8 
 The vertical height of primary injection point (m). 

<<Hidden>>   logical 
 Make T, if hidden parameters are to be written. 

<<IncMuonPolari>>   logical 
 if T, consider muon polarization 

<<InitRN>>   integer 
 Initial random number seed. 2 integers. 

<<IntModel>>   character*6 
 interaction model. If 'int1', Lund Fritiof/Nucrin/adhoc models are used 
 for hadronic interactions. If 'int2' Gheisha code is used for hadronic interactions, 
 at low energies. (See, Elund?). 

<<Job>>   character*10 
 What kind of job you are going to do. 
 =' ' (default). nothing special. 
 ='skeleton'. Makes skeleton. Later you can add AS, follow particles 
 to a lower energy and/or to a deeper depth, by giving Job='flesh'. 
 'skeletonParm' is automatically created. 
 ='flesh'. Using information created by Job='skeleton' and the 
 initial random number for each event, you can add AS, 
 follow particles to a lower energy and/or to a deeper depth. 
 The namelist parameters you give should be 'skeletonParam' 
 with some necessary modification. 

<<KEminObs>>   real*8 
 The min kinetic energy of particles for observation. 

<<LatitOfSite>>   real*8 
 Latitude of the deepest observation level in degree. East is positive. 

<<LongitOfSite>>   real*8 
 Longitude of the deepest observation level in degree. North is positive. 

<<LpmEffect>>   logical 
 If t, the LPM effect is considered when Ee > LpmBremEmin for electrons and 
 Eg > LpmPairEmin for gamma rays. 

<<MinPhotoProdE>>   real*8 
 Below this energy, no photo-prod of hadron. See also PhotoProd. 

<<ObsPlane>>   integer 
 How to observe particles. 
 0==> no detector plane is used for observation. BorderHeightL is used to detect particles. This is 
 for, say, neutrino observation. See BorderHeightL. (HeigtList and DepthList are neglected. 
 ASHeightList and ASDepthList are still used to observe Air Shower.) 
 1==> detector at the observation place is horizontal. 
 2==> detector is perpendicular to the primary. 
 3==> detector is user defined (not supported) 

<<OneDim>>   integer 
 If 0, 3 dimensional simulation. if 1, one 
 dimensional simulation is performed. For flux calculation, 
 sometimes one dimensional one is better. 

<<PhotoProd>>   logical 
 Switch. if .false., no photo prod. of hadron is considered at all. 
 See also MinPhotoProdE. 

<<PrimaryFile>>   character*70 
 Primary Spectrum data file (full or relative path) 

<<SeedFile>>   character*70 
 Files containing the initial random numbers for those events to which you 
 want to flesh. You can create the file by calling cwriteSeed in a 
 user hook routine (say, in chookEnEvent) at skeleton making time. 
 default is '../Data/seed'. 
 If you want to create the file not at the time of making skeleton, 
 you have to write two integer initial random numbers, 
 and the event number (and more optional ones) 
 This file is opened automatically, if needed. 

<<SkeletonFile>>   character*70 
 Skeleton information file (full or relative path). created if Job = 
 'skeleton'. Default is ../Data/skeletonParam 
 This is the Namelist data which should be used as input at the time of 
 Job='flesh' after modifying parameters in &Param there. 

<<SourceDec>>   real*8 
 Source declination of point source.(deg) 

<<ThinSampling>>   logical 
 if F, thinsampling is not tried. if T, alla Hillas thinning. Don't use with 
 the skeleton/flesh method 

<<TimeStructure>>   logical 
 If T, time information is computed 

<<Trace>>   integer 
 Flag for trace information output. 
 0 ==> no trace information is output. 
 <10==> x, y, z in the primary system (say, 1) 
 <20 ==> x, y, in the primary system. z in kg/m2. (say,11) 
 <30 ==> x, y, z in the detector system 
 <40 ==> x, y, in the detector system. z in kg/m2 
 <50 ==> x, y, z in 'xyz' system 
 <60 ==> x,y, in 'xyz' and z in kg/m2 
 61-100 ==> for Cherenkov observation. For Coord system subtract 60. 
 if the value is even, binary output is made on TraceDev. 
 > 101 ==> The user hook for trace is used. 
 Primary system: Origin is the deepest detector. Z-axis is the primary direction. 
 X-axis is Z x Vertical axis. X-Y plane is orthogonal to the primary. 
 Detector system: origin is the deepest detector. Z-axis is the vertical one. X-axis is 
 directed to the magnetic east. X-Y plane is horizontal. 
 z in kg/m2 : Vertical depth in kg/m2 above the deepest detector to the particle. 

<<TraceDir>>   character*70 
 Directory. Default Trace information is put TraceDir/trace1, 2,.. 
 for event 1, 2, ... The directory should exist. Default is ' ' and in this case 
 /tmp/YourLoginName/ is employed. The last "/" should not be given. 
 *** NOTE that default Cherenkov output is made only using TraceDev, 
 TraceDir is not used. You have to open the disk file at chookbgRun 
 It can by binary or ascii file depending on Trace value. 

<<WaitRatio>>   real*8 
 Wait A.S generation until the electron energy, Ee, becomes < WaitRatio* E0. 
 This many be 1.0 for hadron origin case. But for gamma/electron primary, 
 this should be as low as 0.01 to enjoy full fluctuation. 

<<Within>>   integer 
 The job should end within this minutes from now. Default is 99999. 
 If 0 is given, DeadLine is used. 

<<YearOfGeomag>>   real*8 
 Like 1994.5. Year when Geomagnetic field is to be calculated. 

<<Za1ry>>   character*4 
 Specify the primary angle sampling method by one of 'is', 'ps' or 'aps'. 
 "is" is isotropic. The range is by CosZenith. 
 "ps" is for point source (See also SourceDec) 
 "aps" is around point source (See also SourceDec and Ddelta) 

