c
c Routine directly callable from C/C++ (it could be necessary to
c append an underscore, "_", to its name depending on the C/C++
c compiler used).
c
c
subroutine ghfpars(nmax, xmax, x0, lambda, vrb, irc)
c
c Setting the internal quantities needed to work with the
c Gaisser-Hillas function related routines.
c The Gaisser-Hillas function for the profile of charged particles
c (T. K. Gaisser and A. M. Hillas, in Proc. 15th ICRC (Plovdiv),
c vol. 8, p. 353 (1977)) is given by:
c
c / X - X0 \ [(Xmax-X0)/lambda]
c Ngh(X) = Nmax | ----------- | *
c \ Xmax - X0 /
c
c / Xmax - X \
c * exp | ---------- |
c \ lambda /
c
c for X > X0 (Ngh(X) = 0 for X < X0).
c
c Written by: S. J. Sciutto, La Plata 1999, 2000.
c
c
c Arguments:
c =========
c
c nmax, xmax, x0,
c lambda.......... (input, double precision) The four parameters of
c the Gaisser-Hillas function. Xmax, X0, and
c lambda are expressed in g/cm2.
c vrb............. (input, integer) Verbosity control. If vrb is
c zero or negative then no error/informative
c messages are printed; error conditions are
c communicated to the calling program via the
c return code. If vrb is positive error messages
c will be printed: vrb = 1 means that messages
c will be printed even with successful operations.
c vrb = 2,3 means that only error messages will
c be printed. vrb > 3 is similar to vrb = 3, but
c with the additional action of stopping the
c program if a fatal error takes place.
c irc............. (output, integer) Return code. 0 means
c successful return.
c
c
c
c<--->
|
These pages are maintained by Sergio J. Sciutto |
Back to AIRES library index |
| Back to AIRES home |