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/) Powered by Discuz! X3.1