Initial commit from Polytechnique Montreal

1 files changed, 119 insertions, 0 deletions

diff --git a/Dragon/data/tmatxs2_proc/TMATXS24.c2m b/Dragon/data/tmatxs2_proc/TMATXS24.c2m
new file mode 100644
index 0000000..f4ba048
--- /dev/null
+++ b/Dragon/data/tmatxs2_proc/TMATXS24.c2m
@@ -0,0 +1,119 @@
+*----
+*  TEST CASE TMATXS24
+*  Electron Beam on an intraoperative breast cancer lumpectomy.
+*  NJOY-[CEPXSR]: electron 3MeV - 100 Groups - P2 non-coupled library.
+*  This is a highly simplified Mobetron Benchmark. Please, note that
+*  deterministic parameters are not physical here.
+*  Author: Ahmed Naceur
+*----
+*  Define STRUCTURES and MODULES used
+*----
+MODULE       GEO: LIB: MAC: SNT: ASM: FLU: HEAT: OUT: DELETE: 
+             END: EDI: PSOUR: ;
+LINKED_LIST  MACRO DEPOS MICRO OLDMACRO GEOM TRACK SYSTEM FLUX 
+             ESOUR EEDIT EDITION EEDITc GWT0 ;
+SEQ_ASCII    edep_mob.txt   :: FILE 'edep_mob.txt' ;
+PROCEDURE    assertS ;
+
+* --- Geometry 
+
+GEOM := GEO: :: CAR1D 4
+     X- VOID  X+ VOID
+     MIX      1   2   3   1
+     MESHX 0.0 .6000 1.2000 1.4250 1.5000
+     SPLITX 200 150 150 50
+     ; 
+
+* --- Tracking       
+
+TRACK := SNT: GEOM ::
+      TITLE 'Intraoperative Mobetron Benchmark'
+      EDIT 0 MAXR 1000000  SCHM 1 DIAM 1 NLIVO SN 8 SCAT 8 EPSI 1.E-5
+      MAXI 300 BFPG ;
+
+* --- Library access
+
+MICRO  := LIB: :: 
+          EDIT 1
+          NMIX 3 
+          ANIS 3
+          STERN 1
+          MIXS LIB: MATXS2 FIL: _mobetron
+          * -- water  -- * 
+          MIX  1  300.0 1.0 NOGAS
+          H  = h 0.111894 
+          O  = o 0.888106 
+          * -- aluminium   -- * 
+          MIX  2  300.0 2.69890 NOGAS
+          Al = al 1.00  
+          * -- steel   -- *  
+          MIX  3  300.0 8.060 NOGAS 
+          Fe =  fe  0.641922
+          Cr =  cr  0.1837
+          Ni =  ni  0.1235
+          Mo =  mo  0.02386
+          Mn =  mn  0.01619
+          Si =  si  0.005093
+          Cu =  cu  0.001734
+          W  =  w   0.0011
+          Co =  co  0.001052
+          V  =  v   0.000508
+          C  =  c   0.000445
+          N  =  n   0.0004
+          P  =  p   0.000200
+          S  =  s   0.000175
+          As =  as  0.0001207  ;
+
+* --- Library structuring 
+
+MACRO := MICRO :: STEP UP MACROLIB ;
+
+OLDMACRO := MACRO ;
+MACRO    := MAC: MACRO OLDMACRO :: 
+            MIX 1 1 OLDL 
+            MIX 2 2 OLDL
+            MIX 3 3 OLDL ;
+
+OLDMACRO := DELETE: OLDMACRO ;
+
+* ---- Electron beam definition
+
+ESOUR := PSOUR: MACRO TRACK GEOM :: EDIT 1 MONO 1 X- INTG 1 1.0 
+                                    DIR 1.0 0.0 0.0 ; 
+
+* ---- Electron flux
+
+SYSTEM := ASM: MACRO TRACK :: EDIT 1 ARM ;
+FLUX   := FLU: MACRO TRACK SYSTEM ESOUR :: EDIT 1 TYPE S EXTE 100 ;
+
+* --- Recover flux and reaction rates
+
+EEDIT := OUT: FLUX TRACK MACRO GEOM :: EDIT 1 COND INTG NONE ;
+
+* --- Energy deposition from the e- source
+
+EEDITc := EEDIT ;
+DEPOS  := HEAT: EEDIT EEDITc :: EDIT 1 NORM 1.00 2.69890 8.060 ;
+
+* --- Save data structure for post-treatement 
+
+edep_mob.txt   := DEPOS ; 
+
+* --- AssertS
+
+GWT0 := EEDIT :: STEP UP GROUP STEP AT 1 ;
+assertS GWT0 :: 'FLUX-INTG' 1  3.311044E-03 ;
+assertS GWT0 :: 'FLUX-INTG' 30 3.885913E-03 ;
+assertS GWT0 :: 'FLUX-INTG' 60 4.248134E-03 ;
+assertS GWT0 :: 'H-FACTOR'  1  1.504317E+00 ;
+assertS GWT0 :: 'H-FACTOR'  30 1.609681E+00 ;
+assertS GWT0 :: 'H-FACTOR'  60 1.639779E+00 ;
+
+* --- Delete data structure 
+
+DEPOS EEDIT FLUX SYSTEM MACRO MICRO TRACK GEOM GWT0 := DELETE: 
+  DEPOS EEDIT FLUX SYSTEM MACRO MICRO TRACK GEOM GWT0 ; 
+
+*
+ECHO "test TMATXS24 completed" ;
+END: ;