Compare results of using mctal_merge and ZWelch's version

[makeclean]

Compare the results produced by mctal_merge and @zwelchWI version of this implemented using python. To ensure that we get the same results, (averages and errors) in both cases.

• Compare the averages
• Repeat the same for merge meshtal

[zwelchWI]

There are a few differences between mctal_merge and mctal_combine.py that I am not sure why they differ. Lines 12, 13, and 15 have an extra 0 at the end of the line in merge_mctal that isn't in any of the input files. Other than those lines they produce very similar results. I ran into a bug when testing the meshtal_combine.py, hopefully will have that fixed later today.

merge_mctal

     1  mcnp    5         04/12/13 11:44:29    2        22258819    104995695800
     2   c     Phoenix Nuclear Labs - Neutron Shielding                                 
     3  ntal     1
     4     12
     5  tally   12    1    0
     6  f        1
     7      4.2
     8  d        1
     9  u        0
    10  s        0
    11  m        0
    12  c        0   0
    13  et       3   0
    14    1.00000E-08  3.00000E-07
    15  t        0   0
    16  vals
    17    7.45113E+04 0.0314  1.35557E+06 0.0076  1.43009E+06 0.0073
    18  tfc   14       1       1       1       1       1       1       3       1
    19        64000  1.50329E+06  1.09204E-01  5.62148E+01
    20       128000  1.35252E+06  8.05827E-02  5.14733E+01
    21       192000  1.33087E+06  6.70195E-02  4.96333E+01
    22       256000  1.49274E+06  7.89228E-02  2.68712E+01
    23       320000  1.44718E+06  6.85464E-02  2.85268E+01
    24       384000  1.46089E+06  5.95659E-02  3.14530E+01
    25       448000  1.46250E+06  5.34746E-02  3.34329E+01
    26       512000  1.46124E+06  4.89703E-02  3.48573E+01
    27       576000  1.46943E+06  4.53193E-02  3.61633E+01
    28       640000  1.42125E+06  4.32978E-02  3.56249E+01
    29       704000  1.41815E+06  4.06640E-02  3.67036E+01
    30       768000  1.44388E+06  3.84701E-02  3.75932E+01
    31       832000  1.45455E+06  3.64805E-02  3.86153E+01
    32       855496  1.44711E+06  3.59495E-02  3.86778E+01

mctal_combine.py

     1  mcnp            5 04/12/2013 11:44:02        2 22258819 104995695800
     2  c     Phoenix Nuclear Labs - Neutron Shielding
     3  ntal      1
     4      12   
     5  tally    12     1     0
     6  f         1
     7      4.2            
     8  d         1
     9  u         0
    10  s         0
    11  m         0
    12  c         0
    13  et        3
    14    1.00000e-08 3.00000e-07
    15  t         0
    16  vals
    17    7.45113e+04 3.11816e-02 1.35557e+06 7.60226e-03 1.43008e+06 7.39022e-03
    18  tfc 14  1   1   1   1   1   1   3   1
    19        64000  1.50329e+06  1.09204e-01  5.62148e+01
    20       128000  1.35252e+06  8.05827e-02  5.14733e+01
    21       192000  1.33087e+06  6.70195e-02  4.96333e+01
    22       256000  1.49274e+06  7.89228e-02  2.68712e+01
    23       320000  1.44718e+06  6.85464e-02  2.85268e+01
    24       384000  1.46089e+06  5.95659e-02  3.14530e+01
    25       448000  1.46250e+06  5.34746e-02  3.34329e+01
    26       512000  1.46124e+06  4.89703e-02  3.48573e+01
    27       576000  1.46943e+06  4.53193e-02  3.61633e+01
    28       640000  1.42125e+06  4.32978e-02  3.56249e+01
    29       704000  1.41815e+06  4.06640e-02  3.67036e+01
    30       768000  1.44388e+06  3.84701e-02  3.75932e+01
    31       832000  1.45455e+06  3.64805e-02  3.86153e+01
    32       855496  1.44711e+06  3.59495e-02  3.86778e+01

[makeclean]

Is there another one coming, or is this somehow both?

[zwelchWI]

sorry, I accidentally hit comment before I was done. It should be correct now

[zwelchWI]

I've mostly got meshtal_combine.py working, the only difference I am seeing is in the first line with the ld parameter. Im looking through the mcnp documentation but havent found anything helpful yet.

The only other difference is one prints 1.0 out as 1.00000e-00 and the other prints it out as 1.00000e+00. I can't imagine that being a problem but I figured I should mention it just in case

merge_mctal

zwelch@kermit:/filespace/groups/cnerg/users/zwelch/condor_MCNP > head -n 30 MESHTALMRG 
mcnp   version 5     ld=00002010  probid =  04/12/2013 13:07:27
 c     Phoenix Nuclear Labs - Neutron Shielding                                  
 Number of histories used for normalizing tallies =      2.22588e+07

 Mesh Tally Number   114
 This is a neutron mesh tally.
 This mesh tally is modified by a dose response function.

 Tally bin boundaries:
    X direction:   -500.00   -480.00   -460.00   -440.00   -420.00   -400.00   -380.00   -360.00   -340.00   -320.00   -300.00   -280.00   -260.00   -240.00   -220.00   -200.00   -180.00   -160.00   -140.00   -120.00   -100.00    -80.00    -60.00    -40.00    -20.00      0.00     20.00     40.00     60.00     80.00    100.00    120.00    140.00    160.00    180.00    200.00    220.00    240.00    260.00    280.00    300.00    320.00    340.00    360.00    380.00    400.00    420.00    440.00    460.00    480.00    500.00    520.00    540.00    560.00    580.00    600.00    620.00    640.00    660.00    680.00    700.00    720.00    740.00    760.00    780.00    800.00    820.00    840.00    860.00    880.00    900.00    920.00    940.00    960.00    980.00   1000.00
    Y direction:  -1000.00   -980.00   -960.00   -940.00   -920.00   -900.00   -880.00   -860.00   -840.00   -820.00   -800.00   -780.00   -760.00   -740.00   -720.00   -700.00   -680.00   -660.00   -640.00   -620.00   -600.00   -580.00   -560.00   -540.00   -520.00   -500.00   -480.00   -460.00   -440.00   -420.00   -400.00   -380.00   -360.00   -340.00   -320.00   -300.00   -280.00   -260.00   -240.00   -220.00   -200.00   -180.00   -160.00   -140.00   -120.00   -100.00    -80.00    -60.00    -40.00    -20.00      0.00     20.00     40.00     60.00     80.00    100.00    120.00    140.00    160.00    180.00    200.00    220.00    240.00    260.00    280.00    300.00    320.00    340.00    360.00    380.00    400.00    420.00    440.00    460.00    480.00    500.00    520.00    540.00    560.00    580.00    600.00    620.00    640.00    660.00    680.00    700.00    720.00    740.00    760.00    780.00    800.00    820.00    840.00    860.00    880.00    900.00    920.00    940.00    960.00    980.00   1000.00
    Z direction:    -10.00     10.00
    Energy bin boundaries: 0.00e+00 1.00e+36

        X         Y         Z     Result     Rel Error
   -490.000  -990.000     0.000 2.12560e-01 4.46549e-01
   -490.000  -970.000     0.000 1.68887e-01 4.09402e-01
   -490.000  -950.000     0.000 1.04797e-01 5.32169e-01
   -490.000  -930.000     0.000 2.44580e-01 6.04800e-01
   -490.000  -910.000     0.000 6.70120e-02 7.16681e-01
   -490.000  -890.000     0.000 1.33087e+00 8.93379e-01
   -490.000  -870.000     0.000 1.76535e+00 9.57541e-01
   -490.000  -850.000     0.000 1.12476e+00 9.04308e-01
   -490.000  -830.000     0.000 1.52138e-01 5.07835e-01
   -490.000  -810.000     0.000 0.00000e+00 0.00000e+00
   -490.000  -790.000     0.000 1.86568e-01 5.23251e-01
   -490.000  -770.000     0.000 2.20780e-01 5.70877e-01
   -490.000  -750.000     0.000 6.44030e-02 7.20003e-01
   -490.000  -730.000     0.000 1.35785e-01 6.33157e-01
   -490.000  -710.000     0.000 4.92995e-02 7.27811e-01

mine

zwelch@kermit:/filespace/groups/cnerg/users/zwelch/condor_MCNP > head -n 30 meshtest
mcnp   version 5    ld=10072009  probid = 04/12/2013 14:19:37
 c     Phoenix Nuclear Labs - Neutron Shielding                                  
 Number of histories used for normalizing tallies = 2.22588e+07

 Mesh Tally Number   114
 This is a neutron mesh tally.
 This mesh tally is modified by a dose response function.

 Tally bin boundaries:
    X direction:   -500.00   -480.00   -460.00   -440.00   -420.00   -400.00   -380.00   -360.00   -340.00   -320.00   -300.00   -280.00   -260.00   -240.00   -220.00   -200.00   -180.00   -160.00   -140.00   -120.00   -100.00    -80.00    -60.00    -40.00    -20.00      0.00     20.00     40.00     60.00     80.00    100.00    120.00    140.00    160.00    180.00    200.00    220.00    240.00    260.00    280.00    300.00    320.00    340.00    360.00    380.00    400.00    420.00    440.00    460.00    480.00    500.00    520.00    540.00    560.00    580.00    600.00    620.00    640.00    660.00    680.00    700.00    720.00    740.00    760.00    780.00    800.00    820.00    840.00    860.00    880.00    900.00    920.00    940.00    960.00    980.00   1000.00
    Y direction:  -1000.00   -980.00   -960.00   -940.00   -920.00   -900.00   -880.00   -860.00   -840.00   -820.00   -800.00   -780.00   -760.00   -740.00   -720.00   -700.00   -680.00   -660.00   -640.00   -620.00   -600.00   -580.00   -560.00   -540.00   -520.00   -500.00   -480.00   -460.00   -440.00   -420.00   -400.00   -380.00   -360.00   -340.00   -320.00   -300.00   -280.00   -260.00   -240.00   -220.00   -200.00   -180.00   -160.00   -140.00   -120.00   -100.00    -80.00    -60.00    -40.00    -20.00      0.00     20.00     40.00     60.00     80.00    100.00    120.00    140.00    160.00    180.00    200.00    220.00    240.00    260.00    280.00    300.00    320.00    340.00    360.00    380.00    400.00    420.00    440.00    460.00    480.00    500.00    520.00    540.00    560.00    580.00    600.00    620.00    640.00    660.00    680.00    700.00    720.00    740.00    760.00    780.00    800.00    820.00    840.00    860.00    880.00    900.00    920.00    940.00    960.00    980.00   1000.00
    Z direction:    -10.00     10.00
    Energy bin boundaries: 0.00e+00 1.00e+36

        X         Y         Z     Result     Rel Error
   -490.000  -990.000     0.000 2.12560e-01 4.46549e-01
   -490.000  -970.000     0.000 1.68887e-01 4.09402e-01
   -490.000  -950.000     0.000 1.04797e-01 5.32169e-01
   -490.000  -930.000     0.000 2.44580e-01 6.04800e-01
   -490.000  -910.000     0.000 6.70120e-02 7.16681e-01
   -490.000  -890.000     0.000 1.33087e+00 8.93379e-01
   -490.000  -870.000     0.000 1.76535e+00 9.57541e-01
   -490.000  -850.000     0.000 1.12476e+00 9.04308e-01
   -490.000  -830.000     0.000 1.52138e-01 5.07835e-01
   -490.000  -810.000     0.000 0.00000e+00 0.00000e+00
   -490.000  -790.000     0.000 1.86568e-01 5.23251e-01
   -490.000  -770.000     0.000 2.20780e-01 5.70877e-01
   -490.000  -750.000     0.000 6.44030e-02 7.20003e-01
   -490.000  -730.000     0.000 1.35785e-01 6.33157e-01
   -490.000  -710.000     0.000 4.92995e-02 7.27811e-01

[makeclean]

Need to define expected functionality of all operation modes of the mesh_merge script, see #12