ChinaNet

标题: 求助mcnp连续计算的办法 [打印本页]

作者: Evan 时间: 2015-4-8 10:27
标题: 求助mcnp连续计算的办法
遇到一个情况，需要计算不同能量光子沉积能量情况，需要连续计算几十个能量点，当然我可以每计算一个后手动修改输入文件接着计算完成，
但我想mcnp应该可以解决，即在没有我的干预下完成上述过程，这个应该如何实现？Thanks in advance.

可以批处理啊。教你一个经验的方法。先新建一个txt 文档，在里面写上
                           copy  test .inp inp
                                                         del  outp
                                                            del runtpe
                                                            mcnp
                                                            copy  outp  test.out 这里面test是文件名，你可以将test改成你想用的文件名，如果你想运行一些输入文件的话，就接在第一行后面写输入文件的名字，象第一行那样，当然了，也要在最后面加上象最后一行语句的语句。懂我的意思不？

批处理我会的，已经解决，我是问能否mcnp自己解决。
我的批处理是这样的：（需要配合按键精灵，因为需要第三方程序终止mcnp的一次运行，而dos单程序串行执行）
可以从10keV计算到25keV（将输入文件直接通过批处理输入及更改）
cd c:\mcnp\mcnp4c-t
set e=5
:mclp
set /a e+=5
echo ArCO2 70:30 5mm 10*10cm2 5:5:100keV >EKK.i
echo 1 1 -1.8391E-3 1 -2 3 -4 5 -6 >>EKK.i
echo 2 0          #1 -7 >>EKK.i
echo 3 0          7 >>EKK.i
echo.>>EKK.i
echo 1 PZ -0.25 >>EKK.i
echo 2 PZ  0.25 >>EKK.i
echo 3 PX  -5 >>EKK.i
echo 4 PX 5 >>EKK.i
echo 5 PY  -5 >>EKK.i
echo 6 PY 5 >>EKK.i
echo 7 SO  100 >>EKK.i
echo.>>EKK.i
echo mode p e >>EKK.i
echo imp:p 1 1 0 >>EKK.i
echo imp:e 1 0 0 >>EKK.i
echo sdef pos 0 -6 0 vec 0 1 0 dir 1 erg 0.0%e% par 2 >>EKK.i
echo C +f8:e 1 >>EKK.i
echo f6:p 1 >>EKK.i
echo m1  18000.42c  0.7  6000.60c 0.3  8016.60c 0.6 >>EKK.i
echo CTME  0.25 >>EKK.i
echo NPS  1E8 >>EKK.i
mcnp inp=EKK.i ixr
taskkill /im mcnp.exe
if %e% lss 25 goto mclp
pause

转自52mc

欢迎光临 ChinaNet (http://www.nuclear.net.cn/)