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 )
ILOFLK (NPFLKA) = IJHION
....
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
carefully.
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