FLUKA: problem with dose scoring, comscw


To fluka-discuss@listbox.cern.ch
From Lev Chekhtman <Lev.Chekhtman@cern.ch>
Date Tue, 23 Oct 2001 06:59:22 +0200 (CEST)
Reply-To Lev Chekhtman <Lev.Chekhtman@cern.ch>
Sender owner-fluka-discuss@listbox.cern.ch

Hello, I am trying to use comscw to have proper dose scoring. I changed
comscw according to the comments (see attached), compiled and linked
new executable but got no changes. Please comment. Thanks, Lev.
 

*$ CREATE COMSCW.FOR
*COPY COMSCW
*
*=== comscw  ===========================================================*
*
      DOUBLE PRECISION FUNCTION COMSCW(IJ,XA,YA,ZA,MREG,RULL,LLO,ICALL)

      INCLUDE '(DBLPRC)'
      INCLUDE '(DIMPAR)'
      INCLUDE '(IOUNIT)'
*
*----------------------------------------------------------------------*
*                                                                      *
*     Input variables:                                                 *
*                                                                      *
*           Ij = (generalized) particle code                           *
*     Xa,Ya,Za = position                                              *
*         Mreg = region number                                         *
*         Rull = amount to be deposited                                *
*          Llo = particle generation                                   *
*        Icall = call id                                               *
*                                                                      *
*     Output variables:                                                *
*                                                                      *
*       Comscw = factor the scored amount will be multiplied by        *
*       Lsczer = logical flag, if true no amount will be scored        *
*                regardless of Comscw                                  *
*                                                                      *
*     Useful variables (common SCOHLP):                                *
*                                                                      *
*     Energy/Star binnings/scorings (Comscw):                          *
*          ISCRNG = 1 --> Energy density  binning                      *
*          ISCRNG = 2 --> Star   density  binning                      *
*          ISCRNG = 3 --> Residual nuclei scoring                      *
*          JSCRNG = # of the binning                                   *
*                                                                      *
*     Useful variables (common SOUEVT):                                *
*                                                                      *
*          X,Y,Zsoevt(i) = position    of the i_th source particle     *
*          TX,Y,Zsoev(i) = direction   of the i_th source particle     *
*              Wtsoev(i) = weight      of the i_th source particle     *
*              Pmsoev(i) = momentum    of the i_th source particle     *
*              Tksoev(i) = kin. energy of the i_th source particle     *
*              Agsoev(i) = age         of the i_th source particle     *
*              Aksoev(i) = Kaon ampl.  of the i_th source particle     *
*              Ussoev(i) = user var.   of the i_th source particle     *
*              Ijsoev(i) = identity    of the i_th source particle     *
*              Nrsoev(i) = region      of the i_th source particle     *
*              Nlsoev(i) = lattice     of the i_th source particle     *
*                Npsoev  = number of the source particles              *
*                                                                      *
*                                                                      *
*----------------------------------------------------------------------*
*
      INCLUDE '(MAPA)'
      INCLUDE '(SCOHLP)'
      INCLUDE '(SOUEVT)'
*
* ======== In order to compute doses ========= *
*     (Medium(n) is the material number of region n
*      Rho(m) is the density of material m)
      COMSCW = ONEONE
      IF ( ISCRNG .EQ. 1 )  COMSCW = ONEONE / RHO (MEDIUM(MREG))
        LSCZER=.FALSE.
      RETURN
      END



Your name :
Your email :
Subject :
Body :
 

Partial thread listing: