Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

MCNP error about dose claculation of phantom

  1. May 23, 2013 #1
    Hello,everyone.Lastly,I use MCNP to claculate the dose distribution of CT phantom,but when I run my code,I get a message that says
    ' bad trouble in subroutine newcel of mcrun
    source particle no. 18859
    starting random number = 112983690142369
    the surface crossed is not a surface of this cell.
    1problem summary
    run terminated because of bad trouble.. '

    What exactly does that mean and how can i correct it? Help,please~
    And another question is when I change the 'mode p e' to 'mode p',this question is solved,but that makes me cannot claculate the dose distribution valid. SOS~

    Part of my code :
    c *-1----*----2----*----3----*----4----*----5----*----6----*----7----*----8----*
    1 0 +1:-2:+3:-4:+5:-6 imp:p=0 $out space;
    2 0 (-1 +2 -3 +4 -5 +6)#(+7 -8 +9 -10 +11 -12) imp:p=1 $
    3 0 +7 -8 +9 -10 +11 -12 imp:p=1 fill=1 $phantom
    c *-1----*----2----*----3----*----4----*----5----*----6----*----7----*----8----*
    100 10 -0.00127 +13 -14 +15 -16 +17 -18 imp:p=1 lat=1 u=1
    fill=-63:64 -63:64 -19:20
    1 3497R 2 1R 1 122R 2 4R 1 121R 2 3R
    (omit other...)
    101 1 -0.195 +19 -20 +21 -22 +23 -24 imp:p=1 u=2 $lung
    ...
    109 9 -1.609 +19 -20 +21 -22 +23 -24 imp:p=1 u=10

    c * surface section:
    1 px 500
    2 px -500
    3 py 500
    4 py -500
    5 pz 500
    6 pz -500
    7 px -11.61
    8 px 11.43 9 py -11.61
    10 py 11.43
    11 pz -6.15
    12 pz 5.85
    c *-1----*----2----*----3----*----4----*----5----*----6----*----7----*----8----*
    c * (1) container:
    13 px -0.09
    14 px 0.09
    15 py -0.09
    16 py 0.09
    17 pz -0.15
    18 pz 0.15
    c * (2) cell:
    19 px -0.1
    20 px 0.1
    21 py -0.1
    22 py 0.1
    23 pz -0.16
    24 pz 0.16

    c * data card section:
    mode p e
    sdef erg 6 par 2 x d1 y 94.49 z d2 vec 0 -1 0 dir 1 ara=100
    si1 -5 5
    sp1 0 1
    si2 -5 5
    sp2 0 1
    f15:p -0.13 12.77 -0.37 0.01 -0.13 11.77 -0.37 0.01 -0.13 9.77 -0.37 0.01
    de15 0.01 0.015 0.02 0.03 0.04 0.05 0.06 0.08 0.1 0.15 0.2 0.3 0.4
    0.5 0.6 0.8 1. 1.5 2. 3. 4. 5. 6. 8. 10. $(MeV)
    df15 lin 2.78E-06 1.11E-06 5.883-07 2.563-07 1.56E-07 1.20E-07 1.11E-07
    1.20E-07 1.47E-07 2.38E-07 3.45E-07 5.56E-07 7.69E-07 9.09E-07 1.14E-06
    1.47E-06 1.79E-06 2.44E-06 3.03E-06 4.00E-06 4.76E-06 5.56E-06 6.25E-06
    7.69E-06 9.09E-06 $ (rem/hr)/(photons/cm2.s)
    nps 1000000
    (omit material...)
     
  2. jcsd
  3. May 24, 2013 #2
    Can anybody solve my questions? Please ~~
     
  4. May 25, 2013 #3
    Any sugestions?
     
  5. Jul 4, 2013 #4
    may I ask something about your lattice; what is fill=-63:64 -63:64 -19:20///// from where did you get that numbers?? I am also working on phantoms, but I am a begginer.. :(((
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: MCNP error about dose claculation of phantom
  1. Mcnp error (Replies: 7)

Loading...