Discountinuities on surface net current and flux

Report any suspected bugs and unexpected behavior here
Post Reply
yunlin
Posts: 4
Joined: Thu Dec 12, 2019 4:58 pm
Security question 1: No
Security question 2: 7

Discountinuities on surface net current and flux

Post by yunlin » Fri Jun 12, 2020 5:45 am

The surface fluxes and currents on same surface shared by two neighboring GCU have different values.
Is it possible to fix this problem?

Antti Rintala
Posts: 85
Joined: Tue Jul 03, 2012 1:25 pm
Security question 1: No
Security question 2: 92

Re: Discountinuities on surface net current and flux

Post by Antti Rintala » Fri Jun 12, 2020 8:03 am

Could you provide more information on your problem? It is impossible to say anything based on this report. What values did you compare? Detector values or values provided by set adf card? What are the boundary surfaces of the to universes to be homogenized? Are you comparing heterogeneous or homogeneous surface fluxes? The surface currents and heterogeneous surface fluxes should have same values on a shared boundary.

yunlin
Posts: 4
Joined: Thu Dec 12, 2019 4:58 pm
Security question 1: No
Security question 2: 7

Re: Discountinuities on surface net current and flux

Post by yunlin » Sat Jun 13, 2020 1:24 am

Values provided by set adf card.
here are lines copied from output,
DF_SURFACE (idx, [1: 3]) = '102' ;
DF_SURF_NET_CURR (idx, [1: 24]) = [ 0.00000E+00 0.0E+00 0.00000E+00 0.0E+00 0.00000E+00 0.0E+00 0.00000E+00 0.0E+00 2.06949E-05 0.40325 6.80538E-05 0.34425 6.11223E-06 1.00000 -2.91043E-05 0.72128 -6.21219E-04 0.00484 -5.26589E-03 0.00252 3.00178E-03 0.00315 4.64663E-03 0.00358 ];
DF_SURFACE (idx, [1: 3]) = '103' ;
DF_SURF_NET_CURR (idx, [1: 24]) = [ 0.00000E+00 0.0E+00 0.00000E+00 0.0E+00 0.00000E+00 0.0E+00 0.00000E+00 0.0E+00 4.70006E-04 0.03288 3.04060E-04 0.09167 2.95974E-04 0.05989 1.48080E-04 0.17504 -3.00177E-03 0.00315 -4.64661E-03 0.00358 1.02640E-02 0.00161 -1.35603E-03 0.01398 ];

The last two net currents of '102', 3.00178E-03 and 4.64663E-03 are for top of cell.
the corresponding curents of '103', -3.00177E-03 and -4.64661E-03.

there are not too much difference, but can be more significant for other cases.

I will provide my input here
set title "3D Simplified Fuel Assembly for 1 by 1 Mesh XS Generation"
%----------------------------------------------------------------------------------------
% --- SURFACES

% 1 = Fuel
% 2 = Fuel Can Inner
% 3 = Fuel Can/End Tube Outer
surf 1 octa 0.0 0.0 4.82600 6.04775
surf 2 octa 0.0 0.0 4.96570 6.25511
surf 3 octa 0.0 0.0 5.02920 6.31861

% 4 = Short Plug Hole
% 5 = Short Plug Narrow (Square)
% 6 = Short Plug/Long Plug
% 7 = End Tube Inner
surf 4 cyl 0.0 0.0 0.95250
surf 5 sqc 0.0 0.0 4.16560 0.00000
surf 6 octa 0.0 0.0 4.80060 6.02072
surf 7 octa 0.0 0.0 4.90220 6.19161

% 8 = CR hole inner
% 9 = CR hole outer for bottom ref and fuel section
% 10 = Bearing tube outer
% 11 = Bushing Inner
% 12 = Lower Plug CR Inner
% 13 = FUEL block Inner
% 14 = Upper Plug CR Inner
% 15 = CR FUEL outgas tube
% 16 = CR rod follower inner
% 17 = CR rod follower outer
%
surf 8 cyl 0.0 0.0 2.54000
surf 9 cyl 0.0 0.0 2.85750
surf 10 cyl 0.0 0.0 3.17500
surf 11 cyl 0.0 0.0 2.25552
surf 12 cyl 0.0 0.0 2.87734
surf 13 cyl 0.0 0.0 2.95910
surf 14 cyl 0.0 0.0 3.25755
surf 15 cyl 3.14309 3.14309 0.68937
surf 16 cyl 0.0 0.0 1.90500
surf 17 cyl 0.0 0.0 2.22250
%
surf 99 inf
%
%basic blocks:
%
% fuel block
%
cell FA14-FB 91 FUEL -1
cell FA14-Gap 91 AIR 1 -2
cell FA14-CAN 91 ZR-3 2 -3
cell FA14-AIRG 91 AIR 3
% reflector long plug (Graphite in Al tube)
cell FA11-LSP 92 GRPH-CP2 -1
cell FA11-Gap 92 AIR 1 -2
cell FA11-ENDT 92 AL-6063 2 -3
cell FA11-AIRG 92 AIR 3
%
% follower in Reflector
%
cell CS13-FOL 93 GRPH-CP2 -16
cell CS13-ZRC 93 ZR-2 16 -17
cell CS13-CAIR 93 AIR 17 -8
cell CS13-CTUB 93 ZR-3 8 -9
cell CS13-TAIR 93 AIR 9 -12
cell CS13-LSP 93 GRPH-CP2 12 -6
cell CS13-Gap 93 AIR 6 -7
cell CS13-ENDT 93 AL-6063 7 -3
cell CS13-AIRG 93 AIR 3
%
% Flower in Fuel Block
%
cell CS14-FOL 94 GRPH-CP2 -16
cell CS14-ZRC 94 ZR-2 16 -17
cell CS14-CAIR 94 AIR 17 -8
cell CS14-CTUB 94 ZR-3 8 -9
cell CS14-TAIR 94 AIR 9 -13
cell CS14-FB 94 FUEL 13 -1
cell CS14-Gap 94 AIR 1 -2
cell CS14-CAN 94 ZR-3 2 -3
cell CS14-AIRG 94 AIR 3
%
% Control rod in Fuel Block
%
cell CS15-FOL 95 B4C -16
cell CS15-ZRC 95 M-STEEL 16 -17
cell CS15-CAIR 95 AIR 17 -8
cell CS15-CTUB 95 ZR-3 8 -9
cell CS15-TAIR 95 AIR 9 -13
cell CS15-FB 95 FUEL 13 -1
cell CS15-Gap 95 AIR 1 -2
cell CS15-CAN 95 ZR-3 2 -3
cell CS15-AIRG 95 AIR 3
%
% Control rod in Reflector
%
cell CS16-FOL 96 B4C -16
cell CS16-ZRC 96 M-STEEL 16 -17
cell CS16-CAIR 96 AIR 17 -8
cell CS16-CTUB 96 ZR-3 8 -9
cell CS16-TAIR 96 AIR 9 -12
cell CS16-LSP 96 GRPH-CP2 12 -6
cell CS16-Gap 96 AIR 6 -7
cell CS16-ENDT 96 AL-6063 7 -3
cell CS16-AIRG 96 AIR 3
%
% air block
%
cell FA54-Gap 97 AIR -99
%
%

% Define assemblies

% assembly 1
cell 101 101 fill 92 -99
cell 102 102 fill 92 -99
cell 103 103 fill 92 -99
cell 104 104 fill 91 -99
cell 105 105 fill 91 -99
cell 106 106 fill 91 -99
cell 107 107 fill 91 -99
cell 108 108 fill 91 -99
cell 109 109 fill 91 -99
cell 110 110 fill 92 -99
cell 111 111 fill 92 -99
cell 112 112 fill 92 -99
lat 1 9 0 0 12
0.00000 101
20.66396 102
41.32792 103
61.99188 104
82.15313 105
102.31438 106
122.47563 107
142.63688 108
162.79813 109
182.95938 110
204.15251 111
225.34564 112
% assembly 2
cell 201 201 fill 92 -99
cell 202 202 fill 92 -99
cell 203 203 fill 92 -99
cell 204 204 fill 91 -99
cell 205 205 fill 91 -99
cell 206 206 fill 91 -99
cell 207 207 fill 91 -99
cell 208 208 fill 91 -99
cell 209 209 fill 91 -99
cell 210 210 fill 92 -99
cell 211 211 fill 92 -99
cell 212 212 fill 92 -99
lat 2 9 0 0 12
0.00000 201
20.66396 202
41.32792 203
61.99188 204
82.15313 205
102.31438 206
122.47563 207
142.63688 208
162.79813 209
182.95938 210
204.15251 211
225.34564 212
% assembly 3
cell 301 301 fill 92 -99
cell 302 302 fill 92 -99
cell 303 303 fill 92 -99
cell 304 304 fill 91 -99
cell 305 305 fill 91 -99
cell 306 306 fill 91 -99
cell 307 307 fill 91 -99
cell 308 308 fill 91 -99
cell 309 309 fill 91 -99
cell 310 310 fill 92 -99
cell 311 311 fill 92 -99
cell 312 312 fill 92 -99
lat 3 9 0 0 12
0.00000 301
20.66396 302
41.32792 303
61.99188 304
82.15313 305
102.31438 306
122.47563 307
142.63688 308
162.79813 309
182.95938 310
204.15251 311
225.34564 312
% assembly 4
cell 401 401 fill 92 -99
cell 402 402 fill 92 -99
cell 403 403 fill 92 -99
cell 404 404 fill 91 -99
cell 405 405 fill 91 -99
cell 406 406 fill 91 -99
cell 407 407 fill 91 -99
cell 408 408 fill 91 -99
cell 409 409 fill 91 -99
cell 410 410 fill 92 -99
cell 411 411 fill 92 -99
cell 412 412 fill 92 -99
lat 4 9 0 0 12
0.00000 401
20.66396 402
41.32792 403
61.99188 404
82.15313 405
102.31438 406
122.47563 407
142.63688 408
162.79813 409
182.95938 410
204.15251 411
225.34564 412
% assembly 5
cell 501 501 fill 92 -99
cell 502 502 fill 92 -99
cell 503 503 fill 92 -99
cell 504 504 fill 97 -99
cell 505 505 fill 97 -99
cell 506 506 fill 97 -99
cell 507 507 fill 97 -99
cell 508 508 fill 97 -99
cell 509 509 fill 97 -99
cell 510 510 fill 92 -99
cell 511 511 fill 92 -99
cell 512 512 fill 92 -99
lat 5 9 0 0 12
0.00000 501
20.66396 502
41.32792 503
61.99188 504
82.15313 505
102.31438 506
122.47563 507
142.63688 508
162.79813 509
182.95938 510
204.15251 511
225.34564 512
% assembly 6
cell 601 601 fill 92 -99
cell 602 602 fill 92 -99
cell 603 603 fill 92 -99
cell 604 604 fill 97 -99
cell 605 605 fill 97 -99
cell 606 606 fill 97 -99
cell 607 607 fill 97 -99
cell 608 608 fill 97 -99
cell 609 609 fill 97 -99
cell 610 610 fill 92 -99
cell 611 611 fill 92 -99
cell 612 612 fill 92 -99
lat 6 9 0 0 12
0.00000 601
20.66396 602
41.32792 603
61.99188 604
82.15313 605
102.31438 606
122.47563 607
142.63688 608
162.79813 609
182.95938 610
204.15251 611
225.34564 612
% assembly 7
cell 701 701 fill 92 -99
cell 702 702 fill 92 -99
cell 703 703 fill 92 -99
cell 704 704 fill 91 -99
cell 705 705 fill 91 -99
cell 706 706 fill 91 -99
cell 707 707 fill 91 -99
cell 708 708 fill 91 -99
cell 709 709 fill 91 -99
cell 710 710 fill 92 -99
cell 711 711 fill 92 -99
cell 712 712 fill 92 -99
lat 7 9 0 0 12
0.00000 701
20.66396 702
41.32792 703
61.99188 704
82.15313 705
102.31438 706
122.47563 707
142.63688 708
162.79813 709
182.95938 710
204.15251 711
225.34564 712
% assembly 8
cell 801 801 fill 92 -99
cell 802 802 fill 92 -99
cell 803 803 fill 92 -99
cell 804 804 fill 91 -99
cell 805 805 fill 91 -99
cell 806 806 fill 91 -99
cell 807 807 fill 91 -99
cell 808 808 fill 91 -99
cell 809 809 fill 91 -99
cell 810 810 fill 92 -99
cell 811 811 fill 92 -99
cell 812 812 fill 92 -99
lat 8 9 0 0 12
0.00000 801
20.66396 802
41.32792 803
61.99188 804
82.15313 805
102.31438 806
122.47563 807
142.63688 808
162.79813 809
182.95938 810
204.15251 811
225.34564 812
% assembly 9
cell 901 901 fill 92 -99
cell 902 902 fill 92 -99
cell 903 903 fill 92 -99
cell 904 904 fill 91 -99
cell 905 905 fill 91 -99
cell 906 906 fill 91 -99
cell 907 907 fill 91 -99
cell 908 908 fill 91 -99
cell 909 909 fill 91 -99
cell 910 910 fill 92 -99
cell 911 911 fill 92 -99
cell 912 912 fill 92 -99
lat 9 9 0 0 12
0.00000 901
20.66396 902
41.32792 903
61.99188 904
82.15313 905
102.31438 906
122.47563 907
142.63688 908
162.79813 909
182.95938 910
204.15251 911
225.34564 912

% core configuration
lat 100 1 0.0 0.0 3 3 10.16000
1 2 3
4 5 6
7 8 9
surf 100 cuboid -15.24000 15.24000 -15.24000 15.24000 0.00000 246.53877
cell a1 0 fill 100 -100
cell outrf 0 outside 100
%
% define gcu
set gcu
101 201 301 401 501 601 701 801 901
102 202 302 402 502 602 702 802 902
103 203 303 403 503 603 703 803 903
104 204 304 404 504 604 704 804 904
105 205 305 405 505 605 705 805 905
106 206 306 406 506 606 706 806 906
107 207 307 407 507 607 707 807 907
108 208 308 408 508 608 708 808 908
109 209 309 409 509 609 709 809 909
110 210 310 410 510 610 710 810 910
111 211 311 411 511 611 711 811 911
112 212 312 412 512 612 712 812 912

%define surfaces for ADF
% surfaces for DF of assembly 1
surf 101 cuboid -15.24000 -5.08000 -15.24000 -5.08000 0.00000 20.66396
surf 102 cuboid -15.24000 -5.08000 -15.24000 -5.08000 20.66396 41.32792
surf 103 cuboid -15.24000 -5.08000 -15.24000 -5.08000 41.32792 61.99188
surf 104 cuboid -15.24000 -5.08000 -15.24000 -5.08000 61.99188 82.15313
surf 105 cuboid -15.24000 -5.08000 -15.24000 -5.08000 82.15313 102.31438
surf 106 cuboid -15.24000 -5.08000 -15.24000 -5.08000 102.31438 122.47563
surf 107 cuboid -15.24000 -5.08000 -15.24000 -5.08000 122.47563 142.63688
surf 108 cuboid -15.24000 -5.08000 -15.24000 -5.08000 142.63688 162.79813
surf 109 cuboid -15.24000 -5.08000 -15.24000 -5.08000 162.79813 182.95938
surf 110 cuboid -15.24000 -5.08000 -15.24000 -5.08000 182.95938 204.15251
surf 111 cuboid -15.24000 -5.08000 -15.24000 -5.08000 204.15251 225.34564
surf 112 cuboid -15.24000 -5.08000 -15.24000 -5.08000 225.34564 246.53877
% surfaces for DF of assembly 2
surf 201 cuboid -5.08000 5.08000 -15.24000 -5.08000 0.00000 20.66396
surf 202 cuboid -5.08000 5.08000 -15.24000 -5.08000 20.66396 41.32792
surf 203 cuboid -5.08000 5.08000 -15.24000 -5.08000 41.32792 61.99188
surf 204 cuboid -5.08000 5.08000 -15.24000 -5.08000 61.99188 82.15313
surf 205 cuboid -5.08000 5.08000 -15.24000 -5.08000 82.15313 102.31438
surf 206 cuboid -5.08000 5.08000 -15.24000 -5.08000 102.31438 122.47563
surf 207 cuboid -5.08000 5.08000 -15.24000 -5.08000 122.47563 142.63688
surf 208 cuboid -5.08000 5.08000 -15.24000 -5.08000 142.63688 162.79813
surf 209 cuboid -5.08000 5.08000 -15.24000 -5.08000 162.79813 182.95938
surf 210 cuboid -5.08000 5.08000 -15.24000 -5.08000 182.95938 204.15251
surf 211 cuboid -5.08000 5.08000 -15.24000 -5.08000 204.15251 225.34564
surf 212 cuboid -5.08000 5.08000 -15.24000 -5.08000 225.34564 246.53877
% surfaces for DF of assembly 3
surf 301 cuboid 5.08000 15.24000 -15.24000 -5.08000 0.00000 20.66396
surf 302 cuboid 5.08000 15.24000 -15.24000 -5.08000 20.66396 41.32792
surf 303 cuboid 5.08000 15.24000 -15.24000 -5.08000 41.32792 61.99188
surf 304 cuboid 5.08000 15.24000 -15.24000 -5.08000 61.99188 82.15313
surf 305 cuboid 5.08000 15.24000 -15.24000 -5.08000 82.15313 102.31438
surf 306 cuboid 5.08000 15.24000 -15.24000 -5.08000 102.31438 122.47563
surf 307 cuboid 5.08000 15.24000 -15.24000 -5.08000 122.47563 142.63688
surf 308 cuboid 5.08000 15.24000 -15.24000 -5.08000 142.63688 162.79813
surf 309 cuboid 5.08000 15.24000 -15.24000 -5.08000 162.79813 182.95938
surf 310 cuboid 5.08000 15.24000 -15.24000 -5.08000 182.95938 204.15251
surf 311 cuboid 5.08000 15.24000 -15.24000 -5.08000 204.15251 225.34564
surf 312 cuboid 5.08000 15.24000 -15.24000 -5.08000 225.34564 246.53877
% surfaces for DF of assembly 4
surf 401 cuboid -15.24000 -5.08000 -5.08000 5.08000 0.00000 20.66396
surf 402 cuboid -15.24000 -5.08000 -5.08000 5.08000 20.66396 41.32792
surf 403 cuboid -15.24000 -5.08000 -5.08000 5.08000 41.32792 61.99188
surf 404 cuboid -15.24000 -5.08000 -5.08000 5.08000 61.99188 82.15313
surf 405 cuboid -15.24000 -5.08000 -5.08000 5.08000 82.15313 102.31438
surf 406 cuboid -15.24000 -5.08000 -5.08000 5.08000 102.31438 122.47563
surf 407 cuboid -15.24000 -5.08000 -5.08000 5.08000 122.47563 142.63688
surf 408 cuboid -15.24000 -5.08000 -5.08000 5.08000 142.63688 162.79813
surf 409 cuboid -15.24000 -5.08000 -5.08000 5.08000 162.79813 182.95938
surf 410 cuboid -15.24000 -5.08000 -5.08000 5.08000 182.95938 204.15251
surf 411 cuboid -15.24000 -5.08000 -5.08000 5.08000 204.15251 225.34564
surf 412 cuboid -15.24000 -5.08000 -5.08000 5.08000 225.34564 246.53877
% surfaces for DF of assembly 5
surf 501 cuboid -5.08000 5.08000 -5.08000 5.08000 0.00000 20.66396
surf 502 cuboid -5.08000 5.08000 -5.08000 5.08000 20.66396 41.32792
surf 503 cuboid -5.08000 5.08000 -5.08000 5.08000 41.32792 61.99188
surf 504 cuboid -5.08000 5.08000 -5.08000 5.08000 61.99188 82.15313
surf 505 cuboid -5.08000 5.08000 -5.08000 5.08000 82.15313 102.31438
surf 506 cuboid -5.08000 5.08000 -5.08000 5.08000 102.31438 122.47563
surf 507 cuboid -5.08000 5.08000 -5.08000 5.08000 122.47563 142.63688
surf 508 cuboid -5.08000 5.08000 -5.08000 5.08000 142.63688 162.79813
surf 509 cuboid -5.08000 5.08000 -5.08000 5.08000 162.79813 182.95938
surf 510 cuboid -5.08000 5.08000 -5.08000 5.08000 182.95938 204.15251
surf 511 cuboid -5.08000 5.08000 -5.08000 5.08000 204.15251 225.34564
surf 512 cuboid -5.08000 5.08000 -5.08000 5.08000 225.34564 246.53877
% surfaces for DF of assembly 6
surf 601 cuboid 5.08000 15.24000 -5.08000 5.08000 0.00000 20.66396
surf 602 cuboid 5.08000 15.24000 -5.08000 5.08000 20.66396 41.32792
surf 603 cuboid 5.08000 15.24000 -5.08000 5.08000 41.32792 61.99188
surf 604 cuboid 5.08000 15.24000 -5.08000 5.08000 61.99188 82.15313
surf 605 cuboid 5.08000 15.24000 -5.08000 5.08000 82.15313 102.31438
surf 606 cuboid 5.08000 15.24000 -5.08000 5.08000 102.31438 122.47563
surf 607 cuboid 5.08000 15.24000 -5.08000 5.08000 122.47563 142.63688
surf 608 cuboid 5.08000 15.24000 -5.08000 5.08000 142.63688 162.79813
surf 609 cuboid 5.08000 15.24000 -5.08000 5.08000 162.79813 182.95938
surf 610 cuboid 5.08000 15.24000 -5.08000 5.08000 182.95938 204.15251
surf 611 cuboid 5.08000 15.24000 -5.08000 5.08000 204.15251 225.34564
surf 612 cuboid 5.08000 15.24000 -5.08000 5.08000 225.34564 246.53877
% surfaces for DF of assembly 7
surf 701 cuboid -15.24000 -5.08000 5.08000 15.24000 0.00000 20.66396
surf 702 cuboid -15.24000 -5.08000 5.08000 15.24000 20.66396 41.32792
surf 703 cuboid -15.24000 -5.08000 5.08000 15.24000 41.32792 61.99188
surf 704 cuboid -15.24000 -5.08000 5.08000 15.24000 61.99188 82.15313
surf 705 cuboid -15.24000 -5.08000 5.08000 15.24000 82.15313 102.31438
surf 706 cuboid -15.24000 -5.08000 5.08000 15.24000 102.31438 122.47563
surf 707 cuboid -15.24000 -5.08000 5.08000 15.24000 122.47563 142.63688
surf 708 cuboid -15.24000 -5.08000 5.08000 15.24000 142.63688 162.79813
surf 709 cuboid -15.24000 -5.08000 5.08000 15.24000 162.79813 182.95938
surf 710 cuboid -15.24000 -5.08000 5.08000 15.24000 182.95938 204.15251
surf 711 cuboid -15.24000 -5.08000 5.08000 15.24000 204.15251 225.34564
surf 712 cuboid -15.24000 -5.08000 5.08000 15.24000 225.34564 246.53877
% surfaces for DF of assembly 8
surf 801 cuboid -5.08000 5.08000 5.08000 15.24000 0.00000 20.66396
surf 802 cuboid -5.08000 5.08000 5.08000 15.24000 20.66396 41.32792
surf 803 cuboid -5.08000 5.08000 5.08000 15.24000 41.32792 61.99188
surf 804 cuboid -5.08000 5.08000 5.08000 15.24000 61.99188 82.15313
surf 805 cuboid -5.08000 5.08000 5.08000 15.24000 82.15313 102.31438
surf 806 cuboid -5.08000 5.08000 5.08000 15.24000 102.31438 122.47563
surf 807 cuboid -5.08000 5.08000 5.08000 15.24000 122.47563 142.63688
surf 808 cuboid -5.08000 5.08000 5.08000 15.24000 142.63688 162.79813
surf 809 cuboid -5.08000 5.08000 5.08000 15.24000 162.79813 182.95938
surf 810 cuboid -5.08000 5.08000 5.08000 15.24000 182.95938 204.15251
surf 811 cuboid -5.08000 5.08000 5.08000 15.24000 204.15251 225.34564
surf 812 cuboid -5.08000 5.08000 5.08000 15.24000 225.34564 246.53877
% surfaces for DF of assembly 9
surf 901 cuboid 5.08000 15.24000 5.08000 15.24000 0.00000 20.66396
surf 902 cuboid 5.08000 15.24000 5.08000 15.24000 20.66396 41.32792
surf 903 cuboid 5.08000 15.24000 5.08000 15.24000 41.32792 61.99188
surf 904 cuboid 5.08000 15.24000 5.08000 15.24000 61.99188 82.15313
surf 905 cuboid 5.08000 15.24000 5.08000 15.24000 82.15313 102.31438
surf 906 cuboid 5.08000 15.24000 5.08000 15.24000 102.31438 122.47563
surf 907 cuboid 5.08000 15.24000 5.08000 15.24000 122.47563 142.63688
surf 908 cuboid 5.08000 15.24000 5.08000 15.24000 142.63688 162.79813
surf 909 cuboid 5.08000 15.24000 5.08000 15.24000 162.79813 182.95938
surf 910 cuboid 5.08000 15.24000 5.08000 15.24000 182.95938 204.15251
surf 911 cuboid 5.08000 15.24000 5.08000 15.24000 204.15251 225.34564
surf 912 cuboid 5.08000 15.24000 5.08000 15.24000 225.34564 246.53877

%define ADF for Surface quantity output
% assembly 1
set adf 101 101 0
set adf 102 102 0
set adf 103 103 0
set adf 104 104 0
set adf 105 105 0
set adf 106 106 0
set adf 107 107 0
set adf 108 108 0
set adf 109 109 0
set adf 110 110 0
set adf 111 111 0
set adf 112 112 0
% assembly 2
set adf 201 201 0
set adf 202 202 0
set adf 203 203 0
set adf 204 204 0
set adf 205 205 0
set adf 206 206 0
set adf 207 207 0
set adf 208 208 0
set adf 209 209 0
set adf 210 210 0
set adf 211 211 0
set adf 212 212 0
% assembly 3
set adf 301 301 0
set adf 302 302 0
set adf 303 303 0
set adf 304 304 0
set adf 305 305 0
set adf 306 306 0
set adf 307 307 0
set adf 308 308 0
set adf 309 309 0
set adf 310 310 0
set adf 311 311 0
set adf 312 312 0
% assembly 4
set adf 401 401 0
set adf 402 402 0
set adf 403 403 0
set adf 404 404 0
set adf 405 405 0
set adf 406 406 0
set adf 407 407 0
set adf 408 408 0
set adf 409 409 0
set adf 410 410 0
set adf 411 411 0
set adf 412 412 0
% assembly 5
set adf 501 501 0
set adf 502 502 0
set adf 503 503 0
set adf 504 504 0
set adf 505 505 0
set adf 506 506 0
set adf 507 507 0
set adf 508 508 0
set adf 509 509 0
set adf 510 510 0
set adf 511 511 0
set adf 512 512 0
% assembly 6
set adf 601 601 0
set adf 602 602 0
set adf 603 603 0
set adf 604 604 0
set adf 605 605 0
set adf 606 606 0
set adf 607 607 0
set adf 608 608 0
set adf 609 609 0
set adf 610 610 0
set adf 611 611 0
set adf 612 612 0
% assembly 7
set adf 701 701 0
set adf 702 702 0
set adf 703 703 0
set adf 704 704 0
set adf 705 705 0
set adf 706 706 0
set adf 707 707 0
set adf 708 708 0
set adf 709 709 0
set adf 710 710 0
set adf 711 711 0
set adf 712 712 0
% assembly 8
set adf 801 801 0
set adf 802 802 0
set adf 803 803 0
set adf 804 804 0
set adf 805 805 0
set adf 806 806 0
set adf 807 807 0
set adf 808 808 0
set adf 809 809 0
set adf 810 810 0
set adf 811 811 0
set adf 812 812 0
% assembly 9
set adf 901 901 0
set adf 902 902 0
set adf 903 903 0
set adf 904 904 0
set adf 905 905 0
set adf 906 906 0
set adf 907 907 0
set adf 908 908 0
set adf 909 909 0
set adf 910 910 0
set adf 911 911 0
set adf 912 912 0

% --- MATERIALS
% --- Thermal Scattering Data
therm gre7 gre7.00t
therm lwe7 lwe7.00t
therm hwe7 hwe7.00t

% Graphite HEU Fuel
mat FUEL sum moder gre7 6000
6000.80c 5.1039E-02
6002.80c 3.5468E-02
5010.80c 1.4440E-07
5011.80c 5.8124E-07
8016.80c 1.8685E-05
23050.80c 1.5339E-09
23051.80c 6.1201E-07
26054.80c 2.9114E-07
26056.80c 4.5703E-06
26057.80c 1.0555E-07
26058.80c 1.4047E-08
92234.80c 8.5472E-08
92235.80c 8.7202E-06
92236.80c 4.0790E-08
92238.80c 4.9986E-07

% Zircaloy-3: Fuel can
mat ZR-3 sum
1001.80c 8.5822E-05
1002.80c 9.8707E-09
3006.80c 2.1501E-08
3007.80c 2.6177E-07
4009.80c 2.1817E-07
5010.80c 3.6192E-08
5011.80c 1.4568E-07
6000.80c 3.9289E-05
7014.80c 1.4266E-05
7015.80c 5.2119E-08
8016.80c 2.2604E-04
8017.80c 8.5927E-08
11023.80c 8.5526E-07
12024.80c 6.3901E-08
12025.80c 8.0898E-09
12026.80c 8.9068E-09
13027.80c 7.2873E-07
14028.80c 6.4564E-07
14029.80c 3.2799E-08
14030.80c 2.1647E-08
15031.80c 6.3480E-07
19039.80c 4.6899E-07
19040.80c 5.8838E-11
19041.80c 3.3846E-08
20040.80c 4.7559E-07
20042.80c 3.1742E-09
20043.80c 6.6231E-10
20044.80c 1.0234E-08
20046.80c 1.9624E-11
20048.80c 9.1742E-10
21045.80c 4.3737E-06
22046.80c 5.4221E-08
22047.80c 4.8898E-08
22048.80c 4.8451E-07
22049.80c 3.5556E-08
22050.80c 3.4044E-08
23050.80c 9.6494E-11
23051.80c 3.8501E-08
24050.80c 1.6430E-09
24052.80c 3.1685E-08
24053.80c 3.5928E-09
24054.80c 8.9432E-10
25055.80c 7.1579E-08
26054.80c 1.0701E-05
26056.80c 1.6799E-04
26057.80c 3.8795E-06
26058.80c 5.1630E-07
27059.80c 3.3363E-08
28058.80c 2.2806E-08
28060.80c 8.7847E-09
28061.80c 3.8186E-10
28062.80c 1.2175E-09
28064.80c 3.1007E-10
29063.80c 2.1396E-08
29065.80c 9.5455E-09
30064.80c 1.4510E-08
30066.80c 8.4094E-09
30067.80c 1.2331E-09
30068.80c 5.7187E-09
30070.80c 1.8968E-10
31069.80c 1.6951E-08
31071.80c 1.1250E-08
32070.80c 5.5164E-09
32072.80c 7.3923E-09
32073.80c 2.1005E-09
32074.80c 9.9394E-09
32076.80c 2.1194E-09
33075.80c 2.6244E-08
34074.80c 2.2162E-10
34076.80c 2.3333E-09
34077.80c 1.9000E-09
34078.80c 5.9191E-09
34080.80c 1.2354E-08
34082.80c 2.1739E-09
37085.80c 1.6603E-08
37087.80c 6.4024E-09
38084.80c 1.2567E-10
38086.80c 2.2126E-09
38087.80c 1.5708E-09
38088.80c 1.8531E-08
39089.80c 4.4231E-07
41093.80c 4.2327E-07
42092.80c 3.0270E-09
42094.80c 1.8916E-09
42095.80c 3.2586E-09
42096.80c 3.4184E-09
42097.80c 1.9592E-09
42098.80c 4.9576E-09
42100.80c 1.9818E-09
44096.80c 1.0778E-09
44098.80c 3.6379E-10
44099.80c 2.4823E-09
44100.80c 2.4512E-09
44101.80c 3.3189E-09
44102.80c 6.1377E-09
44104.80c 3.6223E-09
45103.80c 1.9107E-08
46102.80c 1.8846E-10
46104.80c 2.0582E-09
46105.80c 4.1257E-09
46106.80c 5.0495E-09
46108.80c 4.8887E-09
46110.80c 2.1654E-09
47107.80c 2.8348E-06
47109.80c 2.6336E-06
48106.80c 2.1864E-10
48108.80c 1.5567E-10
48110.80c 2.1847E-09
48111.80c 2.2389E-09
48112.80c 4.2207E-09
48113.80c 2.1374E-09
48114.80c 5.0253E-09
48116.80c 1.3101E-09
49113.80c 1.0285E-08
49115.80c 2.2946E-07
50112.80c 9.6398E-07
50114.80c 6.5590E-07
50115.80c 3.3789E-07
50116.80c 1.4450E-05
50117.80c 7.6323E-06
50118.80c 2.4070E-05
50119.80c 8.5367E-06
50120.80c 3.2378E-05
50122.80c 4.6013E-06
50124.80c 5.7541E-06
51121.80c 9.2384E-09
51123.80c 6.9099E-09
52120.80c 1.3868E-11
52122.80c 3.9293E-10
52123.80c 1.3714E-10
52124.80c 7.3040E-10
52125.80c 1.0894E-09
52126.80c 2.9031E-09
52128.80c 4.8909E-09
52130.80c 5.2515E-09
55133.80c 1.4794E-08
56130.80c 1.5177E-11
56132.80c 1.4461E-11
56134.80c 3.4606E-10
56135.80c 9.4383E-10
56136.80c 1.1245E-09
56137.80c 1.6082E-09
56138.80c 1.0266E-08
57138.80c 1.2740E-11
57139.80c 1.4142E-08
58136.80c 2.5961E-11
58138.80c 3.5222E-11
58140.80c 1.2412E-08
58142.80c 1.5596E-09
59141.80c 1.3954E-08
60142.80c 3.7077E-09
60143.80c 1.6630E-09
60144.80c 3.2443E-09
60145.80c 1.1314E-09
60146.80c 2.3446E-09
60148.80c 7.7699E-10
60150.80c 7.6336E-10
62144.80c 4.0146E-10
62147.80c 1.9602E-09
62148.80c 1.4698E-09
62149.80c 1.8072E-09
62150.80c 9.6506E-10
62152.80c 3.4980E-09
62154.80c 2.9750E-09
63151.80c 6.1860E-09
63153.80c 6.7527E-09
64152.80c 2.5008E-11
64154.80c 2.7258E-10
64155.80c 1.8506E-09
64156.80c 2.5595E-09
64157.80c 1.9568E-09
64158.80c 3.1059E-09
64160.80c 2.7333E-09
65159.80c 1.2372E-08
66156.80c 6.7759E-12
66158.80c 1.1495E-11
66160.80c 2.8180E-10
66161.80c 2.2855E-09
66162.80c 3.0824E-09
66163.80c 3.0124E-09
66164.80c 3.4194E-09
67165.80c 1.1921E-08
68162.80c 1.6340E-11
68164.80c 1.8821E-10
68166.80c 3.9385E-09
68167.80c 2.6884E-09
68168.80c 3.1714E-09
68170.80c 1.7527E-09
69169.80c 1.1639E-08
71175.80c 1.0947E-08
71176.80c 2.9105E-10
72174.80c 2.7848E-09
72176.80c 9.1550E-08
72177.80c 3.2373E-07
72178.80c 4.7481E-07
72179.80c 2.3706E-07
72180.80c 6.1057E-07
73180.80c 1.3039E-12
73181.80c 1.0865E-08
74180.80c 1.2834E-11
74182.80c 2.8342E-09
74183.80c 1.5305E-09
74184.80c 3.2770E-09
74186.80c 3.0407E-09
75185.80c 3.9492E-09
75187.80c 6.6101E-09
77191.80c 3.8155E-09
77193.80c 6.4137E-09
79197.80c 9.9825E-09
80196.80c 1.4703E-11
80198.80c 9.7728E-10
80199.80c 1.6536E-09
80200.80c 2.2643E-09
80201.80c 1.2919E-09
80202.80c 2.9269E-09
80204.80c 6.7341E-10
81203.80c 2.8399E-09
81205.80c 6.7803E-09
82204.80c 1.5942E-09
82206.80c 2.7444E-08
82207.80c 2.5166E-08
82208.80c 5.9670E-08
83209.80c 9.4086E-09
90232.80c 8.4737E-09
92234.80c 4.4606E-13
92235.80c 5.9508E-11
92238.80c 8.2005E-09
40090.80c 2.2022E-02
40091.80c 4.8024E-03
40092.80c 7.3405E-03
40094.80c 7.4390E-03
40096.80c 1.1985E-03

% AL-6063: End Tube
mat AL-6063 sum
12024.80c 3.5471E-04
12025.80c 4.4906E-05
12026.80c 4.9441E-05
13027.80c 5.8984E-02
14028.80c 2.1238E-04
14029.80c 1.0789E-05
14030.80c 7.1205E-06
22046.80c 1.3649E-06
22047.80c 1.2309E-06
22048.80c 1.2196E-05
22049.80c 8.9504E-07
22050.80c 8.5699E-07
24050.80c 6.7559E-07
24052.80c 1.3028E-05
24053.80c 1.4773E-06
24054.80c 3.6772E-07
25055.80c 1.4716E-05
26054.80c 2.9616E-06
26056.80c 4.6491E-05
26057.80c 1.0737E-06
26058.80c 1.4289E-07
27059.80c 6.8592E-06
28058.80c 4.6886E-06
28060.80c 1.8060E-06
28061.80c 7.8507E-08
28062.80c 2.5032E-07
28064.80c 6.3748E-08
29063.80c 8.7976E-06
29065.80c 3.9249E-06
30064.80c 5.9660E-06
30066.80c 3.4578E-06
30067.80c 5.0701E-07
30068.80c 2.3514E-06
30070.80c 7.7993E-08
50112.80c 3.3031E-08
50114.80c 2.2474E-08
50115.80c 1.1578E-08
50116.80c 4.9512E-07
50117.80c 2.6152E-07
50118.80c 8.2474E-07
50119.80c 2.9251E-07
50120.80c 1.1094E-06
50122.80c 1.5766E-07
50124.80c 1.9716E-07

% CP2 Graphite: Graphite plugs
mat GRPH-CP2 sum moder gre7 6000 moder lwe7 1001 moder hwe7 1002
1001.80c 2.2327E-05
1002.80c 2.5679E-09
3006.80c 7.3682E-09
3007.80c 8.9709E-08
4009.80c 2.7898E-09
5010.80c 6.4792E-09
5011.80c 2.6079E-08
6000.80c 8.3677E-02
7014.80c 9.7812E-06
7015.80c 3.5734E-08
8016.80c 1.1161E-05
8017.80c 4.2426E-09
9019.80c 2.6468E-09
11023.80c 1.0936E-09
12024.80c 1.9611E-09
12025.80c 2.4827E-10
12026.80c 2.7334E-10
13027.80c 8.2001E-08
14028.80c 2.3116E-06
14029.80c 1.1743E-07
14030.80c 7.7503E-08
15031.80c 2.1105E-08
16032.80c 1.3407E-06
16033.80c 1.0585E-08
16034.80c 5.9984E-08
16036.80c 1.4114E-10
17035.80c 9.6708E-09
17037.80c 3.0943E-09
19039.80c 1.1994E-09
19040.80c 1.5047E-13
19041.80c 8.6557E-11
20040.80c 3.8921E-06
20042.80c 2.5977E-08
20043.80c 5.4202E-09
20044.80c 8.3752E-08
20046.80c 1.6060E-10
20048.80c 7.5079E-09
21045.80c 1.1185E-10
22046.80c 2.0800E-08
22047.80c 1.8758E-08
22048.80c 1.8586E-07
22049.80c 1.3640E-08
22050.80c 1.3060E-08
23050.80c 5.9226E-09
23051.80c 2.3631E-06
24050.80c 2.1010E-10
24052.80c 4.0515E-09
24053.80c 4.5941E-10
24054.80c 1.1436E-10
25055.80c 3.6612E-10
26054.80c 6.5261E-09
26056.80c 1.0245E-07
26057.80c 2.3659E-09
26058.80c 3.1486E-10
27059.80c 8.5325E-11
28058.80c 1.6331E-08
28060.80c 6.2905E-09
28061.80c 2.7344E-10
28062.80c 8.7186E-10
28064.80c 2.2204E-10
29063.80c 1.6416E-09
29065.80c 7.3236E-10
30064.80c 1.8554E-10
30066.80c 1.0753E-10
30067.80c 1.5767E-11
30068.80c 7.3125E-11
30070.80c 2.4255E-12
31069.80c 4.3350E-11
31071.80c 2.8770E-11
32070.80c 7.0540E-11
32072.80c 9.4526E-11
32073.80c 2.6859E-11
32074.80c 1.2710E-10
32076.80c 2.7101E-11
33075.80c 3.3558E-10
34074.80c 2.8339E-12
34076.80c 2.9836E-11
34077.80c 2.4295E-11
34078.80c 7.5688E-11
34080.80c 1.5797E-10
34082.80c 2.7798E-11
35079.80c 3.1900E-10
35081.80c 3.1031E-10
37085.80c 2.1230E-10
37087.80c 8.1868E-11
38084.80c 7.7131E-11
38086.80c 1.3581E-09
38087.80c 9.6414E-10
38088.80c 1.1374E-08
39089.80c 5.6559E-10
40090.80c 1.9285E-09
40091.80c 4.2056E-10
40092.80c 6.4283E-10
40094.80c 6.5145E-10
40096.80c 1.0495E-10
41093.80c 2.7062E-10
42092.80c 3.8707E-11
42094.80c 2.4188E-11
42095.80c 4.1668E-11
42096.80c 4.3712E-11
42097.80c 2.5053E-11
42098.80c 6.3393E-11
42100.80c 2.5341E-11
44096.80c 1.3781E-11
44098.80c 4.6518E-12
44099.80c 3.1742E-11
44100.80c 3.1344E-11
44101.80c 4.2439E-11
44102.80c 7.8484E-11
44104.80c 4.6319E-11
45103.80c 4.8865E-11
46102.80c 2.4098E-12
46104.80c 2.6319E-11
46105.80c 5.2756E-11
46106.80c 6.4569E-11
46108.80c 6.2513E-11
46110.80c 2.7689E-11
47107.80c 1.2083E-10
47109.80c 1.1226E-10
48106.80c 2.7958E-12
48108.80c 1.9906E-12
48110.80c 2.7936E-11
48111.80c 2.8629E-11
48112.80c 5.3970E-11
48113.80c 2.7332E-11
48114.80c 6.4259E-11
48116.80c 1.6752E-11
49113.80c 9.3940E-12
49115.80c 2.0958E-10
50112.80c 2.0544E-12
50114.80c 1.3979E-12
50115.80c 7.2011E-13
50116.80c 3.0795E-11
50117.80c 1.6266E-11
50118.80c 5.1297E-11
50119.80c 1.8193E-11
50120.80c 6.9003E-11
50122.80c 9.8061E-12
50124.80c 1.2263E-11
51121.80c 1.1813E-10
51123.80c 8.8357E-11
52120.80c 1.7734E-13
52122.80c 5.0245E-12
52123.80c 1.7537E-12
52124.80c 9.3397E-12
52125.80c 1.3931E-11
52126.80c 3.7122E-11
52128.80c 6.2540E-11
52130.80c 6.7151E-11
53127.80c 3.9624E-11
55133.80c 1.8917E-10
56130.80c 1.9407E-13
56132.80c 1.8491E-13
56134.80c 4.4251E-12
56135.80c 1.2069E-11
56136.80c 1.4379E-11
56137.80c 2.0564E-11
56138.80c 1.3127E-10
57138.80c 1.6290E-13
57139.80c 1.8084E-10
58136.80c 3.3196E-13
58138.80c 4.5039E-13
58140.80c 1.5871E-10
58142.80c 1.9943E-11
59141.80c 1.7843E-10
60142.80c 4.7411E-11
60143.80c 2.1265E-11
60144.80c 4.1485E-11
60145.80c 1.4467E-11
60146.80c 2.9981E-11
60148.80c 9.9355E-12
60150.80c 9.7611E-12
62144.80c 1.0267E-12
62147.80c 5.0131E-12
62148.80c 3.7590E-12
62149.80c 4.6218E-12
62150.80c 2.4681E-12
62152.80c 8.9459E-12
62154.80c 7.6082E-12
63151.80c 1.5820E-11
63153.80c 1.7270E-11
64152.80c 6.3955E-14
64154.80c 6.9711E-13
64155.80c 4.7327E-12
64156.80c 6.5458E-12
64157.80c 5.0045E-12
64158.80c 7.9432E-12
64160.80c 6.9903E-12
65159.80c 3.1640E-11
66156.80c 1.7329E-14
66158.80c 2.9397E-14
66160.80c 7.2069E-13
66161.80c 5.8451E-12
66162.80c 7.8831E-12
66163.80c 7.7039E-12
66164.80c 8.7449E-12
67165.80c 3.0488E-11
68162.80c 4.1789E-14
68164.80c 4.8132E-13
68166.80c 1.0072E-11
68167.80c 6.8753E-12
68168.80c 8.1106E-12
68170.80c 4.4825E-12
71175.80c 2.7995E-11
71176.80c 7.4435E-13
72174.80c 4.5075E-14
72176.80c 1.4819E-12
72177.80c 5.2400E-12
72178.80c 7.6854E-12
72179.80c 3.8371E-12
72180.80c 9.8828E-12
73181.80c 2.7789E-07
74180.80c 1.6411E-13
74182.80c 3.6242E-11
74183.80c 1.9571E-11
74184.80c 4.1904E-11
74186.80c 3.8881E-11
75185.80c 1.0100E-11
75187.80c 1.6905E-11
77191.80c 9.7578E-12
77193.80c 1.6403E-11
79197.80c 2.5529E-10
80196.80c 3.7602E-13
80198.80c 2.4993E-11
80199.80c 4.2290E-11
80200.80c 5.7908E-11
80201.80c 3.3040E-11
80202.80c 7.4854E-11
80204.80c 1.7222E-11
81203.80c 3.6314E-11
81205.80c 8.6701E-11
82204.80c 1.6988E-12
82206.80c 2.9244E-11
82207.80c 2.6817E-11
82208.80c 6.3584E-11
83209.80c 1.2031E-10

% Air in IDAHO FALLS,ID :)
mat AIR -0.0012 moder lwe7 1001 moder hwe7 1002
1001.80c 5.0303E-07
1002.80c 5.7856E-11
2003.80c 1.7656E-16
2004.80c 1.3176E-10
6000.80c 9.5977E-09
7014.80c 3.9125E-05
7015.80c 1.4294E-07
8016.80c 1.0800E-05
8017.80c 4.1056E-09
18036.80c 7.9027E-10
18038.80c 1.4843E-10
18040.80c 2.3391E-07
36078.80c 1.0176E-13
36080.80c 6.5528E-13
36082.80c 3.3231E-12
36083.80c 3.2965E-12
36084.80c 1.6335E-11
36086.80c 4.9530E-12

% AL-1100: CR Bearing Tube
% mat AL-1100 -2.71
% 13027.03c 6.0032E-02
% 14028.03c 1.2727E-04
% 14029.03c 6.4656E-06
% 14030.03c 4.2672E-06
% 25055.03c 7.4265E-06
% 26054.03c 4.0568E-06
% 26056.03c 6.3683E-05
% 26057.03c 1.4707E-06
% 26058.03c 1.9572E-07
% 27059.03c 6.9230E-06
% 28058.03c 4.7323E-06
% 28060.03c 1.8229E-06
% 28061.03c 7.9238E-08
% 28062.03c 2.5265E-07
% 28064.03c 6.4342E-08
% 29063.03c 2.2199E-05
% 29065.03c 9.9036E-06
% 30000.03c 1.2475E-05
% 50112.03c 3.3338E-08
% 50114.03c 2.2684E-08
% 50115.03c 1.1686E-08
% 50116.03c 4.9973E-07
% 50117.03c 2.6396E-07
% 50118.03c 8.3242E-07
% 50119.03c 2.9523E-07
% 50120.03c 1.1197E-06
% 50122.03c 1.5913E-07
% 50124.03c 1.9900E-07

% Graphitar: Bushing
mat GRPHTAR -1.8 moder gre7 6000
6000.80c 1.0

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%% Control Information %%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Neutron population and criticality cycles

set pop 100000 1000 200
%set pop 1000 100 100

% Cross section library file path:

set acelib "/scratch/downar_fluxoe/vseker/ACE/endf71x/sss_endf71.xsdata"

% Relective BC
set bc 2 2 1
set sym 0

% 23 group structure from dragon article

%set nfg 23 5.0E-8 1.0E-7 1.4E-7 1.8E-7 2.5E-7 3.0E-7 3.5E-7 4.0E-7
% 5.0E-7 8.5E-7 9.72E-7 1.045E-6 1.097E-6 1.5E-6 4.0E-6
% 3.673E-4 9.118E-3 6.738E-2 1.11E-1 5.0E-1 1.353 3.679
%ene egrid 1 1E-11 5.0E-8 1.0E-7 1.4E-7 1.8E-7 2.5E-7 3.0E-7 3.5E-7 4.0E-7
% 5.0E-7 8.5E-7 9.72E-7 1.045E-6 1.097E-6 1.5E-6 4.0E-6
% 3.673E-4 9.118E-3 6.738E-2 1.11E-1 5.0E-1 1.353 3.679 20

set nfg 2 1E-4
ene egrid 1 1E-11 1E-4 20

% 14 group structure from dragon article

%set nfg 14 5E-8 1.4E-7 3.5E-7 5E-7 8.5E-7 9.72E-7 1.097E-6 4E-6
%3.673E-4 9.118E-3 6.738E-2 1.110E-1 5E-1
%
%ene egrid 1 1E-11 5E-8 1.4E-7 3.5E-7 5E-7 8.5E-7 9.72E-7 1.097E-6 4E-6
%3.673E-4 9.118E-3 6.738E-2 1.110E-1 5E-1 20

set micro egrid
% --- Plotters
% --- Geometry and mesh plots
plot 3 1000 1000
plot 2 500 4000
mesh 3 500 500

User avatar
Jaakko Leppänen
Site Admin
Posts: 2441
Joined: Thu Mar 18, 2010 10:43 pm
Security question 2: 0
Location: Espoo, Finland
Contact:

Re: Discountinuities on surface net current and flux

Post by Jaakko Leppänen » Sat Jun 13, 2020 1:12 pm

How large are the maximum differences that you see? The routine that scores the ADF currents checks the position where the boundary is crossed. This is not 100% accurate near the edges and corners, so there maybe some random inconsistencies between the corresponding inward and outward currents.
- Jaakko

yunlin
Posts: 4
Joined: Thu Dec 12, 2019 4:58 pm
Security question 1: No
Security question 2: 7

Re: Discountinuities on surface net current and flux

Post by yunlin » Tue Jun 23, 2020 9:25 pm

sorry for late response. The maximum difference I saw was about 1.0E-4, not very large, but large enough to be noticeable.

Post Reply