Dear Yang Tao,
I refer to the attachments in your original post
http://www.fluka.org/web_archive/earchive/new-fluka-discuss/8373.html
Please check them...they are as I described.
> Dear Paola,
>
> Thanks for your reply. I donnot exactly understand your suggestion. I the
> attachment "comscw.f", I haven't commented the default assignment
> comscw=oneone, it is just under the <include> statements as shown below,
> do you mean comscw=oneone should be placed in another place? And there is
> only one USRBIN card in the dose.inp.
>
> INCLUDE '(FLKMAT)'
> INCLUDE '(SCOHLP)'
> INCLUDE '(SOUEVT)'
> INCLUDE '(TRACKR)'
>
> LSCZER = .FALSE.
> COMSCW = ONEONE
>
> IF ( ISCRNG .EQ. 1 .AND. JSCRNG .EQ. 1 ) THEN
> COMSCW = 1.6 D-07 / RHO (MEDFLK(4,1))
> END IF
> RETURN
> END
>
> regards,
> Tao Yang
>
>> -----原始邮件-----
>> 发件人: "Paola Sala" <paola.sala_at_mi.infn.it>
>> 发送时间: 2015年9月3日 星期四
>> 收件人: "YANG Tao" <yangt_at_ihep.ac.cn>
>> 抄送: fluka-discuss_at_fluka.org
>> 主题: Re: [fluka-discuss]: PROBLEM about comscw.f to obtain absorption
>> dose !!!
>>
>> 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
>>
>
>
> --
> 杨涛
> 中科院高能所东莞分部（东莞中子科学中心）加速器技术部
> 地址：东莞市大朗镇中子源路1号中国散裂中子源A2栋606室
> E-mail：yangt_at_ihep.ac.cn
> 电话：0769-38944239
> 邮编：523803
>
>
>
> __________________________________________________________________________
> You can manage unsubscription from this mailing list at
> https://www.fluka.org/fluka.php?id琧_info
>
>
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 Thu Sep 03 2015 - 17:55:26 CEST