IMPORTANT: Scoring of Neutron Yield

Dear Jing,

Thank you for finding out this bug. The problem is that the usysuw.f was
not modified to accept the new neutron library that has now 260 groups.
The parameter defining the max size of neutron groups is the NMXGRP,
which was set to 100, enough for the 72 groups for the old library.
Please modify the respective line in usysuw.f to
        PARAMETER ( NMXGRP =3D 300 )
And recompile the program by typing make
$ make
Then it should work with no problem.

Best Regards
Vasilis

-----Original Message-----
From: owner-fluka-discuss_at_mi.infn.it
[mailto:owner-fluka-discuss_at_mi.infn.it] On Behalf Of Jing Hantao
Sent: Monday, 27 October 2008 09:53
To: fluka-discuss_at_fluka.org
Subject:=20

Hello everyone,
When I calculated the neutron yield from a target, I used the following
cards,=20
*
..+....1....+....2....+....3....+....4....+....5....+....6....+....7....
+....8
USRYIELD 2399.0 NEUTRON -47.0 FE DE1 1.
neutron0
USRYIELD 3.0 1.0D-14 300.0 180.0 90.0 3.0&
*
..+....1....+....2....+....3....+....4....+....5....+....6....+....7....
+....8
USRYIELD 2399.0 NEUTRON -47.0 DE1 DE2 1.
neutron1
USRYIELD 3.0 1.0D-14 300.0 180.0 90.0 3.0&

After finishing the calculation, I compiled the usysuw.f to usysuw and
the following warning happened,=20

[root_at_localhost flutil]# $FLUPRO/flutil/lfluka usysuw.f -o usysuw
$FLUPRO =3D /root/fluka
g77 -O3 -g -mcpu=3Dpentiumpro -fexpensive-optimizations -funroll-loops
-fstrength-reduce -Wall -fno-automatic -fno-silent -m32
-ffortran-bounds-check -I/root/fluka/flukapro -v -o usysuw -Xlinker -Map
-Xlinker usysuw.map usysuw.f -L/root/fluka -lflukahp
Driving: g77 -O3 -g -mcpu=3Dpentiumpro -fexpensive-optimizations
-funroll-loops -fstrength-reduce -Wall -fno-automatic -fno-silent -m32
-ffortran-bounds-check -I/root/fluka/flukapro -v -o usysuw -Xlinker -Map
-Xlinker usysuw.map usysuw.f -L/root/fluka -lflukahp -lfrtbegin -lg2c
-lm -shared-libgcc
Reading specs from /usr/lib/gcc/i386-redhat-linux/3.4.5/specs
Configured with: ../configure --prefix=3D/usr --mandir=3D/usr/share/man
--infodir=3D/usr/share/info --enable-shared --enable-threads=3Dposix
--disable-checking --with-system-zlib --enable-__cxa_atexit
--disable-libunwind-exceptions --enable-java-awt=3Dgtk
--host=3Di386-redhat-linux
Thread model: posix
gcc version 3.4.5 20051201 (Red Hat 3.4.5-2)
 /usr/libexec/gcc/i386-redhat-linux/3.4.5/f771 usysuw.f
-mtune=3Dpentiumpro -quiet -dumpbase usysuw.f -m32 -auxbase usysuw -g =
-O3
-Wall -version -fexpensive-optimizations -funroll-loops
-fstrength-reduce -fno-automatic -fno-silent -ffortran-bounds-check
-I/root/fluka/flukapro -o /tmp/cc5GFYdf.s
GNU F77 version 3.4.5 20051201 (Red Hat 3.4.5-2) (i386-redhat-linux)
        compiled by GNU C version 3.4.5 20051201 (Red Hat 3.4.5-2).
GGC heuristics: --param ggc-min-expand=3D64 --param =
ggc-min-heapsize=3D64265
usysuw:
usysuw.f: In program `usysuw':
In file included from usysuw.f:0:
usysuw.f:184: warning: unused variable 'accumt'
usysuw.f:183: warning: unused variable 'domega'
usysuw.f:182: warning: unused variable 'costhe'
usysuw.f:183: warning: 'deltae' might be used uninitialized in this
function
usysuw.f:184: warning: 'efact' might be used uninitialized in this
function
usysuw.f:316: warning: 'kmax' might be used uninitialized in this
function
usysuw.f:465: warning: 'energe' might be used uninitialized in this
function
  lnnbln:
 as -V -Qy -o /tmp/cc7Jk7Kj.o /tmp/cc5GFYdf.s
GNU assembler version 2.15.92.0.2 (i386-redhat-linux) using BFD version
2.15.92.0.2 20040927
 /usr/libexec/gcc/i386-redhat-linux/3.4.5/collect2 --eh-frame-hdr -m
elf_i386 -dynamic-linker /lib/ld-linux.so.2 -o usysuw
/usr/lib/gcc/i386-redhat-linux/3.4.5/../../../crt1.o
/usr/lib/gcc/i386-redhat-linux/3.4.5/../../../crti.o
/usr/lib/gcc/i386-redhat-linux/3.4.5/crtbegin.o -L/root/fluka
-L/usr/lib/gcc/i386-redhat-linux/3.4.5
-L/usr/lib/gcc/i386-redhat-linux/3.4.5
-L/usr/lib/gcc/i386-redhat-linux/3.4.5/../../.. -Map usysuw.map
/tmp/cc7Jk7Kj.o -lflukahp -lfrtbegin -lg2c -lm -lgcc_s -lgcc -lc -lgcc_s
-lgcc /usr/lib/gcc/i386-redhat-linux/3.4.5/crtend.o
/usr/lib/gcc/i386-redhat-linux/3.4.5/../../../crtn.o

Then, I neglected the warning and extracted data from the output files
using the following command:

[root_at_localhost CSNS]# $FLUPRO/flutil/usysuw
 Type the input file: CSNS001_fort.47
 MC-CAD Test=20
  DATE: 10/27/ 8, TIME: 11:29:38
  100000.
 100000
 0
 1
Subscript out of range on file line 334, procedure usysuw.f/MAIN.
Attempt to access the 102-th element of variable engmax.
Aborted

I don't know why the errors happen. The version for fluka is 2008.3(I
can't afirm wether the version is 2008.3.2). but I didn't meet such
errors before by using fluka 2006.3b. Thank you very much!
                         =09

Jing Hantao
jinght_at_mail.ihep.ac.cn
2008-10-27
Received on Thu Oct 30 2008 - 12:20:23 CET