Dear FLUKA experts,
Here I am again. :/
This solution (making a two steps simulation) was suggested at the topic "Estimating hardware configurations for bremsstrahlung simulations",
In my last e-mail, I present a big mistake trying to score a lot of values (energy, type, position, etc) using a combination of USRBDX, USRBIN and FORTRAN programming with a wrong Gaussian function.
The right solution for doing this two-steps simulation is using USERDUMP card according to some publications/discussions:
https://www.fluka.org/free_download/course/triumf2012/Workshop/Laxdal.pdf
http://www.fluka.org/web_archive/earchive/new-fluka-discuss/6806.html
http://www.fluka.org/web_archive/earchive/prova/0395.html
So, I introduced USERDUMP, edited mgdraw.f (attached) and I typed " $FLUPRO/flutil/fff $FLUPRO/usermvax/mgdraw.f" to create mgdraw.o.
It works. So, the next step here is understand why I can not generate the executable using any of the codes below (found on internet) and understand how to read and concatenate such output files. Nevertheless, the existence of this BXDRAW function left me quite optimistic.
$FLUPRO/flutil/lflukac -m fluka -C $FLUPRO/usermvax/usrini.f $FLUPRO/usermvax/usrout.f $FLUPRO/usermvax/mgdraw.f -o flukaseamu
$FLUPRO/flutil/lfluka -m fluka -o flukamy $FLUPRO/usermvax/mgdraw.o
$FLUPRO/flutil/lfluka -o mgdraw $FLUPRO/usermvax/mgdraw.f
$FLUPRO/flutil/lfluka -o $FLUPRO/MyFluka -m fluka $FLUPRO/flutil/s160302.o $FLUPRO/usermvax/mgdraw.o
Thanks again,
Marlon
Some errors showed:
/opt/fluka/libflukahp.a(mgdraw.o): In function `mgdraw_':
/home/alfredo/fluprogfor/usermvax/mgdraw.f:6: multiple definition of `mgdraw_'
/tmp/cceb4mg7.o:/opt/fluka/usermvax/mgdraw.f:6: first defined here
/opt/fluka/libflukahp.a(mgdraw.o): In function `bxdraw_':
/home/alfredo/fluprogfor/usermvax/mgdraw.f:105: multiple definition of `bxdraw_'
/tmp/cceb4mg7.o:/opt/fluka/usermvax/mgdraw.f:67: first defined here
/usr/bin/ld: cannot find -lmathlib
/usr/bin/ld: cannot find -lpawlib
/usr/bin/ld: cannot find -lgraflib
/usr/bin/ld: cannot find -lgrafX11
/usr/bin/ld: cannot find -lpacklib
/usr/bin/ld: cannot find -lkernlib
collect2: error: ld returned 1 exit status
------------------------------------------------------
gfortran: error: -m: No such file or directory
gfortran: error: fluka: No such file or directory
gfortran: error: -o: No such file or directory
gfortran: error: flukamy: No such file or directory
------------------------------------------------------
collect2: error: ld returned 1 exit status
------------------------------------------------------
/opt/fluka/libflukahp.a(mgdraw.o): In function `mgdraw_':
/home/alfredo/fluprogfor/usermvax/mgdraw.f:6: multiple definition of `mgdraw_'
/opt/fluka/usermvax/mgdraw.o:/opt/fluka/usermvax/mgdraw.f:6: first defined here
/opt/fluka/libflukahp.a(mgdraw.o): In function `bxdraw_':
/home/alfredo/fluprogfor/usermvax/mgdraw.f:105: multiple definition of `bxdraw_'
/opt/fluka/usermvax/mgdraw.o:/opt/fluka/usermvax/mgdraw.f:67: first defined here
/usr/bin/ld: i386 architecture of input file `/opt/fluka/flutil/s160302.o' is incompatible with i386:x86-64 output
/opt/fluka/flutil/s160302.o: In function `source_':
/opt/fluka/flutil/s160302.f:133: undefined reference to `s_wsle'
/opt/fluka/flutil/s160302.f:133: undefined reference to `do_lio'
/opt/fluka/flutil/s160302.f:133: undefined reference to `do_lio'
/opt/fluka/flutil/s160302.f:133: undefined reference to `e_wsle'
/opt/fluka/flutil/s160302.f:159: undefined reference to `s_rnge'
/opt/fluka/flutil/s160302.f:159: undefined reference to `s_rnge'
/opt/fluka/flutil/s160302.f:159: undefined reference to `s_rnge'
/opt/fluka/flutil/s160302.f:159: undefined reference to `s_rnge'
/opt/fluka/flutil/s160302.f:159: undefined reference to `s_rnge'
=================================================================
USERDUMP configurations: Type Dump, What Complete, Unit 96, Score Local Losses
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....
USERDUMP 100. 96. 3. Dump
-----Mensagem original-----
De: Mao, Stan [mailto:mao_at_slac.stanford.edu]
Enviada em: sexta-feira, 4 de março de 2016 15:06
Para: 'fluka-discuss_at_fluka.org' <fluka-discuss_at_fluka.org>
Cc: Marlon Saveri Silva <marlon.saveri_at_lnls.br>
Assunto: RE: Splitting FLUKA runs
Please help Marlon on his question.
Thanks
Stan
-----Original Message-----
From: Mao, Stan
Sent: Friday, March 04, 2016 10:05 AM
To: 'fluka-discuss_at_fluka.org'
Cc: 'marlon.saveri_at_lnls.br'
Subject: FW: Splitting FLUKA runs
-----Original Message-----
From: Marlon Saveri Silva [mailto:marlon.saveri_at_lnls.br]
Sent: Friday, March 04, 2016 5:27 AM
To: hvincke_at_slac.stanford.edu; Mao, Stan; Rokni, Sayed H.
Cc: Roberto Madacki; MARILIA LISBOA DE OLIVEIRA; Beamlines Common Systems
Subject: Splitting FLUKA runs
Dear Mr. Vincke, Mao and Rokni,
In December 2003 there was a very interesting publication from SLAC (SLAC-PUB-10010 - FLUKA Calculations for the Shielding Design of the SPPS Project at SLAC), in which FLUKA simulation were made in order to determine shielding.
A greatly important point is quoted in the 4th page: "the calculations were split into two runs. In a first set of runs (called RUN 1) information of particles reaching a virtual plane outside the FFTB north wall was stored on a file" and such information included: "number of total primaries, total weight of the primary praticles, energy, particle type, weight, co-ordinates [and] direction cosines".
I've been working on FLUKA simulations in order to estimate life of some equipment subject to high dose inside a beamline and, recently, aiding the Radiological Protection Group with other calculations for Sirius beamlines.
Perhaps you have seen my topic at FLUKA-discuss forum some days ago about the difficulty in getting faster simulations and how to relate time and hardware configuration. One of the indicated solution is splitting the simulation in two or more parts as you did in your work thirteen years ago.
In a perfect world, there would be a card in FLUKA (RUN1) wherewith we just draw a plane and it would give us a file containing all information to be exported as an input file for RUN2 (it would make me very happy).
I've started trying an alternative way using USRBDX card to get spectrum of photons, using USRBIN card to get the position distribution in this transversal section of the pipe and, finally, trying to generate a FORTRAN code inside source.f (based on random numbers) to create an input for RUN2 similar to those outputs of RUN1.
However, I'm having troubles programming it and the new input for RUN2 is different from RUN1 output (as you can check in the attached example). Moreover, this solution don't give me information about particle type, weight, direction cosines and I'm not sure about the best way to select energies (random vs import from a discrete table).
Therefore, if you may share, I'd like to understand how have you done such splitting.
Thanks,
Marlon Saveri Silva
Mechanical Engineer
Beamlines Instrumentation and Support Group - SIL Brazilian Synchrotron Light Laboratory- LNLS - CNPEM
+ 55 (19) 3512-2490
+ 55 (19) 99478-8097
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
- application/octet-stream attachment: mgdraw.f
Received on Fri Mar 04 2016 - 23:48:49 CET