source subroutine

From: mir inamul haq jeri <mirinamulhaq_at_gmail.com>
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