Re: source subroutine

From: Paola Sala <paola.sala_at_mi.infn.it>
Date: Thu, 19 May 2011 08:58:48 +0200 (CEST)

hello
the attachment cannot be read, it went corrupted
 From the error messages, it seems that tere are formatting errors in the
source.f code . Did you put 7 spaces before starting each line? BPlease
keep the same formatting as in the example source.f. I would also suggest
to leave the fff script in the $FLUPRO/flutil directory,
If you cannot solve your ptoblem, please try to send the source code again
Paola

> --001517440fd0b1e7f304a362deea
> Content-Type: multipart/alternative; boundary=001517440fd0b1e7ea04a362dee8
>
> --001517440fd0b1e7ea04a362dee8
> Content-Type: text/plain; charset=ISO-8859-1
>
> Hello,
> Thank you for your previous responces because of which i have been able to
> atleast write a source subroutine(without angular distrribution for the
> time
> being), however when i am trying to compile it i am getting following
> error,
>
> *inam_at_inam:~$ cd Desktop/test
> inam_at_inam:~/Desktop/test$ ls
> fff SOURCE.f test1.flair test1.inp test1.inp~
> inam_at_inam:~/Desktop/test$ ./fff SOURCE.f
> -e g77 -o SOURCE.o -c -O3 -g -mcpu=pentiumpro -fexpensive-optimizations
> -funroll-loops -fstrength-reduce -Wall -fno-automatic -fno-silent -m32
> -ffortran-bounds-check -I/home/inam/FLUKA_CODE2011/flukapro SOURCE.f
> `-mcpu=' is deprecated. Use `-mtune=' or '-march=' instead.
> SOURCE.f:1:
> SUBROUTINE SOURCE (NOMORE)
> ^
> Invalid first character at (^) [info -f g77 M LEX]
> SOURCE.f:2:
> INCLUDE '(DBLPRC)' !CONTAINING PHYSICAL CONSTANTS ,THEIR VALUES AND
> DEFINITION FOR DOUBLE AND SINGLE PRECISION.
> ^*
> Invalid first character at (^) [info -f g77 M LEX]
> SOURCE.f:3:
> INCLUDE '(IOUNIT)'
> ^
> Invalid first character at (^) [info -f g77 M LEX]
> SOURCE.f:4:
> INCLUDE '(DIMPAR)'
> ^
> Invalid first character at (^) [info -f g77 M LEX]
> SOURCE.f:5:
> INCLUDE '(BEAMCM)' !CONTAINING BEAM INFORMATION LIKE MOMENTUM, ENERGY,
> SPECTRUM SHAPE, RADIUS ETC
> ^
> Invalid first character at (^) [info -f g77 M LEX]
> SOURCE.f:6:
> INCLUDE '(FHEAVY)'
> ^
> Invalid first character at (^) [info -f g77 M LEX]
> SOURCE.f:7:
> INCLUDE '(FLKSTK)'
> ^
> Invalid first character at (^) [info -f g77 M LEX]
> SOURCE.f:8:
> INCLUDE '(SOURCM)'
> ^
> Invalid first character at (^) [info -f g77 M LEX]
> SOURCE.f:9:
> INCLUDE '(LTCLCM)'
> ^
> Invalid first character at (^) [info -f g77 M LEX]
> SOURCE.f:10:
> INCLUDE '(PAPROP)'
> ^
> Invalid first character at (^) [info -f g77 M LEX]
> SOURCE.f:11:
> INCLUDE '(SUMCOU)'
> ^
> Invalid first character at (^) [info -f g77 M LEX]
> SOURCE.f:12:
> INCLUDE '(IOIOCM)'
> ^
> Invalid first character at (^) [info -f g77 M LEX]
> SOURCE.f:13:
> LOGICAL LFIRST
> ^
> Invalid first character at (^) [info -f g77 M LEX]
> SOURCE.f:14:
> PARAMETER (NDP=24) !NUMBER OF DATA POINTS=24+1
> ^
> Invalid first character at (^) [info -f g77 M LEX]
> SOURCE.f:17:
> SAVE LFIRST
> ^
> Invalid first character at (^) [info -f g77 M LEX]
> SOURCE.f:20:
> & 5.0235D-3, 5.9617D-3, 7.1081E-3, 8.0485D-3,
> ^
> Continuation indicator at (^) invalid on first non-comment line of file or
> following END or INCLUDE [info -f g77 M LEX]
> SOURCE.f:21:
> & 9.1061D-3, 1.0169D-2, 1.1016D-2, 1.2070D-2,
> ^
> Continuation indicator at (^) invalid on first non-comment line of file or
> following END or INCLUDE [info -f g77 M LEX]
> SOURCE.f:22:
> & 1.3140D-2, 1.4133D-2, 1.5118D-2, 1.6082D-2,
> ^
> Continuation indicator at (^) invalid on first non-comment line of file or
> following END or INCLUDE [info -f g77 M LEX]
> SOURCE.f:23:
> & 1.7154D-2, 1.8130D-2, 1.9105D-2, 2.0069D-2,
> ^
> Continuation indicator at (^) invalid on first non-comment line of file or
> following END or INCLUDE [info -f g77 M LEX]
> SOURCE.f:24:
> & 2.1136D-2, 2.1989D-2, 2.3057D-2, 2.4128D-2,
> ^
> Continuation indicator at (^) invalid on first non-comment line of file or
> following END or INCLUDE [info -f g77 M LEX]
> SOURCE.f:25:
> & 2.5094D-2, 2.6056D-2, 2.7017D-2, 2.7976D-2 /
> ^
> Continuation indicator at (^) invalid on first non-comment line of file or
> following END or INCLUDE [info -f g77 M LEX]
> SOURCE.f:28:
> & 0.0419133236D0, 0.0863376028D0, 0.1532320378D0, 0.2179548878D0,
> ^
> Continuation indicator at (^) invalid on first non-comment line of file or
> following END or INCLUDE [info -f g77 M LEX]
> SOURCE.f:29:
> & 0.2970625522D0, 0.3785352978D0, 0.4445988641D0, 0.5304797924D0,
> ^
> Continuation indicator at (^) invalid on first non-comment line of file or
> following END or INCLUDE [info -f g77 M LEX]
> SOURCE.f:30:
> & 0.6185002645D0, 0.6896345035D0, 0.7443335649D0, 0.7893005009D0,
> ^
> and the same thing for the whole programme file.
>
> *please find attached the programme file also*
> cuold you please let me know the problem
> Thanku
>
> On Tue, May 10, 2011 at 11:40 AM, mir inamul haq jeri<
> mirinamulhaq_at_gmail.com> wrote:
>
>> 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.
>>
>>
>
>
> --
> MIR INAMUL HAQ JERI
> INSTN SACLAY.
> UNIVERSITY OF PARIS SUD 11.
>

Paola Sala
INFN Milano
tel. Milano +39-0250317374
tel. CERN +41-227679148
Received on Thu May 19 2011 - 09:53:19 CEST

This archive was generated by hypermail 2.2.0 : Thu May 19 2011 - 09:53:20 CEST