COSMOS v7.655  COSMOSv7655
(AirShowerMC)
cstdatmos0.fNew2.h File Reference
#include "Zatmos.h"
Include dependency graph for cstdatmos0.fNew2.h:

Go to the source code of this file.

Typedefs

using ans = atmos%rho(atmos%nodes) **exp(-(z-atmos%z(atmos%nodes))/Hinf) elseif(z .lt. atmos%z(1)) then ans=atmos%rho(1) **exp((atmos%z(1) -z)/atmos%H(1)) else call kdwhereis(z, atmos%nodes, atmos%z, 1, i) a=atmos%a(i) if(a .ne. 0.d0) then ans=atmos%rho(i) **(1+a *(z-atmos%z(i))/atmos%H(i)) **(-1.0d0-1.d0/a) else ans=*atmos%rho(i) *exp(-(z-atmos%z(i))/atmos%H(i)) endif endif ! zsave=z ! endif cvh2den=ans end ! ---------------------------------- real *8 function cvh2temp(z) implicit none ! vettical height to temperatur(Kelvin) real *8 z ! input. vertical height in m ! output is temperature of the atmospher in Kelvin real *8 ans integer i if(z .gt. atmos%z(atmos%nodes)) then ans=atmos%T(atmos%nodes) elseif(z .lt. atmos%z(1)) then ans=atmos%T(1)+atmos%b(1) *(z - atmos%z(1)) else call kdwhereis(z, atmos%nodes, atmos%z, 1, i) ans=atmos%T(i)+atmos%b(i) *(z-atmos%z(i)) endif cvh2temp=ans end !--------------------------------------------- real *8 function cthick2h(t) implicit none real *8 t ! input. air thickness in kg/m^2 real *8 logt, ans integer i real *8 dod0, fd, a logt=log(t) if(t .ge. atmos%cumd(1)) then ans=atmos%z(1) - *(logt - atmos%logcumd(1)) *atmos%H(1) elseif(t .le. atmos%cumd(atmos%nodes)) then ans=atmos%z(atmos%nodes) - *Hinf *log(t/atmos%cumd(atmos%nodes)) else call kdwhereis(t, atmos%nodes, atmos%cumd, 1, i) ! i is such that X(i) > x >=x(i+1)
 

Functions

nodes a ! write (0, *) ' atmos.H
 
nodes a atmos H (i)
 
nodes a atmos atmos z (i) ! endif cthick2h
 

Variables

 ! i =1
 
nodes dod0
 
nodes t
 
nodes z ='
 
nodes i =',i, ' dod0='
 
nodes a ='
 
nodes a atmos atmos temp realcthick2h
 

Typedef Documentation

◆ ans

nodes a atmos atmos ans

Definition at line 82 of file cstdatmos0.fNew2.h.

Function Documentation

◆ ! write()

nodes a ! write ( ,
 
)

◆ H()

nodes a atmos H ( i  )

◆ z()

nodes a atmos atmos z ( i  )

Variable Documentation

◆ ! i

! i =1

Definition at line 83 of file cstdatmos0.fNew2.h.

◆ a

dE dx *! Nuc Int sampling table a ='

Definition at line 93 of file cstdatmos0.fNew2.h.

Referenced by __det2exyz.f__(), __xyz2det.f__(), atmncobptcl(), cbremsenergy(), cgpxs1(), cinippx(), cinippxn(), ckmpelaxs(), ckmptotxs(), ckpptotxs(), cmydecay(), cnptotxs(), cpimpelaxs(), cpimptotxs(), cpippelaxs(), cpipptotxs(), cppelaxs(), cpptotxs(), crossc_kk(), cutbyangle(), dpmjet(), epgencolbyphoton(), gamfun(), gamfun_kk(), gaucr(), gencol(), ixxdef(), ixxson(), k4ptdis(), kadbp(), kadth(), kadthi(), kb50j2(), kcelei(), kcosd(), kctoq(), kdcmjd(), kdhtoh(), kdifva(), kdtoa(), kdtoe(), kdtoh(), kdzth2(), kdztoh(), keqtog(), ketod(), ketoh(), kgcrc(), kgcttc(), kgdted(), kgtoeq(), khtad(), khtoe(), kj2b50(), kj2j90(), kj2tox(), kj90j2(), kjxjy(), kmjd(), kmjdst(), kmjdym(), kmobec(), kmoon(), kmoont(), koangl(), komeg0(), kpmtrx(), kqtoc(), krkfehl(), krkfhstep(), krkflstep(), krungekutfds(), krungekutfs(), kside0(), ksided(), ksidet(), ksind(), ksun(), ksuneq(), ktu(), kvtoa(), kxtoj2(), mfnems(), minuit(), mnamin(), mnbins(), mncalf(), mncler(), mncntr(), mncomd(), mncont(), mncrck(), mncros(), mncuve(), mnderi(), mndxdi(), mneig(), mnemat(), mnerrs(), mneval(), mnexcm(), mnexin(), mnfixp(), mnfree(), mngrad(), mnhes1(), mnhess(), mnimpr(), mninex(), mninit(), mninpu(), mnintr(), mnlims(), mnline(), mnmatu(), mnmigr(), mnmnos(), mnmnot(), mnparm(), mnpars(), mnpfit(), mnpint(), mnplot(), mnpout(), mnprin(), mnpsdf(), mnrazz(), mnread(), mnrn15(), mnrset(), mnsave(), mnscan(), mnseek(), mnset(), mnseti(), mnsimp(), mnstat(), mnstin(), mntiny(), mnvers(), mnvert(), mnwarn(), mnwerr(), outresul(), outresula(), outresulb(), outresulc(), outresulx(), outtrace(), psaini(), psapint(), psasetc(), psbint(), psborn(), pscajet(), psconf(), pscs(), psdefrot(), psdeftr(), psdr(), psfap(), psfaz(), psfborn(), psfsh(), psftild(), psgea(), psgint(), pshard(), pshot(), psjdef(), psjet(), psjet1(), psjint(), psjint0(), psjint1(), pslam(), psnorm(), psqint(), psrec(), psrejs(), psrejv(), psrjint(), psroot(), psrotat(), psshar(), pstrans(), pstrans1(), psudint(), psuds(), psudt(), psv(), psvdef(), pszsim(), puttrace(), qgs01init(), sectnu(), sortbyke(), stand(), translate_aline(), xxaini(), xxaset(), xxddfr(), xxdec2(), xxdec3(), xxdpr(), xxdtg(), xxfau(), xxfrag(), xxfragm(), xxfz(), xxgau(), xxgau1(), xxgener(), xxjetsim(), xxreg(), xxrot(), xxstr(), xxt(), and xxtwdec().

◆ cthick2h

nodes a atmos atmos temp real* cthick2h

◆ dod0

nodes dod0
Initial value:
=( atmos%cumd(i) - t )/ atmos%d0(i)
a = atmos%a(i)
fd = 1. - dod0
if(a .ne. 0.) then
ans = (fd**(-a)- 1.0d0)*atmos%H(i)/a + atmos%z(i)
else
ans = -log(fd)* atmos%H(i) + atmos%z(i)
endif
!
! write(0,*) ' t='
nodes z
atmos%rho(atmos%nodes) **exp(-(z-atmos%z(atmos%nodes))/Hinf) elseif(z .lt. atmos%z(1)) then ans=atmos%rho(1) **exp((atmos%z(1) -z)/atmos%H(1)) else call kdwhereis(z, atmos%nodes, atmos%z, 1, i) a=atmos%a(i) if(a .ne. 0.d0) then ans=atmos%rho(i) **(1+a *(z-atmos%z(i))/atmos%H(i)) **(-1.0d0-1.d0/a) else ans=*atmos%rho(i) *exp(-(z-atmos%z(i))/atmos%H(i)) endif endif ! zsave=z ! endif cvh2den=ans end ! ---------------------------------- real *8 function cvh2temp(z) implicit none ! vettical height to temperatur(Kelvin) real *8 z ! input. vertical height in m ! output is temperature of the atmospher in Kelvin real *8 ans integer i if(z .gt. atmos%z(atmos%nodes)) then ans=atmos%T(atmos%nodes) elseif(z .lt. atmos%z(1)) then ans=atmos%T(1)+atmos%b(1) *(z - atmos%z(1)) else call kdwhereis(z, atmos%nodes, atmos%z, 1, i) ans=atmos%T(i)+atmos%b(i) *(z-atmos%z(i)) endif cvh2temp=ans end !--------------------------------------------- real *8 function cthick2h(t) implicit none real *8 t ! input. air thickness in kg/m^2 real *8 logt, ans integer i real *8 dod0, fd, a logt=log(t) if(t .ge. atmos%cumd(1)) then ans=atmos%z(1) - *(logt - atmos%logcumd(1)) *atmos%H(1) elseif(t .le. atmos%cumd(atmos%nodes)) then ans=atmos%z(atmos%nodes) - *Hinf *log(t/atmos%cumd(atmos%nodes)) else call kdwhereis(t, atmos%nodes, atmos%cumd, 1, i) ! i is such that X(i) > x >=x(i+1) ans
nodes i
nodes dod0
int ne[nl][nth]
Definition: Zprivate.h:11
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
Definition: cblkEvhnp.h:5
nodes a
nodes t
nodes a atmos H(i)

Definition at line 84 of file cstdatmos0.fNew2.h.

◆ i

dE dx i =',i, ' dod0='

Definition at line 93 of file cstdatmos0.fNew2.h.

Referenced by __add2hyb.f__(), __add2nrfai.f__(), __ascii2bin.f__(), __atmosd.f__(), __atmosd2.f__(), __atmosd3.f__(), __atmosd4.f__(), __bin2ascii.f__(), __bin2bin.f__(), __ccalcstdatmos.f__(), __checknrfai.f__(), __getacceptance.f__(), __getst.f__(), __getst2.f__(), __getst3.f__(), __gpinte.f__(), __mknrfai.f__(), __mknrfaifromdat.f__(), __randomsel.f__(), __ranseeascii.f__(), __read_show_skel.f__(), __readbin.f__(), __reanal.f__(), __reducebinsize.f__(), __reduceeachsize.bin.f__(), __reduceeachsize.binbyl.f__(), __reduceeachsize.f__(), __reducesize.f__(), __seeascii.f__(), __select.f__(), __showspec.f__(), __smashskel.f__(), __smashskelsmallfileno.f__(), __smashskeltrial.f__(), __test.f__(), __testkl3.f__(), __testsampprim.f__(), __timeslice.f__(), __vh2den.f__(), atmncobptcl(), crossc_kk(), cscalerprod(), cwhatjob(), gamfun(), gamfun_kk(), idtrafo(), ixxson(), k4ptdis(), kadbp(), kcfg(), kcfgb(), kcfp(), kcfpb(), kfrges(), kintp3s(), kmoon(), knrml(), knrml2(), kqsort(), kqsortc(), kqsortd(), kqsorth(), kqsorti(), kqsortr(), krkfehl(), krkfhstep(), krungekutfs(), ksun(), mfrgas(), minuit(), mnamin(), mncalf(), mncler(), mncntr(), mncomd(), mncont(), mncros(), mncuve(), mnderi(), mndxdi(), mneig(), mnemat(), mnerrs(), mneval(), mnexcm(), mnexin(), mnfixp(), mnfree(), mngrad(), mnhes1(), mnhess(), mnimpr(), mninex(), mninit(), mninpu(), mnintr(), mnlims(), mnline(), mnmatu(), mnmigr(), mnmnos(), mnmnot(), mnparm(), mnpars(), mnpfit(), mnpint(), mnplot(), mnpout(), mnprin(), mnpsdf(), mnrazz(), mnread(), mnrset(), mnsave(), mnscan(), mnseek(), mnset(), mnseti(), mnsimp(), mnstat(), mnstin(), mnunpt(), mnvers(), mnvert(), mnwarn(), mnwerr(), psaini(), psbint(), psborn(), pscajet(), psconf(), psdeftr(), psfaz(), psfsh(), psftild(), psgea(), psgint(), pshard(), pshot(), psjdef(), psjet(), psjet1(), psjint(), psjint1(), psnorm(), psqint(), psrec(), psrejs(), psrjint(), psshar(), pstrans(), pstrans1(), psudt(), rndcsng(), sectnu(), xxaini(), xxddfr(), xxdec2(), xxdec3(), xxdpr(), xxdtg(), xxfrag(), xxfragm(), xxgau(), xxgau1(), xxgener(), xxjetsim(), xxreg(), xxrot(), xxstr(), and xxt().

◆ t

◆ z

nodes a z ='

Definition at line 93 of file cstdatmos0.fNew2.h.

Referenced by __ccalcstdatmos.f__(), __getcoef.f__(), __test.f__(), atmncobptcl(), chookbgevent(), chookobs(), cmintime2websec(), crossc_kk(), cvthick2h(), dpmjet(), gamfun(), gamfun_kk(), gaucr(), ixxdef(), ixxson(), kadbp(), kadth(), kadthi(), kb50j2(), kcelei(), kcosd(), kctoq(), kdcmjd(), kdhtoh(), kdifva(), kdtoa(), kdtoe(), kdtoh(), kdzth2(), kdztoh(), keqtog(), ketod(), ketoh(), kgcrc(), kgcttc(), kgdted(), kgtoeq(), khtad(), khtoe(), kide(), kideb(), kj2b50(), kj2j90(), kj2tox(), kj90j2(), kjxjy(), kmjd(), kmjdst(), kmjdym(), kmobec(), kmoon(), kmoont(), koangl(), komeg0(), kpmtrx(), kqsort(), kqtoc(), krkfehl(), krkfhstep(), krkflstep(), krungekutfds(), krungekutfs(), kside0(), ksided(), ksidet(), ksind(), ksun(), ksuneq(), ktu(), kvtoa(), kxtoj2(), mfgmf0(), mflcra(), minuit(), mkafile(), mkaskelfile(), mkfiles(), mkskelfiles(), mnamin(), mnbins(), mncalf(), mncler(), mncntr(), mncomd(), mncont(), mncrck(), mncros(), mncuve(), mnderi(), mndxdi(), mneig(), mnemat(), mnerrs(), mneval(), mnexcm(), mnexin(), mnfixp(), mnfree(), mngrad(), mnhes1(), mnhess(), mnimpr(), mninex(), mninit(), mninpu(), mnintr(), mnlims(), mnline(), mnmatu(), mnmigr(), mnmnos(), mnmnot(), mnparm(), mnpars(), mnpfit(), mnpint(), mnplot(), mnpout(), mnprin(), mnpsdf(), mnrazz(), mnread(), mnrn15(), mnrset(), mnsave(), mnscan(), mnseek(), mnset(), mnseti(), mnsimp(), mnstat(), mnstin(), mntiny(), mnvers(), mnvert(), mnwarn(), mnwerr(), psaini(), psapint(), psasetc(), psbint(), psborn(), pscajet(), psconf(), pscs(), psdefrot(), psdeftr(), psdr(), psfap(), psfaz(), psfborn(), psfsh(), psftild(), psgea(), psgint(), pshard(), pshot(), psjdef(), psjet(), psjet1(), psjint(), psjint0(), psjint1(), pslam(), psnorm(), psqint(), psrec(), psrejs(), psrejv(), psrjint(), psroot(), psrotat(), psshar(), pstrans(), pstrans1(), psudint(), psuds(), psudt(), psv(), psvdef(), pszsim(), sectnu(), slice(), stand(), stndg(), stzss(), thinning(), translate_aline(), xxaini(), xxaset(), xxddfr(), xxdec2(), xxdec3(), xxdpr(), xxdtg(), xxfau(), xxfrag(), xxfragm(), xxfz(), xxgau(), xxgau1(), xxgener(), xxjetsim(), xxreg(), xxrot(), xxstr(), xxt(), and xxtwdec().