*-------------------------------------------------------------------------------------------------
*  Name          : TDCM69.x2m
*  Author        : B. Cui
*  MODEL TO TEST GLOW GEOMETRIES WITH MACRO-GEOMETRIES
*-------------------------------------------------------------------------------------------------
*  Define STRUCTURES and MODULES used
*-------------------------------------------------------------------------------------------------
SEQ_ASCII TDT_LAT_DAT :: FILE './irena_tiso_s30.dat' ;
SEQ_ASCII TDT_LAT_PS :: FILE './irena_tiso_s30.eps' ;
SEQ_ASCII TDT_LAT_M :: FILE './irena_tiso_s30.m' ;
MODULE
    G2S: MAC: SALT: TLM: ASM: FLU: EDI: DELETE: END: ;
LINKED_LIST
    MACRO4G_0 MACRO4G TRACK CP FLUX OUT ;
SEQ_BINARY
    TRACK_f ;
PROCEDURE assertS ;
*----
*  Input Macrolib
*----
MACRO4G_0 := MAC: ::
 EDIT 1 NGRO 4 NMIX 9 NIFI 1 ANIS 2
 READ INPUT
 MIX     1 (* NA *)
    TOTAL  5.316E-01 9.352E-01 9.931E-01 1.5931E+00
    TRANC  2.880E-01 5.529E-01 4.936E-01 4.4340E-01
     SCAT 1 1 (*1->1*) 4.973E-01
          2 2 (*2->2*) 6.570E-01 (*1->2*) 3.134E-02
          4 4 (*4->3*) 1.163E-02 (*3->3*) 7.176E-01 (*2->3*) 2.528E-01
                                                    (*1->3*) 2.304E-03
          4 4 (*4->4*) 1.565E+00 (*3->4*) 2.707E-01 (*2->4*) 2.412E-02
                                                    (*1->4*) 2.228E-04
          1 1 (*1->1*) 2.735E-01
          2 2 (*2->2*) 4.253E-01 (*1->2*) 1.410E-02
          4 4 (*4->3*) 5.064E-03 (*3->3*) 4.286E-01 (*2->3*) 1.203E-01
                                                    (*1->3*) 3.050E-04
          4 4 (*4->4*) 4.384E-01 (*3->4*) 6.426E-02 (*2->4*) 3.373E-03
                                                    (*1->4*) 1.245E-05
 MIX     2 (* CLAD *)
    TOTAL  4.240E-01 1.0476E+00 2.7503E+00 1.1194E+01
    TRANC  0.000E+00 0.000E+00  0.000E+00  0.000E+00
     SCAT 1 1 (*1->1*) 3.738E-01
          2 2 (*2->2*) 7.102E-01 (*1->2*) 2.439E-04
          3 4 (*4->3*) 4.985E-03 (*3->3*) 8.956E-01 (*2->3*) 7.370E-04
          2 4 (*4->4*) 8.996E-01 (*3->4*) 3.248E-03
          1 1 (*1->1*) 0.000E+00
          1 2 (*2->2*) -6.161E-02
          1 3 (*3->3*) -2.326E-01
          1 4 (*4->4*) 0.000E+00
 MIX     3 (* VOID *)
    TOTAL  4.240E-02 1.0476E-01 2.7503E-01 1.1194E+00
    TRANC  0.000E+00 0.000E+00  0.000E+00  0.000E+00
     SCAT 1 1 (*1->1*) 3.738E-02
          2 2 (*2->2*) 9.102E-02 (*1->2*) 2.439E-04
          3 4 (*4->3*) 4.985E-04 (*3->3*) 9.956E-02 (*2->3*) 7.370E-04
          2 4 (*4->4*) 9.996E-02 (*3->4*) 3.248E-04
          1 1 (*1->1*) 0.000E+00
          1 2 (*2->2*) 0.000E+00
          1 3 (*3->3*) 0.000E+00
          1 4 (*4->4*) 0.000E+00
 MIX     4 (* FUEL *)
    TOTAL  4.097E-01 6.286E-01 4.594E-01 7.508E-01
    TRANC  4.841E-02 7.939E-03 8.080E-03 8.038E-03
    NUSIGF 1.931E-02 1.259E-02 2.129E-01 2.508E-01
    CHI    1.000E+00 0.000E+00 0.000E+00 0.000E+00
     SCAT 1 1 (*1->1*) 3.908E-01
          2 2 (*2->2*) 4.407E-01 (*1->2*) 8.546E-04
          3 4 (*4->3*) 5.318E-03 (*3->3*) 3.785E-01 (*2->3*) 8.057E-03
          2 4 (*4->4*) 3.983E-01 (*3->4*) 1.203E-02
          1 1 (*1->1*) 4.793E-02
          2 2 (*2->2*) 5.419E-03  (*1->2*) -2.640E-04
          3 4 (*4->3*) -7.338E-04 (*3->3*) 9.014E-03 (*2->3*) -2.453E-03
          2 4 (*4->4*) 8.769E-03  (*3->4*) -2.572E-03
;
MACRO4G := MAC: MACRO4G_0 ::
 EDIT 1 NMIX 9
 MIX     5 4
 MIX     6 4
 MIX     7 4
 MIX     8 3
 MIX     9 2
;
*-----------------------------------------------------
*  Tracking of the imported geometry
*-----------------------------------------------------
TDT_LAT_PS := G2S: TDT_LAT_DAT :: DRAWMIX ;

TRACK TRACK_f := SALT: TDT_LAT_DAT ::
    TITL 'TDCM69 imported'
    EDIT 1
    ANIS 1
    TISO 8 10.0
    REND LONG
    IC EPSJ 1.0E-4
    ;
*-------------------------------------------------------------------------------------------------
*  Flux calculation IC
*-------------------------------------------------------------------------------------------------
CP := ASM: MACRO4G TRACK TRACK_f :: PIJ ;
FLUX := FLU: MACRO4G CP TRACK TRACK_f :: TYPE K ;
assertS FLUX :: 'K-EFFECTIVE' 1 1.008920 ;

TRACK TRACK_f CP FLUX := DELETE: TRACK TRACK_f CP FLUX ;

TRACK TRACK_f := SALT: TDT_LAT_DAT ::
    TITL 'TDCM69 imported TISO'
    EDIT 1
    ANIS 1
    TISO 8 10.0
    REND LONG
    NOIC
    ;

TDT_LAT_M := TLM: TRACK TRACK_f ::
    EDIT 1
    NTPO 1
    DIRECTIONS NoPause DIR 0
    ;
*-------------------------------------------------------------------------------------------------
*  Flux calculation TISO
*-------------------------------------------------------------------------------------------------
CP := ASM: MACRO4G TRACK TRACK_f :: PIJ ;
FLUX := FLU: MACRO4G CP TRACK TRACK_f :: TYPE K ;
assertS FLUX :: 'K-EFFECTIVE' 1 1.008918 ;

*-------------------------------------------------------------------------------------------------
*  Editing
*-------------------------------------------------------------------------------------------------
OUT := EDI: MACRO4G TRACK FLUX ::
    EDIT 4
    MERG COMP
    COND
    ;

ECHO "test TDCM69 completed" ;
END: ;
QUIT "LIST" .