Re: SOURCE composed of several heavy ions

From: Giuseppe Battistoni <>
Date: Sat, 14 Mar 2009 08:21:15 +0100 (CET)

Dear Sharmalee
from the technical point of you, one you have sample which
which is your projectile nucleus, you will define in two integer variables
of you choice Z and A. For instance IZPRI and IAPRI
Then, after setting IJBEAM = -2 you must also set
        IPROZ = IZPRI
        IPROA = IAPRI
Then you can see already in the templase source.f that at some place
there are the following lines which manage the situation:
       ELSE IF ( IJBEAM .EQ. -2 ) THEN
         IJHION = IPROZ * 1000 + IPROA
         IJHION = IJHION * 100 + KXHEAV
         IONID = IJHION
         CALL DCDION ( IONID )
         CALL SETION ( IONID )
For protons you just put IJBEAM=1

As far as physics is concerned, the real attention in these cosmic ray
problems is to remind that the sampling of projectile species is
correlated to energy sampling and of course you need to treat this

        Giuseppe Battistoni

On Thu, 12 Mar 2009, Sharmalee Randeniya wrote:

> Dear Experts,
> I am working on a problem on Galactic cosmic rays. The source particles are
> composed of about 28 heavy isotopes ranging from proton to Ni(A=58, Z=28), and
> each has its own fluence, dphi/dE, distribution.
> I thought of using the source.f routine define a source containing all
> isotopes. I was able to read in the fluence spectra and cumulative fluence
> spectra for each. But I don't know how to define and choose between the
> isotopes during the simulation. I see that if I set IJBEAM=-2, means a heavy
> ion. I also see that IJHION refers to heavy ion type. But I am not sure how to
> proceed further. please advice me.
> An alternative would be to run 28 simulations, one for each isotope using BEAM
> and HI-PROPE cards in inp file, and post process data. if possible, I would
> like to avoid it.
> Thank you !!!
> Sharmalee Randeniya
> Rice University
Received on Sat Mar 14 2009 - 09:00:49 CET

This archive was generated by hypermail 2.2.0 : Sat Mar 14 2009 - 09:00:49 CET