Re: Problems in Fluka first running

From: Alberto Fasso' (fasso@SLAC.Stanford.EDU)
Date: Wed Sep 13 2006 - 23:51:31 CEST

  • Next message: Alfredo Ferrari: "Release 2006.3: reminder for those who were using a pre-release"

    Dear Ruiguang Wang,

    your run never stops because you have asked for too many particles!
    When you test a new input, try it first with only a few particles
    to see how long it takes, and then run it again with a number of particles
    for which a run will be completed in a reasonable time.
    I have run your problem for 1000 particles without any problem.
    It has taken 405 seconds.
    You had asked for 1 million particles: on my machine the run would
    take 405000 sec, i.e. more than 4 days! Of course you have the impression
    that the run never stops!

    By the way, the previous run aborted because there were 3 errors:

    - CALSIUM instead of CALCIUM (already pointed out by G. Battistoni)

    - LOW-MAT 7.0 8.0 -2.0 293.0 NITROGEN
    Should have been
      LOW-MAT 7.0 7.0 -2.0 293.0 NITROGEN

    - LOW-MAT 20.0 8.0 -2.0 293.0 ARGON
    Should have been
      LOW-MAT 20.0 18.0 -2.0 293.0 ARGON

    Check the low-energy neutron material identifiers in the manual!

    In any case, as Battistoni again correctly pointed out, command
    LOW-MAT should be used ONLY IN A VERY FEW RARE CASES. In most normal cases,
    avoiding it saves you from many occasions of error.
    This is written very clearly in the manual. I report it here below because
    many users seem to love this command (I have never understood why).

       Alberto

    >From the manual:

       Option LOW-MAT should be avoided if it is not really necessary
       (experience has shown that it is often misinterpreted by beginner
       users). The option is NOT REQUIRED if the following 3 conditions
       are all true:
       1) the low-energy neutron material desired is unique or is
          listed before any other material with the same name in
          list 13}
                                and
       2) that name is the same as one in the FLUKA list 6} or as
          given by a MATERIAL option
                                and
       3) there is only one FLUKA material associated with that
          low-energy neutron material

       On the other hand, the option IS REQUIRED in any one of the
       following cases:

       1) there is more than one low-energy neutron material with
          that name (this can happen because of data sets coming
          from different sources, or corresponding to different
          neutron temperatures, or concerning different isotopes, or
          weighted on different spectra, etc), and the one desired
          is not coming first in the list. In this case it is
          sufficient to provide just as many identifiers as required
          to remove ambiguity
                                or
       2) The FLUKA name is different from the name of the
          low-energy neutron material
                                or
       3) There is more than one FLUKA material associated with the
          given low-energy neutron material. This can happen for
          instance when the same material is present with different
          densities in different regions. In reality this is a
          special case of 2) above, since to distinguish between the
          different densities, different names must be used and one
          at least will not be equal to the name of the
          low-energy neutron material.

    On Wed, 13 Sep 2006, Wang RuiGuang wrote:

    > Dear Fluka users and authors,
    >
    > Thanks Dr. Giuseppe Battistoni for test for me.
    > According to his/her advice, I comment the three LOW-MAT lines of NITROGEN,
    > ARGON and CALCIUM.
    > The running is never stop. I don't know why. Could any person please test and
    > check it for me?
    >
    > Thanks and best regards,
    > R.G. Wang
    >
    > ---------- Forwarded Message -----------
    > From: Giuseppe Battistoni <Giuseppe.Battistoni@mi.infn.it>
    > To: Wang Rui-Guang <wangrg@mail.ihep.ac.cn>
    > Cc: support <support@fluka.org>, physics <physics@fluka.org>
    > Sent: Mon, 11 Sep 2006 08:46:34 +0200 (CEST)
    > Subject: Re: Problems in Fluka first running
    >
    > Dear Ruiguang Wang,
    > The mistake is in the use of LOW-MAT cards for NITROGEN, ARGON and CALCIUM
    > (by the way you should write CALCIUM and not CALSIUM)
    > You do not need to put the LOW-MAT card for these elements (and for many
    > others: please read carefully the manual for the LOW-MAT use): just take
    > them away, or comment them. Everything is working as far as the rest is
    > concerned (I tested it).
    > Next time please send your requests to fluka-discuss@fluka.org
    > Good luck
    > Giuseppe Battistoni for the FLUKA team
    >
    > On Mon, 11 Sep 2006, Wang Rui-Guang wrote:
    >
    > > Date: Mon, 11 Sep 2006 10:26:11 +0800
    > > From: Wang Rui-Guang <wangrg@mail.ihep.ac.cn>
    > > To: support <support@fluka.org>, physics <physics@fluka.org>
    > > Cc: wangrg <wangrg@mail.ihep.ac.cn>
    > > Subject: Problems in Fluka first running
    > >
    > > Dear Fluka users and authors,
    > >
    > > I am a new Fluka user. About one month ago I installed the current version of
    > > Fluka in Linux.
    > > I can normally run fluka with example input file. The simulation with Fluka is
    > > neutron trasportation and production of underground muons in our detector.
    > > Based on example.inp and a similar input file of neutron detection provided by
    > > a enthusiastic friend, I wrote a test input file (named test1.inp). However,
    > > it does,t work.
    > >
    > > The information is as bellow.
    > >
    > > When run fluka, we can see " No rantest1002 generated!, No test1001.err
    > > generated! ......" on the terminal lines of the first cycle, and see " No
    > > rantest1003 generated!, ......" on the terminal lines of the second cycle,
    > > the same for other cycles.
    > >
    > > The test1.inp file is
    > > "
    > > *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    > > TITLE
    > > Muon-induced neutrons in underground rock
    > > DEFAULTS NEW-DEFA
    > > *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    > > BEAM 100.E+00 MUON-
    > > BEAMPOS 0.0 0.0 -1000.0
    > > *
    > > GEOBEGIN COMBINAT
    > > Detector in Daya Bay exp.
    > > RPP 1 -10000.0 +10000.0 -10000.0 +10000.0 -10000.0 +10000.0
    > > RPP 2 -1500.0 +1500.0 -1500.0 +1500.0 -1000.0 0.0
    > > RPP 3 -1000.0 +1000.0 -1000.0 +1000.0 0.0 +1500.0
    > > RPP 4 -500.0 +500.0 -500.0 +500.0 +1500.0 +2400.0
    > > RCC 5 0.0 0.0 +1700.0 0.0 0.0 +500.0
    > > 250.0
    > > RCC 6 0.0 0.0 +1740.0 0.0 0.0 +420.0
    > > 210.0
    > > RCC 7 0.0 0.0 +1780.0 0.0 0.0 +340.0
    > > 100.0
    > > END
    > > * black hole
    > > 1 5 1 -2
    > > * rock over the exp. hall
    > > 2 5 2 -3
    > > * air in the exp. hall
    > > 3 5 3 -4
    > > * water of the buffer
    > > 4 5 4 -5
    > > * oil
    > > 5 5 5 -6
    > > * liquid scintillator
    > > 6 5 6 -7
    > > * LS+Gd
    > > 7 5 7 -8
    > > END
    > > GEOEND
    > > *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    > > MATERIAL 1.0 1.0079 0.0000899 3.0 HYDROGEN
    > > LOW-MAT 3.0 1.0 -3.0 293.0 HYDROGEN
    > > MATERIAL 6.0 12.01 2.25 6.0 CARBON
    > > LOW-MAT 6.0 6.0 -2.0 293.0 CARBON
    > > MATERIAL 7.0 14.0067 0.00117 7.0 NITROGEN
    > > LOW-MAT 7.0 8.0 -2.0 293.0 NITROGEN
    > > MATERIAL 8.0 15.9994 0.00133 8.0 OXYGEN
    > > LOW-MAT 8.0 8.0 +16.0 293.0 OXYGEN
    > > MATERIAL 11.0 22.989770 0.971 19.0 SODIUM
    > > LOW-MAT 19.0 11.0 -2.0 293.0 SODIUM
    > > MATERIAL 12.0 24.3050 1.740 9.0 MAGNESIU
    > > LOW-MAT 9.0 12.0 -2.0 293.0 MAGNESIU
    > > MATERIAL 13.0 26.981538 2.699 10.0 ALUMINUM
    > > LOW-MAT 10.0 13.0 -2.0 293.0 ALUMINUM
    > > MATERIAL 14.0 28.0855 2.329 14.0 SILICON
    > > LOW-MAT 14.0 14.0 -2.0 293.0 SILICON
    > > MATERIAL 18.0 39.948 0.00166 20.0 ARGON
    > > LOW-MAT 20.0 8.0 -2.0 293.0 ARGON
    > > MATERIAL 20.0 40.078 1.550 21.0 CALSIUM
    > > LOW-MAT 21.0 20.0 -2.0 293.0 CALSIUM
    > > MATERIAL 26.0 55.845 7.874 11.0 IRON
    > > LOW-MAT 11.0 26.0 -2.0 293.0 IRON
    > > MATERIAL 28.0 58.6934 8.902 25.0 NICKEL
    > > LOW-MAT 19.0 11.0 -2.0 293.0 SODIUM
    > > *
    > > MATERIAL 15.0 30.97376 1.82 26.0 PHOSPHO
    > > LOW-MAT 26.0 15.0 -2.0 293.0 PHOSPHO
    > > MATERIAL 16.0 32.06 2.07 27.0 SULFUR
    > > LOW-MAT 27.0 16.0 -2.0 293.0 SULFUR
    > > MATERIAL 19.0 39.0983 0.862 28.0 POTASSIU
    > > LOW-MAT 28.0 19.0 -2.0 293.0 POTASSIU
    > > MATERIAL 24.0 51.996 7.19 29.0 CHROMIUM
    > > LOW-MAT 29.0 24.0 -2.0 293.0 CHROMIUM
    > > MATERIAL 25.0 54.9380 7.43 30.0 MANGANES
    > > LOW-MAT 30.0 25.0 55.0 293.0 MANGANES
    > > MATERIAL 64.0 157.25 7.895 31.0 GADOLINI
    > > LOW-MAT 31.0 64.0 -2.0 293.0 GADOLINI
    > > * No Xsec of Gd for low energy netron
    > > *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    > > * granite ROCK
    > > MATERIAL 2.70 32.0 GRANITE
    > > COMPOUND -0.4850 8.0 -0.3430 14.0 -0.0800 10.0 GRANITE
    > > COMPOUND -0.02 11.0 -0.001 9.0 -0.0020 21.0 GRANITE
    > > COMPOUND -0.0240 19.0 -0.0450 28.0 GRANITE
    > > * air
    > > MATERIAL 0.001205 33.0 AIR
    > > COMPOUND -0.7550 7.0 -0.2321 8.0 -0.0128 20.0 AIR
    > > COMPOUND -0.0001 6.0 AIR
    > > * water of the buffer
    > > MATERIAL 1.00 34.0 WATER
    > > COMPOUND -0.1119 3.0 -0.8881 8.0 WATER
    > > * stainless steel
    > > MATERIAL 8.1 35.0 STEEL
    > > COMPOUND -0.70845 11.0 -0.0008 6.0 -0.02 30.0 STEEL
    > > COMPOUND -0.00045 26.0 -0.00030 27.0 -0.01 14.0 STEEL
    > > COMPOUND -0.18 29.0 -0.08 25.0 STEEL
    > > * mineral oil
    > > MATERIAL 0.838 36.0 OIL
    > > COMPOUND -0.8514 6.0 -0.1486 3.0 OIL
    > > * acrylic
    > > MATERIAL 1.18 37.0 ACRYLIC
    > > COMPOUND -0.59984 6.0 -0.08055 3.0 -0.31961 8.0 ACRYLIC
    > > * liquid scintillator
    > > MATERIAL 0.855 38.0 LISCIN
    > > COMPOUND -0.8544 6.0 -0.1289 3.0 -0.0003 7.0 LISCIN
    > > COMPOUND -0.0164 8.0 LISCIN
    > > * liquid scintillator; 0.1% Gd loaded
    > > MATERIAL 1.032 39.0 LISCINGD
    > > COMPOUND -0.8535 6.0 -0.1288 3.0 -0.0003 7.0 LISCINGD
    > > *COMPOUND -0.0174 8.0 LISCINGD
    > > COMPOUND -0.0164 8.0 -0.0010 31.0 LISCINGD
    > > *
    > > * External black hole
    > > ASSIGNMAT 1.0 1.0
    > > * granite rock over the exp. hall
    > > ASSIGNMAT 32.0 2.0
    > > * air in the exp. hall
    > > ASSIGNMAT 33.0 3.0
    > > * water of the buffer
    > > ASSIGNMAT 34.0 4.0
    > > * mineral oil
    > > ASSIGNMAT 36.0 5.0
    > > * liquid scintillator
    > > ASSIGNMAT 38.0 6.0
    > > * LS+Gd
    > > ASSIGNMAT 39.0 7.0
    > > *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    > > STEPSIZE 1000.0 1.0 7.0
    > > DELTARAY 1e-5 39.0
    > > *EMF EMF-OFF
    > > EMF
    > > EMFCUT 0.01 0.01 0.0 1.0 7.0 1.0
    > > *FLUKAFIX 0.001 39.0
    > > * LANDAU 1.0 1.0 26.0
    > > LOW-NEUT
    > > MUPHOTON 1.0 39.0
    > > PAIRBREM 3.0 0.0 0.0001 39.0 KELNER
    > > PHOTONUC 1.0 39.0
    > > *OUTLEVEL 1.0 0.0
    > > THRESHOL 0.000001
    > > *PART-THR -0.001 1.0 64.0 0.0 1.0
    > > *DISCARD 1. 2. 3. 4. 7. 8.
    > > *DISCARD 9. 12. 13. 14. 15. 16.
    > > *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    > > USRTRACK 1.0 8.0 47.0 7.0 .00E8 40. centdet
    > > USRTRACK 2. 0.0 &
    > > USRTRACK 1.0 212.0 47.0 7.0 1.00E8 24. Muons
    > > USRTRACK 12. 0.0 &
    > > *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    > > USRBIN 10.0 8.0 1000.0 1000.0 4000.0 BIN20
    > > USRBIN -1000.0 -1000.0 0.0 1.0 1.0 20.0 &
    > > USRBIN 0.0 8.0 1000.0 1000.0 4000.0 BIN20
    > > USRBIN -1000.0 -1000.0 0.0 1.0 1.0 20.0 &
    > > RANDOMIZE 1.0
    > >
    > *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
    > > *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789 *23456789
    > > START 1000000.0 99999999. 0.0 0.0
    > > STOP
    > > "
    > >
    > > The end statment of test1001.out is
    > > "***** Neutron Xsec: ****
    > > **** Group limits, average energies, velocities and momenta start at
    > location 1124211, end at location 1124834 (I*4 addr.) *****
    > > *** Recoil proton production activated for Xsec mat. # 1 ***
    > > **** Low energy neutron xsec not found for some media 9 12 ****
    > > NITROGEN
    > > ARGON
    > > CALSIUM
    > > "
    > > The end statment of test1002.log is "STOP NO RANDOM FILE"
    > >
    > > The test1002.err is
    > > "Abort called from FLRM64 reason NO RANDOM FILE Run stopped!
    > > STOP NO RANDOM FILE"
    > >
    > > The end statment of test1002.out is
    > > "
    > > **** No Random file available !!!!!! ****
    > > Abort called from FLRM64 reason NO RANDOM FILE Run stopped!
    > > STOP NO RANDOM FILE
    > > "
    > >
    > > I checked again and again, but have not found any bugs until now. This problem
    > > has stop me three weeks.
    > > So I sincerely hope warmhearted persons to give me a hand.
    > >
    > > Thanks in advance!
    > > With best regards,
    > >
    > > Ruiguang Wang
    > >
    > >
    > >
    > ------- End of Forwarded Message -------
    >
    >
    > --
    > Open WebMail Project (http://openwebmail.org)
    >
    >

    -- 
    Alberto FassÚ
    SLAC-RP, MS 48, 2575 Sand Hill Road, Menlo Park CA 94025
    Phone: (1 650) 926 4762   Fax: (1 650) 926 3569
    fasso@slac.stanford.edu
    

  • Next message: Alfredo Ferrari: "Release 2006.3: reminder for those who were using a pre-release"

    This archive was generated by hypermail 2.1.6 : Thu Sep 14 2006 - 08:57:28 CEST