COSMOS v7.655  COSMOSv7655
(AirShowerMC)
interface2.h
Go to the documentation of this file.
1  integer leng, lengn, lengid
2  integer lengdir
3  character*120 dir
4  integer kgetenv2
5  character*8 numb
6  character*64 execid
7  character*128 msg
8  logical takehist
9  save
10 
11  do i= 1, nsites
12  histdep(i) = 0
13  indivdep(i) = 0
14  enddo
15 
16 
17 
18  leng = kgetenv2("INDIVDEP", input)
19  read(input, *) indivdep
20 
21  leng = kgetenv2("OUTPUT", input)
22  read(input, *) tkrtspec, tkarspec
23 
24  leng = kgetenv2("HISTDEP", input)
25  read(input, *) histdep
26  takehist = .false.
27  do i= 1, nsites
28  if(histdep(i) .gt. 0 ) then
29  takehist=.true.
30  exit
31  endif
32  enddo
33  takehist = takehist .and. (tkrtspec .or. tkarspec)
34 
35  leng = kgetenv2("BINW", input)
36  read(input, *) binw
37 
38 
39 
40  input = ' '
41 #if defined (KEKB) || (KEKA)
42  numb="rank-#"
43  lengn=6
44 #else
45  lengn = kgetenv2("NUMB", numb)
46 #endif
47  lengdir = kgetenv2("OUTDIR", input)
48  dir = input
49  lengid = kgetenv2("EXECID", execid)
50 
51 
52 ! msg =input(1:lengdir)//"/"//execid(1:lengid)//
53 ! * "-@."//numb(1:lengn)//".nrfai"
54 ! call copenfw2(fnonrfai, msg, 1, icon)
55 ! if(icon .gt. 1) then
56 ! write(0,*) ' icon=', icon
57 ! call cerrorMsg(msg, 1)
58 ! call cerrorMsg('could not be opened', 0)
59 ! endif
60 #if FNODATDEF > 0 && BUFSIZE < 200000
61 ! if(fnodat .gt. 0) then
62  msg =input(1:lengdir)//"/"//execid(1:lengid)//
63  * "-@."//numb(1:lengn)//".dat"
64  call copenfw2(fnodat, msg, 2, icon)
65  if(icon .gt. 1) then
66  write(0,*) ' icon=', icon
67  call cerrorMsg(msg, 1)
68  call cerrorMsg('could not be opened', 0)
69  endif
70 ! endif
71 #endif
72  if(takehist) then
73  msg =input(1:lengdir)//
74  * "/"//execid(1:lengid)//"-@."//numb(1:lengn)//".hist"
75  call copenfw2(fno, msg, binw, icon) ! binary write
76  if(icon .gt. 1) then
77  write(0,*) ' icon=', icon
78  call cerrorMsg(msg, 1)
79  call cerrorMsg('could not be opened', 0)
80  endif
81  endif
read
do nsites input read(input, *) histdep call cqUhookc(2
histdep
integer lengid integer lengdir character *dir integer kgetenv2 character *numb character *execid character *msg logical takehist save do nsites histdep(i)=0indivdep(i)=0enddo leng
integer
int integer
Definition: Zdef.h:3
kgetenv2
integer function kgetenv2(envname, envresult)
Definition: cgetLoginN.f:77
i
integer lengid integer lengdir character *dir integer kgetenv2 character *numb character *execid character *msg logical takehist save do i
Definition: interface2.h:11
lengn
integer lengn
Definition: interface2.h:1
true
block data cblkElemag data *AnihiE ! Eposi< 1 TeV, anihilation considered *X0/365.667/, ! radiation length of air in kg/m2 *Ecrit/81.e-3/, ! critical energy of air in GeV *MaxComptonE/1./, ! compton is considered below 1 GeV *MaxPhotoE/1.e-3/, ! above this, PhotoElectric effect neg. *MinPhotoProdE/153.e-3/, ! below 153 MeV, no gp --> hadrons ! scattering const not MeV *Knockon true
Definition: cblkElemag.h:7
leng
integer leng
Definition: interface2.h:1
tkrtspec
logical KeepWeight ! see setupenv sh logical tkarspec ! get lateral histo in a web fai bin logical tkrtspec ! get time histo in each web bin logical SeeLowdE common Zprivatec2 tkrtspec
Definition: Zprivate4.h:7
copenfw2
subroutine copenfw2(io, fnin, form, icon)
Definition: copenf.f:205
cerrorMsg
integer lengid integer lengdir character *dir integer kgetenv2 character *numb character *execid character *msg logical takehist save do nsites tkarspec nsites icon ! call cerrorMsg(msg, 1) ! call cerrorMsg('could not be opened'
tkarspec
logical KeepWeight ! see setupenv sh logical tkarspec ! get lateral histo in a web fai bin logical tkrtspec ! get time histo in each web bin logical SeeLowdE common Zprivatec2 tkarspec
Definition: Zprivate4.h:7
logical
int logical
Definition: Zdef.h:1
if
integer lengid integer lengdir character *dir integer kgetenv2 character *numb character *execid character *msg logical takehist save do nsites tkarspec nsites if(histdep(i) .gt. 0) then takehist
false
block data cblkElemag data *AnihiE ! Eposi< 1 TeV, anihilation considered *X0/365.667/, ! radiation length of air in kg/m2 *Ecrit/81.e-3/, ! critical energy of air in GeV *MaxComptonE/1./, ! compton is considered below 1 GeV *MaxPhotoE/1.e-3/, ! above this, PhotoElectric effect neg. *MinPhotoProdE/153.e-3/, ! below 153 MeV, no gp --> hadrons ! scattering const not MeV *Knockon ! knockon is considered Obsolete *PhotoProd false
Definition: cblkElemag.h:7
modassize::fno
integer, parameter fno
Definition: chook.f:43