*----
*  TEST CASE TDCM58
*  MACROSCOPIC CROSS SECTIONS
*  4-group R120 lozenge Alamos geometry
*  Author: A. Hebert
*----
*  Define STRUCTURES and MODULES used
*----
MODULE G2S: SALT: MCCGT: MAC: ASM: FLU: TLM: DELETE: END: ;
SEQ_ASCII GT_ALAMOS :: FILE './30AV5_SSH.ap3' ;
SEQ_ASCII GT_SAL ;
LINKED_LIST TRACKF MACRO4G SYS FLUX ;
SEQ_BINARY INTLINS ;
SEQ_ASCII alamos1.ps ;
SEQ_ASCII alamos2.ps :: FILE './30AV5_SSH.ps' ;
SEQ_ASCII TDCM58.m :: FILE './TDCM58.m' ;
SEQ_BINARY INTLINS ;
PROCEDURE assertS ;
INTEGER an2d := 6 ;
REAL densur := 20.0 ;
*----
*  Input Macrolib
*----
MACRO4G := MAC: ::
 EDIT 1 NGRO 4 NMIX 6 NIFI 1 ANIS 2
 READ INPUT
 MIX     1 (* E110 *)
    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*) 4.102E-02 (*1->2*) 2.439E-05
          3 4 (*4->3*) 4.985E-04 (*3->3*) 4.956E-02 (*2->3*) 7.370E-05
          2 4 (*4->4*) 2.996E-02 (*3->4*) 3.248E-04
          1 1 (*1->1*) 0.000E+00
          1 2 (*2->2*) -6.161E-03
          1 3 (*3->3*) -2.326E-02
          1 4 (*4->4*) 0.000E+00
 MIX     2 (* E635 *)
    TOTAL  8.240E-02 2.0476E-01 3.7503E-01 2.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*) 4.102E-02 (*1->2*) 2.439E-05
          3 4 (*4->3*) 4.985E-04 (*3->3*) 4.956E-02 (*2->3*) 7.370E-05
          2 4 (*4->4*) 2.996E-02 (*3->4*) 3.248E-04
          1 1 (*1->1*) 0.000E+00
          1 2 (*2->2*) -6.161E-03
          1 3 (*3->3*) -2.326E-02
          1 4 (*4->4*) 0.000E+00
 MIX     3 (* G2.4 *)
    TOTAL  4.097E-00 6.286E+00 4.594E+00 7.508E+00
    TRANC  4.841E-02 7.939E-03 8.080E-03 8.038E-03
    NUSIGF 1.532E-02 9.259E-02 5.009E-01 9.008E-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
 MIX     4 (* U3.0 *)
    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 3.632E-02 1.259E-01 8.529E-01 9.898E-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
 MIX     5 (* void *)
    TOTAL  4.240E-05 1.0476E-04 2.7503E-03 1.1194E-04
    TRANC  0.000E+00 0.000E+00  0.000E+00  0.000E+00
     SCAT 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
          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     6 (* water *)
    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
;
*----
* Surfacic file conversion for RA120 rotation and isotropic reflection
*----
GT_SAL alamos1.ps := G2S: GT_ALAMOS :: ALAMOS -11 ;
*----
*  Tracking of surfacic file for equilateral triangular geometry
*----
TRACKF INTLINS := SALT: GT_SAL ::
  EDIT 1
  ALLG
  TISO <<an2d>> <<densur>> REND LONG
;

TRACKF := MCCGT: TRACKF INTLINS ::
  EDIT 2
  LCMD 3 AAC 150 TMT SCR 0 EPSI 1E-5
  MAXI 1 KRYL 0 HDD 0.0 ;
*----
*  Compute flux
*----
SYS := ASM: MACRO4G TRACKF INTLINS :: ARM EDIT 1 ;    

FLUX := FLU: MACRO4G SYS TRACKF INTLINS :: 
     EDIT 1 TYPE K EXTE 100 1E-5 ;

assertS FLUX :: 'K-EFFECTIVE' 1 1.033788 ;
GT_SAL TRACKF INTLINS SYS FLUX := DELETE: GT_SAL TRACKF INTLINS SYS FLUX ;

*----
* Surfacic file conversion for RA120 rotation and specular reflection
*----
GT_SAL alamos2.ps := G2S: GT_ALAMOS :: ALAMOS 11 ;
*----
*  Tracking of surfacic file for equilateral triangular geometry
*----
TRACKF INTLINS := SALT: GT_SAL ::
  EDIT 1
  ALLG
  TSPC EQW2 <<an2d>> <<densur>> REND LONG
;

TRACKF := MCCGT: TRACKF INTLINS ::
  EDIT 2
  LCMD 3 AAC 150 TMT SCR 0 EPSI 1E-5
  MAXI 1 KRYL 0 HDD 0.0 ;
*----
*  Compute flux
*----
SYS := ASM: MACRO4G TRACKF INTLINS :: ARM EDIT 1 ;    

FLUX := FLU: MACRO4G SYS TRACKF INTLINS :: 
     EDIT 1 TYPE K EXTE 100 1E-5 ;
*----
* Plot tracks
*----
TDCM58.m := TLM: TRACKF INTLINS ::
  EDIT 1
  NTPO 1
  DIRECTIONS NoPause DIR 0 ;

assertS FLUX :: 'K-EFFECTIVE' 1 1.041901 ;
ECHO "test TDCM58 completed" ;
END: ;
QUIT .