AW: USERDUMP and User written SOURCE routine

From: Chris Theis (
Date: Tue Nov 28 2006 - 01:19:24 CET

  • Next message: Alfredo Ferrari: "Re: USERDUMP and User written SOURCE routine"

    Hi Martin,
    I think that some time ago a colleague and I have already solved the problem that you are tackling at the moment. I have a bunch of routines (modified mgdraw, source + parser extensions) that allow you to flag geometry regions as particle sinks in the input file. All particles entering those regions are stopped and their type, weight, energy, position & direction is saved to a file. In a second step these files can be loaded and the particles are re-started for further transport and also the normalization with respect to the original run is taken care of. In case you're interested send me a private message and I will dig those routines up.
    Chris Theis
    CERN/SC-RP - European Organization for Nuclear Research
    1211 Geneva 23, Switzerland
    Phone: +41 22 767 8069 Office: 892-2A-015


    Von: im Auftrag von Holbourn, MP (Martin)
    Gesendet: Mo 27.11.2006 18:27
    An: Fluka-Discuss (E-mail)
    Betreff: USERDUMP and User written SOURCE routine

    My geometry is similar to that found at most storage rings - a concrete ratchet wall with a beampipe through the end wall. My job runs successfully and I do indeed get high energy photons streaming through the hole in the concrete from electrons lost upstream on a collimator. What I want to do is use these photons as the source to a separate phased run as described in the section of the manual relating to collision tapes. The manual seems to suggest that one can record all particles crossing a given boundary but how do you specify which boundary crossing you are interested in? Option USERDUMP does not seem to contain any info on how to specify a boundary so am I missing something? Should WHAT(3) of USERDUMP be = 1?

    Also how does one write subroutine SOURCE to read the collision tape? Page 416 of the manual says for case 3 (source particles) the first record contains
    and the next record
    (ILO(I), ETOT(I), WT(I), XA(I), YA(I), ZA(I), TX(I), TY(I), TZ(I), I=1,LSTACK)

    BUT P438 of the manual suggests reading the collision tape by the statement

    which is a different order of variables than written to the tape - again I might have misunderstood what is going on here.

    Can anyone point me in the right direction or perhaps share a working source.f routine that does what I am trying to do.

    Thanks for any help,


    Martin Holbourn
    Radiation Protection Adviser

    Daresbury Laboratory
    WA4 4AD

    Tel: 01925 603266
    Fax: 01925 603381

  • Next message: Alfredo Ferrari: "Re: USERDUMP and User written SOURCE routine"

    This archive was generated by hypermail 2.1.6 : Tue Nov 28 2006 - 14:50:59 CET