# MCNP error about dose claculation of phantom

1. May 23, 2013

### Frank Yang

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. May 24, 2013

### Frank Yang

Can anybody solve my questions? Please ~~

3. May 25, 2013

### Frank Yang

Any sugestions?

4. Jul 4, 2013

### crnjak

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