USERDUMP and User written SOURCE routine

From: Holbourn, MP \(Martin\) (m.p.holbourn@dl.ac.uk)
Date: Mon Nov 27 2006 - 18:27:17 CET

  • Next message: Francesco Cerutti: "Re: your mail"

    Experts,
    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
    NCASE, LSTACK, LSTMAX, TKESUM, WEIPRI
    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
    READ(21,*) IPART,X ,Y ,Z, COSX, COSY, COSZ, ENERGY, WEIGHT

    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,

    Regards

    Martin Holbourn
    Radiation Protection Adviser
     
    Daresbury Laboratory
    Daresbury
    Warrington
    Cheshire
    WA4 4AD
     
    Tel: 01925 603266
    Fax: 01925 603381
    mailto:m.p.holbourn@dl.ac.uk

  • Next message: Francesco Cerutti: "Re: your mail"

    This archive was generated by hypermail 2.1.6 : Mon Nov 27 2006 - 20:41:52 CET