The MCNP code never ends, what went wrong?

[momomo]

I'm using the fill card to repeat the 60Co source to simulate the radiation source in the irradiation chamber, the code works, and the vised shows the graph. But it just doesn't end. Thank you very much for your request.

 

[Alex A]

Hi @momomo,
Welcome to Physics Forums. I can't fully read your SDEF line. Are you able to share your input file? If you are then changing the file name to add .txt will let you add it as an attachment.

 

[momomo]

Hi @momomo,
Welcome to Physics Forums. I can't fully read your SDEF line. Are you able to share your input file? If you are then changing the file name to add .txt will let you add it as an attachment.

25  0   -34  33  -36  35  -38  37  imp:p=1 fill=1 (-831.95 0 289)
26  2   -0.001293  -400(-33:34:-35:36:-37:38)#34#35#36#37
                 #38#39#40#41#43#44#45#46#48#49#50#51
                  #52#53#54#55#56#57#58  imp:p=1 
27  0   -26  25 -27 28 -30 29  imp:p=1 lat=1 u=1 fill=0:0 -24:25  0:0
             2 3 3 3 3 2 3 3 3 3 2 3 3 3 3 2 3 3 3 3 2 3 3 3 3 
             2 3 3 3 3 2 3 3 3 3 2 3 3 3 3 2 3 3 3 3 2 3 3 3 3 
28  4   -7.8  -32    imp:p=1  u=2 $60 Co
29  3   -7.83  -31  32   imp:p=1 u=2  
30  2   -0.001293  31   imp:p=1 u=2 
31  2   -0.001293 -32   imp:p=1 u=3 
32  3   -7.83  -31  32   imp:p=1 u=3 
33  2   -0.001293  31  imp:p=1 u=3 
34 like 25 but trcl=(-60.9 0 0) fill=1  (-892.85  0 289) 
35 like 25 but trcl=(-248.6 0 0) fill=1  (-1080.55 0 289) 
36 like 25 but trcl=(-309.5 0 0) fill=1  (-1141.45 0 289) 
37 like 25 but trcl=(0  0  -150.45) fill=1  (-831.95  0 138.55) 
38 like 25 but trcl=(-60.9 0 -150.45) fill=1  (-892.85 0 138.55)
39 like 25 but trcl=(-248.6  0 -150.45) fill=1  (-1080.55 0 138.55) 
40 like 25 but trcl=(-309.5  0 -150.45) fill=1  (-1141.45 0 138.55) 
41 0  -40  39  -36 35 -38 37   imp:p=1 fill=4  (-953.75 0 289) 
42 0  -26  25 -27 28 -30 29  imp:p=1 lat=1 u=4  fill=0:0 -24:25 0:0
             2 3 3 3 3 2 3 3 3 3 2 3 3 3 3 2 3 3 3 3 2 3 3 3 3 
             2 3 3 3 3 2 3 3 3 3 2 3 3 3 3 2 3 3 3 3 2 3 3 3 3  
43 like 41 but trcl=(-65.9  0  0) fill=4  (-1019.65  0 289)
44 like 41 but trcl=(0  0 -150.45)  fill=4  (-953.75 0 138.55)  
45 like 41 but trcl=(-65.9 0 -150.45) fill=4  (-1019.65 0 138.55) 
46 0  -34 33 -36 35 -42  41   imp:p=1 fill=5  (-831.95 0 238.85) 
47  0  400  imp:p=0

25  py -0.6
26  py  0.6
27  px  0
28  px  -1.118
29  pz  -23
30  pz   23
31  rcc   -0.559  0  -22.575 0   0  45.15  0.555
32  rcc   -0.559  0  -20.65   0   0   41.3   0.3
33  px  -859.9
34  px  -804
35  py  -0.6
36  py   0.6
37  pz   266
38  pz   312
39  px   -981.7
40  px  -925.8
41  pz   215.85
42  pz   261.85
400  so  100000

mode  P
m1   1000    -0.0056 
         8000    -0.4993 
         11000  -0.0171  
          12000  -0.0024 
          13000    -0.0454
           14000     -0.315 
           16000   -0.0012
            19000   -0.0192
           20000   -0.0826 
           26000   -0.0122
m2        6000  -0.000124          
              7014  -0.755268                 
              8016  -0.231781         
              18000 -0.012827 
m3        26000    -0.99 
              6000     -0.01
m4        27059     1      
SDEF  ERG=D1  par=2  axs=0 0 1  rad=d2  ext=d3   cel=d4 pos=-0.559  0 0
si1  L 1.17 1.33
sp1   0.5 0.5 
si2    0  0.3
sp2   -21  1
si3  -20.65  20.65 
sp3   0  1
si4  s  d5   d6   
sp4    0.4  0.6  
si5 L 25:27:28 34:27:28 35:27:28 36:27:28 37:27:28 38:27:28 39:27:28 40:27:28
sp5   1 1 1 1 1 1 1 1
si6  L  41:42:28   43:42:28 44:42:28 45:42:28 
sp6   1 1 1 1
f15：p 1000 550 236.35 5
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.0 1.5 2.0
DF15 林 2.78E-06 1.11E-06 5.88E-07 2.56E-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
NPS 500 系列
打印

 

[Alex A]

These versions seem unfinished, but when I removed the bits that were complaining it did just hang the program.

My quick answer is "I don't know". You seem to be building a 2D lattice by translating a 1D lattice. That may be not the best way to do it. I suggest you define a 2D lattice. If nothing else this should be easier to debug.

I will look again when I have more time.

 

[momomo]

这些版本似乎尚未完成，但是当我删除抱怨的部分时，它确实只是挂起了程序。

我的快速回答是“我不知道”。您似乎正在通过平移 1D 晶格来构建 2D 晶格。这可能不是最好的方法。我建议你定义一个 2D 晶格。如果不出意外，这应该更容易调试。

当我有更多时间时，我会再看一遍。

 

[momomo]

这些版本似乎尚未完成，但是当我删除抱怨的部分时，它确实只是挂起了程序。

我的快速回答是“我不知道”。您似乎正在通过平移 1D 晶格来构建 2D 晶格。这可能不是最好的方法。我建议你定义一个 2D 晶格。如果不出意外，这应该更容易调试。

当我有更多时间时，我会再看一遍。

非常感谢您的回答。

Translated by a Mentor using Google Translate:

Thank you very much for your answer.