Dear Carmona,
This is due to the fact that there is now an additional output variable
(IDCDAU, a one-dimensional array of 4 integers) in the subroutine. Here
is the relevant part of the header for completeness:
DOUBLE PRECISION FUNCTION TIM1O2 ( IA, IZ, KA, KZ, T12DAU,
& BRDAUG, IDCDAU, ISOMER )
* *
* 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 *
* 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) *
* Idcdau(i) = index of the decay channels *
* 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 *
Hope this helps,
Cesc
On Thu, Jun 07 2018, at 15:19 +0200, JUAN CARLOS ARMENTEROS CARMONA wrote:
>
>
>Dear FLUKA users,
>
>These days, I’m facing some issues because the function TIM1O2 is not
>working properly in FLUKA.2x.2 version. It was previously used in
>version FLUKA.2c.6 to discard nuclides, filtering the half-decay time
>in a personalized routine, but now it crashes. Please, could you give
>me any advice of how to deal with this problem?
>
>I attach the example I run and the output from the debugger. Hope that
>anybody can help.
>
>Best regards,
>
>Carmona
>
>---------------------------------------------------------------------------------------
>
>info locals
> __g77_cilist_1 = ( 0, 6, 0, 0x0, 0 )
> lfluka = 1
> lgbias = 0
> lgbana = 1
> lflgeo = 1
> loflts = 0
> lusrin = 0
> lusrgl = 0
> lnmgeo = 1
> lnminp = 1
> lfrfmt = 0
> ldmpcr = 0
> lldflr = 1
> ldpbdy = 0
> kflgeo = 0
> kfldnr = 0
> crvrfl = 'FLUKA011'
> ka = (0, 0, 0, 0)
> kz = (0, 0, 0, 0)
> brdaug = (0, 0, 0, 0)
> t12dau = (0, 0, 0, 0)
> thalf = 0
>x/f thalf
> 0x0: Cannot access memory at address 0x0
>crash message in lines 6, 65, 66, 67, 68, 69, and 73
> tim1o2_ (ia=0x9eeb008 <iares.7>, iz=0x9eeb00c <izres.8>,
>ka=0x97abb20 <ka.0>,
> kz=0x97abb30 <kz.1>, t12dau=0x97abb60 <t12dau.3>,
> brdaug=0x97abb40 <brdaug.2>, idcdau=0x9eeb010 <iisres.9>,
>isomer=0x0)
> at tim1o2.f:73
>* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
>GLOBAL 1.
>DEFAULTS PRECISIO
>GEOBEGIN COMBNAME
> 0 0
>SPH blkbody 0. 0. 0. 100000.
>SPH void 0. 0. 0. 10000.
>RCC a 0. 0. 0. 0. 0. 1. 5.
>RCC b 0. 0. 2. 0. 0. 3. 5.
>TRC c 0. 0. 10. 0. 0. -4. 5. 2.
>YZP center 0.
>END
>BLKBODY 5 +blkbody -void
>VOID 5 +void -a -b -c
>a1 5 +a +center
>a2 5 +a -center
>b1 5 +b +center
>b2 5 +b -center
>c1 5 +c +center
>c2 5 +c -center
>END
>GEOEND
>* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
>ASSIGNMA BLCKHOLE BLKBODY
>ASSIGNMA VACUUM VOID
>ASSIGNMA GOLD a1 a2
>ASSIGNMA GOLD b1 b2
>ASSIGNMA GOLD c1 c2
>* Set Peanut for All Energies
>PHYSICS 20000. 20000. 20000. 20000. 20000. 20000.PEATHRES
>* Maximum momentum for DPMJET
>PHYSICS 20000.0 LIMITS
>* For activation studies
>PHYSICS 1. COALESCE
>PHYSICS 3. EVAPORAT
>IONTRANS HEAVYION
>* For Sesame, kill prompt shower (no replicas in first step)
>RADDECAY 1. 1. 1.0 99999
>#if 0
>* Open a new unit to save a summary file
>OPEN 99. NEW
>tsresnuc_prompt_summary
>#endif
>#if 0
>OPEN 98. NEW
>tsresnuc_prompt_nuclides
>#endif
>* weight sccording to personalized USRRNC routine
>USERWEIG 1.
>*USROCALL
>RANDOMIZ 1. 3679.
>START 1.E1
>STOP
>*$ CREATE USRRNC.FOR
>*COPY USRRNC
>*
>*=== Usrrnc ===========================================================*
>*
> SUBROUTINE USRRNC ( IZ, IA, IS, X, Y, Z, MREG, WEE, ICALL )
>
> INCLUDE '(DBLPRC)'
> INCLUDE '(DIMPAR)'
> INCLUDE '(IOUNIT)'
>*
>*----------------------------------------------------------------------*
>* *
>* Copyright (C) 2005-2005 by Alfredo Ferrari & Paola Sala *
>* All Rights Reserved. *
>* *
>* *
>* USeR Residual NuClei: *
>* *
>* Created on 06 april 2005 by Alfredo Ferrari & Paola Sala *
>* Infn - Milan *
>* *
>* Last change on 06-apr-05 by Alfredo Ferrari *
>* *
>* *
>*----------------------------------------------------------------------*
>*
>* Needed to check the half-life of nuclides
> EXTERNAL TIM1O2
> DIMENSION KA (4), KZ (4), BRDAUG (4), T12DAU (4)
> INTEGER IZ, IA, IS, MREG
> DOUBLE PRECISION X, Y, Z, WEE, THALF, TIM1O2
>
> WRITE(*,*) 'GDB$0:usrrnc ','Z=',IZ,', A=',IA,', S=',IS,', x=',X,
> & ', y=',Y,', z=',Z,', reg=',MREG,'vwee=',WEE,', icall=',ICALL
> THALF = TIM1O2(IA, IZ, KA, KZ, T12DAU, BRDAUG, IS)
> WRITE(*,*) 'GDB$0:usrrnc ','THALF=',THALF
>
> RETURN
>*=== End of subroutine Usrrnc =========================================*
> END
--
Francesc Salvat Pujol
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
Tel: +41 22 76 64011
Fax: +41 22 76 69474
__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Fri Jun 08 2018 - 11:11:10 CEST