MCDancoff.html 42.1 KB
Newer Older
Batson Iii's avatar
Batson Iii committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622


<!DOCTYPE html>
<html class="writer-html5" lang="en" >
<head>
  <meta charset="utf-8">
  
  <meta name="viewport" content="width=device-width, initial-scale=1.0">
  
  <title>MCDancoff Data Guide &mdash; SCALE Manual 0.0.1 documentation</title>
  

  
  <link rel="stylesheet" href="_static/css/theme.css" type="text/css" />
  <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
  <link rel="stylesheet" href="_static/custom.css" type="text/css" />

  
  
  
  

  
  <!--[if lt IE 9]>
    <script src="_static/js/html5shiv.min.js"></script>
  <![endif]-->
  
    
      <script type="text/javascript" id="documentation_options" data-url_root="./" src="_static/documentation_options.js"></script>
        <script src="_static/jquery.js"></script>
        <script src="_static/underscore.js"></script>
        <script src="_static/doctools.js"></script>
        <script src="_static/language_data.js"></script>
        <script async="async" src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.7/latest.js?config=TeX-AMS-MML_HTMLorMML"></script>
    
    <script type="text/javascript" src="_static/js/theme.js"></script>

    
    <link rel="index" title="Index" href="genindex.html" />
    <link rel="search" title="Search" href="search.html" />
    <link rel="next" title="CAJUN: Module for Combining and Manipulating CENTRM Continuous-Energy Libraries" href="CAJUN.html" />
    <link rel="prev" title="CRAWDAD: Module to Produce CENTRM-Formatted Continuous-Energy Nuclear Data Libraries" href="CRAWDAD.html" /> 
</head>

<body class="wy-body-for-nav">

   
  <div class="wy-grid-for-nav">
    
    <nav data-toggle="wy-nav-shift" class="wy-nav-side">
      <div class="wy-side-scroll">
        <div class="wy-side-nav-search" >
          

          
            <a href="index.html" class="icon icon-home" alt="Documentation Home"> SCALE Manual
          

          
          </a>

          
            
            
          

          
<div role="search">
  <form id="rtd-search-form" class="wy-form" action="search.html" method="get">
    <input type="text" name="q" placeholder="Search docs" />
    <input type="hidden" name="check_keywords" value="yes" />
    <input type="hidden" name="area" value="default" />
  </form>
</div>

          
        </div>

        
        <div class="wy-menu wy-menu-vertical" data-spy="affix" role="navigation" aria-label="main navigation">
          
            
            
              
            
            
              <p class="caption"><span class="caption-text">Reactor Physics</span></p>
<ul>
<li class="toctree-l1"><a class="reference internal" href="Polaris.html">Polaris: 2D Light Water Reactor Lattice Physics Module</a></li>
<li class="toctree-l1"><a class="reference internal" href="PolarisA.html">SCALE 6.3 Polaris Input Format</a></li>
</ul>
<p class="caption"><span class="caption-text">Criticality Safety</span></p>
<ul>
<li class="toctree-l1"><a class="reference internal" href="Criticality%20Safety%20Overview.html">Criticality Safety Overview</a></li>
<li class="toctree-l1"><a class="reference internal" href="CSAS5.html">CSAS5:  Control Module For Enhanced Criticality Safety Analysis Sequences With KENO V.a</a></li>
<li class="toctree-l1"><a class="reference internal" href="CSAS5App.html">Additional Example Applications of CSAS5</a></li>
<li class="toctree-l1"><a class="reference internal" href="CSAS6.html">CSAS6:  Control Module for Enhanced Criticality Safety Analysis with KENO-VI</a></li>
<li class="toctree-l1"><a class="reference internal" href="CSAS6App.html">Additional Example Applications of CSAS6</a></li>
<li class="toctree-l1"><a class="reference internal" href="STARBUCS.html">STARBUCS: A Scale Control Module for Automated Criticality Safety Analyses Using Burnup Credit</a></li>
<li class="toctree-l1"><a class="reference internal" href="Sourcerer.html">Sourcerer: Deterministic Starting Source for Criticality Calculations</a></li>
<li class="toctree-l1"><a class="reference internal" href="DEVC.html">DEVC: Denovo EigenValue Calculation</a></li>
<li class="toctree-l1"><a class="reference internal" href="KMART.html">KMART5 and KMART6: Postprocessors for KENO V.A and KENO-VI</a></li>
<li class="toctree-l1"><a class="reference internal" href="K5C5.html">K5toK6 and C5toC6: Input File Conversion Programs for KENO and CSAS</a></li>
</ul>
<p class="caption"><span class="caption-text">Material Specification and Cross Section Processing</span></p>
<ul class="current">
<li class="toctree-l1"><a class="reference internal" href="Material%20Specification%20and%20Cross%20Section%20Processing%20Overview.html">Material Specification and Cross Section Processing Overview</a></li>
<li class="toctree-l1"><a class="reference internal" href="XSProc.html">XSPROC: The Material and Cross Section Processing Module for SCALE</a></li>
<li class="toctree-l1"><a class="reference internal" href="XSProcAppA.html">XSProc: Standard Composition Examples</a></li>
<li class="toctree-l1"><a class="reference internal" href="XSProcAppB.html">XSProc Standard Composition Examples</a></li>
<li class="toctree-l1"><a class="reference internal" href="XSProcAppC.html">Examples of Complete XSProc Input Data</a></li>
<li class="toctree-l1"><a class="reference internal" href="stdcmp.html">Standard Composition Library</a></li>
<li class="toctree-l1"><a class="reference internal" href="BONAMI.html">BONAMI: Resonance Self-Shielding by the Bondarenko Method</a></li>
<li class="toctree-l1"><a class="reference internal" href="CENTRM.html">CENTRM: A Neutron Transport Code for Computing Continuous-Energy Spectra in General One-Dimensional Geometries and Two-Dimensional Lattice Cells</a></li>
<li class="toctree-l1"><a class="reference internal" href="PMC.html">PMC: A Program to Produce Multigroup Cross Sections Using Pointwise Energy Spectra from CENTRM</a></li>
<li class="toctree-l1"><a class="reference internal" href="PMCAppAB.html">PMC Appendices A and B</a></li>
<li class="toctree-l1"><a class="reference internal" href="CHOPS.html">CHOPS: Module to Compute Pointwise Disadvantage Factors and Produce a Cell-Homogenized CENTRM Library</a></li>
<li class="toctree-l1"><a class="reference internal" href="CRAWDAD.html">CRAWDAD: Module to Produce CENTRM-Formatted Continuous-Energy Nuclear Data Libraries</a></li>
<li class="toctree-l1 current"><a class="current reference internal" href="#">MCDancoff Data Guide</a><ul>
<li class="toctree-l2"><a class="reference internal" href="#introduction">Introduction</a></li>
<li class="toctree-l2"><a class="reference internal" href="#input-data-description">Input data description</a></li>
<li class="toctree-l2"><a class="reference internal" href="#calculation-and-use-of-3d-dancoff-factors">Calculation and use of 3D Dancoff factors</a></li>
<li class="toctree-l2"><a class="reference internal" href="#example-case">Example Case</a></li>
</ul>
</li>
<li class="toctree-l1"><a class="reference internal" href="CAJUN.html">CAJUN: Module for Combining and Manipulating CENTRM Continuous-Energy Libraries</a></li>
</ul>
<p class="caption"><span class="caption-text">Monte Carlo Transport</span></p>
<ul>
<li class="toctree-l1"><a class="reference internal" href="Monte%20Carlo%20Transport%20Overview.html">Monte Carlo Transport Overview</a></li>
<li class="toctree-l1"><a class="reference internal" href="Keno.html">Keno: A Monte Carlo Criticality Program</a></li>
<li class="toctree-l1"><a class="reference internal" href="KenoA.html">Keno Appendix A: KENO V.a Shape Descriptions</a></li>
<li class="toctree-l1"><a class="reference internal" href="KenoB.html">Keno Appendix B: KENO VI Shape Descriptions</a></li>
<li class="toctree-l1"><a class="reference internal" href="KenoC.html">Keno Appendix C: Sample problems</a></li>
<li class="toctree-l1"><a class="reference internal" href="Monaco.html">Monaco: A Fixed-Source Monte Carlo Transport Code for Shielding Applications</a></li>
</ul>
<p class="caption"><span class="caption-text">Radiation Shielding</span></p>
<ul>
<li class="toctree-l1"><a class="reference internal" href="MAVRIC.html">MAVRIC: Monaco with Automated Variance Reduction using Importance Calculations</a></li>
<li class="toctree-l1"><a class="reference internal" href="CAAScapability.html">MAVRIC Appendix A: CAAS Capability</a></li>
<li class="toctree-l1"><a class="reference internal" href="appendixb.html">MAVRIC Appendix B: MAVRIC Utilities</a></li>
<li class="toctree-l1"><a class="reference internal" href="appendixc.html">MAVRIC Appendix C: Advanced Features</a></li>
</ul>

            
          
        </div>
        
      </div>
    </nav>

    <section data-toggle="wy-nav-shift" class="wy-nav-content-wrap">

      
      <nav class="wy-nav-top" aria-label="top navigation">
        
          <i data-toggle="wy-nav-top" class="fa fa-bars"></i>
          <a href="index.html">SCALE Manual</a>
        
      </nav>


      <div class="wy-nav-content">
        
        <div class="rst-content">
        
          















<div role="navigation" aria-label="breadcrumbs navigation">

  <ul class="wy-breadcrumbs">
    
      <li><a href="index.html" class="icon icon-home"></a> &raquo;</li>
        
      <li>MCDancoff Data Guide</li>
    
    
      <li class="wy-breadcrumbs-aside">
        
            
            <a href="_sources/MCDancoff.rst.txt" rel="nofollow"> View page source</a>
          
        
      </li>
    
  </ul>

  
  <hr/>
</div>
          <div role="main" class="document" itemscope="itemscope" itemtype="http://schema.org/Article">
           <div itemprop="articleBody">
            
  <div class="section" id="mcdancoff-data-guide">
<span id="id1"></span><h1>MCDancoff Data Guide<a class="headerlink" href="#mcdancoff-data-guide" title="Permalink to this headline"></a></h1>
<p><em>L. M. Petrie, B. T. Rearden</em></p>
<p>ABSTRACT</p>
<p>The MCDancoff program is used to calculate Dancoff factors in
complicated, three-dimensional (3-D) geometries using Monte Carlo
integrations. The geometries are standard SCALE geometry descriptions,
with the current restriction that Dancoff factors can only be calculated
for regions bounded by cuboids, spheres, or cylinders. Multiple Dancoff
factors can be calculated with one input file.</p>
<p>ACKNOWLEDGMENT</p>
<p>This work was sponsored in part by Atomic Energy of Canada, Ltd. The
contribution of S. J. Poarch in preparing this document is gratefully
acknowledged.</p>
<div class="section" id="introduction">
<span id="id2"></span><h2>Introduction<a class="headerlink" href="#introduction" title="Permalink to this headline"></a></h2>
<p>MCDancoff (Monte Carlo Dancoff) is a program that calculates Dancoff
factors for complicated, three-dimensional geometries. Its input is a
slight modification of a CSAS6 input file which uses the standard SCALE
geometry as detailed for KENO-VI. The modifications to the input involve
different input in the START data block describing which Dancoff factors
are to be calculated. The calculation involves starting histories
isotropically on the surface of the region for which the Dancoff factor
is to be calculated and following the path of each history until it has
encountered all the elements of the material in the region, or until it
has exited the system. A one group cross-section library is used to
determine the total cross sections of the mixtures in the problem.</p>
<p>A current restriction of MCDancoff is that it can only calculate Dancoff
factors for regions bounded by cylinders, spheres, or cuboids. Other
simple bodies could be added in the future, but a general bounding
surface would be impractical.</p>
<p>The Dancoff factors are used in SCALE to correctly self-shield
multigroup cross sections for a given problem; either as input to BONAMI
or to determine an equivalent cell for CENTRM. This is most typically
accomplished through the MORE DATA and CENTRM DATA blocks.</p>
<p>The Dancoff factors are actually calculated by a modified version of the
KENO-VI code called KENO_Dancoff. All printed output from these
calculations is suppressed by default. If there is a need to see this
output (for example, to find an error message), it can be turned on by
setting an environment variable <strong>print_dancoff=yes</strong>.</p>
</div>
<div class="section" id="input-data-description">
<span id="id3"></span><h2>Input data description<a class="headerlink" href="#input-data-description" title="Permalink to this headline"></a></h2>
<p>MCDancoff input data is the same as CSAS6 input data with the following
exceptions. A special one group cross-section library will be used. It
can be specified as <strong>xn01</strong> in the input but will be set to this if
anything else is entered for the library. Because MCDancoff is running a
fixed source problem, and the Dancoff factor doesn’t need to be
calculated with the same accuracy as an eigenvalue, there are useful
changes that can be made to the parameters in the PARAMETER data block.
<a class="reference internal" href="#id4"><span class="std std-ref">Calculation and use of 3D Dancoff factors</span></a> discusses this in more detail. Finally, the START data
block is used to define which Dancoff factors will be calculated. This
data block is defined below.</p>
<p><strong>READ START</strong> Begins the data block</p>
<ol class="arabic simple">
<li><p><strong>dancoff</strong></p></li>
</ol>
<p>begins defining a new Dancoff factor. Always start
relative to the global unit in the geometry.</p>
<ol class="arabic simple" start="2">
<li><p><strong>array</strong></p></li>
</ol>
<p>step into an array contained in the current unit – followed
by <strong>karray</strong>, <strong>nbx</strong>, <strong>nby</strong>, <strong>nbz</strong> where <strong>karray</strong> is the
region containing the array in the current unit, <strong>nbx</strong> is the x
position in the array of the next unit, <strong>nby</strong> is the y position
in the array of the next unit, and <strong>nbz</strong> is the z position in
the array of the next unit.</p>
<ol class="arabic simple" start="3">
<li><p><strong>hole</strong></p></li>
</ol>
<p>step into a hole contained in the current unit – followed by
<strong>nhole</strong>, the hole number relative to the current unit.</p>
<ol class="arabic simple" start="4">
<li><p><strong>unit</strong></p></li>
</ol>
<p>final unit in the nesting chain – followed by <strong>nn</strong>, the
unit number</p>
<ol class="arabic simple" start="5">
<li><p><strong>region</strong></p></li>
</ol>
<p>region to calculate the Dancoff factor for – followed by
<strong>k</strong>, the relative geometry word in unit <strong>nn</strong> defining the
outer bound of the region.</p>
<ol class="arabic simple" start="6">
<li><p><strong>nst</strong></p></li>
</ol>
<p>if input, must be 0 (defaults to 0).</p>
<p>Repeat 2 and 3 to get from the global unit to the final unit <strong>nn</strong>.</p>
<p>Repeat 1–5 for each Dancoff factor to be calculated.</p>
<p><strong>END START</strong> Ends the data block</p>
</div>
<div class="section" id="calculation-and-use-of-3d-dancoff-factors">
<span id="id4"></span><h2>Calculation and use of 3D Dancoff factors<a class="headerlink" href="#calculation-and-use-of-3d-dancoff-factors" title="Permalink to this headline"></a></h2>
<ol class="arabic simple">
<li><p>The 3-D Dancoff factors are computed with KENO-VI geometry. If
beginning with CSAS5 model, use C5TOC6 to convert to CSAS6.</p></li>
<li><p>Change sequence name from CSAS6 to MCDancoff and change cross-section
library to <strong>xn01</strong>.</p></li>
<li><p>Input appropriate parameter data.</p></li>
</ol>
<blockquote>
<div><p>Since the Dancoff calculation is fixed source integration, there is
no need to skip generations, and <strong>nsk</strong> should be set to 0. Since
small changes to the Dancoff have very minor effects on the cross
sections, fewer histories are probably needed for calculating the
Dancoff than for calculating <em>keff</em>. Thirty thousand histories
divided as 100 generations of 300 histories per generation has
produced Dancoff factors with deviations of less than 1 percent. It
may be advantageous to turn off plots at this point. Since the same
parameters can be entered more than once, with the final entry being
the one used, adding a separate record with these values immediately
before the <strong>end parameter</strong> keywords would override the original
KENO-VI parameters.</p>
<p>Example:</p>
</div></blockquote>
<ol class="arabic simple" start="4">
<li><p>Identify the region for which Dancoff factors are desired in START
data.</p></li>
</ol>
<blockquote>
<div><p>The start type needs to be set to <strong>0</strong> for the Dancoff calculation
(this is the default). All KENO-VI START data should be removed or
commented out by placing an apostrophe in column 1. Each region for
which a Dancoff calculation is desired then starts with the keyword
<strong>dancoff</strong>. This is followed by data that specify the relationship
of the global unit to the specific geometry description of the
region. If the region is nested inside an array, then the keyword,
<strong>array</strong>, is specified, followed by four integers. The first integer
is the indices of the media record specifying the array relative to
the current unit. The next 3 integers are the X, Y, and Z indexes of
the position of the next unit in the array. If the region is nested
in a hole, then the keyword, <strong>hole</strong>, is specified, followed by the
relative count of the correct hole in the unit. The preceding data
are repeated (in the correct nesting order starting with the global
unit) until reaching the unit where the region is located. Then the
keyword, <strong>unit</strong>, followed by the unit number is given, followed by
the keyword, <strong>region</strong>, followed by the relative index of the
geometry keyword describing the desired region with respect to that
unit. Currently, only cylinders, spheres, and cuboids are programmed
for calculating Dancoff factors.</p>
<p>Examples:</p>
<div class="highlight-python notranslate"><div class="highlight"><pre><span></span><span class="n">read</span> <span class="n">start</span>
  <span class="n">nst</span><span class="o">=</span><span class="mi">0</span>
  <span class="n">dancoff</span>    <span class="n">hole</span> <span class="mi">1</span>    <span class="n">unit</span><span class="o">=</span><span class="mi">1</span>  <span class="n">reg</span><span class="o">=</span><span class="mi">1</span>
<span class="n">end</span> <span class="n">start</span>

<span class="n">read</span> <span class="n">start</span>
   <span class="n">dancoff</span>  <span class="n">array</span>  <span class="mi">1</span> <span class="mi">1</span> <span class="mi">1</span> <span class="mi">1</span>  <span class="n">array</span>  <span class="mi">1</span> <span class="mi">17</span> <span class="mi">17</span> <span class="mi">2</span>  <span class="n">unit</span> <span class="mi">10</span>  <span class="n">region</span> <span class="mi">1</span>
<span class="n">end</span> <span class="n">start</span>
</pre></div>
</div>
</div></blockquote>
<ol class="arabic simple" start="5">
<li><p>Execute MCDANCOFF <em>filename.</em>input file like any other SCALE input
file.</p></li>
<li><p>Examine <em>filename</em>.dancoff file, which will contain Dancoff factors
for each nuclide in the specified region</p></li>
</ol>
<div class="highlight-python notranslate"><div class="highlight"><pre><span></span><span class="n">index</span>        <span class="n">nuclide</span>        <span class="n">dancoff</span>      <span class="n">deviation</span>
                   <span class="mi">1</span>          <span class="mi">92234</span>    <span class="mf">3.36340E-01</span>    <span class="mf">1.81134E-03</span>
                   <span class="mi">2</span>          <span class="mi">92235</span>    <span class="mf">3.36340E-01</span>    <span class="mf">1.81134E-03</span>
                   <span class="mi">3</span>          <span class="mi">92236</span>    <span class="mf">3.36340E-01</span>    <span class="mf">1.81134E-03</span>
                   <span class="mi">4</span>          <span class="mi">92238</span>    <span class="mf">3.36340E-01</span>    <span class="mf">1.81134E-03</span>
                   <span class="mi">5</span>           <span class="mi">8016</span>    <span class="mf">1.00000E+00</span>    <span class="mf">0.00000E+00</span>
</pre></div>
</div>
<ol class="arabic simple" start="7">
<li><p>Once all desired Dancoff factors are obtained, return to original model and</p></li>
</ol>
<blockquote>
<div><p>enter CENTRM DATA for each cell with dan2pitch(mix) specified.</p>
</div></blockquote>
<div class="highlight-python notranslate"><div class="highlight"><pre><span></span><span class="n">read</span> <span class="n">celldata</span>
  <span class="n">latticecell</span> <span class="n">triangpitch</span> <span class="n">fuelr</span><span class="o">=</span><span class="mf">0.633</span>  <span class="mi">1</span> <span class="n">gapr</span><span class="o">=</span><span class="mf">0.637</span> <span class="mi">0</span> <span class="n">cladr</span><span class="o">=</span><span class="mf">0.675</span> <span class="mi">10</span> <span class="n">hpitch</span><span class="o">=</span><span class="mf">0.867</span>  <span class="mi">14</span> <span class="n">end</span>
<span class="n">centrm</span> <span class="n">data</span>
  <span class="n">dan2pitch</span><span class="p">(</span><span class="mi">1</span><span class="p">)</span><span class="o">=</span><span class="mf">0.336</span>
<span class="n">end</span> <span class="n">centrm</span>
</pre></div>
</div>
<ol class="arabic simple" start="8">
<li><p>If executing TSUNAMI-3D, additional steps are necessary because
TSUNAMI‑3D does not treat the dan2pitch input parameter.</p></li>
</ol>
<blockquote>
<div><p>Return to the original TSUNAMI-3D input file and replace the sequence
name to “CSAS-MG PARM=CHECK” and delete all data after the unit cell
data to quickly obtain revised pitch values. (Note: CSAS will not
modify cell dimensions to more than 20 cm, so a revised moderator
density may need to be entered to obtain the desired Dancoff factor.)
Search for the word “<em>desired”</em> in output file to find new pitch
values for each cell.</p>
<div class="highlight-python notranslate"><div class="highlight"><pre><span></span><span class="n">unit</span> <span class="n">cell</span>  <span class="o">=</span>    <span class="mi">1</span>
  <span class="n">original</span> <span class="n">pitch</span>                <span class="o">=</span> <span class="mf">1.7340E+00</span>
  <span class="n">Dancoff</span> <span class="k">for</span> <span class="n">orig</span> <span class="n">pitch</span>        <span class="o">=</span> <span class="mf">2.9728E-01</span>
  <span class="n">desired</span> <span class="n">Dancoff</span>               <span class="o">=</span> <span class="mf">3.3600E-01</span>

<span class="n">pitch</span> <span class="n">to</span> <span class="n">produce</span> <span class="n">desired</span> <span class="n">Dancoff</span><span class="o">=</span> <span class="mf">1.6845E+00</span>
</pre></div>
</div>
</div></blockquote>
<ol class="arabic simple" start="9">
<li><p>Enter revised pitch and revised moderator density (for cell
calculation only, not for geometry model) in TSUNAMI model.</p></li>
</ol>
</div>
<div class="section" id="example-case">
<span id="id5"></span><h2>Example Case<a class="headerlink" href="#example-case" title="Permalink to this headline"></a></h2>
<p>The following is a contrived case to illustrate an input file using both
holes, arrays, and multiple sets of Dancoff factors (although both
factors apply to the same pin, so only one set can be used). The case
represents two fuel assemblies in a cylindrical tank, each assembly
having a poisoned central pin, and four water holes. The Dancoff factors
are calculated for each central pin. The input file is listed in
<a class="reference internal" href="#list7-8-1"><span class="std std-numref">Listing 15</span></a>.</p>
<div class="literal-block-wrapper docutils container" id="list7-8-1">
<div class="code-block-caption"><span class="caption-number">Listing 15 </span><span class="caption-text">Example input file (continued below).</span><a class="headerlink" href="#list7-8-1" title="Permalink to this code"></a></div>
<div class="highlight-scale notranslate"><div class="highlight"><pre><span></span><span class="nf">=mcdancoff</span>
<span class="err">sample</span> <span class="err">case</span> <span class="err">demonstrating</span> <span class="err">calculating</span> <span class="err">Dancoff</span> <span class="err">factors</span>
<span class="err">xn</span><span class="m">01</span>
<span class="n">read</span><span class="err"> </span><span class="n">composition</span>
  <span class="err">uo</span><span class="m">2</span>     <span class="m">1</span>  <span class="err">den=</span><span class="m">10</span><span class="err">.</span><span class="m">38</span>  <span class="m">1</span>  <span class="m">294</span>  <span class="m">92234</span> <span class="err">.</span><span class="m">0303</span>  <span class="m">92235</span> <span class="m">4</span><span class="err">.</span><span class="m">7378</span>  <span class="m">92236</span> <span class="err">.</span><span class="m">1364</span>  <span class="m">92238</span> <span class="m">95</span><span class="err">.</span><span class="m">0955</span>  <span class="n">end</span><span class="err"> </span><span class="n">uo</span><span class="err">2</span>
  <span class="err">zirc</span><span class="m">4</span>   <span class="m">2</span>             <span class="m">1</span>  <span class="m">294</span>  <span class="n">end</span><span class="err"> </span><span class="n">zirc</span><span class="err">4</span>
  <span class="err">h</span><span class="m">2</span><span class="err">o</span>     <span class="m">3</span>             <span class="m">1</span>  <span class="m">294</span>  <span class="n">end</span><span class="err"> </span><span class="n">h</span><span class="err">2</span><span class="n">o</span>
  <span class="err">uo</span><span class="m">2</span>     <span class="m">4</span>  <span class="err">den=</span><span class="m">10</span><span class="err">.</span><span class="m">08</span>  <span class="m">1</span>  <span class="m">294</span>  <span class="m">92234</span> <span class="err">.</span><span class="m">0303</span>  <span class="m">92235</span> <span class="m">4</span><span class="err">.</span><span class="m">7378</span>  <span class="m">92236</span> <span class="err">.</span><span class="m">1364</span>  <span class="m">92238</span> <span class="m">95</span><span class="err">.</span><span class="m">0955</span>  <span class="n">end</span><span class="err"> </span><span class="n">uo</span><span class="err">2</span>
  <span class="err">gd</span>      <span class="m">4</span>  <span class="err">den=</span> <span class="m">0</span><span class="err">.</span><span class="m">3</span>   <span class="m">1</span>  <span class="m">294</span>  <span class="n">end</span><span class="err"> </span><span class="n">gd</span>
<span class="n">end</span><span class="err"> </span><span class="n">composition</span>
<span class="n">read</span><span class="err"> </span><span class="n">param</span>
  <span class="err">nsk=</span><span class="m">0</span> <span class="err">gen=</span><span class="m">100</span> <span class="err">npg=</span><span class="m">300</span>
<span class="n">end</span><span class="err"> </span><span class="n">param</span>
<span class="n">read</span><span class="err"> </span><span class="n">geometry</span>
  <span class="err">unit</span> <span class="m">1</span>
    <span class="err">com=!fuel</span> <span class="err">pin!</span>
    <span class="err">cylinder</span>   <span class="m">10</span>  <span class="m">0</span><span class="err">.</span><span class="m">395</span>  <span class="m">40</span><span class="err">.</span><span class="m">0</span>  <span class="err">-</span><span class="m">40</span><span class="err">.</span><span class="m">0</span>
    <span class="err">cylinder</span>   <span class="m">20</span>  <span class="m">0</span><span class="err">.</span><span class="m">410</span>  <span class="m">40</span><span class="err">.</span><span class="m">0</span>  <span class="err">-</span><span class="m">40</span><span class="err">.</span><span class="m">0</span>
    <span class="err">cylinder</span>   <span class="m">30</span>  <span class="m">0</span><span class="err">.</span><span class="m">470</span>  <span class="m">40</span><span class="err">.</span><span class="m">0</span>  <span class="err">-</span><span class="m">40</span><span class="err">.</span><span class="m">0</span>
    <span class="err">cuboid</span>     <span class="m">40</span>  <span class="m">4</span><span class="err">p</span><span class="m">0</span><span class="err">.</span><span class="m">65</span>  <span class="m">2</span><span class="err">p</span><span class="m">40</span><span class="err">.</span><span class="m">0</span>
    <span class="err">media</span>      <span class="m">1</span>  <span class="m">1</span>  <span class="m">10</span>
    <span class="err">media</span>      <span class="m">0</span>  <span class="m">1</span>  <span class="m">20</span> <span class="err">-</span><span class="m">10</span>
    <span class="err">media</span>      <span class="m">2</span>  <span class="m">1</span>  <span class="m">30</span> <span class="err">-</span><span class="m">20</span>
    <span class="err">media</span>      <span class="m">3</span>  <span class="m">1</span>  <span class="m">40</span> <span class="err">-</span><span class="m">30</span>
    <span class="err">boundary</span>   <span class="m">40</span>
  <span class="err">unit</span> <span class="m">2</span>
    <span class="err">com=!water</span> <span class="err">hole!</span>
    <span class="err">cuboid</span>     <span class="m">40</span>  <span class="m">4</span><span class="err">p</span><span class="m">0</span><span class="err">.</span><span class="m">65</span>  <span class="m">2</span><span class="err">p</span><span class="m">40</span><span class="err">.</span><span class="m">0</span>
    <span class="err">media</span>      <span class="m">3</span>  <span class="m">1</span>  <span class="m">40</span>
    <span class="err">boundary</span>   <span class="m">40</span>
  <span class="err">unit</span> <span class="m">3</span>
    <span class="err">com=!unit</span> <span class="err">containing</span> <span class="err">a</span> <span class="m">2</span><span class="err">x</span><span class="m">2</span> <span class="err">array</span> <span class="err">of</span> <span class="err">fuel</span> <span class="err">pins!</span>
    <span class="err">cuboid</span>     <span class="m">10</span>  <span class="m">4</span><span class="err">p</span><span class="m">1</span><span class="err">.</span><span class="m">30</span>  <span class="m">2</span><span class="err">p</span><span class="m">40</span><span class="err">.</span><span class="m">0</span>
    <span class="err">array</span>      <span class="m">1</span>  <span class="m">10</span>  <span class="err">place</span> <span class="m">1</span> <span class="m">1</span> <span class="m">1</span> <span class="err">-</span><span class="m">0</span><span class="err">.</span><span class="m">65</span> <span class="err">-</span><span class="m">0</span><span class="err">.</span><span class="m">65</span> <span class="m">0</span><span class="err">.</span><span class="m">0</span>
    <span class="err">boundary</span>   <span class="m">10</span>
  <span class="err">unit</span> <span class="m">4</span>
    <span class="err">com=!unit</span> <span class="err">containing</span> <span class="err">a</span> <span class="m">1</span><span class="err">x</span><span class="m">2</span> <span class="err">array</span> <span class="err">of</span> <span class="err">fuel</span> <span class="err">pins!</span>
    <span class="err">cuboid</span>     <span class="m">10</span>  <span class="m">2</span><span class="err">p</span><span class="m">0</span><span class="err">.</span><span class="m">65</span>  <span class="m">2</span><span class="err">p</span><span class="m">1</span><span class="err">.</span><span class="m">30</span>  <span class="m">2</span><span class="err">p</span><span class="m">40</span><span class="err">.</span><span class="m">0</span>
    <span class="err">array</span>      <span class="m">2</span>  <span class="m">10</span>  <span class="err">place</span> <span class="m">1</span> <span class="m">1</span> <span class="m">1</span> <span class="m">0</span><span class="err">.</span><span class="m">0</span> <span class="err">-</span><span class="m">0</span><span class="err">.</span><span class="m">65</span> <span class="m">0</span><span class="err">.</span><span class="m">0</span>
    <span class="err">boundary</span>   <span class="m">10</span>
  <span class="err">unit</span> <span class="m">5</span>
    <span class="err">com=!unit</span> <span class="err">containing</span> <span class="err">a</span> <span class="m">2</span><span class="err">x</span><span class="m">1</span> <span class="err">array</span> <span class="err">of</span> <span class="err">fuel</span> <span class="err">pins!</span>
    <span class="err">cuboid</span>     <span class="m">10</span>  <span class="m">2</span><span class="err">p</span><span class="m">1</span><span class="err">.</span><span class="m">30</span>  <span class="m">2</span><span class="err">p</span><span class="m">0</span><span class="err">.</span><span class="m">65</span>  <span class="m">2</span><span class="err">p</span><span class="m">40</span><span class="err">.</span><span class="m">0</span>
    <span class="err">array</span>      <span class="m">3</span>  <span class="m">10</span>  <span class="err">place</span> <span class="m">1</span> <span class="m">1</span> <span class="m">1</span> <span class="err">-</span><span class="m">0</span><span class="err">.</span><span class="m">65</span> <span class="m">0</span><span class="err">.</span><span class="m">0</span> <span class="m">0</span><span class="err">.</span><span class="m">0</span>
    <span class="err">boundary</span>   <span class="m">10</span>
  <span class="err">unit</span> <span class="m">6</span>
    <span class="err">com=!unit</span> <span class="err">containing</span> <span class="err">a</span> <span class="m">5</span><span class="err">x</span><span class="m">5</span> <span class="err">array</span> <span class="err">of</span> <span class="err">fuel</span> <span class="err">pins!</span>
    <span class="err">cuboid</span>     <span class="m">10</span>  <span class="m">4</span><span class="err">p</span><span class="m">3</span><span class="err">.</span><span class="m">25</span>  <span class="m">2</span><span class="err">p</span><span class="m">40</span><span class="err">.</span><span class="m">0</span>
    <span class="err">array</span>      <span class="m">4</span>  <span class="m">10</span>  <span class="err">place</span> <span class="m">2</span> <span class="m">2</span> <span class="m">1</span> <span class="m">0</span><span class="err">.</span><span class="m">0</span> <span class="m">0</span><span class="err">.</span><span class="m">0</span> <span class="m">0</span><span class="err">.</span><span class="m">0</span>
    <span class="err">boundary</span>   <span class="m">10</span>
  <span class="err">unit</span> <span class="m">7</span>
    <span class="err">com=!unit</span> <span class="err">containing</span> <span class="err">a</span> <span class="m">5</span><span class="err">x</span><span class="m">5</span> <span class="err">array</span> <span class="err">of</span> <span class="err">fuel</span> <span class="err">pins</span> <span class="err">-</span> <span class="err">water</span> <span class="err">hole</span> <span class="err">in</span> <span class="err">the</span> <span class="err">middle!</span>
    <span class="err">cuboid</span>     <span class="m">10</span>  <span class="m">4</span><span class="err">p</span><span class="m">3</span><span class="err">.</span><span class="m">25</span>  <span class="m">2</span><span class="err">p</span><span class="m">40</span><span class="err">.</span><span class="m">0</span>
    <span class="err">array</span>      <span class="m">5</span>  <span class="m">10</span>  <span class="err">place</span> <span class="m">2</span> <span class="m">2</span> <span class="m">1</span> <span class="m">0</span><span class="err">.</span><span class="m">0</span> <span class="m">0</span><span class="err">.</span><span class="m">0</span> <span class="m">0</span><span class="err">.</span><span class="m">0</span>
    <span class="err">boundary</span>   <span class="m">10</span>
</pre></div>
</div>
</div>
<div class="highlight-python notranslate"><div class="highlight"><pre><span></span>unit 8
  com=!unit containing a 5x5 array of fuel pins - poisoned pin in the middle!
  cuboid     10  4p3.25  2p40.0
  array      6  10  place 2 2 1 0.0 0.0 0.0
  boundary   10
unit 9
  com=!poisoned fuel pin!
  cylinder   10  0.395  40.0  -40.0
  cylinder   20  0.410  40.0  -40.0
  cylinder   30  0.470  40.0  -40.0
  cuboid     40  4p0.65  2p40.0
  media      4  1  10
  media      0  1  20 -10
  media      2  1  30 -20
  media      3  1  40 -30
  boundary   40
unit 10
  com=!unit containing a 15x15 fuel assembly!
  cuboid     10  4p9.75  2p40.0
  array      7  10  place  2 2 1 0.0 0.0 0.0
  boundary   10
global
unit 11
  com=!global unit with 2 fuel assemblies!
  cylinder   10  25.0  60.0  -60.0
  hole       10  origin x=-10.0
  hole       10  origin x= 10.0
  media      3  1  10
  boundary   10
end geometry
read array
ara=1 typ=square nux=2 nuy=2 nuz=1  fill f1 end fill
ara=2 typ=square nux=1 nuy=2 nuz=1  fill f1 end fill
ara=3 typ=square nux=2 nuy=1 nuz=1  fill f1 end fill
ara=4 typ=square nux=3 nuy=3 nuz=1  fill 3 4 3 5 1 5 3 4 3 end fill
ara=5 typ=square nux=3 nuy=3 nuz=1  fill 3 4 3 5 2 5 3 4 3 end fill
ara=6 typ=square nux=3 nuy=3 nuz=1  fill 3 4 3 5 9 5 3 4 3 end fill
ara=7 typ=square nux=3 nuy=3 nuz=1  fill 7 6 7 6 8 6 7 6 7 end fill
end array
read start
&#39; first Dancoff - calculate for the poisoned fuel pin in unit 9 for the x=-10 assembly
dancoff
&#39; hole 1 is unit 10 at x=-10
  hole 1
&#39; array in first region of unit 10 is array 7 - 2 2 1 position is unit 8
  array 1 2 2 1
&#39; array in first region of unit 8 is array 6 - 2 2 1 position is unit 9
  array 1 2 2 1
&#39;  cylinder labeled 10 in unit 9 is the first region
  unit 9   region 1
&#39; second Dancoff - calculate for the poisoned fuel pin in unit 9 for the x=+10 assembly
dancoff
&#39; hole 2 is unit 10 at x=+10
  hole 2
&#39; array in first region of unit 10 is array 7 - 2 2 1 position is unit 8
  array 1 2 2 1
&#39; array in first region of unit 8 is array 6 - 2 2 1 position is unit 9
  array 1 2 2 1
&#39;  cylinder labeled 10 in unit 9 is the first region
  unit 9   region 1
end start
end data
end
</pre></div>
</div>
<p>This input file creates two files of Dancoff factors.  The first such file is listed in <a class="reference internal" href="#list7-8-2"><span class="std std-numref">Listing 16</span></a>.</p>
<div class="literal-block-wrapper docutils container" id="list7-8-2">
<div class="code-block-caption"><span class="caption-number">Listing 16 </span><span class="caption-text">Output file of Dancoff factors.</span><a class="headerlink" href="#list7-8-2" title="Permalink to this code"></a></div>
<div class="highlight-scale notranslate"><div class="highlight"><pre><span></span><span class="err">Unit</span> <span class="m">9</span> <span class="err">at</span> <span class="err">global</span> <span class="err">x</span> <span class="err">-</span><span class="m">1</span><span class="err">.</span><span class="m">00000</span><span class="err">E+</span><span class="m">01</span>  <span class="err">y</span>  <span class="m">0</span><span class="err">.</span><span class="m">00000</span><span class="err">E+</span><span class="m">00</span>  <span class="err">z</span><span class="m">  0.00000E+00</span>
          <span class="err">index</span>        <span class="err">nuclide</span>        <span class="err">dancoff</span>      <span class="err">deviation</span>
              <span class="m">1</span>          <span class="m">92234</span>    <span class="m">2</span><span class="err">.</span><span class="m">20873</span><span class="err">E-</span><span class="m">01    1.03436E-03</span>
              <span class="m">2</span>          <span class="m">92235</span>    <span class="m">2</span><span class="err">.</span><span class="m">20873</span><span class="err">E-</span><span class="m">01    1.03436E-03</span>
              <span class="m">3</span>          <span class="m">92236</span>    <span class="m">2</span><span class="err">.</span><span class="m">20873</span><span class="err">E-</span><span class="m">01    1.03436E-03</span>
              <span class="m">4</span>          <span class="m">92238</span>    <span class="m">2</span><span class="err">.</span><span class="m">20873</span><span class="err">E-</span><span class="m">01    1.03436E-03</span>
              <span class="m">5</span>           <span class="m">8016</span>    <span class="m">9</span><span class="err">.</span><span class="m">64748</span><span class="err">E-</span><span class="m">01    4.28121E-04</span>
              <span class="m">6</span>          <span class="m">64000</span>    <span class="m">2</span><span class="err">.</span><span class="m">82254</span><span class="err">E-</span><span class="m">10    3.61320E-11</span>
</pre></div>
</div>
</div>
<p>The second file is statistically the same, as it solved for the mirror image Dancoff factor.</p>
</div>
</div>


           </div>
           
          </div>
          <footer>
  
    <div class="rst-footer-buttons" role="navigation" aria-label="footer navigation">
      
        <a href="CAJUN.html" class="btn btn-neutral float-right" title="CAJUN: Module for Combining and Manipulating CENTRM Continuous-Energy Libraries" accesskey="n" rel="next">Next <span class="fa fa-arrow-circle-right"></span></a>
      
      
        <a href="CRAWDAD.html" class="btn btn-neutral float-left" title="CRAWDAD: Module to Produce CENTRM-Formatted Continuous-Energy Nuclear Data Libraries" accesskey="p" rel="prev"><span class="fa fa-arrow-circle-left"></span> Previous</a>
      
    </div>
  

  <hr/>

  <div role="contentinfo">
    <p>
        
        &copy; Copyright 2020, SCALE developers

    </p>
  </div>
    
    
    
    Built with <a href="http://sphinx-doc.org/">Sphinx</a> using a
    
    <a href="https://github.com/rtfd/sphinx_rtd_theme">theme</a>
    
    provided by <a href="https://readthedocs.org">Read the Docs</a>. 

</footer>

        </div>
      </div>

    </section>

  </div>
  

  <script type="text/javascript">
      jQuery(function () {
          SphinxRtdTheme.Navigation.enable(true);
      });
  </script>

  
  
    
   

</body>
</html>