summaryrefslogtreecommitdiff
path: root/Dragon/data/tmacro_proc/TCM53.c2m
diff options
context:
space:
mode:
Diffstat (limited to 'Dragon/data/tmacro_proc/TCM53.c2m')
-rw-r--r--Dragon/data/tmacro_proc/TCM53.c2m190
1 files changed, 190 insertions, 0 deletions
diff --git a/Dragon/data/tmacro_proc/TCM53.c2m b/Dragon/data/tmacro_proc/TCM53.c2m
new file mode 100644
index 0000000..b1dd7f5
--- /dev/null
+++ b/Dragon/data/tmacro_proc/TCM53.c2m
@@ -0,0 +1,190 @@
+*----
+* TEST CASE TCM53
+* MACROSCOPIC CROSS SECTIONS
+* MONTE-CARLO MULTIGROUP METHOD
+* CANDU-NG UNIT CELL WITH NXT:
+*
+*----
+* Define STRUCTURES and MODULES used
+*----
+MODULE
+ NXT: MAC: DELETE: MC: GREP: ABORT: SYBILT: ASM: FLU: GEO: END:
+ ;
+REAL VALUE REFVALUE DELTA ;
+PROCEDURE
+ CHB_NXT_MC
+ ;
+LINKED_LIST
+ GEOMF
+ MACRO
+ TRAC2
+ OUTMC
+ MACRO2 GEOM2 TRACK2 SYS FLUX
+ ;
+SEQ_ASCII
+ plot.ps
+ ;
+STRING
+ state SolveurSS SolveurFL polquad
+ Library
+ ;
+INTEGER
+ splitc splitm nban2d nban3d nban2dss
+ ;
+REAL
+ density densityss
+ ;
+INTEGER
+ CF11 CF12 CF13 CF14
+ CF21i CF22i CF23i CF24i
+ CF31i CF32i CF33i CF34i
+ CF41i CF42i CF43i CF44i
+ CF21e CF22e CF23e CF24e
+ CF31e CF32e CF33e CF34e
+ CF41e CF42e CF43e CF44e
+ VF11 VF12 VF13 VF14
+ VF21i VF22i VF23i VF24i
+ VF31i VF32i VF33i VF34i
+ VF41i VF42i VF43i VF44i
+ VF21e VF22e VF23e VF24e
+ VF31e VF32e VF33e VF34e
+ VF41e VF42e VF43e VF44e
+ COOL PT HEG CT
+ MOD SH VOID :=
+ 1 1 1 1
+ 2 2 2 2
+ 3 3 3 3
+ 4 4 4 4
+ 2 2 2 2
+ 3 3 3 3
+ 4 4 4 4
+ 1 1 1 1
+ 2 2 2 2
+ 3 3 3 3
+ 4 4 4 4
+ 2 2 2 2
+ 3 3 3 3
+ 4 4 4 4
+ 6 6 5 5
+ 6 5 5
+ ;
+*
+* geometrical parameters
+*
+EVALUATE state := "Voided" ; ! Cooled/Voided/Chboard
+EVALUATE splitc := 1 ;
+EVALUATE splitm := 1 ;
+
+EVALUATE nban2dss := 8 ;
+EVALUATE densityss := 20.0 ;
+
+EVALUATE nban2d := 12 ;
+EVALUATE density := 50.0 ;
+EVALUATE polquad := "LCMD" ;
+EVALUATE nban3d := 2 ;
+*
+* geometries
+*
+GEOMF := CHB_NXT_MC ::
+ <<state>>
+ <<splitc>> <<splitm>>
+ <<CF11>> <<CF12>> <<CF13>> <<CF14>>
+ <<CF21i>> <<CF22i>> <<CF23i>> <<CF24i>>
+ <<CF31i>> <<CF32i>> <<CF33i>> <<CF34i>>
+ <<CF41i>> <<CF42i>> <<CF43i>> <<CF44i>>
+ <<CF21e>> <<CF22e>> <<CF23e>> <<CF24e>>
+ <<CF31e>> <<CF32e>> <<CF33e>> <<CF34e>>
+ <<CF41e>> <<CF42e>> <<CF43e>> <<CF44e>>
+ <<VF11>> <<VF12>> <<VF13>> <<VF14>>
+ <<VF21i>> <<VF22i>> <<VF23i>> <<VF24i>>
+ <<VF31i>> <<VF32i>> <<VF33i>> <<VF34i>>
+ <<VF41i>> <<VF42i>> <<VF43i>> <<VF44i>>
+ <<VF21e>> <<VF22e>> <<VF23e>> <<VF24e>>
+ <<VF31e>> <<VF32e>> <<VF33e>> <<VF34e>>
+ <<VF41e>> <<VF42e>> <<VF43e>> <<VF44e>>
+ <<COOL>> <<PT>> <<HEG>> <<CT>>
+ <<MOD>> <<SH>> <<VOID>>
+ ;
+*
+* geometry tracking
+*
+ TRAC2 := NXT: GEOMF ::
+ EDIT 1 ANIS 1 MC PSPC
+ ;
+*
+* macrolib
+*
+MACRO := MAC: ::
+ EDIT 0 NGRO 2 NMIX 6 NIFI 1
+ READ INPUT
+ MIX 1
+ TOTAL 2.222E-01 8.3333E-01
+ NUSIGF 0.000E+00 2.3600E-01
+ CHI 1.000E+00 0.0000E+00
+ SCAT 1 1 0.1921 2 2 0.7533 0.02
+ MIX 2
+ TOTAL 2.222E-01 8.3333E-01
+ NUSIGF 0.000E+00 2.3600E-01
+ CHI 1.000E+00 0.0000E+00
+ SCAT 1 1 0.1921 2 2 0.7483 0.02
+ MIX 3
+ TOTAL 2.222E-01 7.5333E-01
+ NUSIGF 0.000E+00 2.36000E-01
+ CHI 1.000E+00 0.00000E+00
+ SCAT 1 1 0.1921 2 2 0.7033 0.02
+ MIX 4
+ TOTAL 2.322E-01 7.3333E-01
+ NUSIGF 0.000E+00 2.36000E-01
+ CHI 1.000E+00 0.00000E+00
+ SCAT 1 1 0.2121 2 2 0.7033 0.02
+ MIX 5
+ TOTAL 0.0 0.0
+ SCAT 1 1 0.0 2 2 0.0 0.0
+ MIX 6
+ TOTAL 0.1677 1.1211
+ SCAT 1 1 0.1265 2 2 1.1011 0.04
+ ;
+*
+* Monte-Carlo calculation
+*
+OUTMC := MC: TRAC2 MACRO ::
+ KCODE 300 400 6000
+ SEED 23412
+ EDIT 1
+ TALLY
+ MERG COMP
+ COND NONE
+ ENDT
+ ;
+MACRO2 := OUTMC :: STEP UP 'MACROLIB' ;
+GEOM2 := GEO: :: CARCEL 0
+ X- REFL X+ REFL Y- REFL Y+ REFL
+ MIX 1
+ MESHX 0.0 1.0
+ MESHY 0.0 1.0
+ ;
+TRACK2 := SYBILT: GEOM2 ::
+ TITLE 'INFINITE MEDIUM'
+ EDIT 1 MAXR 5 QUA1 5 ;
+SYS := ASM: MACRO2 TRACK2 ;
+FLUX := FLU: SYS MACRO2 TRACK2 :: TYPE K ;
+*
+* Non regression assert
+*
+EVALUATE REFVALUE := 1.0088 ;
+*
+GREP: FLUX :: GETVAL K-EFFECTIVE 1 >>VALUE<< ;
+EVALUATE DELTA := VALUE REFVALUE - REFVALUE / ABS ;
+IF DELTA 5.0E-3 < THEN
+ PRINT "TEST SUCCESSFUL; DELTA=" DELTA ;
+ELSE
+ PRINT "------------" ;
+ PRINT "TEST FAILURE" ;
+ PRINT "------------" ;
+ PRINT "REFERENCE=" REFVALUE " CALCULATED=" VALUE ;
+ ABORT: ;
+ENDIF ;
+
+ECHO "test TCM53 completed" ;
+END: ;
+QUIT "LIST" .