Help to run input file (MCNPX)

[zakaria]

Hello,

I'm working on an MCNPX simulation related to internal radiotherapy. Unfortunately, I'm currently unable to run the code on my machine due to technical limitations.

If possible, could someone kindly run the input file and share the output or let me know whether the simulation runs successfully?

I'm not asking for result analysis just a quick check to see if my input file is correct and executable.

Thank you very much in advance

 

[berkeman]

Welcome to PF.

Unfortunately, I'm currently unable to run the code on my machine due to technical limitations.

What kind of limitations? We don't usually do folks' work for them.

 

[zakaria]

Welcome to PF.


What kind of limitations? We don't usually do folks' work for them.

Thank you very much for your reply.
I am not asking anyone to do my work for me, but I have tried many times to run the simulation myself without success. MCNPX stops immediately when the simulation starts.

 

[berkeman]

Does it give any error message? What version of MCNP are you using?

 

[Alex A]

There seems to be at least one mistake in the input file.

  cora1 0.0 0.023 2.3
  corb1 0 0.25 2.5
  corc1 0 1 90

Is this intended to produce a mesh containing 8 or 90'000 elements (voxels)?

For cylindrical geometry corc is the angle, it has to start with a 0, so this is implied. In MCNP this is really common, energy bins, histogram distributions etc. Putting an extra 0 on this line is asking for mesh elements (voxels) between angle 0 and angle 0, which have zero volume. This triggers a fatal error. That may be a cause of this input file not working. I also don't know if corc is supposed to end with a 360, some implementations of mesh need this.

So "corc1 1 90" would fix that error, but if the intention is 1 division per degree up to 90 degrees, then "corc1 89i 90" would do it and the other lines might also need fixing. i{n} interpolates n values and saves having long lists in the input file, it would be the same as corc1 [0] 1 2 3 4 5 .... 87 88 89 90.