FLUKA: problem with dose scoring, comscw
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
Partial thread listing: