From: mir inamul haq jeri <mirinamulhaq_at_gmail.com>

Date: Tue, 10 May 2011 11:40:47 +0200

Date: Tue, 10 May 2011 11:40:47 +0200

Dear Fluka Experts

I am a beginer in Fluka and i need to create a source subroutine.the

details are as under

28 Mev protons are incident on carbon target. the neutrons are created and

the neutrons are detected(number of neutrons per steradian per Mev per

projectile) at 0 ,5 ,10 and 20 degree at a distance of say about 1cm, the

energy range of detected neutrons is from 5 Mev to 28 Mev with height of the

neutron distribution peak(as function of energy) decreasing with increasing

angle.I am doing an external fitting to this data with a polynomial of

degree 6.however i also have the data in tabulated form in a text file.

can you provide some guidance in setting up the source subroutine.like

1.whether i should input the polynomial equation in the fluka and make it

reconstruct the polynomial or i should use the tabulated data as input

without the need of inputting the polynomial.in any case what is the

detailed procedure to be followed.(i went through some examples available on

the web however i had some problems in understanding the meaning of some

statements which i am including at the end of this mail)

2.how will i invoke the angular dependence of this spectrum in the code.(the

neutron source is not a siple beam.indeed the neutrons will come out of

target in any direction with some preference to forward direction as seen

from the spectrum).

the following are the key points that i was not able to understabd in an

example provided by Mr. Alberto in respeonce to a question by Zhang Weihua

on Sat, 14 Nov 2009

* | First call initializations:

IF ( LFIRST ) THEN

* | *** The following 3 cards are mandatory ***

TKESUM = ZERZER

LFIRST = .FALSE.

LUSSRC = .TRUE.

* | *** User initialization ***

WRITE(LUNOUT,*)

WRITE(LUNOUT,'(A,132A)') ("*",I=1,132)

WRITE(LUNOUT,*)

WRITE(LUNOUT,*)

IF(NINT(WHASOU(1)) .EQ. 1) THEN

WRITE(LUNOUT,'(A)') "@ Source: Cf-252"

ELSE IF(NINT(WHASOU(1)) .EQ. 2) THEN

WRITE(LUNOUT,'(A)') "@ Source: Am-Be"

ELSE

STOP "Incorrect WHASOU(1)"

END IF

WRITE(LUNOUT,*)

WRITE(LUNOUT,*)

WRITE(LUNOUT,'(A,132A)') ("*",I=1,132)

WRITE(LUNOUT,*)

then later in the same example there is

if(kount .gt. 100) nomore = 1

Thank you

-- MIR INAMUL HAQ JERI INSTN SACLAY. UNIVERSITY OF PARIS SUD 11.Received on Tue May 10 2011 - 14:25:03 CEST

*
This archive was generated by hypermail 2.2.0
: Tue May 10 2011 - 14:25:04 CEST
*