"No tally, no plot" error MCNP6

[khary23]

I am using the TMESH tally to try and get the flux of a 1 MeV point source through a spherical MESH a distance of 1cm away from the origin and covering all polar and azimuthal angles. When I try to plot using the RUNTPE function I get the error "no tally no plot". I also do not see tally results in the MDATA file. I am teaching myself MCNP6 and find the manual to quite difficult as it seems to presuppose knowledge of older versions of MCNP. I would appreciate it if someone could look at my output file that has the tally cards and help me figure out what I am doing wrong.

 

[PSRB191921]

Hi,
in cora, corb, corc you must put the number of step.
For example :
TMESH
SMESH1:p
CORA1 0 10i 1
CORB1 0.01 10i 180
CORC1 0.01 10i 360
ENDMD

note that you have a cukkie cutter.
Put
120 0 2 IMP:P=0 $Graveyard
instead of 120 0 -12 2 IMP:P=0 $Graveyard

it gives Something like that (with vised)

 

[khary23]

Thank you for the response and help! Wjats do you mean by cookie cutter?

 

[PSRB191921]

a cookie cutter is a geometry error :
Each region must be uniquely defined. If not it gives a cookie cutter

 

[khary23]

Ahh that makes sense. What software are you using to visualize?

 

[PSRB191921]

I used vised or xming

 

[khary23]

Thanks!

 

[khary23]

I have changed the deck with your suggestions and I am still getting the no tally no plot error. The mdata file does contain information so I don't know why its not plotting. Using Xming.

 

[PSRB191921]

When you are in xming I suppose you are in "MCPLOT"
if yes, type "PLOT" and try again.
Another suggestion, instead of "CORA1 0 10i 10" try with "CORA1 0 10i 1" because the radius of your sphere is 1 cm not 10 cm.

 

[khary23]

Tried what you advised and still nothing. I attached a screen shot of what the plot window shows

 

[PSRB191921]

I will try to write a tutorial

 

[khary23]

Oh wow Thank you!

 

[khary23]

One more question. How does one get the tally data from the TMESH in a way that is similar to the output of a normal F4/F6.. tally. What I need is flux in each mesh bin for example what is the flux at each cm from 0.5 to 10 cm or the angular flux in 1 degree increments.

 

[PSRB191921]

Try with this tuto

 

[khary23]

Awesome thank you so much! So is there a way to get the angular bin data in a table form?

 

[PSRB191921]

Awesome thank you so much! So is there a way to get the angular bin data in a table form?

yes, you can use "Gridconv" for that

 

[khary23]

I used it, but the output is not in a form that makes sense. No descriptors or labels on the data just unformatted numbers.

 

[PSRB191921]

Yes ! Mesh tally files are not well documented
To understand the file it is necessary to work in 1D (and not in 3D)
For that you can put in your MCNP file :
SMESH1:p
CORA1 0 100i 1
CORB1 .0001
CORC1 .0001
ENDMD

In this MESH tally you will have 100 intervals for the radius and only one in phi and theta.
with "gridconv" the ASCII file is essak.txt

the file essakcomment.txt is commented.
It gives :

0​	2,44E+03​
9,90E-03​	3,48E+02​
1,98E-02​	1,28E+02​
2,97E-02​	6,58E+01​
3,96E-02​	3,99E+01​
4,95E-02​	2,68E+01​
5,94E-02​	1,92E+01​
6,93E-02​	1,44E+01​
7,92E-02​	1,12E+01​
8,91E-02​	8,99E+00​
9,90E-02​	7,36E+00​
0,108911​	6,13E+00​
0,118812​	5,19E+00​
0,128713​	4,45E+00​
0,138614​	3,86E+00​
0,148515​	3,38E+00​
0,158416​	2,98E+00​
0,168317​	2,65E+00​
0,178218​	2,37E+00​
0,188119​	2,13E+00​
0,19802​	1,93E+00​
0,207921​	1,76E+00​
0,217822​	1,60E+00​
0,227723​	1,47E+00​
0,237624​	1,35E+00​
0,247525​	1,25E+00​
0,257426​	1,16E+00​
0,267327​	1,07E+00​
0,277228​	9,99E-01​
0,287129​	9,33E-01​
0,29703​	8,73E-01​
0,306931​	8,18E-01​
0,316832​	7,69E-01​
0,326733​	7,23E-01​
0,336634​	6,82E-01​
0,346535​	6,44E-01​
0,356436​	6,09E-01​
0,366337​	5,77E-01​
0,376238​	5,48E-01​
0,386139​	5,20E-01​
0,39604​	4,95E-01​
0,405941​	4,71E-01​
0,415842​	4,49E-01​
0,425743​	4,29E-01​
0,435644​	4,10E-01​
0,445545​	3,92E-01​
0,455446​	3,75E-01​
0,465347​	3,60E-01​
0,475248​	3,45E-01​
0,485149​	3,31E-01​
0,49505​	3,18E-01​
0,50495​	3,06E-01​
0,514851​	2,95E-01​
0,524752​	2,84E-01​
0,534653​	2,73E-01​
0,544554​	2,64E-01​
0,554455​	2,54E-01​
0,564356​	2,46E-01​
0,574257​	2,37E-01​
0,584158​	2,29E-01​
0,594059​	2,22E-01​
0,60396​	2,15E-01​
0,613861​	2,08E-01​
0,623762​	2,01E-01​
0,633663​	1,95E-01​
0,643564​	1,89E-01​
0,653465​	1,84E-01​
0,663366​	1,78E-01​
0,673267​	1,73E-01​
0,683168​	1,68E-01​
0,693069​	1,63E-01​
0,70297​	1,59E-01​
0,712871​	1,54E-01​
0,722772​	1,50E-01​
0,732673​	1,46E-01​
0,742574​	1,42E-01​
0,752475​	1,39E-01​
0,762376​	1,35E-01​
0,772277​	1,32E-01​
0,782178​	1,28E-01​
0,792079​	1,25E-01​
0,80198​	1,22E-01​
0,811881​	1,19E-01​
0,821782​	1,16E-01​
0,831683​	1,14E-01​
0,841584​	1,11E-01​
0,851485​	1,08E-01​
0,861386​	1,06E-01​
0,871287​	1,04E-01​
0,881188​	1,01E-01​
0,891089​	9,91E-02​
0,90099​	9,70E-02​
0,910891​	9,49E-02​
0,920792​	9,29E-02​
0,930693​	9,09E-02​
0,940594​	8,90E-02​
0,950495​	8,72E-02​
0,960396​	8,54E-02​
0,970297​	8,37E-02​
0,980198​	8,20E-02​
0,990099​	8,04E-02​
1​

 

[khary23]

Wow this was incredibly helpful! I gave it a try for a grid with radii of 1cm and 5cm and polar angles running from 10-180* . If my understanding is correct after the list of angles there is the flux associated with that grid and after that the error ?

 

[PSRB191921]

Yes first you have the fluence value and after the error.
You have the value for the radius between 1 to 5 cm every 10 °

 

[khary23]

Thank you for the clarification.
I am also trynig to print a MTCAL file using the PRDMP card and wanted to ask if this looks correct. The syntax I am using is
PRDMP NDP NDM MCT NDMP DMMP
PRDMP 100000 -1 1 100000 0
Where NDP= 100000 is to print tallies after 100000 histories
NDM = -1 dump after every minute
MCT = 1 create MCTAL file
NDMP = 100 keep a hundred data dumps
DMMP = 100000 write fluctuations for every 100000 histories

 

[PSRB191921]

Sorry I'm not an user of MCTAL or PRDMP. I cannot help you !

 

[khary23]

Yes ! Mesh tally files are not well documented
To understand the file it is necessary to work in 1D (and not in 3D)
For that you can put in your MCNP file :
SMESH1:p
CORA1 0 100i 1
CORB1 .0001
CORC1 .0001
ENDMD

In this MESH tally you will have 100 intervals for the radius and only one in phi and theta.
with "gridconv" the ASCII file is essak.txt

the file essakcomment.txt is commented.
It gives :

0​	2,44E+03​
9,90E-03​	3,48E+02​
1,98E-02​	1,28E+02​
2,97E-02​	6,58E+01​
3,96E-02​	3,99E+01​
4,95E-02​	2,68E+01​
5,94E-02​	1,92E+01​
6,93E-02​	1,44E+01​
7,92E-02​	1,12E+01​
8,91E-02​	8,99E+00​
9,90E-02​	7,36E+00​
0,108911​	6,13E+00​
0,118812​	5,19E+00​
0,128713​	4,45E+00​
0,138614​	3,86E+00​
0,148515​	3,38E+00​
0,158416​	2,98E+00​
0,168317​	2,65E+00​
0,178218​	2,37E+00​
0,188119​	2,13E+00​
0,19802​	1,93E+00​
0,207921​	1,76E+00​
0,217822​	1,60E+00​
0,227723​	1,47E+00​
0,237624​	1,35E+00​
0,247525​	1,25E+00​
0,257426​	1,16E+00​
0,267327​	1,07E+00​
0,277228​	9,99E-01​
0,287129​	9,33E-01​
0,29703​	8,73E-01​
0,306931​	8,18E-01​
0,316832​	7,69E-01​
0,326733​	7,23E-01​
0,336634​	6,82E-01​
0,346535​	6,44E-01​
0,356436​	6,09E-01​
0,366337​	5,77E-01​
0,376238​	5,48E-01​
0,386139​	5,20E-01​
0,39604​	4,95E-01​
0,405941​	4,71E-01​
0,415842​	4,49E-01​
0,425743​	4,29E-01​
0,435644​	4,10E-01​
0,445545​	3,92E-01​
0,455446​	3,75E-01​
0,465347​	3,60E-01​
0,475248​	3,45E-01​
0,485149​	3,31E-01​
0,49505​	3,18E-01​
0,50495​	3,06E-01​
0,514851​	2,95E-01​
0,524752​	2,84E-01​
0,534653​	2,73E-01​
0,544554​	2,64E-01​
0,554455​	2,54E-01​
0,564356​	2,46E-01​
0,574257​	2,37E-01​
0,584158​	2,29E-01​
0,594059​	2,22E-01​
0,60396​	2,15E-01​
0,613861​	2,08E-01​
0,623762​	2,01E-01​
0,633663​	1,95E-01​
0,643564​	1,89E-01​
0,653465​	1,84E-01​
0,663366​	1,78E-01​
0,673267​	1,73E-01​
0,683168​	1,68E-01​
0,693069​	1,63E-01​
0,70297​	1,59E-01​
0,712871​	1,54E-01​
0,722772​	1,50E-01​
0,732673​	1,46E-01​
0,742574​	1,42E-01​
0,752475​	1,39E-01​
0,762376​	1,35E-01​
0,772277​	1,32E-01​
0,782178​	1,28E-01​
0,792079​	1,25E-01​
0,80198​	1,22E-01​
0,811881​	1,19E-01​
0,821782​	1,16E-01​
0,831683​	1,14E-01​
0,841584​	1,11E-01​
0,851485​	1,08E-01​
0,861386​	1,06E-01​
0,871287​	1,04E-01​
0,881188​	1,01E-01​
0,891089​	9,91E-02​
0,90099​	9,70E-02​
0,910891​	9,49E-02​
0,920792​	9,29E-02​
0,930693​	9,09E-02​
0,940594​	8,90E-02​
0,950495​	8,72E-02​
0,960396​	8,54E-02​
0,970297​	8,37E-02​
0,980198​	8,20E-02​
0,990099​	8,04E-02​
1​

Is the end of the tally output the error or uncertainty?

 

[PSRB191921]

it is the error

 

[khary23]

Is there way to get the relitive uncertainty?