From: Sebastien WURTH (wurth@ipno.in2p3.fr)
Date: Fri Mar 17 2006 - 12:35:45 CET
Hello,
In your USERWEIGH card you apply a 4.0 value to WHAT(3).
According to the manual this calls routine FLUSCW and routine FLDSCP and
then you apply a 1.0 value to WHAT(6) that calls routine COMSCW with
which you wanna do your calculations. I think that the problem comes
from here, USERWEIGH card is read and FLUSCW is first called. Then,
maybe routine COMSCW is not called at all and all your changes in COMSCW
do not modify anything at all. This should be confirmed by a member of
the FLUKA development team, I'm not quite sure of what is really
happening here !
My advice is try it again with a 0.0 value in WHAT(3) of USERWEIGH card.
Regards.
Sebastien
Fan, Ph.D., Jiajin a écrit :
> Hi All,
>
>
>
> I want to calculate proton dose using USRBIN and weight deposited
> energy by using COMSCW.f. But I can not get results as I expected.
>
>
>
> First, according to FLUKA manual, parameter ISCRNG in COMSCW.f should
> be 1 when we calculate energy density bin (208.), but in my
> calculation it always be 2, I can not find out why it's not consist
> with manual.
>
>
>
> Second, I change COMSCW to different values but always got same output
> in USRBIN.
>
>
>
> Is there anything wrong in my input file or COMSCW.f file?
>
>
>
> Thank you.
>
>
>
>
>
> In my input file, define USERWEIG and USRBIN as:
>
> USERWEIG 0. 0. 4. 0. 0. 1.0
>
>
>
> USRBIN 10.0 208. 60.0 .5 .5
> 100. ENERGY
>
> USRBIN -.5 -.5 0.0 1. 1.
> 1000. &
>
>
>
> And the contents in COMSCW.f are:
>
>
>
> * ======== In order to compute RBE included dose ======== *
>
> if ( ij .eq. 1 ) then
>
>
>
> enerProton = etrack - 0.93827231
>
> write(*,*) iscrng
>
> if ( iscrng .eq. 1 ) then
>
> if ( enerProton .lt. 0.1) COMSCW = 12.0
>
> end if
>
>
>
> end if
>
>
>
>
>
>
>
> *******************************************
>
> Jiajin Fan, Ph.D.
>
> Research Associate, Radiation Oncology
>
> Fox Chase Cancer Center
>
> 432 Rhawn Street, Philadelphia, PA 19111
>
> Tel: (215) 214-1711
>
> Fax:(215) 728-1148
>
> Email: jiajin.fan@fccc.edu
>
> *******************************************
>
>
>
>------------------------------------------------------------------------
>
>******** McAfee GroupShield for Microsoft Exchange **********
>**********************************************************************
>
>Alert generated on: Thursday, March 16, 2006 06:14:08 PM Eastern Standard Time
>
>
>The file weighted.inp has been replaced. Reason: The file met the blocking options set in the anti-virus system.
>Please consult your administrator for further help quoting your ticket number: OA97_1142550848_EXCHSERVER_3
>
This archive was generated by hypermail 2.1.6 : Fri Mar 17 2006 - 13:09:27 CET