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. 

<<Azimuth>>   complex*16 
 Range of azimuthal angle in deg. Say, (0, 45). Defaulti is (0,360). 
 Can be such as (300., 390.) Used when Za1ry is 'is' 

<<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 "ContInfo". This will be created when job 
 is finished normally. 

<<CosZenith>>   complex*16 
 Range of cos(zenith angle). Say, (0.5, 1.0). Used when Za2ry is 'is' 

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

<<Ddelta>>   real*8 
 SourceDec $\pm$ 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/m$^2$. If $< 0$, HeightList has priority. 
 (See HeightList) 

<<DestEventNo>>   integer 
 2 integers: Final event no. to be generated and events to be generated 
 in the current run. 

<<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). 
 If this is $<$ deepest obs. level and zeinth angle of primary is $< 0$, 
 the primary is assumed to be upgoing even if Reverse =0. 
 NOTE: BorderHeightH must be given explicitly in this case. 

<<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*16 
 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?). 
 If 'inclusive', an inclusive treatment is employed as given in Contrib/Inclusive. 
 To make this work actually, INCLUSIVE in cosmos/Zcondc.h must be given a value of 1 
 when making the cosmos library. At execution time, InclusiveFile must be given 
 the path to the xdist.d file. (normally in Contrib/Inclusive/xdist.d) 
 The treatmen is one dimesional. May be used for neutrino calculations. 

<<Job>>   character*10 
 What kind of job you are going to do.\newline 
 =' ' (default). nothing special.\newline 
 ='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.\newline 
 ='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. \newline 
 0$ \Rightarrow $ 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.)\newline 
 1$ \Rightarrow $ detector at the observation place is horizontal. \newline 
 2$ \Rightarrow $ detector is perpendicular to the primary. \newline 
 3$ \Rightarrow $ 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. 
 1: onedim without use of table. 
 2: table is used for thickness <--> length conversion. if cos < .5 
 3: table is always used for any angle. 
 ( for height > 30 km, table is not used in any case). 

<<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 
 File to contain 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 'Seed'. For a normal run with 
 Job=' ', if SeedFile is not ' ', two integer initial random numbers 
 and the event number are automatically output on the speicfied disk file. 

<<SkeletonFile>>   character*70 
 Skeleton information file (full or relative path). created if Job = 
 'skeleton'. Default is 'skeletonParam'. This is the Namelist data 
 referred by Cosmos automatically if Job='flesh' is specified. For 
 Job='flesh', you have to modify some part of this file. 

<<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 sys. 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 ==> subtract 100 and apply the above, but chookTrace or chookCeren["",S,E] are 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 1999.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) 

