*----
*  TEST CASE TCDR12
*  172-GROUP JEF2.2 DRAGLIB
*  GENERATE A LIBRARY USING INFO AND OTHER OPTIONS
*
*  REF: None
*
*----
*  Define variables and initialize
*    Coolant properties
*    a) Input
*       TempCool    = Coolant temperature (K)
*       Purity      = D2/(D2+H1) Weight % ratio in Coolant
*    b) Output
*       DensCool    = Coolant Density (g/cm**3)
*       WH1C        = H1 Weight % in Coolant
*       WD2C        = D2 Weight % in Coolant
*       WO16C       = O16 Weight % in Coolant
*    Fuel properties
*    a) Input
*       TempFuel    = Fuel temperature  (K)
*       Enrichment  = U235/(U235+U238) Weight % ratio in Fuel
*       DensFuel    = Fuel Density (g/cm**3)
*    b) Output
*       WU235F      = U235 Weight % in Fuel
*       WU238F      = U238 Weight % in Fuel
*       WO16F       = O16 Weight % in Fuel
*----
REAL
  TempCool Purity TempFuel Enrichment DensFuel   :=
  560.66   99.95  560.66   0.72       10.437501 ;
REAL
 WH1C   WD2C   WO16C  DensCool
 WU235F WU238F WO16F ;
*----
*  Define STRUCTURES and MODULES used
*----
LINKED_LIST
  LIBRARY ISOT ;
MODULE
  LIB: INFO: END: ;
PROCEDURE assertS ;
*----
*  Get Coolant properties
*----
ECHO
  "Input - Coolant temperature (K)             " TempCool ;
ECHO
  "Input - D2/(D2+H1) Weight % ratio in Coolant" Purity ;
INFO: ::
  TMP: <<TempCool>> K
  PUR: <<Purity>>   WGT%
  CALC DENS WATER >>DensCool<<
  LIB: DRAGON FIL: DLIB_J2
  ISO: 3 H1_H2O    H2_D2O   O16
  CALC WGT% D2O >>WH1C<< >>WD2C<< >>WO16C<<
  ;
ECHO
  "Output - Coolant Density (g/cm**3)" DensCool ;
ECHO
  "Output - H1 Weight % in Coolant   " WH1C ;
ECHO
  "Output - D2 Weight % in Coolant   " WD2C ;
ECHO
  "Output - O16 Weight % in Coolant  " WO16C ;
*----
*  Get Fuel properties
*----
ECHO
  "Input - Fuel temperature (K)                   " TempFuel ;
ECHO
  "Input - U235/(U235+U238) Weight % ratio in Fuel" Enrichment ;
ECHO
  "Input - Fuel Density (g/cm**3)                 " DensFuel ;
INFO: ::
  ENR: <<Enrichment>> WGT%
  LIB: DRAGON FIL: DLIB_J2
  ISO: 3 U235    U238      O16
  CALC WGT% UO2 >>WU235F<< >>WU238F<< >>WO16F<< 
  ;
ECHO
  "Output - U235 Weight % in Fuel" WU235F ;
ECHO
  "Output - U238 Weight % in Fuel" WU238F ;
ECHO
  "Output - O16 Weight % in Fuel " WO16F ;
*----
*  Microscopic cross sections from Jef 2.2 Draglib
*----
LIBRARY := LIB: ::
  NMIX 8 CTRA WIMS
  MIXS LIB: DRAGON FIL: DLIB_J2
  MIX  1  <<TempCool>> <<DensCool>>  O16      = O16      <<WO16C>>
    D2D2O    = H2_D2O   <<WD2C>>     H1H2O    = H1_H2O   <<WH1C>>
  MIX  2  <<TempFuel>> <<DensFuel>>  O16      = O16      <<WO16F>>
    U235     = U235     <<WU235F>> 1 
    U238     = U238     <<WU238F>> 1 
  MIX  3 COMB 1 0.5 0 0.5
  MIX  4 COMB 1 0.1 2 0.9
  ;
  
ISOT := LIBRARY :: STEP UP ISOTOPESLIST STEP AT 14 ; (*U238    0004*)
assertS ISOT :: NG 35 2.219074E-01 ;
ECHO "test TCDR12 completed" ;
END: ;
QUIT "LIST" .