Hello
in your comscw.f routine you commented out the default assignment
comscw=oneone
BUT, comscw is called for ALL binnings, therefore a meaningful value must
be provided for all binnings. This is why the defaut comscw=oneone is
there.
By commenting this out, you will get for binning 2 ta potentially
ininitialized value, indeed you get he same value that was calculated for
binning 1, so the same final result in the two binnings.
Uncommenting the default everything works.
Please note that you can also use the DOSE generalized particle to get
automatically results in GeV/g
Regards
Paola
> Hi, everyone,
> I use the USERWEIG and comscw.f to score the dose(Gy) according to the
> USRBIN card which obtains the energy deposition density. But nothing has
> changed, it seems comscw.f has no influence to the results. Attachment is
> the input file and the user routine comscw.f. Below is the USERWEIG
> setting and comscw.f input. Could anyone give me some suggestions?
>
>
> USRBIN 11. ENERGY -21. 13.0 0.
> 70.endep
> USRBIN 11.0 0.0 -50. 5.0 10.0 120. &
> USERWEIG 0.0 0.0 1.
>
>
> ********comscw.f
> INCLUDE '(FLKMAT)'
> INCLUDE '(SCOHLP)'
> INCLUDE '(SOUEVT)'
> INCLUDE '(TRACKR)'
> *
> LSCZER = .FALSE.
> COMSCW = ONEONE
> * ======== In order to compute doses ========= *
> * (Medflk(n,iprodc) is the material number of region n
> * Rho(m) is the density of material m)
> IF ( ISCRNG .EQ. 1 .AND. JSCRNG .EQ. 1 ) THEN
> COMSCW = 1.6 D-07 / RHO (MEDFLK(4,1))
> END IF
> RETURN
> *=== End of function Comscw ===========================================*
> END
>
>
>
>
>
Paola Sala
INFN Milano
tel. Milano +39-0250317374
tel. CERN +41-227679148
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Wed Sep 02 2015 - 20:48:31 CEST