Re: [fluka-discuss]: Checking if isotope is stable in usrrnc.f

From: Alfredo Ferrari <alfredo.ferrari_at_cern.ch>
Date: Fri, 14 Mar 2014 15:55:57 +0100

Dear Moritz

yes you can call the TIM1O2 function which returns an isotope half-life.
A half-life of 1E38 s means "stable". How to call the function and
the meaning of the various arguments are explained below.

                Alfredo

       DOUBLE PRECISION FUNCTION TIM1O2 ( IA, IZ, KA, KZ, T12DAU,
      & BRDAUG, ISOMER )

....

* *
* TIMe 1/2: *
* *

....

* Input variables: *
* *
* Ia = Mass number of the isotope *
* Iz = Atomic number of the isotope *
* Isomer = flag for isomers: -1 = maybe there is an isomer, *
* please check *
* 0 = no isomer or anyway no isomer *
* to be considered *
* +n = this is an isomer (the n_th of *
* the given A,Z) *
* *
* Output variables: *
* *
* Tim1o2 = half-life for the considered transition (s) *
* Ka(i) = Mass number of the decay product(s)(-1 -> no decay) *
* If the decay product is the n_th isomer of Ka, Kz *
* n x 1000 is added to the mass number *
* Kz(i) = Atomic number of the decay product(s)(-1 -> no decay) *
* T12dau(i) = half-life of the daughter isotope(s) *
* Brdaug(i) = branching of the daughter isotope(s) *
* Isomer = flag for isomers: -n = there are n isomers, present *
* output refers to the ground *
* state *
* 0 = no isomer requested or found *
* to be considered *
* +n = the output refers to the nth *
* isomer of the given A,Z *

....

       DIMENSION KA (4), KZ (4), BRDAUG (4), T12DAU (4)


+----------------------------------------------------------------------+
| Alfredo Ferrari || Tel.: +41.22.76.76119 |
| CERN-EN/STI || Fax.: +41.22.76.69474 |
| 1211 Geneva 23 || e-mail: Alfredo.Ferrari_at_cern.ch |
| Switzerland || |
+----------------------------------------------------------------------+

On Fri, 14 Mar 2014, Moritz Guthoff wrote:

> Dear FLUKA experts,
> I am using the usrrnc.f user routine to write out all residual nuclei
> created in my simulation into a file. This file gets quite big and I am
> looking to reduce the size. Since I am only interested in radioactive
> isotopes, I would like to omit all stable isotopes.
>
> My question: Is the a built in function that would help me to check wether
> an isotope (identified by the Z and A parameters) is stable or not?
>
> Thanks a lot,
> Moritz
>
>
Received on Fri Mar 14 2014 - 16:50:13 CET