"No tally, no plot" error MCNP6

  • Thread starter khary23
  • Start date
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.
 

Attachments

109
32
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)

242841
 
Thank you for the response and help! Wjats do you mean by cookie cutter?
 
109
32
a cookie cutter is a geometry error :
Each region must be uniquely defined. If not it gives a cookie cutter
 
Ahh that makes sense. What software are you using to visualize?
 
109
32
I used vised or xming
 
Thanks!
 
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.
 

Attachments

109
32
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.
 
Tried what you advised and still nothing. I attached a screen shot of what the plot window shows
 

Attachments

109
32
I will try to write a tutorial
 
Oh wow Thank you!
 
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.
 
Awesome thank you so much! So is there a way to get the angular bin data in a table form?
 
109
32
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
 
I used it, but the output is not in a form that makes sense. No descriptors or labels on the data just unformatted numbers.
 
109
32
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
 

Attachments

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 ?
 

Attachments

109
32
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 °
 
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
 
109
32
Sorry I'm not an user of MCTAL or PRDMP. I cannot help you !
 
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?
 
109
32
it is the error
 
Is there way to get the relitive uncertainty?
 

Want to reply to this thread?

""No tally, no plot" error MCNP6" You must log in or register to reply here.

Physics Forums Values

We Value Quality
• Topics based on mainstream science
• Proper English grammar and spelling
We Value Civility
• Positive and compassionate attitudes
• Patience while debating
We Value Productivity
• Disciplined to remain on-topic
• Recognition of own weaknesses
• Solo and co-op problem solving
Top